# 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,046 questions
Filter by
Sorted by
Tagged with
58 views

### Solving Laplace equation with constraint on boundary

I have found the following PDE problem in a paper: Essentially, we have a rectangular domain where there is a unknown interface $z=\xi(x)$ (liquid-air interface) separating the domain into two medium ...
93 views

### 1D Poisson equation and quadratic basis functions assembly

I'm solving the simple Poisson problem $$-u''(x)=1$$ in the interval $[0,1]$ with $u(0)=u(1)=0$. I discretised my domain as done here, i.e. with ...
84 views

### Discretization of a nonlinear boundary value problem

I am trying to use finite element method to discretize the following problem \begin{align} \min_{u \in H^1_0(\Omega)} \int \| \Delta u(x) - 0.5*[u(x) + \langle e, x \rangle + 1]^3 \|^2_2 \ d\Omega, \...
82 views

### Renumbering the nodes for quadratic basis functions for a 2D domain

I have a simple triangulation for a 2D domain, described by the connectivity matrix $T$ and by the point matrix $P$. For didactic purposes, I assembled the stiffness matrix for $-\Delta u = f$ by ...
30 views

### Differences in using Clausius-Duhem inequality vs Principle of Virtual Work/Power in derriving constitutive equations?

I am a novice getting my toes wet in continuum mechanics and nonlinear elasticity. I have seen papers that use both approaches in developing constitutive connections to compliment balance equations of ...
214 views

### Who uses finite elements with higher continuity?

Lagrange elements of any polynomial describe piecewise continuous functions. Typically, those functions are differentiable. Mixed finite element methods use vector fields of even less continuity, such ...
52 views

### Spline interpolation for vector-valued data in 3D space

