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
2answers
415 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 ...
12
votes
1answer
422 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 ...
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
459 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. ...
7
votes
3answers
485 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
1answer
3k views

Meaning of CFL condition on parabolic problems

I've been studying this FEM theory and for the parabolic problems, there's the analysis of stability of the $\theta$-method. I followed the analysis and they get this CFL (Courant-Friedrich-Lewy) ...
12
votes
1answer
339 views

How to integrate polynomial expression over 3D 4-node element?

I want to integrate a polynomial expression over a 4-node element in 3D. Several books on FEA cover the case where integrating is performed over an arbitrary flat 4-noned element. The usual procedure ...
9
votes
1answer
267 views

Resources on mesh generation for finite element methods

I know that this is not really apart of the rules as this is a recommendation question and these don't really have an answer per say. But, like this forum posting: https://stackoverflow.com/questions/...
9
votes
2answers
744 views

Space-time finite element discretization for time-dependent PDEs

In FEM literature, semi-variational methods are typically used in the solution of time-dependent PDEs. I have not seen a fully-variational approach i.e. where space and time are discretised by FEM, ...
6
votes
2answers
1k views

Projecting Finite Element solution onto new mesh

I am implementing a finite element solver in MATLAB and I have the following problem. Let's say I have a mesh $\mathcal{T}_1$ with triangular elements on a rectangular domain $\Omega\subset\mathbb{R}^...
5
votes
2answers
621 views

Absorbing boundary conditions for acoustics in Discontinuous Galerkin

