การหาผลเฉลยเชิงตัวเลขของสมการเชิงอนุพันธ์ย่อย

สาขาวิชานี้พัฒนาวิธีการที่ใช้การทำให้เป็นส่วนย่อย (discretize) สมการเชิงอนุพันธ์ย่อยในเชิงพื้นที่และเวลา โดยแทนที่ตัวดำเนินการต่อเนื่องด้วยระบบพีชคณิต ซึ่งผลเฉลยของระบบดังกล่าวจะประมาณพฤติกรรมของฟิลด์ที่ควบคุมโดยกฎทางฟิสิกส์

Definition

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

Scope

ครอบคลุมกรอบการทำให้เป็นส่วนย่อยหลักสามแบบ ได้แก่ วิธีผลต่างจำกัด (finite difference), วิธีไฟไนต์เอลิเมนต์ (finite element) และวิธีไฟไนต์วอลุ่ม (finite volume) ซึ่งนำไปใช้กับสมการเชิงวงรี (elliptic), พาราโบลา (parabolic) และไฮเพอร์โบลา (hyperbolic); การวิเคราะห์ความสอดคล้อง (consistency), เสถียรภาพ (stability) และการลู่เข้า (convergence) (รวมถึงทฤษฎีบทสมมูลของ Lax และเงื่อนไข CFL); และระบบเชิงเส้นและไม่เชิงเส้นขนาดใหญ่แบบสปาร์ส (sparse) ที่เกิดจากการทำให้เป็นส่วนย่อย

Sub-topics

• Finite Difference Methods
• Finite Element Methods
• Finite Volume Methods

Core questions

• ตัวดำเนินการเชิงอนุพันธ์ในเชิงพื้นที่และเวลาถูกทำให้เป็นส่วนย่อยเป็นระบบพีชคณิตที่มีเสถียรภาพและลู่เข้าได้อย่างไร?
• ความสอดคล้องและเสถียรภาพรวมกันเพื่อรับประกันการลู่เข้าได้อย่างไร ดังเช่นในทฤษฎีบทสมมูลของ Lax?
• ประเภทของ PDE — เชิงวงรี, พาราโบลา หรือไฮเพอร์โบลา — กำหนดวิธีการที่เหมาะสมและข้อจำกัดด้านเสถียรภาพอย่างไร?
• ระบบสปาร์สขนาดใหญ่ที่เกิดขึ้นถูกแก้ปัญหาอย่างมีประสิทธิภาพได้อย่างไร?

Key theories

ทฤษฎีบทสมมูลของ Lax

สำหรับการประมาณผลต่างจำกัดที่สอดคล้องกับปัญหาค่าเริ่มต้นเชิงเส้นที่มีเงื่อนไขดี (well-posed linear initial value problem) เสถียรภาพเป็นสิ่งจำเป็นและเพียงพอสำหรับการลู่เข้า; ทฤษฎีบทนี้เป็นรากฐานที่ลดการพิสูจน์การลู่เข้าให้เหลือเพียงการตรวจสอบความสอดคล้องและเสถียรภาพ

เงื่อนไขเสถียรภาพและเลข CFL

แผนการเชิงประจักษ์ (explicit schemes) สำหรับ PDE ที่ขึ้นกับเวลาจะมีเสถียรภาพภายใต้ข้อจำกัดของขนาดขั้นเท่านั้น; สำหรับปัญหาไฮเพอร์โบลา เงื่อนไข Courant-Friedrichs-Lewy กำหนดให้โดเมนการพึ่งพาเชิงตัวเลขต้องมีโดเมนทางกายภาพอยู่ภายใน ซึ่งจำกัดขนาดขั้นเวลาเมื่อเทียบกับตาข่ายเชิงพื้นที่

หลักการแปรผันและการอนุรักษ์

วิธีการไฟไนต์เอลิเมนต์อาศัยการกำหนดรูปแบบอ่อน (weak (variational) formulations) และการฉายภาพ Galerkin ในขณะที่วิธีการไฟไนต์วอลุ่มบังคับใช้กฎการอนุรักษ์แบบไม่ต่อเนื่อง; แต่ละกรอบการทำงานให้แนวทางในการทำให้เป็นส่วนย่อยที่สอดคล้องกับคุณสมบัติการประมาณค่าที่พิสูจน์ได้

Clinical relevance

วิธีการ PDE เชิงตัวเลขเป็นรากฐานการคำนวณของการจำลองในสาขาวิศวกรรมศาสตร์และวิทยาศาสตร์กายภาพ — กลศาสตร์โครงสร้างและของแข็ง, พลศาสตร์ของไหลและอากาศพลศาสตร์, การถ่ายเทความร้อน, แม่เหล็กไฟฟ้า, ธรณีฟิสิกส์, การสร้างแบบจำลองสภาพอากาศและภูมิอากาศ, และการสร้างภาพทางการแพทย์ — ในทุกกรณีที่ต้องแก้สมการฟิลด์ต่อเนื่องบนรูปทรงเรขาคณิตที่ซับซ้อนซึ่งไม่สามารถหาผลเฉลยแบบปิดได้

History

การวิเคราะห์ PDE ด้วยวิธีผลต่างจำกัดเริ่มต้นด้วยบทความของ Courant-Friedrichs-Lewy ในปี 1928; วิธีไฟไนต์เอลิเมนต์เกิดขึ้นจากการวิศวกรรมโครงสร้างและคณิตศาสตร์เชิงแปรผันในช่วงทศวรรษ 1940-60 และวิธีไฟไนต์วอลุ่มพัฒนาควบคู่ไปกับพลศาสตร์ของไหลเชิงคำนวณ โดยมีทฤษฎีบทสมมูลของ Lax เป็นกรอบการลู่เข้าที่เป็นหนึ่งเดียวกันในช่วงทศวรรษ 1950

Key figures

• Richard Courant
• Peter Lax
• Olga Ladyzhenskaya
• Randall J. LeVeque

Seminal works

• morton2005
• leveque2007

Frequently asked questions

เหตุใดจึงมีกรอบการทำให้เป็นส่วนย่อยที่แตกต่างกันสามแบบ?

ผลต่างจำกัดนั้นง่ายที่สุดบนกริดปกติ, ไฟไนต์เอลิเมนต์จัดการกับรูปทรงเรขาคณิตที่ซับซ้อนและปัญหาเชิงแปรผันได้อย่างเป็นธรรมชาติ, และไฟไนต์วอลุ่มบังคับใช้การอนุรักษ์เฉพาะที่ ทำให้เหมาะสำหรับการไหลของของไหล การเลือกขึ้นอยู่กับรูปทรงเรขาคณิต, ประเภทของสมการ, และคุณสมบัติที่ต้องรักษาไว้

เงื่อนไข CFL หมายความว่าอย่างไร?

สำหรับแผนการเชิงประจักษ์ในปัญหาไฮเพอร์โบลาที่ขึ้นกับเวลา เงื่อนไข Courant-Friedrichs-Lewy จะจำกัดขนาดขั้นเวลาเมื่อเทียบกับระยะห่างของกริดเชิงพื้นที่ เพื่อให้แน่ใจว่าข้อมูลจะไม่เดินทางเกินหนึ่งเซลล์กริดต่อขั้น การละเมิดเงื่อนไขนี้จะทำให้เกิดความไม่เสถียร