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
121 views

Iterative solver for high order DG methods (3D Laplace problem)

I have a 3D Laplace problem on quite a complicated geometry where I am using Discontinuous Galerkin method. My mesh is composed by hexas, hence I am employing classical tensor product basis functions $...
FEGirl's user avatar
  • 311
1 vote
0 answers
42 views

Jacobian for 6-noded triangle in 3D to calculate the area

I would like to calculate the surface area of a 6-noded triangle element, i.e., the face of a 10-noded tetrahedral element in 3D space. A typical solution is to calculate the surface integral of the ...
reox's user avatar
  • 113
2 votes
1 answer
99 views

FEM textbooks recommendation

I would really appreciate if you could recommend me some textbooks that explain FEM process well and have many exercises to solidify the knowledge. So far, Brenner and Zienkiewicz seem to be a bit ...
Nomad's user avatar
  • 45
0 votes
0 answers
23 views

maxwell and kelvin models via prony series in abaqus

abaqus documentation only shows generalized models and their formulae. What are the prony series in abaqus for the simplest maxwell and kelvin models?
feynman's user avatar
  • 269
2 votes
2 answers
31 views

How do we implement the balance of stress on interface in ALE FSI method?

""we aim for a consistent variational-monolithic coupling scheme in which we need all equations defined in the same domain; therefore, $\mathrm{ALE}_{f x}$ was introduced. In variational ...
吴yuer's user avatar
  • 193
1 vote
2 answers
102 views

How to compute overall inertia properties from FE mass matrix?

How should I evaluate the overall mass and moments of inertia (polar and transverse) of a finite elements model, having its mass matrix? Given that the global mass matrix is composed of symmetric mass ...
temporary_pigeon's user avatar
0 votes
1 answer
66 views

mesh size restriction of spatial discretization in FEMs and FDMs

Is there any mesh size restriction of spatial discretization in FEMs and FDMs (finite difference)? If a mesh is very coarse there is still perhaps nothing to stop the program from running and cause ...
feynman's user avatar
  • 269
0 votes
1 answer
72 views

how Petrov-Galerkin FEM different from others and advantages

How is the Petrov-Galerkin FEM different from more common FEMs? Could anyone give a layman description? What's the advantages of using the Petrov-Galerkin FEM?
feynman's user avatar
  • 269
1 vote
1 answer
137 views

Nodal functionals in finite element analysis

I have a quintic Hermite basis functions [-1,1] for FEM applications, I wanted to check if it's nodal functionals are proper. Could someone explain nodal functionals in details and give an example of ...
Nomad's user avatar
  • 45
1 vote
1 answer
126 views

Simple, easy to install and use Python FEM solver (and example) for 2D cylindrical Laplace equation

In my earlier question Interpolating the gradient of a cylindrically symmetric potential field expected obey the Laplace equation, especially near/across r=0 I give a simple example of a Jacobi ...
uhoh's user avatar
  • 1,038
0 votes
1 answer
50 views

How to implement the interface extension of fluid "displacement" in ALE?

