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.
1,213
questions
0
votes
0
answers
29
views
RFEM: Equilibrium over plate section
I'm calculating a simple plate arrangement in RFEM:
The plate has dimensions of 1m X 1m, thickness of 50mm. Elastic modulus of the material is 205000.0 N/mm2, Poisson's ratio is 0.3.
The plate is ...
- 133
3
votes
1
answer
163
views
Projection (or fractional-step) methods Vs coupled method for incompressible Navier-Stokes
My question is in the context of the finite element method.
Incompressible Navier-Stokes equations can be solved using the coupled method or projection/fractional step methods. Each method has its ...
- 875
0
votes
0
answers
40
views
General dimensional solution to fast sweeping quadratic equation
I am reading Hongkai Zhao's paper The Fast Sweeping method. I have implemented the method in 2D and now I want to move onto 3D.
However section 2.6 confuses me.
The article says:
But I have no idea ...
- 263
3
votes
1
answer
290
views
Under what circumstances is parallel scaling of the finite element method not "solved"?
I see things like this deal.II example, (ctrl + f for "superMUC") which seems to show some pretty impressive scaling (nearly twice as fast for twice as many CPU cores for a wide range of ...
- 145
0
votes
1
answer
75
views
Cohesive zone model finite element implementation
I have a 2D vector $\boldsymbol{u}$ and it's norm is $\lambda$. I have this function:
$$\boldsymbol{T}=\dfrac{\boldsymbol{u}}{\lambda} e^\lambda$$
I need to compute $\boldsymbol{T}$ and it's Jacobean ...
- 394
1
vote
0
answers
63
views
Accessing nodal and degree of freedom organisation in the solution vector in Fenics and Firedrake
Context
Let suppose a solution to a finite element problem (associated to a mesh) stored in a vector $X$ of size $(N \cdot d)$ corresponding to $N$ nodes in the mesh and $d$ degrees of freedom for ...
0
votes
1
answer
47
views
RFEM: Tensile force in member in buckling analysis
Here is a simple T-shaped structure in RFEM. The structure is loaded from above, into the intersection of all members.
I performed a buckling analysis, and obtained the buckling modes that I'm ...
- 133
1
vote
1
answer
110
views
Evaluate the local mass matrix over a quadratic curvilinear tetrahedral element
I am evaluating the local mass matrix over a 3D curvilinear 10-node tetrahedral element. According to an index-naming rule, the shape functions are:
\begin{equation}
\begin{bmatrix}
2L^2_0 - L_0\\
...
- 11
2
votes
0
answers
102
views
Schur complement formulation of linear system
Consider a system of the following form:
$$(A+K)x=b$$
where $A$ is symmetric, positive definite and block diagonal (in fact, a block diagonal matrix made of stiffness matrices arising from FEM ...
- 157
4
votes
1
answer
145
views
Estimating forces on a model from the displacements of nodes
In any FEM problem involving mechanics, we try to solve the differential equation for the displacement field, $u$ given the force vector in the nodes, $F$. In industry, we often see our automobiles ...
- 43
0
votes
0
answers
43
views
Are my boundary conditions in my stiffness matrix correct?
I am trying to find a P1 Lagrange finite element solution to the following ODE:
$\begin{cases}-u''+u'+u=f~~~~~~~~\text{in}~~(0,1)\\ u(0)=1, u'(1)=0\end{cases}$
Where $f(x)=-2e^{x}+2\left(1-x\right)e^{...
1
vote
2
answers
194
views
C^1 continuous element for a triangle?
I am looking for an element for FEM that is piecewise $C^1$ continuous across triangles (i.e. $C^1$ continuous on the edge separating 2 triangles of the mesh).
I have heard about the Bell element:
...
- 263
0
votes
0
answers
81
views
It is possible to solve integro-differential equations using in Fenics?
I am interested in solve the following integro-differential equation:
\begin{align}
\frac{\partial{\rho(\theta, t)}}{\partial{t}} = D \frac{\partial{\rho(\theta, t)}}{\partial{\theta^2}} - \beta \...
- 111
0
votes
1
answer
141
views
FreeFEM++ converting equation into code
I am trying to solve the following problem nuemrically:
$$u_t = \Delta u + \sin t$$
To that effect I scanned the documentation of FreeFEM, the closest example I can find to my problem is the Thermal ...
- 263
1
vote
1
answer
65
views
Solving a boundary value problem with variable number of coupled equations
Let's assume the equation
$$
\nabla^2u_n(\vec{r})+a_n(\vec{r})u_n(\vec{r})=\sum_{m=1}^{N}b_{nm}u_m(\vec{r}),\quad n=1,2,\dots,N,\quad \vec{r}\in\varOmega\tag{1}\label{eq1},
$$
is to be solved for $u_n$...
- 194
1
vote
0
answers
84
views
2D FEM element with 6 DOF per node (including rotation)
I'm writing finite element code in 3D, and I wish to connect vertical (Euler-Bernoulli) beams to a horizontal plate. The deformations of a beam are well described, but for the plate it is more ...
- 11
1
vote
0
answers
111
views
Interface condition for 1D Helmholtz equation using finite element method
I want to implement a 1D Helmholtz equation with jump condition. The domain is $x=[0,1]$ and both ends have Dirichlet boundaries($p$=0). The 1D strong formulation is;
$$c^2\nabla^2p + w^2p=0 \qquad \...
- 11
0
votes
0
answers
100
views
Pyton get fem solving simple 2D differential equation?
I need to solve the differential equation:
$$u_t = \Delta u + \sin t$$
On a 2D domain with homogeneous Dirichlet boundary conditions. Tod o this I am trying to use the python package getfem.
I am ...
- 263
3
votes
0
answers
86
views
Is there any book about fundamental FEM theory using similar terminology as Comsol MultiPhysics?
I am considering to solve complexed PDE systems, like in this post, using Comsol MultiPhysics. The PDEs are different from the General Form provided by Comsol. The Weak Form module may be worth a ...
- 31
3
votes
1
answer
160
views
How to compute numerically the $H^{1/2}$ norm of a function
I'm, in the context of FEM. Let's say I have a discrete function $g$ living on the boundary of my domain $D$. I need to compute numerically $||g||_{1/2,\partial D}$.
The definition I know is the ...
- 281
1
vote
1
answer
119
views
Explicit polynomial for quadratic elements? (FEM)
In this resource the linear (barycentric) elements are explicitly given:
The geometry placement of higher order elements is also given but not expression for the polynomial of $P_2$ is given. I am ...
- 263
1
vote
1
answer
130
views
Could you recommend some books on FEM that explain various data-structures in FEM?
I want to understand the data structure of elements, elements around elements, and so on, and various other data structures in FEM, could you please recommend some books?
- 151
1
vote
2
answers
115
views
Elements on a triangle (FEM)
I am trying to learn about 2D FEM methods. I am trying to understand the generalization of Lagrange polynomial basis from 1D into 2 variable polynomials over triangle domains.
The most basic element ...
- 263
3
votes
0
answers
52
views
Is there an obvious n-dimensional elasticity complex?
Finite element exterior calculus helps to furnish nice stable elements for many common elliptic PDE by looking at the de Rham complex
$$\Lambda^0(\Omega) \stackrel{d_1}{\to} \Lambda^1(\Omega) \...
- 10.1k
6
votes
2
answers
292
views
Continuous vs discontinuous space-time FEM
What are some reasons for approximating a (e.g. parabolic) PDE using the space-time method, with continuous finite elements in time, vs discontinuous finite elements in time?
Are there e.g. ...
- 157
5
votes
0
answers
90
views
How to get an "optimal point" for refinement in FEM adaptive mesh refinement?
Consider the following 1D problem
\begin{align*}
\begin{cases}
\displaystyle
-\frac{d^2u}{dx^2} = f(x), \hspace{0.5cm} x\in (a,b) \\[4mm]
u(a) = u_{a}, \ \ u(b) = u_{b}
\end{cases}
\end{align*}
I ...
- 51
1
vote
1
answer
185
views
Discrete Laplacian operator with finite element discretization
Let a function $u \in H^1_0(\Omega)$ defined by its values at the mesh nodes. Can we compute its Laplacian using the matrix resulting from the finite elements discretization of Laplace's equation? I ...
- 25
0
votes
0
answers
48
views
Translating poisson equation to sfepy code
As the title suggests, I am trying to convert the Poisson equation:
Into sfepy code that can solve the differential euqation problem numerically.
I am trying to modify the tutorial on lienar ...
- 263
2
votes
1
answer
115
views
Cell-based vs face-based finite element methods
Notation:
Denote $T_{h} = \left\{K\right\}$ to be a face-conforming triangulation of a domain $\Omega$ such that $K_{i} \cap K_{j} = \emptyset$ for $i \neq j.$ Additionally, denote $\mathcal{V}_{h} = ...
- 21
2
votes
0
answers
102
views
Conceptual doubt regarding 2D conjugate heat transfer modelling (COMSOL and Mathemtica)
I have been dealing with some conceptual flaws in my understanding of modelling, which I will elaborate herein. I am modelling conjugate heat transfer of a reciprocating fluid, which flows with ...
- 41
0
votes
0
answers
39
views
books/paper recommendation on computational thermal-turbulence by using FEM
I have just learned basic FEM for 2D N-S euqation, now my teacher let me to do the following problem, the document of this problem is in large fluid problem, the system of equations is listed in that ...
- 151
3
votes
1
answer
247
views
Is it really necessary to solve a system of linear equations in the Finite Element Method?
When we solve some boundary value problem by Finite Element Method, the appropriate system of linear equations is built, $$Ax=b.$$
Usually we use the solution x just for plugging it into some ...
0
votes
0
answers
57
views
solve a coupled PDE system with some discontinuity by a mixed FEM
$$
\begin{cases}
u_t=f(u,v,\nabla{v}),\\
\Delta{[g(u)v]}=0.
\end{cases}$$
I want to solve the above PDE system by a mixed FEM, that is, $u_t=f(u,v,\nabla{v})$ by discontinuous Galerkin (DG), where $f$...
- 225
0
votes
1
answer
105
views
How can I solve this PDE system by discontinuous Galerkin method?
As is known to all, the discontinuous Galerkin method (DG) was first used to solve the equation $u_t+u_x=0$. Now I have the following system of PDEs:
$$
\begin{cases}
u_t=f(u,v,\nabla{v}),\\
\Delta{v}...
- 225
7
votes
0
answers
251
views
Matrix-free FEM references
I've seen that many people are using matrix-free fem codes in my community (mechanical engineering). I have to admit that I googled a bit and I didn't manage to find a good reference for the subject. ...
- 281
1
vote
0
answers
160
views
Partition of unity in FEM with bubble functions
In FEM with bubble functions, the field ($\boldsymbol{u}$) is approximated as a linear combination of the standard one ($\tilde{\boldsymbol{u}}$) plus the bubble field ($\boldsymbol{u}^b$). That is,
\...
- 875
2
votes
2
answers
226
views
Errors imposing boundary conditions weakly with DG
I am using interior penalty discontinuous Galerkin to solve a simple Laplace problem:
\begin{align*}
\nabla u=0
\end{align*}
with prescribed 0 and 1 Dirichlet boundary conditions on opposite edges of ...
- 394
4
votes
2
answers
331
views
computing higher order derivatives with linear elements
Consider the following equation on $(0,1)$, with Dirichlet boundary conditions on both ends.
$$
\frac{d}{dx}\left(k(x)\frac{du}{dx}\right) = 0
$$
Let us solve this using simple linear finite elements. ...
- 668
0
votes
1
answer
109
views
$u_t+a u_x=0$ solved by Discontinuous Galerkin
If $u_t+a \cdot u_x=0$ under a periodic boundary condition (to mimic an infinite domain) is solved by Discontinuous Galerkin (DG), how to implement periodic boundary condition and the other details in ...
- 225
4
votes
1
answer
339
views
guiding principles of integration by parts in FEMs
Most FEM texts or tutorials apply FEMs on diffusion equations where the 2nd spatial derivative is integrated by parts during weak formulation. For convection diffusion equations, there is also a first ...
- 225
1
vote
0
answers
60
views
(Dynamic FE analysis) Proportionality between time step and displacement
I am currently adopting a Newmark algorithm, for which the system of equation is
$$(M+{\beta}{\Delta}t^2{K}){d_i}=\beta{\Delta}t^2{F_i}+Md_{i-1}+{\Delta}tM\dot{d}_{i-1}+{\Delta}t^2M\left(\frac{1}{2}-\...
- 11
1
vote
0
answers
43
views
FEM Basis for functions of two variables in $\mathbb{R}^2$ (Applied to Linear Full Stokes Equations)
I'm trying to do FEM for a very basic version of the linear full Stokes equations in two dimensions. Say we are working in the grid $[0,1]\times[0,1]$ in the $xy-$plane.
To solve an FEM problem for a ...
- 111
0
votes
1
answer
125
views
time discretization via FEM to solve ODEs
It seems that FEMs are never seen to implement time discretization. Is time discretization via FEMs to solve ODEs ever possible?
- 225
1
vote
1
answer
54
views
How to implement pressure stabilization in matlab while solving steady stokes equation
the equation is $$
\left\{\begin{array}{l}
-\nabla \cdot \mathbb{T}(\mathbf{u}, p)=\mathbf{f} \text { in } \Omega, \\
\nabla \cdot \mathbf{u}=0 \text { in } \Omega, \\
\mathbf{u}=\mathbf{g} \text { on ...
- 151
2
votes
1
answer
465
views
Meaning of Degree of Freedom in FEM
Assume we want to solve the Poisson eq. with the FEM on some Domain $\Omega$, i.e.
$$\begin{cases} -\Delta u = f, \; \Omega\\
u = 0, \; \partial \Omega \end{cases}$$
For the sake of the discussion let ...
- 123
1
vote
0
answers
133
views
How do the current FEM opensource libraries compare?
Almost all FEM libraries are good enough, but I want to start with a FEM package and stick to it for some time. Instead of trying all of them, or going with what everyone else is using, I want to ...
- 394
0
votes
0
answers
110
views
Once a method is implemented in a computer program, then does learning the theory become less useful use of time?
Once a method is implemented in a computer program, then does learning the theory become less useful use of time?
This has confused me. I used to view that all things should be studied from first ...
- 427
1
vote
1
answer
77
views
Geometrically nonlinear finite element problem and mesh distortion
In my research on topology optimization of fluid-structure interaction problems (2D), I am using a geometrically nonlinear model to represent the structure. Body/surface forces are extracted from the ...
2
votes
1
answer
252
views
Finite Element Modelling of Hyperelastic Material under 2D Plane Strain Conditions
I am currently working on writing a MATLAB code for running a finite element simulation of a hyperelastic material in 2D. Since I am building this simulation as a part of a fluid-structure interaction ...
0
votes
1
answer
120
views
Problem with my Octave code (unsteady heat equation with FEM)
I want help with my Octave code regarding the unsteady heat equation.
My geometry and mesh are generated with FreeFEM++, so there is no problem with that (I tried it with the steady problem with no ...
- 3