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
6
votes
0answers
421 views

Understanding Boundary Condition in FEM

I am trying to understand Dirichlet and Neumann boundary conditions in FEM and I wanted to know if my inference is correct. To articulate my understanding, lets consider a simple case of TE and TM ...
0
votes
1answer
200 views

stiffness matrix for 3D regular grid in FEM

I have understood the stiffness matrix for 3D truss, and programmed Ku=f from scratch (in Java) to find the displacements. Then I moved to 3D solid but lost in too many concepts and equations, such ...
0
votes
1answer
514 views

Finite Element Analysis for Laminated Plates with Holes or Patches

As the title says, I am trying to code in FEM a plate structure that either has a hole in one of the layers or one of the layers is made of patches of plates, rather than one whole plate. However, ...
1
vote
1answer
142 views

Parallel calculation in finite elements

I am trying to solve a 1 Dimensional eigenvalue of poisson problem: $$\nabla \phi ^2 +\nabla \phi = k\phi$$ with the boundary condition: $\phi (0)=0 , \nabla \phi(1) = 0 $. I could solve this ...
0
votes
1answer
67 views

Principle of virtual work - extra term needed for deformation dependent loading?

I'm working on a problem in nonlinear elasticity, for which the external forces (loadings) depend on the displacements. Following Klaus-Jürgen Bathe's book "Finite Element Procedures", the virtual ...
3
votes
1answer
236 views

FEM current toy problem

I am solving the Dirichlet problem $$ \begin{cases} \Delta u = 0, \\ u|_{\partial D} = f, \end{cases} $$ in a $2d$ domain $D$ using the finite element method. What I want to get is the ...
8
votes
1answer
246 views

Resources on mesh generation for finite element methods

I know that this is not really apart of the rules as this is a recommendation question and these don't really have an answer per say. But, like this forum posting: https://stackoverflow.com/questions/...
4
votes
1answer
970 views

How to apply non zero Dirichlet boundary condition in finite elements?

I am writing a code for steady state heat transfer on a rectangular domain. I am specifying temperature on the edges - nonzero Dirichlet boundary condition. The equations can be written in form of $$...
0
votes
1answer
94 views

What is $H^{\frac{1}{2}}(\Gamma_g)$

In the second answer here one considers a function $w\in H^{\frac{1}{2}}(\Gamma_g)$. What is this space and what is the finite element functions that one should use that belongs to this space? I have ...
3
votes
2answers
332 views

Finite Difference and Finite Volume as special cases of Finite Element

I have heard this many times that "FD is a special case of FE" and separately, that "FV is a special case of FE". However I have never come across a brief document that substantiates that claim, ...
2
votes
2answers
185 views

Question about motivation of mixed FE method

I'm trying to understand what a mixed FE method is, and I was recommended this pdf: http://www.ima.umn.edu/~arnold/papers/mixed.pdf (Arnold, Mixed Finite Element Methods for Elliptic Problems). On ...
1
vote
1answer
209 views

Raviart Thomas Mixed Finite Element with Mixed boundary conditions reference request

This is perhaps a more focused version of this question. Using standard notation, I have a code that works and solves the following PDE using a Raviart Thomas mixed method. $$\begin{align} 0 &= ...
6
votes
0answers
128 views

Do practice and theory differ substantially when implementing Neumann Boundary Conditions using a Mixed Method?

I have implemented a pretty straightforward finite element solver for the following Poisson equation. For the purposes of this question we can assume the source term and the Dirichlet data both ...
8
votes
2answers
309 views

Lagrange multipliers space is too rich in a mathematical view

Background: Lagrange multiplier method has been employed in numerous fields, such as contact problems, material interfaces, phase transformation, stiff constraints or sliding along interfaces. It is ...
1
vote
0answers
201 views

the augmented global stiffness matrix is not positive semi-definite using Lagrange Multipliers method within FEM

The augmented global stiffness matrix is not positive semi-definite when using Lagrange Multipliers method to enforce boundary constraints on a simple square domain of integral form: I am considering ...
0
votes
1answer
855 views

Stiffness matrix computation for 4 node quadrilateral element

I am writing a finite element code for heat transfer (scalar field problem) and starting from simple 4 node quadrilateral element. I tried computing conductance (stiffness) matrix in the physical ...
2
votes
1answer
185 views

Galerkin method for a system of nonlinear PDEs

