สภาพ CFL ในรูปแบบไม่ต่อเนื่อง Galerkin


9

ฉันได้ใช้ระบบ Galerkin ของ ADER-Discontinuous สำหรับการแก้ปัญหาระบบเชิงเส้นของกฎหมายการอนุรักษ์ประเภทและสังเกตว่าเงื่อนไข CFL นั้นเข้มงวดมาก ในบรรณานุกรมขอบเขตบนของขั้นตอนเวลาสามารถพบได้โดยที่คือขนาดเซลล์คือจำนวนของ ขนาดและคือระดับสูงสุดของพหุนามtU+AxU+ByU=0Δthd(2N+1)λmaxhdN

มีวิธีใดที่จะหลีกเลี่ยงปัญหานี้หรือไม่? ฉันทำงานร่วมกับรูปแบบปริมาณ จำกัด ของ WENO-ADER และข้อ จำกัด CFL นั้นผ่อนคลายมากขึ้น ตัวอย่างเช่นสำหรับชุดรูปแบบลำดับที่ 5 จะต้องกำหนด CFL ที่ต่ำกว่า 0.04 เมื่อใช้ DG ในขณะที่ CFL = 0.4 ยังคงสามารถใช้ในรูปแบบ WENO-ADER FV

เหตุใดการใช้รูปแบบ DG แทนที่จะเป็น ADER-FV เช่นในการคำนวณ aeroacustics (สมการเชิงเส้นออยเลอร์เชิงเส้น) หรือแอปพลิเคชันที่คล้ายกัน (พลศาสตร์ของแก๊ส, น้ำตื้น, ค่าใช้จ่ายในการคำนวณโดยรวมของโครงการคล้ายกับของ ADER-FV ทั้งๆที่ขั้นตอนเวลาที่ต่ำกว่ามากหรือไม่?

ความคิดและข้อเสนอแนะสำหรับเรื่องนี้ยินดีต้อนรับ

คำตอบ:


6

CFL ที่ จำกัด ของชุดรูปแบบ DG มักจะมาจากการรวมกันของความแม่นยำในการสั่งซื้อสูงและมีขนาดกะทัดรัด (ดูตัวอย่างอ้างอิงนี้ ) CFL ขึ้นอยู่กับขอบเขตของรูปแบบการแปรผันในรูปแบบของ norm ของการแก้ปัญหาซึ่งขึ้นอยู่กับอนุพันธ์และร่องรอยของพหุนาม ขอบเขตสำหรับแต่ละปริมาณเหล่านี้ (หรือใช้สเตนมาร์คอฟพี่น้องความไม่เท่าเทียมและความไม่เท่าเทียมกันร่องรอยที่ไม่ต่อเนื่อง) ให้คงที่ซึ่งขึ้นอยู่ผกผันในและ quadratically คำสั่งซื้อที่ผลใน CFL โดยรวมของ2)L2hNO(h/N2)

FYI - ฉันเคยเห็น CFL ที่คุณพูดถึงการอ้างอิงมาก่อน แต่ฉันจำไม่ได้ว่ามันได้รับการพิสูจน์จากที่ใด ฉันต้องการทราบวิธีที่พวกเขาหลีกเลี่ยงการพึ่งพากำลังสองในในขอบเขตของพวกเขาN

จำกัด ที่แตกต่างกันและ Weno รูปแบบ (เช่นเดียวกับ B-อิสระตามระเบียบวิธีไฟไนต์เอลิเมน ต์ บนตาข่ายเป็นระยะ) มีเงื่อนไขที่โยก CFL เนื่องจากค่าคงที่ในขอบเขตคล้ายเติบโตช้าลงในNนี่คือการเปิดใช้เนื่องจากขนาด stencil มีแนวโน้มเพิ่มขึ้นตามลำดับซึ่งช่วยลดปัญหาเหล่านี้บางอย่างNN

วิธีการของ DG มีราคาแพงกว่า แต่สามารถจัดการกับตาข่ายที่ไม่มีโครงสร้างได้อย่างง่ายดายและสามารถใช้งานได้อย่างมีประสิทธิภาพ มีคำสั่ง WENO ในระดับสูง (หรือไทปันที่คล้ายกัน) สำหรับกริดที่ไม่มีโครงสร้างแม้ว่าสิ่งเหล่านี้สามารถแนะนำภาวะแทรกซ้อนทางคณิตศาสตร์หรือการนำไปใช้เพิ่มเติม


ขอบคุณมากสำหรับคำตอบโดยละเอียดของคุณ Jesse มันทำให้ฉันมีมุมมองที่กว้างขึ้นในเรื่องนี้ ในการทดลองเชิงตัวเลขของฉันกับ DG-ADER ฉันสังเกตเห็นว่าเมื่อใช้ตาข่ายรูปสี่เหลี่ยมขนมเปียกปูนที่มีโครงสร้าง (มีรูปร่างรูปสี่เหลี่ยมขนมเปียกปูนโดยพลการเช่นสี่เหลี่ยมจตุรัสสี่เหลี่ยมคางหมูหรือสี่เหลี่ยมด้านขนาน ... ) โซลูชันเชิงตัวเลขไม่ได้มีการแกว่ง อย่างไรก็ตามเมื่อย้ายไปที่ตาข่ายที่ไม่มีโครงสร้างการแกว่งจะปรากฏขึ้นแม้จะเป็นแบบกึ่งโครงสร้างที่สร้างขึ้นโดยการสุ่มแทนที่โหนดของตาข่ายแบบโครงสร้างเป็นระยะทางสั้น ๆ นี่เป็นพฤติกรรมที่คาดหวังหรือไม่?
Adr

1
@ เอเดรีย - มันค่อนข้างธรรมดาสำหรับความผันผวนที่จะปรากฏขึ้นเมื่อคุณย้ายออกไปจากตาข่ายเครื่องแบบ เมื่อคุณใช้ตาข่ายทั่วไปก็ยังไม่ได้เลยล้างใด ๆ เพิ่มเติมว่าสิ่งที่คุณหมายถึงโดยขนาดตาข่ายชั่วโมงมันอาจเป็นเส้นผ่าศูนย์กลางเซลล์ความยาวของขอบที่สั้นที่สุดสแควร์รูทของพื้นที่ (in 2d) หรือวิธีอื่นใดในการกำหนด "ขนาดตาข่าย" h
Wolfgang Bangerth

ขอบคุณสำหรับการชี้แจง Wolfgang จนถึงตอนนี้ฉันตั้งค่าเป็นความยาวของขอบที่สั้นที่สุด แต่ต่อไปแม้ว่าการลดจำนวน CFL อันดับหนึ่งของขนาดหรือมากกว่าจาก CFL ที่กำหนดโดยสูตรก็ยังคงผันผวน h
Adr
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.