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
8
votes
1answer
193 views

Should we always expect FEM error plots to be straight lines?

The error estimates in FEM are usually of the form $$||u^h-u||\leq Ch.$$ Taking logarithm on both sides, we obtain $$\log ||u^h-u||\leq \log C + \log h.$$ This estimate implies that the error lies ...
4
votes
1answer
104 views

Original paper on the augmented Lagrangian method in FEM

I am writing a paper in which I want to cite the earliest reference to the augmented Lagrangian method in FEM. For the pure Lagrangian method in FEM, the classical work of Babuška [1] is the original ...
3
votes
1answer
174 views

calculation of the right hand side of DG FEM (with code)

I got stuck with Hestaven/Warburton's dG-FEM Matlab code. Starting with the file AdvecRHS1D.m, we see in line 11 ...
5
votes
1answer
245 views

Stability of hyperbolic PDE and DG-FEM

In the book of Hesthaven and Warburton on discontinuous Galerlkin methods in example 2.3 (regarding solutions of the wave equation), the authors regard the following PDE: $$\frac{\partial u }{\...
4
votes
1answer
226 views

DG-FEM integration by parts

I am going through the book of Hesthaven and Warburton on discontinuous Galerkin methods. I have difficulties understanding some basic steps in the calculations. Consider the PDE: $$\frac{\partial u}...
1
vote
1answer
94 views

Abaqus, ANSYS, and FVM solver for thermal expansion problem converges to different values

Is it reasonable for a FEM and FVM code to converge to slightly different solutions for the same physical problem (identical BCs, geometry, properties, etc...), provided stability constraints are ...
0
votes
2answers
499 views

(FEM) 1D time-dependent heat equation convergence problem

I'm simulating a simple 3-node bar with convection BCs at the edges to validate my FEM code. The following data was used: Initial temperature = 25 ºC Temperature surrounding the rod = 10 ºC Thermal ...
0
votes
0answers
76 views

Algebraic multigrid for coupled equations

As far as I understand is algebraic multigrid(AMG) a method that was intentionally developed to solve linear systems where every grid point or node has a single DOF. When AMG should now be used for ...
2
votes
0answers
38 views

How to solve the Poisson equation with KINK aligned with mesh facet

I have a problem that solving the Poisson equation with kink ( discontinuous gradient but solution is continuous ) in the analytical solution, I want to solve this problem with FEM. To approximate ...
3
votes
1answer
115 views

References for the nonlinear reaction-diffusion equation using Finite Element Methods

I want to study how to solve the following PDE \begin{cases} -\nabla \cdot(\ k(x,y) \ \nabla u \ ) + \beta(x,y)\ u^2 = f(x,y), \ (x,y) \in \Omega \subset \mathbb{R^2} \\ \hspace{0.5cm} u = ...
1
vote
2answers
281 views

Order of element vs Degrees of freedom of the element

I have read that the order of the element is the order of the polynomial used to approximate/represent the field variable in that element. If we consider a one-dimensional, 2 degrees of freedom ...
3
votes
1answer
76 views

How to show the stability of $L^2$ projection?

If $\mathcal{T}_h$ is a regular and quasi-uniform triangulation of $\Omega$, and $V_h$ is the $H^1$-conforming linear finite element space. Moreover, let $P_h$ be the $L^2$ projection to $V_h\subset H^...
3
votes
2answers
69 views

Good references for dual-weighted residual (DWR) method for goal-oriented adaptive mesh refinement (AMR)

Can anyone help me with good references (books or papers) where I can learn about dual-weighted residual (DWR) method for goal-oriented adaptive mesh refinement (AMR)?
7
votes
2answers
540 views

Spectral Element vs Finite Element

I am trying to understand the difference between SEM and FEM. If I go by this paper, spectral element methods are a subset of FEM methods and the only difference lies in the choice of basis functions. ...
3
votes
2answers
382 views

Electromagnetism FEM (FEniCS) interpolation - leakage effect

As for the background of what is going on: I'm using FEniCS that is dedicated FEM solver The problem I'm solving is magnetostatic problem where the governing PDE is $$ \bf{\nabla} \times \frac{1}{\mu}...
1
vote
1answer
73 views

Split solution of FEM problem depending on number of DOF

Assume we have a 3D finite element structural problem discretized with hexahedral elements with 8 nodes and 3 degrees of freedom per node. Instead of solving the global stiffness matrix system for all ...
0
votes
0answers
80 views

Trouble with creating correct element matrices in Finite Element Analysis for a cantilever beam

I'm trying to solve for displacements of a cantilever beam numerically with FEA. The beam is modeled as a 3D-solid made up of a set of 8-noded hexahedral elements, which are in their undeformed state, ...
1
vote
1answer
148 views

Is there any fundamental difference between meshing for FEM, FVM and FDM?

I am a novice to the field of computational science and have just started studying the FDM and FEM (haven't started on FVM yet). While trying the understand the subject I got this question and trying ...
0
votes
0answers
36 views

Testing the SUPG method and other methods for hyperbolic equations

I am interesting in integrating the simple equation $$ \frac{\partial \phi}{\partial t} + \mathbf{u}\cdot\nabla \phi = 0 $$ with a Dirichlet boundary condition at the influx boundary ($\mathbf{u} \...
13
votes
3answers
2k views

Mathematically, why does mass matrix / load vector lumping work?

I know that people often replace consistent mass matrices with lumped diagonal matrices. In the past, I've also implemented a code where the load vector is assembled in a lumped fashion rather than ...
4
votes
1answer
573 views

What are all these functions in FEM? Shape function vs Basis Function vs Trial Function vs Test Function vs Interpolation Function

I am a novice in FEM. I have some experience with FDM which was pretty straight forward. Since I have a confusion with a number of concepts, I will try to break them down by writing down what I have ...
0
votes
0answers
52 views

Strange spectral (finite) element results of a solid plate

I tried to implement the spectral element method1 as proposed in [1]. I simulated an aluminum plate and a vertical concentrated force was applied at the middle of the top left quarter of the plate2. I ...
1
vote
0answers
62 views

Metis: how to use and tutorial recommendation

I am new to METIS and trying to use it in my fortran code. I read the manual online. But still, I am not sure about how to implement it my code. I tried the test cases in the graphs directory. For ...
0
votes
2answers
275 views

3d vs 2d finite element method

Is the theory of 3d finite element method just an assembly of 2d finite element analysis by putting planes on top of each other, or, a much more comple and different theory applies for 3d, with ...
2
votes
1answer
240 views

Why do we use hermite interpolation for finite element method in beams?

Why not just Lagrange polynomials basis functions
1
vote
1answer
320 views

(FEM) Nodes reordering for sparse matrix storing techniques

Is it necessary to reorder nodes (using Reverse Cuthill-Mckee algorithm, for example) if I am already using a CSR or CSC storing technique? Because since CSR/CSC stores only non-zero elements I guess ...
5
votes
0answers
82 views

Numerical methods for the continuity equation with Sobolev vector field

Consider the continuity equation $$ \partial_t \rho(x,t) + \operatorname{div} (b(x,t) \rho(x,t)) = 0, \qquad t \in [0,T], \quad x \in \mathbb R^N, $$ with $b \in L^1((0,T), W^{1,p}(\mathbb R^N))$. ...
1
vote
1answer
754 views

Efficiently plot a finite element mesh solution with Matplotlib

I am looking for the most efficient way to plot a mesh using Matplotlib given the following information, coordinates of each node, what nodes belong to each element, and the value each node has. Below ...
4
votes
1answer
273 views

Is a symmetric bilinear form necessary to ensure a weak formulation has a solution?

Problem I want to convert the general second order linear PDE problem \begin{align} \begin{cases} a(x,y)\frac{\partial^2 u}{\partial x^2}+b(x,y) \frac{\partial^2 u}{\partial y^2} +c(x,y)\frac{\...
5
votes
1answer
241 views

Advantages and disadvantages of space-time finite element methods

I have heard of space-time finite element methods. Although I was able to find some articles that describe the different possible methods from a mathematical point of view (thanks to Space-time finite ...
2
votes
0answers
73 views

Numerical solution to N-dimensional diffusion on simplex?

Assume I have a system of at least (but generally only) $N+1$ points in an $N$-dimensional space ($N > 3$ is possible). At each of these points $x_i, i=1,...,N+1$ I know an initial potential/...
4
votes
1answer
65 views

Solver for generalized eigenvalue problem with multipoint constraints

We have the following generalized eigenvalue (set of) problem(s) $$[K_R(\kappa)]\{u_R\} = \omega^2[M_R(\kappa)]\{u_R\}\quad \forall \kappa \in [\kappa_0, \kappa_1]$$ with \begin{align} &K_R(\...
0
votes
0answers
78 views

Constraining the total volume in Finite Element Methods

I have a diffusion problem which can be broken down to be: $-\Delta u = f(u) $ on $\Omega ~/~ \Omega_{int}$ $u = 1$ on $\Omega_{int}$ Note that this is an internal Dirichlet constraint to the ...
-1
votes
1answer
226 views

Finite Element - Flux Calculation

I am solving an advection-diffusion equation using the FEM and am having trouble calculating my fluxes. I start with the equation, $$\frac{\partial n}{\partial t} = \frac{\partial j_{n}}{\partial x}\...
3
votes
1answer
72 views

Parallelizing FEM for elliptical PDEs with n >1

For a little personal project, I am picking up my FEM skills again. I learned a lot about the theory back in university and I am able to implement a simple FEM solver for specific problems but I was ...
6
votes
3answers
356 views

Why is the FVM traditionally used in CFD, and FEM in computational structures?

Most CFD codes use FVM. Most computational structures codes use FEM. Why is the FEM not frequently used in CFD, and why is FVM not frequently used in FEM?
0
votes
0answers
75 views

Finite element method for Surface integrals using polar coordinates

I am trying to solve a 2D elliptic PDE (see complete electrode model for electrical impedance tomography) using the finite element method (FEM) over a circular region $\Omega$. I have discretized the ...
1
vote
1answer
70 views

Area and volume of P2 elements

Context: For a fluid solver, I need to compute areas and volumes of curved elements. My curved elements are quadratic triangles and quadratic tetrahedra, defined by their Lagrange nodes. Question: ...
3
votes
0answers
68 views

Structural Analysis Library

Can anyone recommend a structural analysis library that satisfies the following requirements: C++ API Simulate both beam elements and shell (slab) elements Both static and dynamic analysis Free and/...
1
vote
1answer
68 views

Mixed formulation in 1D

I have been working on a hybrid dimensional model using the mixed FEM formulation, in which 3D elements and 2D elements are combined by certain relationships between the degrees of freedom (DOFs) ...
6
votes
1answer
201 views

Dirichlet boundary conditions in generalized eigenvalue problem

Let us consider a problem of the form $$(\mathcal{L} + k^2) u(\mathbf{x})=0\, ,\quad \forall \mathbf{x} \in \Omega$$ with Dirichlet boundary conditions $$u(\mathbf{x}) = 0, \quad \forall \mathbf{x} ...
0
votes
1answer
105 views

Imposing zero mean condition in FEM

I wanted to solve a periodic elliptic equation of the form $$-\nabla\cdot(A\nabla u)=-\nabla\cdot F$$ on $Y=[0,2\pi)^d$ using FreeFem++, where $A$ and $F$ are $Y$-periodic. The space of solutions is $...
1
vote
0answers
78 views

Finite Element Model of Euler-Bernoulli Beam Theory with Isoparametric Element

In the formulation of Euler-Bernoulli Beam Theory, there are two degrees of freedom at a point, $w$ and $\frac{dw}{dx}$. Typically, the finite element model of this theory uses cubic polynomial for ...
1
vote
1answer
173 views

Normalization of polynomials for discontinuous Galerkin methods (DGM)

I was curious if someone could share their opinion on this matter. I have noticed that some people in literature normalize their Legendre polynomials, i.e. divide or multiply the polynomial by $$\...
7
votes
1answer
221 views

Size of jump for piecewise discontinuous approximations

If one has a sufficiently smooth function $u$ that is approximated by a piecewise constant function $u_h=\Pi^0_h u$ on a mesh of cell size $h$ (where $\Pi^0_h$ is the $L_2$ projection onto the ...
1
vote
1answer
126 views

Node renumbering in a 2D mesh

I have a 2D domain which is discretized using Q4 elements. I have the nodal positions and the element connectivity matrix. I would now like to renumber the nodes in such a way that all the interior ...
1
vote
1answer
110 views

Relation between conjugate gradient method and finite elements method

What is difference beetwen this two method? Are these methods far from each other or are these methods complement each other? Could you take an example?
3
votes
3answers
217 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=...
0
votes
1answer
46 views

Given co-ordinates of 8 vertices, how to calculate the outward normal and surface area for each face of a irregular hexahedron?

I am working on an FEA mesh of hexahedron elements. The elemental level calculations involve finding the surface normals and area for each surface of a hex element. I preferred the vector cross ...
3
votes
0answers
446 views

Gmsh for 3D volume with inclusions [closed]

In an attempt to create three-dimensional volumes with inclusions in Gmsh I stumble upon a problem which was non-existent in the two-dimensional case. I'm using the OpenCASCADE geometry kernel ...

1 2
3
4 5
18