Questions tagged [finite-element]

A means of solving ordinary and partial differential equations. The domain of the problem is broken up into elements, and the solution in each element is expanded in a basis of functions. The Finite Element Method lends itself well to adaptive refinement, irregular geometry, and good error estimates.

Filter by
Sorted by
Tagged with
3 votes
0 answers
93 views

How is the surface Jacobian determinant calculated in FEM?

I am currently trying to evaluate surface forces on a structure. I came across P356 in Bathe's Finite Element Procedures 2014 (example 5.8) in which he related the edge derivative from the global to ...
0 votes
0 answers
55 views

Derivation of the second Piola-Kirchhoff tensor

I tried many formulas to find the components of the second Piola Kirchhoff. I need help to derive equations 27-29 on reference 1. [![enter image description here][1]][1] we have $$ \mathbf{C}=C_{11} \...
9 votes
3 answers
15k views

The real myth of GPU (specifically CUDA) really speed up FEM/CFD

Now I have been believing that FEM/CFD is supposed to be faster on a GPU unit - here I am using CUDA as solid example. However, I have not been able to find a convincing paper where the benchmark ...
1 vote
1 answer
66 views

A priori estimates in finite elements for inhomogeneous heat equation

Consider the problem $$\partial_t u-\Delta u = f\\ u(\Sigma_1)=f_D\\ \partial_\nu u (\Sigma_2)=f_N\\u(0)=u_0$$ where the sides of the space-time cylinder $\Sigma_i$ are disjoint (one of them could be ...
0 votes
1 answer
42 views

How can I define an equipotential surface/volume in FEniCS?

I want to solve electrostatic problem for potential. Charge density and medium permittivity are known, so is the potential of a grounded surface. I know how I can implement that. But I would like to ...
1 vote
0 answers
80 views

Numerical methods for $u_t+c u_x= \frac{-c}{x}u$?

I am looking for possible numerical methods to solve the PDE $$u_t+c u_x= \frac{-c}{x}u$$ I am particularly interested in a Finite elements method, although I am also curious if you can expose some ...
0 votes
0 answers
42 views

A Finite Element Method for a first order PDE?

I want to develop a finite element method to solve for $u(x,t)$ the PDE: $$u_t+c u_x= \frac{-c}{x}u$$ where $c$ is a constant. so I am trying the following ( as Rothe's method? ) : Letting $k= t_n- ...
0 votes
1 answer
407 views

Global to local transformation matrix in 2D frame structures

In section 3.2 of this paper [1], where 2D planar frame structures are being analyzed, the authors mentioned a transformation matrix to be used in extracting the element displacement vector from the ...
2 votes
1 answer
124 views

Determining the voxels between two boundary surfaces

Issue description I am working on human brain tACS simulations where I have the models of the skin, skull, csf, brain and ventricles in STL format. The shape does not matter and there are no ...
2 votes
1 answer
131 views

Distributed Lagrange multiplier approach to impose constraint in Poisson equation

I'm trying to understand how Lagrange multipliers are applied in order to impose constraints in PDEs. Consider $B \subset \Omega$. For instance, a square inside another square domain $\Omega$. Let's ...
1 vote
1 answer
47 views

Applying Stress Boundary Conditions in Commercial Finite Element Analysis Codes

I am trying to replicate a finite element analysis given in a research paper titled On the Detection of Stress Singularities in Finite Element Analysis 1 by G.B.Sinclair et. al. The geometry of the ...
1 vote
1 answer
91 views

What FEM solver should be used for matrix-valued FE spaces?

I am pretty new to using FE solvers. I am trying to solve a system of (up to) 9 complex equations. We write these as a matrix equation (here), (with the implied sum over $j$, for each component ...
1 vote
0 answers
32 views

positivity preservation for mixed finite element discretization

I'm interested in mixed discretizations of the diffusion (and related) equations: $$\begin{align} \frac{\partial h}{\partial t} + \nabla\cdot \mathbf q & = f \\ k^{-1}\mathbf q + \nabla h & = ...
1 vote
0 answers
57 views

Finding the weak form of a PDE with a tensor argument

I am trying to solve for the order parameter ($A$) in the Ginzburg Landau equations. I had asked on the math SE site but was recommended to ask here. We are trying to solve the following equation, (...
1 vote
1 answer
70 views

How to define a 3D surface from a set of points

Sorry in advance if this question has already been asked, I found nothing to help. I want to buid a 3D box whose top surface is topography. This topography is defined by a DEM i.e. a set of points (x,...
3 votes
2 answers
201 views

How to enforce fluid and solid dynamic coupling in fluid-structure interactions using the finite element method?

I apologize in advance if the question has been posted before or if it sounds a bit naive. I am writing my own code in MATLAB for a staggered finite element solver for fluid-structure interaction ...
4 votes
1 answer
81 views

A priori FEM estimates without $H^2$ regularity

In basic lectures on finite elements theory, people always assume $H^2$ regularity of the solution in order to derive $O(h^k)$ a priori estimates in the norms $H^{2-k}$, $k=1,2$. For simplicity let's ...
1 vote
1 answer
154 views

Question about step in the proof of standard discrete trace inequality

I'm studying from Guermond lecture notes available at https://www.math.tamu.edu/~guermond/M661_FALL_2019/chap12.pdf (see Lemma 12.8( Discrete trace inequality).) Consider the simple case $p=r$, i.e. ...
6 votes
2 answers
181 views

Is their a book/paper about mixed finite element method for engineering students (non-math)?

There are a lot of books about FEM, which are really friendly to engineering students. Through these books, we can know how to use shape/test functions based on the variational principle. But I'd like ...
4 votes
1 answer
94 views

Why "Right" and "Left" Cauchy-Green tensor?

$C=F^TF$ is called the "Right" Cauchy-Green tensor, and $b=FF^T$ is called the "Left" Cauchy-Green tensor. I suppose in $C=F^TF$ the non-transposed $F$ stands on the right, and in $...
1 vote
2 answers
88 views

Taylor-Hood elements for Darcy's equation

I would like to know if Taylor-Hood elements $P_2$-$P_1$ form a stable pair for the mixed approximation of Darcy's equation ( or Poisson's equation) with Dirichlet B.C. In the literature I only find ...
5 votes
1 answer
1k views

Which finite-element framework allows the easy implementation of HDG methods?

I have tinkered around with the implementation of hybridizable Discontinuous Galerkin (HDG) methods in DUNE-PDELab and DUNE-FEM for my PhD but since then I have not had the time to work on these codes ...
3 votes
1 answer
139 views

Adaptive Lagrangian-Eulerian methods and practical benchmark results

Does anyone know of any published study that talks about the practical aspects of running Adaptive Lagrangian-Eulerian techniques for solid and/or fluid mechanics problems? I'm looking for things ...
2 votes
1 answer
197 views

SIPG method for $-\nabla \cdot (\nu \nabla u)=f$

Consider the diffusion equation with a coefficient $\nu$: $$-\nabla \cdot (\nu \nabla u)=f$$ with Dirichlet boundary conditions $u = g_D$ in $\partial \Omega$. If the coefficient would be constant, ...
4 votes
3 answers
440 views

What condition ensures the global continuity of the solution in the FEM?

I know this is trivial but I don't seem to understand it. In which step of the FE formulation do we enforce the global continuity of the solution? Or in other words, how the construction of the local ...
2 votes
0 answers
93 views

FEM applied to heat equation and incompatible conditions

Consider the problem $$u_t - \Delta u = f \text{ on } \Omega\times (0,T)\\u=0 \text{ on } \partial \Omega\times (0,T) \\ u(x,0)=g(x) \text{ on } \Omega$$ with $g$ NOT vanishing on the boundary. If I ...
0 votes
0 answers
52 views

Reference request for finite elements theory

Consider a domain $\Omega \subseteq \mathbb{R}^{2,3}$ which is non convex and with $C^2$ boundary. Could you recommend a good reference where it is explained how: without needing isoparametric ...
0 votes
0 answers
33 views

3D Euler Bernoulli Beam for Nonlinear FEA

Anybody has any experience in coding 3D beam elements? I am trying to write a C++ code for a 3D euler bernoulli beam. For 2D, I used Reddy for coding 2D for non linear FEA. How should I proceed with ...
0 votes
0 answers
111 views

How to implement large rotations in total lagrangian formulation (nonlinear FEM)?

I have developed an Octave script to solve the nonlinear Euler-Bernoulli beam equations with linearized von Karman-strains, i.e. higher-order terms are dropped. The simulation results agree with ...
2 votes
1 answer
256 views

Best preconditioner for mixed-poisson problem (RT0 elements)

For a very large mixed-poisson problem with lowest order Raviart-Thomas elements (RT0), I plan on using an iterative solver. However, this kind of problem is not positive-definite (saddle point ...
13 votes
2 answers
2k views

Which preconditioners (and solver) in PETSc for indefinite symmetric systems should I use?

My system is a symmetric FE problem with lagrange multipliers (e.g. incompressible Stokes' flow): \begin{pmatrix}A & B^T \\ B & C\end{pmatrix} where $C = 0$ is the typical case (I have even ...
0 votes
0 answers
74 views

How to solve spatially discretised PDEs (method of lines) in solve_ivp or ODEint?

I can discretise the spatial domain of a system of PDEs using the method of lines, converting the system of PDEs to a system of ODEs (with a time derivative only). These equations (for context they ...
2 votes
0 answers
131 views

Variational loss of hp-Variational Physics Informed Neural Networks for 2D-Poisson Equation in Tensorflow

I am trying to reproduce the results from the hp-VPINN paper (https://arxiv.org/pdf/2003.05385.pdf) on tensorflow (v1) for Poisson's equation, particularly the two-dimensional Poisson equation. In one ...
1 vote
1 answer
376 views

Lower bound for bilinear form in FEM

I'm searching for lower bounds of bilinear forms arising in FEM for elliptic second order PDEs with mixed boundaries. I did some research and found: $$\max_{v_{h}\in\mathcal{V}_h(\mathcal{\Omega})}a(...
3 votes
3 answers
2k views

Scalar vs. vector potential for magnetostatics

When trying to solve a magneto-static boundary-value problem (BVP) ($\nabla \times \mathbf{H} = \mathbf{J}$ and $\nabla \cdot \mathbf{B} = 0$), one can use either the magnetic vector potential $\...
12 votes
1 answer
518 views

DG local equation, how to interpret mean-averaged test function

In the paper http://www.sciencedirect.com/science/article/pii/S0045782509003521, an HDG element-local equation is described on page 584 equation (4), with one of the equations taking the following ...
1 vote
1 answer
406 views

Proving solution existence and uniqueness of the Helmholtz equation with Robin boundary conditions with complex coefficients

I am trying to solve the Helmholtz equation with Robin boundary conditions with complex coefficients and the weak formulation $$ \iint_\limits\Omega\nabla p_0(x,y)\nabla\left(\overline{v(x,y)}\right)...
1 vote
2 answers
78 views

Can I use Q0 finite elements when there are gradients involved?

Could I apply a Q0-discretization to, say, the Poisson equation $\Delta \phi = f$ (where by Q0 I mean piecewise constant, and thus non-continuous, elements)? Solving this with FEM, at least as I know ...
0 votes
1 answer
172 views

deal.ii - ParaView "warp by scalar" of my output is not continuous

During our finite element course, we've solved the linear elasticity problem in 2D on a square (GridGenerator::hyper_cube) with $Q_1$ bilinear finite elements in ...
0 votes
1 answer
47 views

How can I correctly determine velocity of a point inside a grid after using mixed finite element method to solve Poisson equation?

I am using the mixed finite element method (MFEM) to solve the Poisson equation: $$\Delta h = 0,$$where $h$ denotes hydraulic pressure. The MFEM could determine the normal flux rate, $q_n$, through ...
1 vote
1 answer
189 views

Confusion about bilinear form for elasticity equation in deal.ii tutorial

I'm learning how to solve vector-valued problems with deal.II library. In particular, I'm looking at the following introduction from the official website https://www.dealii.org/current/doxygen/deal.II/...
0 votes
2 answers
70 views

calculating the Laplacian of the field variable in estimating the local residual error in the finite element method

to perform adaptive refinement in the finite element method according to the explicit residual method, the quantity $$\eta_K^2=h_K^2\left\lVert r\right\rVert_{L_2(K)}^2+h_K\left\lVert R\right\rVert_{...
0 votes
1 answer
89 views

Applying flux limiters on vertices/faces instead of quadrature points

I have a couple of questions with regards to the usage of flux/slope limiters, which I was hoping to get some input. I apologise in advance if my questions are trivial, but flux limiters lie way ...
1 vote
1 answer
65 views

Algebraic Multigrid on fine mesh vs. just starting with coarse mesh?

So if I've understood it correctly, the Algebraic Multigrid Method (AMG) basically takes a fine mesh, coarsens it, solves the coarse mesh and projects the solution back on the fine mesh. Wouldn't it ...
3 votes
2 answers
799 views

Initial Value Problem using Finite Element

I am trying to implement a FEM solver for the following initial value problem \begin{align} \frac{\partial u}{\partial t} - \nabla^2 u &= f\quad \text{ in } \Omega\times (0,T)\\ u &= g\quad \...
14 votes
2 answers
12k views

In FEM, why is the stiffness matrix positive definite?

In FEM classes, it's usually taken for granted that the stiffness matrix is positive definite, but I just can't understand why. Could anyone give some explanation? For instance, we can consider the ...
0 votes
0 answers
45 views

Why is the solution in FEM bounded by zero?

Consider the following problem: $$ -\nabla^2 u = f,$$ Referencing to this post: FEM, we write the problem in variational form I'm assuming Dirichlet boundary conditions here): Find $u\in H^1_0(\Omega)$...
0 votes
1 answer
82 views

Assembly of the Isoparametric Quadratic load vector in Matlab [duplicate]

I work to solve PDE using FEM in the case P2 on Matlab. I try to correctly assemble load vector using quadratic Lagrange shape functions $$b_i =(f,\phi_i)=\sum_{q=1}^{nq}f(r_q,s_q)*\phi_{i}(r_q,s_q)*...
1 vote
4 answers
176 views

Thin-plate FEM simulation

I want to simulate a vibration of a thin plate (the Kirchhoff-Love model) on triangular meshes. Can you advise me on an introductory-level review of different elements for thin-plate FEM simulation? I ...
1 vote
3 answers
326 views

Morley element implementation reference

I am looking for a detailed reference on the implementation of the Morley element for FEM, specifically for the biharmonic equation. By detailed, I mean that it should discuss the problems associated ...

1
2 3 4 5
22