การแสดงข้อมูลองค์ประกอบ Galerkin / ไฟไนต์ไม่ต่อเนื่อง

ฉันต้องการเห็นภาพผลการจำลองที่ได้รับโดยใช้วิธี Galerkin (DG) ที่ไม่ต่อเนื่องภายใน ParaView เช่นเดียวกับวิธีปริมาณ จำกัด โดเมนปัญหาแบ่งออกเป็นเซลล์รูปทรงลูกบาศก์ ("องค์ประกอบ") ซึ่งแตกต่างจากวิธีไฟไนต์วอลลุ่มภายในแต่ละเซลล์ไม่มีค่าเพียงค่าเดียวสำหรับเวกเตอร์โซลูชันแต่แต่ละเซลล์มีโซลูชันuที่จุดรวมหลาย Gauss$\mathbf{u}$$\mathbf{u}$

คำถามของฉันคือว่าทุกคนมีประสบการณ์ในการแสดงข้อมูลดังกล่าวอย่างมีประสิทธิภาพด้วย ParaView / VTK หรือไม่และคุณเลือกใช้วิธีใดเพื่อแสดงข้อมูลใน VTK มีหลายวิธีที่เป็นไปได้ในใจของฉัน แต่ฉันไม่รู้ว่าวิธีใดที่มีแนวโน้มมากที่สุด:

(1) ใช้ voxels
ใช้หนึ่ง voxel สำหรับแต่ละจุดรวม
Pro:ปลั๊กอินทั้งหมดที่ทำงานกับเซลล์ที่ไม่มีโครงสร้างแบบ VTK มาตรฐานจะยังคงทำงานต่อไปโดยไม่เปลี่ยนแปลงอะไรเลย
คอนดิชั่น:เนื่องจากจุดรวมไม่ได้กระจายอย่างสม่ำเสมอมันอาจเป็นเรื่องยากที่จะหาตำแหน่งที่ถูกต้องของจุดยอด นอกจากนี้ยังสามารถกำหนดโซลูชันได้สองครั้งบนพื้นผิวเซลล์เนื่องจากเฟรมเวิร์ก DG อนุญาตการแก้ปัญหาที่ไม่ต่อเนื่อง นอกจากนี้ข้อมูลลำดับชั้น (โดเมนแบ่งออกเป็นองค์ประกอบแต่ละองค์ประกอบมีหลายจุด) จะหายไป

(2) ใช้ polyvertices
ใช้หนึ่งจุดยอดต่อจุดรวม
Pro:ใช้งานง่ายที่สุดและง่ายต่อการระบุหลายจุดในตำแหน่งเดียวกันด้วยโซลูชันที่แตกต่างกัน
คอนดิชั่น:ความสามารถในการเห็นภาพข้อมูลเป็น "เซลล์" จะหายไปรวมทั้งข้อเสียเช่นเดียวกับข้างต้น

(3) ใช้โครงร่างการสร้างพื้นที่สี่เหลี่ยมจัตุรัส VTK
ใช้การสนับสนุนในตัวสำหรับโครงร่างการสร้างพื้นที่สี่เหลี่ยมจัตุรัส
Pro:การใช้งานค่อนข้างตรงไปตรงมารักษาความสัมพันธ์และคุณสมบัติทั้งหมดของโซลูชันดั้งเดิม
คอนดิชั่น:เนื่องจากนี่เป็นเซลล์ประเภทใหม่ทั้งหมดปลั๊กอินส่วนใหญ่ (ส่วนใหญ่) ที่มีอยู่จะไม่ทำงานอีกต่อไปและอาจต้องเขียนใหม่

deal.II ใช้ตัวเลือก (2) มาเป็นเวลานานแล้วและประสบความสำเร็จเป็นอย่างดี กล่าวอีกนัยหนึ่งใน 2d ทุกจุดยอดปรากฏ 4 ครั้งในไฟล์เอาต์พุตซึ่งหมายความว่ามีข้อมูลเพิ่มเติมที่จะเขียน แต่ฟิลด์ไม่ต่อเนื่อง ข้อเสียเปรียบอื่น ๆ ที่ฉันทราบคือโปรแกรมสร้างภาพบางโปรแกรมไม่สามารถติดตามส่วนต่อประสานระหว่างเซลล์ที่ไม่แบ่งปันจุดยอดเมื่อทำการปรับปรุงประสิทธิภาพ แต่วันนี้ดูเหมือนว่าจะดีขึ้น

คุณอาจต้องการดูรายการนี้ในคำถามที่พบบ่อยของเรา: http://dealii.sourceforge.net/index.php/Deal.II_Questions_and_Answers#The_graphical_output_files_don.27t_make_sense_to_me_E_80_4

อัปเดต: 31 พฤษภาคม 2018

ลิงก์ถูกย้ายไปที่: https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#the-graphical-output-files-dont-make-sense-to-me----they- ดูเหมือนการมีมากเกินไปหลายองศาของเสรีภาพ