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

learn more… | top users | synonyms (1)

0
votes
1answer
27 views

How does the animation work in eigenvalue problem of FEM

I have used free vibration analysis in FEM. After analysis, we can usually use animation to see the motion of each eigenmode (In Abaqus or Comsol, I would choose either half harmonic or full ...
0
votes
1answer
52 views

openfoam - Programming customized PDEs

I am looking for a method to automatically solve custom PDEs on a custom control volume. Specifically I would like to solve equations similar (but not exactly alike) to: $$\frac{\partial y}{\partial ...
1
vote
0answers
47 views

Arpack and Matlab give different values for eigenvalues

I am solving a generalized eigenvalues problem with inversed complex shift: $$(M-\sigma J)^{-1}J \boldsymbol{x} = \boldsymbol{x} \nu \enspace .$$ My matrices are obtained from a finite element ...
1
vote
0answers
23 views

Integro-differential eigenvalue problem, COMSOL

In my research I encounter an eigenvalue integro-differential equation of the form: $$f_n(x,y)=\lambda_n\iint_D\frac{\nabla'\cdot\big\lbrace ...
2
votes
0answers
39 views

Selecting a Finite Element discretization for subsurface flows

Are there guidelines to follow in selecting the right FE discretization. Specifically, for subsurface flow models like Darcy's equation (both single and multi-phase), Richard's equation, etc. I know ...
4
votes
1answer
80 views

How to formulate lumped mass matrix in FEM

When solving time dependent PDE's using the finite element method, for example say the heat equation, if we use explicit time stepping then we have to solve a linear system because of the mass matrix. ...
1
vote
2answers
75 views

How to deal with nonlinear term in Navier Stokes equations (finite element code)

I am trying to solve the Navier Stokes equations using the finite element method. I plan on using the pressure correction method to deal with the pressure and an implicit time stepping scheme for ...
-2
votes
0answers
198 views

Sequence & convergent series [closed]

I have discovered the following solution (one real root for trinomial equation of bellow form), can you kindly check my solution, Thanks, Define a function f(x) Where $x$ is any real number between 0 ...
2
votes
1answer
101 views

Is my matrix symmetric?

I obtained a mass matrix through Finite Elements discretization. Now, I want to check if it is symmetric. To do that I subtract to my matrix $M$ its transposed $M^T$. The result is another matrix of ...
3
votes
1answer
72 views

Meshless Methods and the Kronecker Delta Property

In texts I often read that meshless methods, as opposed to the FEM, do not exhibit the Kronecker Delta Property [1]: $N_I(x_J)\neq\delta_{IJ} = \begin{cases}1 &\mbox{if } I=J \\0 ...
0
votes
0answers
29 views

Gauss Integration over Zero Order Element

I'm working with the Boundary Element Method and want to integrate an expression over a triangular region. I would like to use Gauss Integration to do this, but I'm having trouble since the triangular ...
0
votes
0answers
49 views

Proof of convergence for Newton method in finite element analysis

The Newton method in FEA is to solve for a non-linear equation where stiffness matrix is a function of displacement: And to assemble the stiffness matrix from constitutive laws, we need the tangent ...
3
votes
2answers
95 views

Is “tangent stiffness matrix” the same as “stiffness matrix”?

I'm trying to implement nonzero Displacement Boundary Conditions in VegaFEM on a non-linear model, using the method outlined in §3.6.2 of University of Colorado's intro to FEM (modify $f = Ku$: set ...
2
votes
2answers
133 views

The meaning of conservative discretization in Galerkin FEM and Discontinuous Galerkin

I do understand the meanning of "conservative discretization" within the FVM/FDM framework, indeed it is well explained in this post. Now, according to the table in this slide (pp.4), it concludes: ...
0
votes
1answer
64 views

Finite element method for odd order DE

What are theoretical hurdles in applying Galerikin method on, say, first order time dependent ODE? Is there no way we can form an inner product??
3
votes
1answer
62 views

Shell vs frame element model stiffness differences

I have a model of a tall, slender structure that I am investigating using both shell and 3D frame elements. The shell elements are type MITC4, 4-node membrane elements. The frame elements are the ...
2
votes
1answer
82 views

Inverse isoparametric mappings for quadrilateral finite elements

I have an isoparametric mapping $F_{E}: \hat{E} \to E$ where $\hat{E}$ is a reference quadrilateral (square) and $E$ is some quadrilateral in the domain $\Omega$. If $E$ is a parallelogram, then the ...
3
votes
1answer
134 views

Boundary conditions in conforming Galerkin method for biharmonic equation

I am trying to solve simple scalar biharmonic equation using bubnov-galerkin finite element method. I am using $H^2$ conforming basis functions. I was wondering that if anyone can give me some ...
4
votes
1answer
89 views

Find the direction of the gradient on a finite element mesh

Suppose we have a triangular mesh of a two dimensional shape $\Omega$, and on this mesh we define a P1 finite element structure. I know that given $u,v$ by their values at the vertices of the ...
0
votes
0answers
44 views

Best way to average surface data or function data in 2d grid

I have pressure data on a 2D triangular surface. I can compute average(data). What is the best way to average the surface data?
4
votes
1answer
41 views

Estimating the local compression/expansion ratio for a transformation on a point cloud

Let's say we have an unorganized point cloud P1 with N points, each with coordinates {x,y,z}. We apply non-rigid transformation to P1 (translation + rotation + warping), to obtain point cloud P2. ...
0
votes
1answer
72 views

FEM - Shape function of a HEX20 - plot in MATLAB

I have a FE model of a simple plate with hole (tension load) with HEX20 mesh. I need to obtain the Shape Function of one of the elements (the one with highest stress) and plot it (with MATLAB). After ...
0
votes
0answers
62 views

Conditions for always positive gradient of heat field in evolutionary thermo-elastic system

I am investigating stability and convergence of series of approximations for coupled thermoelasticity problem yielded by one-step recurrent time-integration scheme. I've managed to show that the ...
0
votes
2answers
73 views

Finite element programming tutorial [closed]

I am new to this site and it might therefore not be the best place for this question. I have been using finite element programs (CFD mainly) for some time and I want to learn more about the basics. ...
6
votes
1answer
138 views

How to project a vector into the H(div) space (in the context of finite elements)?

Say I have a simple elliptic PDE: $$ -\nabla\cdot(K\nabla p) = f \;\;\;\text{in}\;\Omega $$ with the appropriate boundary conditions. I solve for $p$ using a FEM (a discontinuous Galerkin method to ...
6
votes
1answer
73 views

Space-time finite element discretization for time-dependent PDEs

In FEM literature, semi-variational methods are typically used in the solution of time-dependent PDEs. I have not seen a fully-variational approach i.e. where space and time are discretised by FEM, ...
1
vote
0answers
73 views

Determining Youngs Modulus of defined material in FEA (ABAQUS)

I am quite new to FEA but need to determine if I am performing my simulations correctly. I have a cube of material with defined elastic properties (Youngs Modulus and poison ratio). I perform a ...
6
votes
1answer
108 views

How to avoid negative values of numerical solution of transport equation using FEM scheme?

The transport equation is actually an advection-diffussion-reaction equation, which has the form as $$\frac{\partial C}{\partial t} + v_1 \frac{\partial C}{\partial x} + v_2 \frac{\partial ...
1
vote
1answer
51 views

Reaction-Diffusion problem A->B, solving for B

I need to solve a Reaction-Diffusion using Finite Elements, piecewise linear elements. In this problem, a reaction $A \rightarrow B$, with rate law $ r_A = - k_A \cdot u_A $, takes part, where $u_i$ ...
1
vote
0answers
42 views

Can variational formulations be solved using series solutions?

What I specifically mean is, given some functional $F\left[\mathbf{x}\right]$ which is stationary with respect to $\dot{\mathbf{x}}=f(\mathbf{x})$ and some boundary or initial conditions, can one ...
3
votes
1answer
64 views

Solving pure Neumann problem enforcing B.C. with Lagrange Multiplier

I want to solve the Laplace Equation with pure Neumann B.C. using Finite Element Method: $- \Delta u = f \ $ in $ \ \Omega $ $- \partial u/\partial n = g \ $ on $ \ \Gamma = \partial \Omega$ With ...
2
votes
1answer
82 views

Best preconditioner for mixed-poisson problem (RT0 elements)

For a very large mixed-poisson problem with lowest order Raviart-Thomas elements (RT0), I plan on using an iterative solver. However, this kind of problem is not positive-definite (saddle point ...
5
votes
2answers
143 views

How do hexahedral FEM meshes improve approximation quality per degree of freedom, compared to tetrahredal meshes?

From the deal.II FAQ : ...quadrilaterals and hexahedra typically provide a significantly better approximation quality than triangular meshes with the same number of degrees of freedom; you ...
1
vote
1answer
140 views

Step-wise finite element formulations: can this be done?

Given the functional: $$ F[\mathbf{x}]=\frac{1}{2}[\mathbf{x}^{\text{T}} * D(\mathbf{x})]-\frac{1}{2}[\mathbf{x}^{\text{T}} * \mathbf{Ax}]-\frac{1}{2}\mathbf{x}^{\text{T}}(0)\mathbf{x}(t) $$ Where ...
9
votes
3answers
165 views

Finite elements on manifold

I'd like to solve some PDEs on manifolds, say for example an elliptic equation on a sphere. Where do I start? I'd like to find something that use preexisting code/libraries in 2d , nothing so fancy ...
4
votes
1answer
167 views

Why are functional representations of systems important in numerical applications?

I tried asking a similar question in SE.Physics, and I got some information regarding the abstract side of this, but I figured I should post here to get more complete information about the numerical ...
3
votes
0answers
50 views

Creating FEM mesh for image region — what is the most suitable shape function?

I wish to create a FEM mesh to solve an inverse elasticity problem, for an irregular domain. This domain is given by a medical image, so it is discretised and each square on the grid has one scalar ...
0
votes
0answers
29 views

Resources for viscous behavior in simple FEM

I am working on a simple explicit-integration lumped-mass elastic FEM code which implements CST+DKT triangles (plate+shell) and constant-strain tetrahedra ...
3
votes
2answers
140 views

Why do structured and unstructured discretizations give different errors?

It is necessary for me to solve a Poisson problem with a numerical method on a square domain with two types of triangular mesh: uniform triangular mesh (using uniform distribution nodes on square) and ...
9
votes
2answers
106 views

Finite elements $W^{1,\infty}$ error estimates

Are there finite element method setups that provide error estimates in the $W^{1,\infty}$ norm (i.e., bounds on $\|u'_h - u'\|_\infty$)? Which families of elements can be used for implementing them? ...
0
votes
1answer
56 views

Solving Initial Value problem ignoring the time-derivative

I am looking at a heat initial value problem \begin{align} \frac{\partial u}{\partial t}-\nabla^2u = f\quad&\text{in}\quad \Omega\times(0,T)\\ u = g \quad&\text{on}\quad ...
3
votes
2answers
144 views

Initial Value Problem using Finite Element

I am trying to implement a FEM solver for the following initial value problem \begin{align} \frac{\partial u}{\partial t} - \nabla^2 u &= f\quad \text{ in } \Omega\times (0,T)\\ u &= g\quad ...
1
vote
0answers
73 views

Thin plate stiffness: analytical formula to validate FEM model

I tried to compute analitically the stiffness of a cantilever thin plate (shown in picture). The plate is also homogeneous and isotropic. The aim is to compare the result I obtain with the result I ...
0
votes
1answer
160 views

How to solve Energy Balance equation by numerical method

Good Day I am new to heat transfer technique please give me some suggestion on solving energy balance equation $$a \frac{\partial T_p}{\partial t}=\frac{\partial}{\partial x}\left(b\frac{\partial ...
2
votes
2answers
99 views

Model of heat sink problem with fan

I am trying to solve this problem using advection-diffusion model and finite element method for the solution, due to the complex geometry. Basically the problem i'm trying to solve using OpenFOAM is ...
5
votes
1answer
104 views

Effect of subdomain topologies on overlapping additive Schwarz?

Is there a reference on the effect of subdomain topology on performance of the overlapping additive Schwarz method for (high order) finite elements? For example, taking subdomains to be vertex ...
3
votes
2answers
185 views

Projecting Finite Element solution onto new mesh

I am implementing a finite element solver in MATLAB and I have the following problem. Let's say I have a mesh $\mathcal{T}_1$ with triangular elements on a rectangular domain ...
0
votes
1answer
99 views

Importing results of FEM analysis into Matlab

I need to import in Matlab the results (like time histories of diplacement or frequency response at a specific point) obtained from a FEM analysis in Nastran. At the moment I ask Nastran to save the ...
4
votes
2answers
227 views

P versus Q elements

I am currently developing a project that uses finite elements for multi-dimensional PDEs and I'm still wondering if I will use P elements (triangles in 2D and tetra in 3D) or Q elements (squares in 2D ...
0
votes
2answers
85 views

Structural FEM analysis: transiet response vs frequency response

I am running 2 simulations on a cantilever plate in Nastran: one is a transient analysis (time domain) and the other one is a frequency response analysis. The transient analysis computes the response ...