I have output from a 3D linear elasticity finite element simulation which uses linear tetrahedral elements, such that the displacement is continuous over the nodes but the gradient is not ($C_0$ ...
137 views

### FEM solution for Poisson is not exact at nodes

Let's consider the usual Poisson problem for FEM $$-u''(x)=1 \quad x \in [0,1]$$ with homogeneous Dirichlet boundary conditions. The solution is $u(x)=-\frac{1}{2}x(x-1)$ I know that the FEM solution (...
80 views

### C.S.R method in finite element matrix assembly

I have solved the 2D Poisson equation using finite element method with simplex triangular element in MATLAB. First, I generated the triangular mesh using pdetool in ...
117 views

### diffusivity matrix assembly in nonlinear finite element analysis

I want to solve a diffusion analysis using finite elements. According to fick's law, governing equation is $$\frac{\partial h}{\partial t} = D \frac{\partial^{2} h}{\partial x^{2}}$$ . h is relative ...
141 views

### Poisson equation, stiffness matrix positive definiteness, Dirichlet boundary conditions

I have a question regarding the positive definiteness of the stiffness matrix. Specifically, I believe that it should be positive definite only when at least one Dirichlet point is given, so I would ...
105 views

### Finite element interpolation

I have a finite element solver that I implemented in MATLAB. I am calculating a specific potential function that is constant in each tetrahedral element. The question is, I want to be able to ...
30 views

### velocity in CFL condition

I am studying the evolution of the density and velocity field of a core in a molecular cloud in 1 D. I defined the radial grid (let us say x between 100 and 101) and the time grid. I am using the ...
43 views

### Fourier transform in finite element

I have a finite element solver where I am using tetrahedral elements. I am solving for electric potentials and then calculate the current densities in each element, which are constant in each element. ...
69 views

### Sparse linear solver in fortran working with REAL16

I need some (direct) sparse linear solver for fortran, which works with REAL16 data type. Any suggestions? Both Pardiso and MUMPS support only REAL8. (identical question: https://math.stackexchange....
120 views

### Finite elements with CFL condition - How to obtain correct order of convergence

I have discretized a PDE with continuous finite element method in spatial variable and with implicit Euler or Crank-Nicolson in temporal variable. In both cases, I have error estimates in $L_2$ norm ...
141 views

### L2 projection with bounds

In some problems we are currently working on, we are working with discontinuous functions that are defined on a finite element mesh and are established using Lagrangian particles. To obtain them on a ...
79 views

### relation between different tangent stiffness

I need to find a relation between the tangent stiffness $L_1$ of the first Piola-Kirchhoff stress tensor with the tangents stiffness $L_2$ of the second Piola-kirchoff stress tensor. They are defined ...
141 views

### How can I implement second order derivatives of shape functions of a 3D elements?

I am developing an Abaqus UEL with 3D 8 nodes brick elements and I need second order derivatives of the shape functions, I have already mapped the first order derivatives from the element coordinates ...
56 views

### Mesh partitioning with METIS

I am trying to use METIS-5.1.0 edition in order to partion a FE mesh. For demostration purposes I created 2x2 rectangle mesh and tried to partition it. However, I notice a weird behaviour in my code. ...
38 views

### Output solution vector in Petsc

I am using petsc to solve a linear elasticity problem discretized by finite elements.The initial mesh is read by a mesh file and the distribution in each processor is done using METIS.I am using only ...
70 views

### How to solve a Stiffness matrix?

I am quite new in this field. For my university I prepared a stiffness matrix to solve for a project group. This matrix consists of 450 equations with 450 unknows (it's a Matlab script) and I have the ...
90 views

### Fortran compiles for legacy Finite Element Fortran program (1980)

The version of Fortran used come from Montreal Ecole Polytechical in 1980. I need a compiler for Fortran for Windows 7 or Windows 8.
85 views

### Technique to find the CFL condition using the Galerkin method in space and finite-difference in time?

I am using the Galerkin method (Discontinuous to be precise) to discretize in space the scalar linear wave equation and the explicit second order centered finite difference scheme to discretize in ...
151 views

### Solving PDEs in parallel

I have read different approaches on how to solve pdes in parallel which are discretized using finite element method. For example: Non-overlapping domain decomposition approach as mentioned in https://...
47 views

### Finite element analysis software for acoustic and electrostatic

I need to do a simulation for my thesis project involving some piezoelectric nanoparticles in a fluid beamed with ultrasounds. I'm looking for a software for such simulation and for now it seems me ...
182 views

### How to compute gradient of each node in finite element method?

I have a problem with computing gradient of each node in finite element method. I can get the value of each node. But how can I get the gradient? I know $u = \sum u_i \phi_i$ where $\phi_i$ are the ...
143 views

### Finite Element Analysis based on contractions of a subset of edges in enterior of mesh

I am modelling a problem that is "driven", not by typical boundary conditions but, by contractions in its interior. In a finite element analysis, I can specify the new lengths (not ...
80 views

### Transition from 2D to 3D finite element code, what are the inevitable modifications to be implemented?

Imagine we have a simple 2D FEM solver (we are dealing with solid mechanics) and we would like to develop it to a 3D FEM solver (let's say for the same solid mechanics problem) in this case what are ...
180 views

### Problem with solving coupled ODE and DAE equations with mass matrix (Error using daeic12 (line 77) This DAE appears to be of index greater than 1)

I am trying to solve 6 ODE equations coupled with 1 DAE one. The ODE equations have been discritized in space domain and ode15s MATLAB solver is used to solve the equations in time domain. I have ...
509 views

### Gauss-Lobatto quadrature and nodal points for FEM

By using the Legendre-Gauss-Lobatto (LGL) quadrature formula (QF) and LGL nodal points one achives a diagonal mass-matrix for finite element problems. (More specifically, the spectral element method.) ...
655 views

### FEM Python book

Is there any book or site available with Finite element Method for partial differential equations with python code apart from Fenics?
86 views

### inclined/general Dirichlet boundary conditions

For simpilcity, consider a single quad linear elasticity finite element in 2D. The Dirichlet boundary conditions on node 1 and node 2 are easy to implement and can be handled in the standard way. ...
47 views

### Learning the art/science of structural idealization

I am a mechanical engineer working in the field of aerospace structures. During the course of my studies, I have studied a course on structural analysis in which I learned 3D Euler-Bernoulli beam ...
99 views

### Which 2D PDE with an exact solution can I use to test/verify my FEM-PDE code?

I have created a program to solve 2D, time-dependent PDEs with the finite element method and get reasonable looking results for the 2D acoustic wave equation. Now I would like to go further and solve ...
54 views

### How can I improve the accuracy of the calculation of the magnetic field in Gmsh/GetDP?

I need to calculate the magnetic field along a straight line in proximity of an array of 6 magnets. I used the tutorial files "magnets" included in Gmsh and I slightly modified the file in ...
28 views

### Mesh transition between beam and shell element types

I am using NASTRAN solver & FEMAP as preprocessor for reduce modelling of wing using 1D and 2D finite elements. Beside transition of 1D & 2D elements to 3D, I had not found any method/solution ...
86 views

37 views

### How can the choice of coarsening factor affect Multigrid's convergence?

The linear system $Ax=b$ is coming from the discretization of an elliptic PDE. Multigrid method is used in order to solve it. Suppose $c_0$ is the coarsening factor on level 0 and $c_m$ the coarsening ...
151 views

### how to Implement linear tetrahedral elements for finite element computations?

I am trying to implement 3D tetrahedral elements in my finite element code (which works fine for linear triangles and quadrangles in 2D). But my simulations are crashing with tetrahedral elements. My ...
188 views

### Heat equation in non-dimensional form behaving differently than in usual format

Starting from $$c_p \frac{\partial u }{\partial t} = k \nabla^2 u$$ in a one dimensional domain [0,1] where $c_p$ and $k$ are modeling two different materials:  k = \begin{cases} 1 ~\text{if} ~x &...
151 views

### Parallel In Time with Multigrid

I am trying to solve the linear finite element equation $M\ddot{u}+Ku=F(t)$, where $M$ is the mass matrix ,$K$ the stiffness matrix and $F(t)$ the external load vector, parallel in time using XBraid ...
68 views

### Solution predictors for accelerating convergence in nonlinear FEM

I am looking for the details of commonly-used predictors for accelerating the convergence of iterations using Newton-Raphson scheme for nonlinear problems in FEM. I am looking specifically for static ...
I am trying to solve $M\ddot{u}=-Ku+F_\text{ext}$ for a 2D linear elastic model with $M$ be the mass matrix,$K$ the stiffness matrix and $F_\text{ext}$ the external load vector coming from a uniformly ...