คำถามติดแท็ก fenics

2
FeniCS: แสดงองค์ประกอบลำดับสูง
ฉันเพิ่งเริ่มสับสนกับ FEniCS ฉันกำลังแก้ไข Poisson ด้วยองค์ประกอบลำดับที่ 3 และต้องการให้เห็นภาพผลลัพธ์ อย่างไรก็ตามเมื่อฉันใช้ plot (u) การสร้างภาพข้อมูลเป็นเพียงการแก้ไขเชิงเส้นของผลลัพธ์ ฉันได้สิ่งเดียวกันเมื่อฉันส่งออกไปยัง VTK ในรหัสอื่นที่ฉันทำงานด้วยฉันเขียนเอาต์พุต VTK ที่จะเพิ่มองค์ประกอบการสั่งซื้อที่สูงขึ้นเพื่อให้พวกเขาดูลำดับที่สูงขึ้นใน Paraview มีอะไรเช่นนี้ (หรือดีกว่า) ใน FEniCS หรือไม่

3
สมการปัวซอง: กำหนดการไล่ระดับสีแบบเต็มเป็นเงื่อนไขขอบเขตผ่านตัวคูณแบบลากรองจ์
ฉันมีปัญหาทางกายภาพที่ควบคุมโดยสมการปัวซองในสองมิติ ฉันมีการวัดองค์ประกอบการไล่ระดับสีทั้งสองและตามส่วนของขอบเขตดังนั้นต้องการกำหนด และแพร่กระจายไปยังทุ่งไกล∂ u / ∂ x ∂ u / ∂ y Γ m ∂ u−∇2u=f(x,y),inΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y), \; in \; \Omega ∂u/∂x∂u/∂x\partial{u}/\partial{x}∂u/∂y∂u/∂y\partial{u}/\partial{y}ΓmΓm\Gamma_m∂u∂xi0=gm,onΓm∂u∂xi0=gm,onΓm \frac{\partial{u}}{\partial{x_i}}_0 = g_m, \; on \; \Gamma_m องค์ประกอบการไล่ระดับสีแทนเจนต์ , ฉันสามารถรวมเข้าด้วยกันแล้วบังคับใช้ผ่านเงื่อนไข Dirichlet เพื่อกำหนดองค์ประกอบปกติพร้อมกัน , ฉันรวบรวมฉันจะต้องผ่านตัวคูณลากรองจ์∫แกมมาเมตร∂ยู∂u∂x(t,0)∂u∂x(t,0)\frac{\partial{u}}{\partial{x}}_{(t,0)}∂ u∫Γm∂u∂x(t,0)ds=u0∫Γm∂u∂x(t,0)ds=u0 \int_{\Gamma_m}\frac{\partial{u}}{\partial{x}}_{(t,0)} \, ds = u_0 ∂u∂x(n,0)∂u∂x(n,0)\frac{\partial{u}}{\partial{x}}_{(n,0)} ดังนั้นฉันคิดว่ารูปแบบความแปรปรวนคือ ฉันใช้เวลานานในการพยายามรวมเข้าด้วยกันจากข้อมูลเกี่ยวกับปัญหาที่เกี่ยวข้องเช่น https: //answers.launchpad.net/fenics/+question/212434 …

2
Galerkin / Poisson / Fenics ไม่ต่อเนื่อง
ฉันพยายามที่จะแก้สมการปัวซอง 2D โดยใช้วิธีไม่ต่อเนื่อง Galerkin (DG) และการแยกประเภทต่อไปนี้ (ฉันมีไฟล์ png แต่ฉันไม่ได้รับอนุญาตให้อัปโหลดขออภัย): สมการ: ∇⋅(κ∇T)+f=0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 สมการใหม่: q=κ∇T∇⋅q=−fq=κ∇T∇⋅q=−fq = \kappa \nabla T\\\nabla \cdot q = -f รูปแบบที่อ่อนแอพร้อมฟลักซ์ตัวเลขและ : QT^T^\hat{T}q^q^\hat{q} ∫q⋅wdV=−∫T∇⋅(κw)dV+∫κT^n⋅wdS∫q⋅∇vdV=∫vfdV+∫q^⋅nvdS∫q⋅wdV=−∫T∇⋅(κw)dV+∫κT^n⋅wdS∫q⋅∇vdV=∫vfdV+∫q^⋅nvdS\int q \cdot w dV = - \int T \nabla \cdot (\kappa w) dV + \int \kappa \hat{T} …

3
มีแพ็คเกจ FEM“ น้ำหนักเบา” หรือไม่?
โดยทั่วไปแล้ว FEM ดูเหมือนจะเป็นปัญหาที่ "แก้ไข" ได้ค่อนข้างมาก มีเฟรมเวิร์กที่ทรงพลังมากมายเช่น Trilinos, PETSc, FEniCS, Libmesh หรือ MOOSE สิ่งหนึ่งที่พวกเขามีเหมือนกัน: พวกเขาเป็น "น้ำหนักมาก" ก่อนอื่นการติดตั้งปกติจะเจ็บปวดสุด ๆ ประการที่สองอินเทอร์เฟซ / API หนาและหนา - คุณต้องแปลความคิดทั้งหมดของคุณให้เป็นความคิดของห้องสมุดที่เกี่ยวข้อง นั่นหมายความว่าการทำงานร่วมกันและการขยายสำหรับข้อกำหนดพิเศษหรือรหัสที่มีอยู่เป็นเรื่องยาก โครงการอื่น ๆ เช่น (ตัวอย่างแบบสุ่ม) Boost, LibIGL, Aztec (ตัวแก้ปัญหาเชิงเส้น), Eigen หรือ CGAL แสดงให้เห็นว่าเป็นไปได้อย่างแน่นอนที่จะเขียนไลบรารีที่ทรงพลังที่รวมเข้ากับโค้ด C ++ หรือ Python ได้อย่างราบรื่นโดยไม่ต้องติดตั้ง ของกรอบหนักสุด มีแพ็คเกจที่เบามากสำหรับ FEM หรือไม่? ฉันไม่ได้มองหาวิธีแก้ปัญหาที่ง่ายและเป็นระบบอัตโนมัติ - ฉันกำลังมองหาห้องสมุดที่มีฟังก์ชั่นที่ทรงพลังในขณะที่รักษาอินเตอร์เฟสแบบลีนความสามารถในการทำงานร่วมกันกับโครงสร้างข้อมูลทั่วไป (เช่น C …
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.