In ALE, we first set a referenced space for fluid, then we extend the boundary fluid displacement to the whole fluid region, take harmonic extension as an example, we need $$\Delta \left ( \hat{u} \...
吴yuer's user avatar
  • 193
0 votes
1 answer
101 views

Keeping the surface flat for finite element analysis

I ran the FEM analysis for a simple supported beam problem using solid elements in plane strain condition, with linear kinematics and linear elastic material. I'm not happy with the solution since ...
kstn's user avatar
  • 231
0 votes
1 answer
71 views

Any code on ALE method for N-S equations or advection-diffusion equations?

I'm currently reading the book computational fluid-structure interaction, but the book doesn't provide any code, and I can't understand the ALE method, could you please provide some codes written in ...
吴yuer's user avatar
  • 193
1 vote
1 answer
91 views

Reference request: graph Laplacian approximation for domains/manifolds

Is there any reference on solving the heat equation on irregular geometries via creating a mesh and using the graph Laplacian instead of FEM techniques? (Convergence to real solution). That is to say, ...
Aner's user avatar
  • 151
0 votes
0 answers
38 views

Continuous Finite Element vs Space Frames

Finite Element allows one to divide a continuous medium into a mesh and compute structural properties of it. Instead of discretizing the continuous, why not start with a discrete method to begin with? ...
BBSysDyn's user avatar
  • 239
1 vote
1 answer
75 views

How to generate mesh for space-time FEM method in FEniCS?

suppose that I have a region called $\Omega $,suppose that this is a 2D or 3D region, how to generate the corresponding space-time mesh by using FEniCS? and how to extract the boundary and implement ...
吴yuer's user avatar
  • 193
1 vote
1 answer
57 views

Rolling friction settings in COMMERCIAL (Ansys) finite element software

Is there a standard way of implementing rolling friction in COMMERCIAL finite element software? There's no book about such a subject in COMMERCIAL FEM software though there are books on algorithms ...
feynman's user avatar
  • 269
0 votes
1 answer
138 views

Quintic Hermite shape functions

I am trying to use quintic Hermite basis functions for FEM applications, could someone please direct me to the general formula that would help me generate quintic Hermite shape functions? In natural ...
Nomad's user avatar
  • 45
0 votes
1 answer
49 views

FEA with order reduced for a large system with a very localized force application

I have a large structure with many DOFs. But the application of force is very localized say over one or two nodes. If I try to run FEA on this with Guyan Reduction, I think the efficiency of reduction ...
s6292_1997's user avatar
2 votes
1 answer
94 views

Why can the weak forms of the Stokes and continuity equations be combined into a single equation?

Consider this Stokes equations, $$ \left\{ \begin{array}{r} - \mu \Delta \vec{u} + \nabla P = \vec{f} \\ \nabla \cdot \vec{u} = 0 \end{array} \right. $$ Weak form I is: $$ \...
Hao's user avatar
  • 21
2 votes
1 answer
232 views

Solving of KU=F leads to numpy.linalg.LinAlgError: Singular matrix

Solver.py: ...
famatto's user avatar
  • 23
2 votes
1 answer
128 views

How to calculte the contributions of a lift operator in FEM?

This question is largely motivated by the famous paper "UNIFIED ANALYSIS OF DISCONTINUOUS GALERKIN METHODS FOR ELLIPTIC PROBLEMS" The Bassi and Rebay biliniar DG form is given in page 11 in ...
CuteCompute's user avatar
1 vote
2 answers
186 views

Query about FE approximation of a Poisson equation with non-constant coefficients

Consider the standard weak form of the Poisson equation with coefficient $\alpha$: \begin{equation} \int_\Omega \dfrac{\partial v}{\partial x} \dfrac{ \partial \left( \alpha u \right) }{\partial x} = ...
CuteCompute's user avatar
0 votes
0 answers
33 views

How to embed linear elasticity/constrain solver in non-linear soft-body dynamics

I often do simulation of dynamics of mass points connected by strings (e.g. molecular dynamics, soft body dynamics etc.). Typically I do it simply by integration of equations of motion by e.g. verlet ...
Prokop Hapala's user avatar
0 votes
0 answers
51 views

How to get Gauss points for nodal force vector (surface integral) in tetrahedral elements in the isoparametric coordinate system $({\xi \eta \zeta})$?

I am aware of the following question Evaluating the surface integral in an FEM (Finite Elements Method) procedure. But they use the volumetric coordinates while I want to use the cartesian ...
Mr Thomas Anderson's user avatar
2 votes
0 answers
90 views

Why multigrid is inefficient?

I am trying to solve the Stokes equation containing viscosity nonlinearity using the open source finite element software underworld2 with nested PETSc. The resolution is 2000*200. The solution results ...
Darcy's user avatar
  • 21
1 vote
0 answers
82 views

Seeking open-source PDE Solver for inhomogeneous material properties

I'm currently in search of an open-source PDE solver (Finite Element Method is preferred) that can effectively handle the challenge of material properties coefficients associated with each element in ...
Sadjad Abedi's user avatar
3 votes
2 answers
178 views

what preconditioner for incompressible hyperelasticity in 3d (similar to stokes equation?)?

I am working on modeling incompressible elasticity at finite strains. $$ \mathrm{Div} \boldsymbol P = \boldsymbol 0, \quad \boldsymbol X \in \Omega_0 \subset \mathbb R^3, \\ J = 1, \qquad \boldsymbol ...
Simon's user avatar
  • 165
0 votes
0 answers
20 views

uniform refinement is not working in gmsh

I get 0 node and elements when I use -refine command line flag in gmsh. Explicitly writing RefineMesh; in the script produces the original mesh and not a refined one. Am I missing something while ...
Ashb's user avatar
  • 11
1 vote
0 answers
30 views

Recommendations for some new books about computational contact mechanics in solid mechanics

I want to simulate some frictional contact problems, but I'm not familiar with this field, could you please recommend some new books as introductions?thank you
吴yuer's user avatar
  • 193
1 vote
0 answers
56 views

3D Quadrature schemes with points on boundary

In one dimension there are two types of quadrature schemes. asymmetric rules like Newton-Cotes like formulas (Trapezodi, Simpson), and Clenshaw-Curtis place sampling points on boundary of the ...
Prokop Hapala's user avatar
0 votes
0 answers
54 views

eigenvalues of inhomogeneous Helmholtz equation violate superposition using FEM

I am trying to solve the in-homogeneous Helmholtz equation with damping and forcing using finite element method (FEM) with FEniCSx. The equation is; $c^2\nabla^2p - c^2\omega i d \nabla^2p + \omega^2 ...
Ekrem Ekici's user avatar
1 vote
1 answer
190 views

Solving PDEs using FEM using cubic Hermite polynomials

everyone. I am a beginner in Numerical mathematics, I have some idea of how to use Galerkin method to solve PDEs numerically, but so far I had no luck finding an example of how to solve a simple PDE ...
Nomad's user avatar
  • 45
2 votes
1 answer
94 views

references for optimization in the context of parameter identification with finite elements

i am performing parameter identification for a non-linear partial differential equation (elasticity) that I solve with finite elements. My optimization problem is a non-linear least squares data-...
Simon's user avatar
  • 165
1 vote
1 answer
123 views

A confusion about the bubble function in lumped mass FEM

I am studying knowledge related to lumped mass finite elements. As is well known, lumped mass finite element methods higher than 2nd order on simplex mesh require the construction of new function ...
Owen Jun's user avatar
  • 131
0 votes
0 answers
59 views

How to set Neumann BC for coupled transport problem in weak form?

Consider $$\begin{aligned} \partial_t v + b\cdot \nabla \phi &=0 \\ \partial_t \phi + b\cdot \nabla v &= 0 \end{aligned}$$ for $v:(x,y)\mapsto \mathbb{R}$, unknown and time-dependent ($...
l'étudiant's user avatar
3 votes
0 answers
72 views

Confusion between standard finite element and mass-lumping finite element methods

Consider the following equation, subject to homogeneous Neumann boundary condition. $$ u_t = \Delta u + f(u). $$ The weak formulation is as follows: $$ (u_t,w) = (\Delta u, w) + (f(u),w), \quad \...
Owen Jun's user avatar
  • 131
3 votes
0 answers
190 views

Most promising reduced order modeling method

Many players in the field of engineering simulation software are investing on digital twinning and reduced order modeling techniques, meaning that the field bears potential. I was wondering if among ...
Lilla's user avatar
  • 259
1 vote
1 answer
81 views

Is the weak form on this book about a level set FSI problem wrong?

I'm trying to repeat a level set FSI problem on the book <Level Set Methods for Fluid-Structure Interaction>, on the page 89, the provided freefem code define a weak form of the discretized ...
吴yuer's user avatar
  • 193
1 vote
0 answers
140 views

FEM for nonlinear first-order ODE

Currently I am trying to solve nonlinear Ricatti equation using FEM (Matlab language): $$\frac{d r(z)}{dz} = i k(z) r(z) + \frac{i k(z)}{2}(\epsilon(z) - 1)(1 + r(z))^2$$ $z = [-h/2, h/2]$ and $r(-h/2)...
Andrew's user avatar
  • 31
5 votes
1 answer
287 views

Does the weighted residual method not use energy minimization in any form?

I've come across several texts/papers utilizing the concept of a minimum potential energy state corresponding to an equilibrium state, and I know that it is used in FEM formulations that are based on ...
SNIreaPER's user avatar
2 votes
1 answer
193 views

FEM for a biphasic drying problem

I have the following PDE originating from a biphasic drying model, where $\xi \in [0,\Xi]$ is a radial coordinate attached to the dry skeleton of a wet cylindrical body: $$ \frac{\partial u}{\partial ...
IPribec's user avatar
  • 597
0 votes
0 answers
86 views

How to determine whether the symmetric stiffness matrix is positive definite or not? Is it related to the problem?

For two-dimensional or three-dimensional elliptic equations, when will the stiffness matrix be asymmetric and positive definite? This affected the solution efficiency so much that I had to choose an ...
Darcy's user avatar
  • 21
2 votes
2 answers
121 views

Defining the restriction operator in nested multigrid FEM

I'm new in the multigrid approach and I'm trying to implement an algorithm from this paper: https://www.researchgate.net/publication/242913687_A_Multigrid_Algorithm_for_the_p-Laplacian I'm stuck with ...
Timofej's user avatar
  • 21
0 votes
0 answers
87 views

Why aren't mortar domain decomposition techniques used as much as schwarts type DD?

Schwartz type domain decomposition techniques require a transmission condition which can be hard to come by. Mortar type techniques enforce continuity with a Lagrange multiplier across domains. Are ...
CuteCompute's user avatar
1 vote
0 answers
38 views

How does the time needed for force propagation increase with the discretization (using symplectic Euler schemes)?

I am trying to model a physical system by (lets say the system is a long deformable object, on which I can apply forces). It can be described by Cosserat Rod theory and discretized, by modelling it ...
the2second's user avatar
2 votes
1 answer
165 views

Does this second-order implicit Runge-Kutta method have a name?

I am studying the time-integration of the following paper, Young, L. C. (1981). A finite-element method for reservoir simulation. Society of Petroleum Engineers Journal, 21(01), 115-128. A copy (PDF)...
IPribec's user avatar
  • 597
0 votes
0 answers
46 views

recommendation on some papers/books about frontal solver used in FEM

I'm reading a program about computational plasticity, this program use frontal solver to solve the program, but I'm not familiar with frontal solver even after reading some papaers, so could you ...
吴yuer's user avatar
  • 193
1 vote
1 answer
97 views

Formulation of $-\mathrm{div}(k\nabla u)=f$ in $\Omega$ for the Finite Element Method

...
Phillip Russo's user avatar
0 votes
1 answer
76 views

In level set fluid structure interaction method why can we rewrite the elastic force in this form while there is no shear force?

when we consider a Immersed Membrane Without Shear, we can define a regularized elastic energy as $$\mathcal{E}(\varphi)=\int_{\Omega} E(|\nabla \varphi|) \frac{1}{\varepsilon} \zeta\left(\frac{\...
吴yuer's user avatar
  • 193

1
2 3 4 5
25