Pertanyaan yang diberi tag «finite-element»

Saya baru saja mulai bermain-main dengan FEniCS. Saya memecahkan Poisson dengan elemen urutan ke-3 dan ingin memvisualisasikan hasilnya. Namun, ketika saya menggunakan plot (u), visualisasi hanyalah interpolasi linier dari hasilnya. Saya mendapatkan hal yang sama ketika saya output ke VTK. Dalam kode lain yang saya kerjakan, saya menulis keluaran VTK …

14 finite-element  fenics 

Saya belajar teori di balik metode DG-FEM menggunakan buku Hesthaven / Warburton dan saya agak bingung tentang peran 'fluks numerik.' Saya minta maaf jika ini adalah pertanyaan mendasar, tetapi saya telah mencari dan tidak menemukan jawaban yang memuaskan untuk itu. Pertimbangkan persamaan gelombang skalar linier: ∂u∂t+∂f(u)∂x=0∂u∂t+∂f(u)∂x=0\frac{\partial u}{\partial t} + \frac{\partial …

13 finite-element  discontinuous-galerkin 

Dalam persamaan gelombang: c2∇⋅∇u(x,t)−∂2u(x,t)∂t2=f(x,t)c2∇⋅∇kamu(x,t)-∂2kamu(x,t)∂t2=f(x,t)c^2 \nabla \cdot \nabla u(x,t) - \frac{\partial^2 u(x,t)}{\partial t^2} = f(x,t) Mengapa kita pertama kali mengalikan dengan fungsi tes v(x,t)v(x,t)v(x,t) sebelum mengintegrasikan?

13 finite-element 

Saya sedang mengembangkan simulasi FEM. Untuk pengujian awal, saya akan menggunakan mesher sederhana yang ditulis sendiri dan visualisasi grafik mesh. Tapi saya ingin menyiapkan program saya untuk menggunakan data yang dihasilkan oleh mesher yang ada dan output ke alat visualisasi yang ada. Apakah ada standar (kuasi) yang direkomendasikan untuk format …

13 finite-element  mesh  graph-theory 

Mari kita mulai dengan masalah formulir (L+k2)u=0(L+k2)u=0(\mathcal{L} + k^2) u=0 dengan seperangkat kondisi batas yang diberikan ( Dirichlet , Neumann , Robin , Periodic , Bloch-Periodic ). Ini sesuai dengan menemukan nilai eigen dan vektor eigen untuk beberapa operator LL\mathcal{L} , dalam beberapa kondisi geometri, dan batas. Seseorang dapat memperoleh …

13 pde  finite-element  eigensystem  verification 

Pertanyaan: Metode apa yang tersedia untuk secara akurat dan efisien menghitung struktur sparsity dari matriks elemen hingga? Info: Saya sedang mengerjakan pemecah Persamaan Tekanan Poisson, menggunakan metode Galerkin dengan basis Lagrange kuadratik, ditulis dalam C, dan menggunakan PETSc untuk penyimpanan matriks jarang dan rutinitas KSP. Untuk menggunakan PETSc secara efisien, …

13 matrix  finite-element  petsc  performance 

\newcommand{\v}[1]{\boldsymbol{#1}} Misalkan kita memiliki persamaan model aliran Stokes berikut: { - d i v ( ν∇ u ) + ∇ hlmdivu=f=0{−div(ν∇u)+∇p=fdivu=0 \tag{1} \left\{ \begin{aligned} -\mathrm{div}(\nu \nabla \v{u}) + \nabla p &= \v{f} \\ \mathrm{div} \v{u} &= 0 \end{aligned} \right. di mana viskositas ν(x)ν(x)\nu(x) adalah fungsi, untuk elemen hingga campuran standar, …

13 finite-element 

Saya ingin menulis solver saya sendiri untuk persamaan Euler yang dapat dikompres, dan yang paling penting saya ingin itu bekerja dengan kuat dalam semua situasi. Saya ingin ini berbasis FE (DG tidak apa-apa). Apa metode yang mungkin? Saya sadar melakukan urutan ke-0 DG (volume terbatas) dan itu harus bekerja sangat …

13 pde  hyperbolic-pde  finite-element 

Saya tahu bahwa orang sering mengganti matriks massa yang konsisten dengan matriks diagonal yang disatukan. Di masa lalu, saya juga menerapkan kode di mana vektor beban dirakit dalam mode yang disatukan daripada mode yang konsisten dengan FEM. Tapi saya tidak pernah melihat mengapa kita diizinkan melakukan ini sejak awal. Apa …

13 finite-element 

Ini mungkin pertanyaan yang lebih cocok untuk sisi Rekomendasi Perangkat Lunak dari SE, namun saya percaya bahwa orang yang sering menggunakan bagian SE ini cenderung menjawab pertanyaan ini. Saya mencari alternatif gratis (tidak hanya dalam kebebasan) untuk Comsol Multiphysics . Berikut ini sedikit rumit: Saya tidak hanya mencari paket pemodelan …

12 finite-element  comsol  octave  openmodelica  modeling 

Sistem saya adalah masalah FE simetris dengan pengganda lagrange (misalnya aliran Stokes yang tidak dapat dikompresi): ( ABBTC)(SEBUAHBTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} di mana adalah kasus khas (saya bahkan telah memastikan bahwa persamaan diberi nomor sehingga pengali Lagrange muncul terakhir). Sistem ini cukup besar (garis + 100k).C= 0C=0C …

12 finite-element  petsc  krylov-method  preconditioning 

Ketika diskritisasi FEM dan menyelesaikan masalah difusi reaksi, misalnya, dengan (perturbasi singular), solusi dari masalah diskrit biasanya akan menunjukkan lapisan berosilasi dekat dengan batas. Dengan , dan elemen hingga linier, solusi terlihat seperti0 < ε ≪ 1 Ω = ( 0 , 1 ) ε = 10 - 5 u …

12 finite-element  discretization  numerical-analysis  singular-perturbation 

Saya ingin mengintegrasikan ekspresi polinomial melalui elemen 4-simpul dalam 3D. Beberapa buku tentang FEA membahas kasus di mana pengintegrasian dilakukan pada elemen flat 4-noned yang sewenang-wenang. Prosedur yang biasa dalam hal ini adalah menemukan matriks Jacobi dan menggunakannya untuk mengubah basis integrasi menjadi yang dinormalisasi di mana saya memiliki batas …

12 finite-element  quadrature 

Aliran bilangan Reynolds yang tinggi menghasilkan lapisan batas yang sangat tipis. Jika resolusi dinding digunakan dalam Simulasi Eddy Besar, rasio aspek mungkin berada di urutan . Banyak metode menjadi tidak stabil dalam rezim ini karena konstanta inf-sup menurun sebagai akar kuadrat dari rasio aspek atau lebih buruk. Konstanta inf-sup penting …

12 pde  finite-element  fluid-dynamics  stability  incompressible 

Saat membaca literatur tentang metode elemen hingga, istilah "node gantung" sering kali ditemukan. Adakah yang bisa memberi tahu saya apa itu simpul gantung?

11 finite-element  mesh-generation