Suppose I have a nonlinear system of PDEs. I am actually interested in Navier-Stokes, but, for the sake of simplicity and example, suppose I had $$ \frac{\partial f}{\partial t} - f \frac{\partial g}{...
3
votes
0answers
95 views

Piecewise Constant Enrichment of the continuous galerkin method

I am interested to study crack propagation in a hyperelastic material in a variational setting. The crack surface exhibits a jump discontinuity. The function space for displacement field should ...
3
votes
1answer
130 views

Discontinuous Galerkin FEM : Control points are mid-points of edges instead of nodes

I am thinking to use discontinuous galerkin FEM (DGFEM) method to estimate discontinuous displacement field $u: \Omega \rightarrow \mathbb{R}^2$ at the crack surface of a material. The domain is ...
6
votes
1answer
176 views

Pressure definition/convergence issues for the Incompressible Navier-Stokes when using a stabilized P1-P1 finite element formulation

I believe this might be a recurring topic, but i have not found a post that directly related to this issue. I come from a finite volume background and my experience is more with predictor-corrector ...
1
vote
0answers
63 views

Plotting yield surface for geotechnical constitutive model

I have found the following example which is produced as documentation for RS3 geotechnical FE software. This is a validation document that goes through an analytical verification of triaxial testing ...
1
vote
1answer
135 views

Comparison between FEM libraries & languages

Is there any modern online resources which compares the most popular finite element method libraries/packages/languages?
1
vote
1answer
121 views

Correct way to model an embedded reinforcement (non linear FEM)?

I need to add to an existing FEM solver some embedded reinforcement element. This would give me the possibility to model/solve concrete structure (reinforced with steel rebar) taking into account the ...
2
votes
1answer
129 views

Some questions on Trace (operators) on the boundary in the context of PDEs

Background: The solution space of original problem (which requires a fine enough mesh to resolve the microstructure) can be split into a macroscale solution space and microscale solution space. This ...
5
votes
2answers
876 views

What does optimal convergence rate mean? How to prove whether a numerical method, e.g. FEM, has an optimal convergence rate or not?

I searched online about the way optimal is defined mathematically,but without any information acquired? How to prove whether a numerical method, e.g. FEM, has an optimal convergence rate or not? ...
2
votes
2answers
262 views

error estimator VS error indicator in the context of FEM error estimation

Is an error indicator an index for size measurement of something? And how is the relationship between the error estimator and error indicator? Is the error indicator just used for the (e.g. derivative)...
3
votes
1answer
145 views

Bubnov-Galerkin method in 1D: how to handle convective-type nonlinearity?

Consider the BVP: find $u = u(x)$, for $x \in (0,1)$ that satisfies \begin{align} u'' + u u' = f, \\ u'(0) = g_n, u(1) = g_d. \end{align} To derive the weak form for this BVP, we multiply the first ...
1
vote
1answer
78 views

What does Consistent Variation mean?

Consistent variation occurs in papers multiple times when reading . So what does it really mean from a maths view? I searched online and found an answer for it. Someone explained it as ‘Solution is ...
0
votes
1answer
87 views

Is FEM discretization error equivalent to quadrature error of bi-linear term in the weak form?

In FEM, to discrete the weak form, elements are needed to build the FE approximation basis space. FE method involves integration of the weak form. Due to discretization, the solution is an ...
1
vote
3answers
599 views

Create mesh for complicated 3D object for finite element analysis

I see images of steel connections, concrete dams, and other complicated 3D objects in papers which finite element analysis has been performed on them. My questions are: How these objects are created ...
3
votes
1answer
495 views

p-refinement in adaptive methods

Some reference in adaptive techniques say that when the solution $u$ is smooth enough we can use p-refinement instead of h-refinement. And when we have for example singularity, we should use h-...
1
vote
1answer
175 views

Barycentric interpolation equivalent for irregular hexahedra

I have a mesh with irregular hexaedra and I need a fast way to interpolate values at points inside these cells. I know that trilinear interpolation does not work well for large skews. Barycentric ...
6
votes
1answer
226 views

Why is the continuous Galerkin Finite Element Method a poor choice for the inverse problem for the Navier-Lame equation?

I work in a field (elasticity reconstruction) that frequently uses standard Finite Element methods to solve the first order PDE, where we are given u and solve for mu and lambda: $$ (\mu(u_{i,j} + u_{...
0
votes
3answers
206 views

quadratic trial functions for a 2d FEM calculation

I want to solve \begin{align} \nabla^4\psi+\alpha\nabla^2\psi+\beta\psi=F(x,y);\quad \nabla \psi\cdot \hat{n}=\nabla^3\psi\cdot \hat{n}=0\quad \text{on boundaries} \end{align} with a 2d FEM scheme. ...
1
vote
1answer
2k views

Determinant of jacobian matrix

I am using an FEM code written in Fortran which I did not design. For a particular problem, the program complains that the determinant of the Jacobian matrix is inferior to zero. I vaguely ...
4
votes
1answer
128 views

Reference Request: Raviart Thomas with hanging nodes

I am interested in reading about the analysis (existence, uniqueness, error estimates) of elliptic problems solved with a Mixed method that uses the Raviart Thomas elements (so far so good, easy to ...
7
votes
3answers
2k views

Galerkin method: Test functions vs. Basis functions

I'm a novice to finite element and I'm finding quite hard to find the actual difference between Test function(s) and Basis function(s). I would be glad if somone could explain me that and point out ...
1
vote
0answers
424 views

FEM Stiffness and Mass matrices for 2d cubic trial functions

I want to use FEM to solve a 2D eigenvalue problem for the biharmonic equation \begin{align} &K \psi= \lambda M \psi\\ &\psi=\nabla \psi\cdot \hat{n}=0\quad \text{on boundaries} \end{align} ...
0
votes
2answers
292 views

Boundary conditions for streamlines in enclosed flow

I am trying to solve Lid driven square cavity flow problem of Stokes equation using finite element method. I have boundary conditions for velocity as zeros on every boundary but u=1 on top boundary. ...
2
votes
2answers
180 views

How to precondition FEM problems using domain decomposition?

Let's say that I have a FEM code which yields the following problem: $$ \mathbf{A}\mathbf{x} = \mathbf{b}. $$ In order to solve this more efficiently with an iterative method, I would like to ...
1
vote
3answers
770 views

Pressure boundary condition in lid driven cavity using finite element method

Thank you all 1.) I am trying to solve lid driven cavity problem for an incompressible Stokes and Navier Stokes equations using general "Mixed" finite element method. dirchlet boundary conditions are ...
0
votes
2answers
336 views

Imposing total pressure over surface in FEM

I am trying to solve Stokes problem using Finite element method. My question is how to impose that total pressure over the surface is zero to remove the constant pressure mode?
1
vote
0answers
50 views

Computing dilogarithm

I'm measuring the integral of a quantity which, mathematically, requires the computation of a dilogarithm function. $$\operatorname{Li}_2(be^{ax})$$ where $b$ and $a$ (are real and) can be positive ...
1
vote
1answer
163 views

FiPy: Make diffusion coefficient dependent to orientation

I'm trying to solve the heat equation in FiPy and right now it works, but I have one problem: The material I have to simulate has different diffusion coefficients in the x and y direction (due to its ...
1
vote
1answer
975 views

Line integral along the edge of an isoparametrically mapped triangle

I need to integrate the following function on the line segment from $P_{1} = \begin{bmatrix} -2\\-1 \end{bmatrix}$ to $P_{2} = \begin{bmatrix} 1\\2 \end{bmatrix}$: $$\int_{P_{1}}^{P_{2}} 4x + y \ ds$$...
1
vote
0answers
97 views

Do DG methods for the Helmholtz equation always return positive quantities?

Helmholtz Diffusion equation with reaction term: $$ k\Delta u + u = f ~ \text{in} ~\Omega $$ $$ \nabla u \cdot \mathbf{n} = 0 ~ \text{in} ~\partial \Omega $$ For sufficiently small $k$ (relative to ...
1
vote
1answer
714 views

Computation of stiffness matrix with variable coefficient

I am implementing a finite element solver (in 2D) to solve the generic differential equation : $$-\nabla(a(x) \nabla u) = f$$ Brief explanation By integrating and multipling by a test function, the ...
0
votes
1answer
100 views

Need clarification on a piece of book excerpt about spectral element method!

I am reading "Using MPI (3rd edition)" from William Gropp, where in chapter 4 application section 4.13, it introduces an MPI application Nek5000/NekCEM which is based on spectral element method (SEM) ...
1
vote
1answer
76 views

How to deal with multi-region problems (in PETSc)

Consider a problem where space is filled with a liquid and a solid phase, with a large, complicated geometry. On each of the phases, there's an electrical potential/poisson equation. The equations are ...
1
vote
1answer
319 views

Algebraic multigrid in PETSc

Consider a potential/poisson equation on a very large, complicated geometry. Currently, an self-written FEM and linear solvers from NumPy are used. Performance is, of course, not good enough for ...

1 3 4 5 6 7 17