# 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.

973 questions
Filter by
Sorted by
Tagged with
50 views

72 views

### Galerkin Least-Squares stabilization for FEM solution advection (hyperbolic) equations

I am playing with Galerkin Least-Squares stabilization to solve advection diffusion problem in the context of the finite element method. This works very well for steady-state advection-diffusion ...
72 views

### Error $L_{2}$ convergence in Finite Element for Poisson Equation

I have written a Matlab code to solve the equation $-u'' = f$ with conditions $u(0) = u'(1) = 0$ on the domain $x \in [0,1]$. I tested the code with $f(x) = -2, \forall x\in [0,1]$. I check the plot ...
98 views

### Efficient schemes for solving the extended Saddle point problem

I am interested in knowing some efficient techniques for solving the following extended Saddle point problem. \begin{align} \begin{bmatrix} A & B^T & C^T \\ B & 0 & 0 \\ C & ...
122 views

### Finite element method for high-frequency electromagnetics

I am writing a project about the Finite element method for use in high-frequency solutions of Maxwell's equations. This could be for use in antenna design and similar. I have some trouble ...
80 views

### Boundary conditions for an FEM approximation of the Laplace operator

Using FEM, I want to approximate the Laplacian $$u = \nabla \cdot \nabla h \, ,$$ where $h(x,y)$ is an FEM approximated scalar field on the same mesh, i.e. piecewise differentiable. I am using MOOSE ...
87 views

### preconditioner for Laplace “without” boundary values

I'm looking at solving systems with the FEM discretization $$-\int_\Omega (\Delta u) v = \int_\Omega \nabla u \cdot \nabla v - \int_\Gamma (n\cdot\nabla u) v.$$ without applying Dirichlet- or ...
116 views

### Step3 in deal.II - Convergence of the mean

I'm trying to understand the Convergence of the mean part of the Step-3 tutorial in deal.II. The authors say that $\frac{1}{|\Omega|}\int_{\Omega} u_h(x)dx$ converges with $\mathcal{O}(h^2)$, but I ...
50 views

72 views

### Lumped mass matrices for higher-order finite elements for CFD

Given that some of the mass lumping techniques, for example, row-sum lumping does not produce practically viable lumped mass matrices for all the element shapes, what are the techniques used for mass ...
97 views

### How to apply the boundary condition when global stiffness matrix is stored in csr format? [duplicate]

I am solving the poisson equation and I constructed the global stiffness matrix in compressed row storage format. Then I wrote the preconditioned conjugate gradient solver for solving the system of ...
51 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 ...
74 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 ...
79 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, \...
76 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 ...
189 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 ...
37 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$ ...
118 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 (...
73 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 ...
111 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 ...
88 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 ...
97 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 ...
28 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 ...
39 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. ...
64 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....
89 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 ...
103 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 ...
69 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 ...
112 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 ...
53 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. ...
29 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 ...
59 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 ...
83 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.