Podcast #128: We chat with Kent C Dodds about why he loves React and discuss what life was like in the dark days before Git. Listen now.

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
5
votes
1answer
854 views

Meshing options to generate number of the sides of and element (tetgen-triangle)

I wrote a finite element code in fortran 90. This code is really fast, except the meshing process. I used triangle and tetgen for meshing in 2D and 3D, respectively, so this process is fast, of ...
31
votes
9answers
4k views

Modern resources for learning FEM

I need to get started using Finite Element Methods. I am about to start reading Numerical solutions of partial differential equations by the finite element method by Claes Johnson, but it's dated 1987....
9
votes
5answers
941 views

Motivation behind Galerkin method

I have a question about Galerkin method. I don't understand why the Galerkin method weights the residual by the shape functions and sets it equal to zero. I want to know what is reason of this. Why we ...
20
votes
9answers
35k views

Basic explanation of shape function

I just started studying FEM in a more structured basis compared to what I used to do during my undergraduate courses. I am doing this because, despite the fact that I can use the "FEM" in commercial (...
8
votes
2answers
4k views

How to properly apply non-homogeneous Dirichlet boundary conditions with FEM?

In general, Dirichlet boundary conditions won't be satisfied exactly for FEM for non-homogeneous boundary conditions. The FEM codes I've seen set the degrees of freedom to interpolate the Dirichlet ...
46
votes
5answers
8k views

What are criteria to choose between finite-differences and finite-elements

I am used to thinking of finite-differences as a special case of finite-elements, on a very constrained grid. So what are the conditions on how to choose between Finite Difference Method (FDM) and ...
24
votes
3answers
4k views

What is the purpose of using integration by parts in deriving a weak form for FEM discretization?

When going from the strong form of a PDE to the FEM form it seems one should always do this by first stating the variational form. To do this you multiply the strong form by an element in some (...
5
votes
3answers
771 views

Computing accurate fluxes with FEM

I have solved Poisson equation on a 3d domain with neumann and dirichlet boundary condition. I get the potential, take the gradient for each element and integrate on a surface of an element, I do this ...
4
votes
2answers
2k views

The meaning of conservative discretization in Galerkin FEM and Discontinuous Galerkin

I do understand the meanning of "conservative discretization" within the FVM/FDM framework, indeed it is well explained in this post. Now, according to the table in this slide (pp.8), it concludes: ...
8
votes
1answer
2k views

PDE discretization with the method of rothe and the method of lines (Modular implementation)

The Heat equation is discretized in space with FV (or FEM), and a semi-discrete equation is obtained (system of ODEs). This approach, known as the method of lines, allows to easily switch from one ...
16
votes
1answer
6k views

What is the general idea of Nitsche's method in numerical analysis?

I know that the Nitsche's method is a very attractive methods since it allows to take into account Dirichlet type boundary conditions or contact with friction boundary conditions in a weak way without ...
15
votes
1answer
5k views

How to derive the Weak Formulation of a Partial Differential Equation for Finite Element Method?

I have taken a basic introduction to Finite Element Method, which did not emphasize a sophisticated understanding of a 'weak formulation'. I understand that with the galerkin method, we multiply both ...
5
votes
3answers
1k views

Books on mathematical foundation of finite element methods

After reading three books about finite element method, with two of them covering also finite volume and grid generation, I found myself lost when I have to discuss these topics with library developers ...
10
votes
1answer
14k views

How to formulate lumped mass matrix in FEM

When solving time dependent PDE's using the finite element method, for example say the heat equation, if we use explicit time stepping then we have to solve a linear system because of the mass matrix. ...
13
votes
2answers
895 views

Impose the compatibility conditions for mixed finite elements method in Stokes equation

$\newcommand{\v}[1]{\boldsymbol{#1}}$ Suppose we have following Stokes flow model equation: $$ \tag{1} \left\{ \begin{aligned} -\mathrm{div}(\nu \nabla \v{u}) + \nabla p &= \v{f} \\ \mathrm{div} \...
5
votes
2answers
2k views

Does the finite element method impose any restrictions on the Peclet number for numerical stability?

Background on finite volume method When discretising the flux with a central difference stencil of the the advection-diffusion equation restriction $\frac{ah}{d} < 2$ must be observed for the ...
16
votes
3answers
4k views

Finding which triangles points are in

Suppose I have a 2D mesh consisting of nonoverlapping triangles $\{T_k\}_{k=1}^N$, and a set of points $\{p_i\}_{i=1}^M \subset \cup_{k=1}^N T_K$. What is the best way to determine which triangle each ...
14
votes
2answers
2k views

FeniCS: Visualizing high order elements

I've just started messing around with FEniCS. I am solving Poisson with 3rd order elements and would like to visualize the results. However, when I use plot(u), the visualization is just a linear ...
9
votes
2answers
497 views

Which novel data structures are used in adaptive FEM?

A lot of adaptive FEM libraries use more advanced mesh data structures to handle adding/removing nodes, edges, triangles, tetrahedra, etc. For example, the p4est library uses octree data structures ...
7
votes
0answers
500 views

Numerical implementation of the Dirichlet-to-Neumann map

I am solving the Dirichlet problem $$ \begin{cases} \Delta u = 0, \\ u|_{\partial D} = f, \end{cases} $$ in a $2d$ domain $D$ using the finite element method. What I want to get is the ...
4
votes
1answer
2k views

Comparison of Lattice Boltzmann Method vs Traditional Navier-Stokes based Methods

I have a choice of two options, analysing and implementing Lattice Boltzmann methods or traditional Navier Stokes based methods. I'm a CFD newbie and I have a rough idea (though not rigorous enough to ...
2
votes
1answer
153 views

Finite element method for odd order DE

What are theoretical hurdles in applying Galerikin method on, say, first order time dependent ODE? Is there no way we can form an inner product??
-1
votes
1answer
542 views

Pde problem with robin boundary condition

I have my pde 2D problem with robin condition (form: du/dn +ku=g) to solve with matlab. i have the exact function u and I want to find the function g in robin condition. How can i do it? thanks for ...
6
votes
0answers
419 views

How to do FEM in sector elements?

Suppose we have an element in the polar coordinate as $0\leq r\leq 1$, $0\leq\theta\leq \alpha$. What are the correct basis functions in this element? For this kind of mesh, there are a lot of ...
4
votes
1answer
2k views

Is there a mesh generator that will generate zero thickness elements for interfaces?

I have a geometry and a finite element method where there are a couple internal interfaces across which I want to enforce some fairly simple jump conditions. One recommendation I've gotten has been ...
3
votes
1answer
915 views

Implementation of gradient zero boundary conditon in advection-diffusion equation

My question is about Finite Element Method. I want to know how to implement "gradient zero" conditions to advection-diffusion equations in conservative form like, $\frac{\partial \rho}{\partial t} + ...
2
votes
1answer
2k views

3D Solid 8 Node FEM Matlab Code

So this semester, I'm taking a Finite Element Method course at my graduate school. We started out making codes for 1D bars and came all the way to 8 node solid elements. However, I seem to have run ...
1
vote
1answer
317 views

FEM on tet10 element: negetive determinant at the Gauss point

I am trying to implement a fem code on tet10 elements. I follow the lecture notes for tet10 implementation given in http://www.colorado.edu/engineering/CAS/courses.d/AFEM.d/AFEM.Ch10.d/AFEM.Ch10.pdf ...
5
votes
2answers
474 views

What does symmetrize mean? (imposing multifreedom constraints to stiffness matrix)

I have a small FEM implementation program. And I want to add imposing multifreedom constraints (MFC) feature to it. The theory of master-slave method is given here (page 10 for general case). ...
4
votes
1answer
259 views

Interpolating a mathematical function using a Hermite Cubic Finite Element Space

I have a Hermite Cubic Finite Element Space on a computer in the form of Matlab m-files. More specifically, I can evaluate four "shape functions" $N_1, N_2, N_3,$ and $N_4$, for which the following ...
4
votes
3answers
206 views

Is there a minimum angle requirement for cells in the finite volume method?

In his talk "What is a good linear finite element?", Shewchuk states that small dihedral angles in linear tetrahedra elements cause ill-conditioning of the stiffness matrix. Do small dihedral angles ...
3
votes
1answer
299 views

Finite Elements Weak Formulation generalization

I am struggling with an equation that represents the Weak form of Galerkin method: $ \phi^{T}F(\textbf{u})\sim \int_{\Omega}^{ } \phi.f_{0}(\mathit{u},\nabla \mathit{u}) + \nabla\phi:f_{1}(\mathit{u},...
48
votes
3answers
32k views

What are the conceptual differences between the finite element and finite volume method?

There is an obvious difference between finite difference and the finite volume method (moving from point definition of the equations to integral averages over cells). But I find FEM and FVM to be very ...
29
votes
9answers
1k views

What is a good way to run parameter studies in C++

The problem I'm currently working on a Finite Element Navier Stokes simulation and I would like to investigate the effects of a variety of parameters. Some parameters are specified in an input file ...
13
votes
2answers
4k views

What is the purpose of the test function in Finite Element Analysis?

In the wave equation: $$c^2 \nabla \cdot \nabla u(x,t) - \frac{\partial^2 u(x,t)}{\partial t^2} = f(x,t)$$ Why do we first multiply by a test function $v(x,t)$ before integrating?
21
votes
4answers
2k views

How to incorporate the boundary conditions with the Galerkin method?

I've been reading some resources on the web about Galerkin methods to solve PDEs, but I'm not clear about something. The following is my own account of what I have understood. Consider the following ...
10
votes
1answer
2k views

Raviart-Thomas elements on reference square

I'd like to learn how the Raviart-Thomas (RT) element works. To that end I'd like to analytically describe how the basis functions look on the reference square. The goal here is not to implement it ...
9
votes
1answer
373 views

$L^2$-convergence of finite element method when right hand side is only in $H^{-1}$ (Poisson eqn)

I know that the piecewise linear finite element approximation $u_h$ of $$ \Delta u(x)=f(x)\quad\text{in }U\\ u(x)=0\quad\text{on }\partial U $$ satisfies $$ \|u-u_h\|_{H^1_0(U)}\leq Ch\|f\|_{L^2(U)} $...
12
votes
2answers
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 ...
12
votes
2answers
343 views

Oscillations in singularly perturbed reaction-diffusion problems with finite elements

When FEM-discretizing and solving a reaction-diffusion problem, e.g., $$ - \varepsilon \Delta u + u = 1 \text{ on } \Omega\\ u = 0 \text{ on } \partial\Omega $$ with $0 < \varepsilon \ll 1$ (...
9
votes
2answers
4k views

How to efficiently implement Dirichlet boundary conditions in global sparse finite element stiffnes matrices

I am wondering how Dirichlet boundary conditions in global sparse finite element matrices are actually implemented efficiently. For example lets say that our global finite element matrix was: $$K = \...
7
votes
2answers
292 views

Stabilization of convection-dominated flow and turbulence modeling

Are stabilization techniques for convection-dominated flows like SUPG+PSPG, interior penalty methods, etc. able to handle turbulent flows without tubulence model being employed, at least up to some ...
6
votes
1answer
1k views

Find the direction of the gradient on a finite element mesh

Suppose we have a triangular mesh of a two dimensional shape $\Omega$, and on this mesh we define a P1 finite element structure. I know that given $u,v$ by their values at the vertices of the ...
6
votes
1answer
1k views

Appropriate space for weak solutions to an elliptical pde with mixed inhomogeneous boundary conditions

I'm working with the following mixed inhomogeneous boundary value problem: $\nabla(\kappa\nabla u)=f$ in $\Omega$ with $\partial\Omega = \Omega_1 \bigcup\Omega_2$ such that $u=g$ on $\partial\...
12
votes
1answer
401 views

What spatial discretizations work for incompressible flow with anisotropic boundary meshes?

High Reynolds number flows produce very thin boundary layers. If wall resolution is used in Large Eddy Simulation, the aspect ratio may be on the order of $10^6$. Many methods become unstable in this ...
10
votes
3answers
8k views

What is the difference between implicit FEM and explicit FEM?

What is the difference between explicit FEM and implicit FEM exactly? According to the post here, it seems that the only difference is whether implicit or explicit time integration is used. As I ...
9
votes
5answers
1k views

Is discontinuous Galerkin really any more parallelizable than continuous Galerkin?

I've always heard that easy parallelization was one of the advantages of DG methods, but I don't really see why any of those reasons don't also apply to continuous Galerkin.
8
votes
1answer
412 views

Finite element convergence rates for mixed problems

I've coded up a Stokes Flow problem using finite elements and am in the process of verifying that it works. I'm just not sure what convergence rate I should be expecting as I globally refine the mesh. ...
6
votes
3answers
432 views

Mixed Finite Element Method for the Stokes System—Some Implementation Details

I am currently working on my bachelor’s diploma. The research concerns mixed finite element method for the 2D Stokes system $$ - \Delta \boldsymbol u + \nabla p = \boldsymbol f, \quad \boldsymbol x \...
4
votes
2answers
370 views

Solving Poisson equation with current BC using FEM

I have a 3-dimensional domain D, with 3 types of BC which I am trying to solve the Poisson equation on. $$-\nabla \cdot (\sigma(x,y,z)\nabla u)=0$$ Insulator (Neumann BC) Electrode set at some ...