Note: I'm trying to implement a Discontinuous Galerkin method, as kind of a way to learn about these things. As of now, I've taken the acoustic wave equation $c^2 \nabla \cdot \nabla u(x,t) - \frac{\...
5
votes
5answers
862 views

Introductions to hp-FEM

do you know good introductions into or surveys $hp$-adaptive finite elements? In particular I do not know how the heuristics for choosing spatial refinement or increased polynomial degree are ...
4
votes
0answers
83 views

Extrapolation after successive finite element refinement

I wish to compute a certain function $\lambda$ (in my case an eigenvalue of the Laplacian) to a certain accuracy, which I wish to guarantee, if possible. I started with a finite element method. I know ...
3
votes
1answer
133 views

Space-time Galerkin of Burgers changes the convection speed

tldr: Can space-time Galerkin schemes applied to convection-diffusion problems lead to effects on the convection velocity? For time $t\in (0,1)$ and the spatial variable $\xi \in (0,1)$, I am ...
3
votes
1answer
474 views

Jump condition for elliptic equation in standard finite element method

Consider the elliptic PDE $ -\mathrm{div}(k(x)\nabla u) = f(x) \in \Omega$ and $u = 0$ on $\partial \Omega$ If $k(x)$ is piecewise constant across an interface $\Gamma$ in the domain, then we ...
3
votes
2answers
821 views

Dirichlet boundary condition for sparse matrix - Solving Ax=b only for free nodes?

I am solving Biot equation with sparse matrix in MATLAB. I have no problem with global sparse matrices assembly, but when I assign Dirichlet boundary condition, it is so slow. From this topic, How to ...
1
vote
2answers
2k views

Finding rate of convergence by curve fitting in Matlab

I have some data: number of nodes $N$ and error in energy norm corresponing to it. I have seen in some references that the rate of convergence is reported by $$\| u-u_h\| _E=CN^{\alpha} $$ How can ...
0
votes
1answer
281 views

Multi-domain 3D Geometries for MATLAB PDE Toolbox

In principle the PDE Toolbox in MATLAB can handle multi-domain 3D geometries as noted here. This feature and the associated function geometryfromMesh were introduced in MATLAB R2018a. The associated ...
10
votes
2answers
1k views

Are 8 Gauss points required for second order hexahedral finite elements?

Is it possible to get second order accuracy for hexahedral finite elements with fewer than 8 Gauss points without introducing unphysical modes? A single central Gauss point introduces an unphysical ...
7
votes
1answer
483 views

Why are functional representations of systems important in numerical applications?

I tried asking a similar question in SE.Physics, and I got some information regarding the abstract side of this, but I figured I should post here to get more complete information about the numerical ...
6
votes
0answers
131 views

Do practice and theory differ substantially when implementing Neumann Boundary Conditions using a Mixed Method?

I have implemented a pretty straightforward finite element solver for the following Poisson equation. For the purposes of this question we can assume the source term and the Dirichlet data both ...
5
votes
1answer
206 views

Method of Manufactured Solutions for non-differentiable coefficients

The Method of manufactured is commonly used for verification of computational science codes. I want to use the method for verification of Navier-Cauchy (elasticity) equations with periodic and Bloch-...
5
votes
1answer
1k views

Mixed boundary conditions Finite Element Method

I have the following problem in Finite Element Method $$ -(\alpha u')' + \beta u' + \gamma u = f$$ with $ \Omega = (0, 1)$, $ u(0) = 0 $ and $ u'(1) = 3 $ to be able to write the weak formulation ...
4
votes
1answer
1k views

Inverse isoparametric mappings for quadrilateral finite elements

I have an isoparametric mapping $F_{E}: \hat{E} \to E$ where $\hat{E}$ is a reference quadrilateral (square) and $E$ is some quadrilateral in the domain $\Omega$. If $E$ is a parallelogram, then the ...
4
votes
1answer
323 views

$L^2$-error in FEM: how to compute integral over reference element?

I have the following problem. The domain is $(0,1)$ and we consider a uniform triangulation on $\hat{\Omega}$ with elements $K_i = [i/N,(i+1)/N]$ and $X_h^1$ the linear finite element space. I wrote ...
4
votes
2answers
227 views

Preconditioner for the GMRES method in the Uzawa algorithm

I'm trying to solve \begin{equation}\left\{ \begin{split} \frac{\partial u}{\partial t}+(u\cdot\nabla)u-\nu\Delta u+\frac1\rho\nabla p&=f\;\;\;\text{in }\Lambda\\ u&=0\;\;\;\text{on }\partial\...
3
votes
2answers
1k views

Efficently invert tiny matrix in Fortran

I have a piece of code in Fortran90 in which I have to solve both a non-linear (with the Newton-Raphson method, for which I have to invert the Jacobian matrix) and a linear system of equations. When I ...
3
votes
0answers
86 views

numerical instabilities in Fluid Dynamics, Finite Element Method

I'm looking for references to understand where the numerical instabilities come from in hydrodynamics in general, and notably when the Péclet number: $Pe>1$. I'm using the finite element method. ...
2
votes
1answer
237 views

1D analytical solution vs FEM solution for a bar under compression

I simulated the compression problem in ANSYS and compared to the analytical solution and found some discrepancies. The classical solution to the 1-D compression problem is: \begin{align} u(x) = Cx ...
2
votes
1answer
5k views

Calculate Jacobian of triangular element given coordinates of vertices and displacements?

I am trying to determine quality of my mesh elements using the Jacobian determinant as the measure. My algorithm takes vertices of nodes in triangular mesh and moves them around so as to form ...
1
vote
3answers
310 views

Introductory Resources on FEM [duplicate]

I've currently begun studying Finite Element Method (FEM) and I'm finding it a little difficult to find resources that break it down into something comprehensible. All the resources I've found are ...
1
vote
2answers
2k views

Help with basic Triangular Finite Elements, constructing stiffness/mass matrices

I am going through this example problem (picture below) and working it out myself. The problem is that I do not get the same mass matrix (M) as the example problem. I am able to get the same ...
0
votes
1answer
2k views

FEM Stiffness Matrix is always close to Singular or Badly Scaled

I am making a code for an 18-node (3x3x2) 3D element FEM. However, even though I am (pretty) sure that all the shape functions are correct and whatnot, whenever I try and invert the stiffness matrix ...
0
votes
1answer
224 views

How to simulate thermal expansion in a 2D plane using FEA?

I am trying to model 2D thermal expansion of a square area inside another square using FEATool. I have simulated plane strain by incorporating forces pointing along the $[1 \,\,\, -1]^T$ direction ...
9
votes
1answer
661 views

What numerical quadrature to choose to integrate a function with singularities?

For example, I would like to numerically compute the $L^2$-norm of $\displaystyle u = \frac{1}{(x^2+y^2+z^2)^{1/3}}$ in some domain that includes zero, I tried Gauss quadrature and it fails, it is ...
7
votes
1answer
527 views

Which preconditioning for large linear elasticity problem?

The problem I want to solve is the displacement formulation of the linear elasticity : $$ \nabla \cdot \sigma = 0 \quad \text{in} \quad \Omega \\ \sigma = \lambda ( \nabla \cdot u ) I + \mu (\nabla \...
6
votes
2answers
383 views

Implementing the pressure correction method using finite elements

Ok so I am nearing the completion of my finite element Navier-Stokes solver that uses the $\theta$-method for time stepping and the pressure correction method for the pressure. I am following the ...
5
votes
1answer
571 views

How is rigid bodies implemented in finite element codes

I am writing a finite element code for structural analysis, and I want to implement rigid bodies. How is this usually done? Say that I have a square mesh, with one half of the mesh being defined rigid ...
5
votes
1answer
1k views

FEM: which is the correct way to impose Dirichlet B.C

I know the Neumann B.C. is implicit in FEM language. However, I have seen at least two ways to impose Dirichlet B.C. e.g. for the following problem 1D, $$\nabla^2 u + \nabla u= 0, u_{left}= 1, u_{...
4
votes
1answer
1k views

Using Gmsh to create a mesh with zero thickness (quad) interface elements

I acknowledge the following post where a similar question is posed and a very nice answer has been provided: Is there a mesh generator that will generate zero thickness elements for interfaces? ...
3
votes
2answers
212 views

Does length unit in FEM affect numerical condition?

I try to solve a system of coupled PDEs using FEM. Unfortunately, the originating matrix has very poor condition. After days of double checking and thinking, I suspect the following reason: Given a ...
3
votes
3answers
225 views

$H^1$-convergence rate of finite element method for Poisson equation, depending on element order

I wanted to verify my FEM-program by applying the method of manufactured solutions, while solving the Poisson equation in two dimensions using the continuous Galerkin method $$-\nabla^2u=f$$ with $$u=...
3
votes
1answer
160 views

Numerical computation of the velocity in the steady Navier-Stokes equation

I've asked this question on Math.SE too. Let $d\in\left\{1,\ldots,4\right\}$ $\Lambda\subseteq\mathbb R^d$ be bounded, nonempty and open and $\partial\Lambda$ be Lipschitz $V:=\left\{u\in H_0^1(\...
3
votes
2answers
365 views

Finite Difference and Finite Volume as special cases of Finite Element

I have heard this many times that "FD is a special case of FE" and separately, that "FV is a special case of FE". However I have never come across a brief document that substantiates that claim, ...
2
votes
1answer
136 views

Boundary conditions in a finite element eigenvalue problem

I've been reading multiple papers and related posts for a while now, but I can't seem to find a specific answer to the issues I'm having so I hope someone can clarify things here. I'll provide some ...
2
votes
1answer
70 views

Computation of equivalent thermal resistance and thermal capacitance from FE model

I need to compute the equivalent thermal resistance and thermal capacitance of a structure used for heat transfer. For illustration purpose let’s say it’s the 2D problem of the following figure. In ...
2
votes
2answers
670 views

Schrödinger equation with time dependent Hamiltonian

I need to solve the Schrödinger equation with a time dependent Hamiltonian $$i\hbar \frac{\partial}{\partial t} \Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 +\frac{1}{2} k(t)(x^2+y^2) + V(r)\right]\Psi $...
2
votes
2answers
118 views

General Lagrange basis formula (usual problem in finite element context)

It is easy to prove that, $$\{p_1(x,y)=1-x-y\;,\;p_2(x,y)=x\;,\;p_3(x,y)=y\}$$ is a Lagrangian basis of $\mathbb{P}_1(\hat{T})$ (polynomials of total degree less that 1 living on $\hat T$), where $\...
2
votes
1answer
120 views

Prevent single node spikes in a FEM-simulation (using continuous Galerkin)

I am trying to solve a non-linear time-dependent heat equation $$\partial_tT=\nabla \left(k_T(T)\nabla T\right) + f$$ (similar to question Solving a non-linear heat equation with the galerkin method ...
2
votes
0answers
61 views

How can I numerically solve a saddle point problem with repeated constraints?

I am interested in numerically solving the following constrained minimization problem; Find the value of $x\in \mathbb{R}^n$ that minimizes $f$ where $f\colon \mathbb{R}^n\to \mathbb{R}$ is defined ...