«finite-element» etiketlenmiş sorular

FEniCS ile uğraşmaya yeni başladım. Poisson'u 3. dereceden elemanlarla çözüyorum ve sonuçları görselleştirmek istiyorum. Ancak, çizim (u) kullandığımda, görselleştirme sadece sonuçların doğrusal enterpolasyonudur. VTK'ya çıkış yaptığımda da aynı şeyi alıyorum. Birlikte çalıştığım başka bir kodda, aslında Paraview'de daha yüksek sıraya bakmaları için daha yüksek dereceli öğeleri örnekleyen bir VTK çıktısı …

14 finite-element  fenics 

Hesthaven / Warburton kitabını kullanarak DG-FEM yöntemlerinin arkasındaki teoriyi öğreniyorum ve 'sayısal akının' rolü hakkında biraz kafam karıştı. Bu temel bir soru ise özür dilerim, ama baktım ve tatmin edici bir cevap bulamadım. Doğrusal skaler dalga denklemini düşünün: burada lineer akıf(u)=au olarak verilir.∂u∂t+ ∂f( u )∂x= 0∂u∂t+∂f(u)∂x=0\frac{\partial u}{\partial t} + …

13 finite-element  discontinuous-galerkin 

Dalga denkleminde: c2∇⋅∇u(x,t)−∂2u(x,t)∂t2=f(x,t)c2∇⋅∇u(x,t)−∂2u(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) Entegrasyondan önce neden bir test fonksiyonuyla v(x,t)v(x,t)v(x,t) çarpıyoruz?

13 finite-element 

Bir FEM simülasyonu geliştiriyorum. Erken test için, kendi kendine yazılan basit bir yazı ve kafes grafiğinin görselleştirilmesini kullanacağım. Ama programımı mevcut bir mesher tarafından oluşturulan verileri kullanacak şekilde hazırlamak ve mevcut görselleştirme araçlarına göndermek istiyorum. (FEM) kafesleri için dosya biçimi ve dahili veri biçimi için önerilen (yarı) bir standart var …

13 finite-element  mesh  graph-theory 

Formun bir problemiyle başlayalım (L+k2)u=0(L+k2)u=0(\mathcal{L} + k^2) u=0 verilen bir dizi sınır koşulu ile ( Dirichlet , Neumann , Robin , Periyodik , Bloch-Periyodik ). Bu , bazı geometri ve sınır koşulları altında bazı operatör için özdeğerlerin ve özvektörlerin bulunmasına karşılık gelir . Örneğin akustik, elektromanyetizma, elastodinamik, kuantum mekaniğinde böyle …

13 pde  finite-element  eigensystem  verification 

Soru: Sonlu elemanlar matrisinin seyreklik yapısını doğru ve verimli bir şekilde hesaplamak için hangi yöntemler mevcuttur? Bilgi: Ben Poisson Basınç Denklemi çözücü üzerinde çalışıyorum, Galerkin yöntemi kuadratik Lagrange temelli, C ile yazılmış ve seyrek matris depolama ve KSP rutinleri için PETSc kullanarak. PETSc'yi verimli bir şekilde kullanmak için, belleği global …

13 matrix  finite-element  petsc  performance 

\newcommand{\v}[1]{\boldsymbol{#1}} Diyelim ki aşağıdaki Stokes akış modeli denklemine sahibiz: {−div(ν∇u)+∇pdivu=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. Burada viskozite ν(x)ν(x)\nu(x) bir işlevdir, standart karışık sonlu eleman için, sabit çifti kullandığımızı varsayalım: Crouzeix-Raviart alanı \ v {V} _h VhVh\v{V}_hhız için …

13 finite-element 

Sıkıştırılabilir Euler denklemleri için kendi çözücümü yazmak istiyorum ve en önemlisi her durumda sağlam çalışmasını istiyorum. FE tabanlı olmasını istiyorum (DG tamam). Olası yöntemler nelerdir? 0. derece DG (sonlu hacimler) yaptığımın farkındayım ve bu çok sağlam bir şekilde çalışmalı. Temel bir FVM çözücü uyguladım ve harika çalışıyor, ancak yakınsama oldukça …

13 pde  hyperbolic-pde  finite-element 

İnsanların genellikle tutarlı kütle matrislerini topaklı diyagonal matrislerle değiştirdiğini biliyorum. Geçmişte, yük vektörünün FEM tutarlı bir şekilde değil, topaklı bir şekilde birleştirildiği bir kod da uyguladım. Ama neden ilk etapta bunu yapmamıza izin verildiğine hiç bakmadım. Topaklamanın ardındaki kütle ve yük vektörlerine uygulanmasına izin veren sezgi nedir? Bunun matematiksel gerekçesi …

13 finite-element 

Bu, SE'nin Yazılım Önerileri tarafına daha uygun bir soru olabilir, ancak SE'nin bu kısmını sık yapan kişilerin bu soruya cevap verebilme ihtimalinin daha yüksek olduğuna inanıyorum. Comsol Multiphysics'e ücretsiz (sadece özgürlükte değil) alternatif arıyorum . İşte biraz zor: Ben sadece yükleri olan bir modelleme ve simülasyon paketi arıyorum, daha ziyade …

12 finite-element  comsol  octave  openmodelica  modeling 

Sistemim lagrange çarpanları (örn. Sıkıştırılamaz Stokes akışı) ile simetrik bir FE problemidir: ( ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} Burada tipik bir durumdur (denklemlerin Lagrange çarpanlarının en son görünmesini sağlayacak şekilde numaralandığından bile emin oldum). Sistem oldukça büyük (+ 100 bin satır).C= 0C=0C = 0 Bu sorunun cevabını okuduktan …

12 finite-element  petsc  krylov-method  preconditioning 

FEM-diskretize ve bir tepkime ile difüzyon problemi, örneğin çözerken ile 0 < ε « 1 (tekil bozulması), ayrık sorunun çözümü genellikle yakın sınıra salınım tabakaları sergileyecektir. İle Q = ( 0 , 1 ) , ε = 10 - 5 ve sonlu elemanlar doğrusal, çözelti u h oluşacağı- s Æ …

12 finite-element  discretization  numerical-analysis  singular-perturbation 

Ben 3D bir 4-düğüm eleman üzerine bir polinom ifade entegre etmek istiyorum. FEA üzerine birkaç kitap, entegrasyonun keyfi düz 4-elemansız bir eleman üzerinden gerçekleştirildiği durumu kapsamaktadır. Bu durumda olağan prosedür, Jacobi matrisini bulmak ve entegrasyon temelini, daha basit entegrasyon limitlerine [-1; 1] sahip olduğum ve Gauss-Legendre kareleme tekniğinin kolayca kullanıldığı …

12 finite-element  quadrature 

Yüksek Reynolds sayısı akışları çok ince sınır tabakaları üretir. Duvar çözünürlüğü büyük Eddy Simülasyon kullanılırsa, en-boy oranı mertebesinde olabilir . Birçok rejim bu rejimde kararsız hale gelir, çünkü inf-sup sabiti en / boy oranının kare kökü veya daha kötüsü olarak bozulur. İnf-sup sabiti önemlidir, çünkü lineer sistemin durum numarasını ve …

12 pde  finite-element  fluid-dynamics  stability  incompressible 

Sonlu elemanlar yöntemi ile ilgili literatürleri okurken, "asılı düğümler" terimine sıklıkla rastlanmaktadır. Biri bana gerçekten asılı bir düğümün ne olduğunu söyleyebilir mi?

11 finite-element  mesh-generation