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
15 votes
1 answer
2k views

Visualizing discontinuous Galerkin/finite element data

I would like to visualize simulation results, obtained using the discontinuous Galerkin (DG) approach, within ParaView. Similarly to finite volume methods, the problem domain is divided into cube-...
Michael Schlottke-Lakemper's user avatar
10 votes
1 answer
3k views

PDE discretization with the method of rothe and the method of lines (Modular implementation)

The Heat equation is discretized in space with FV (or FEM), and a semi-discrete equation is obtained (system of ODEs). This approach, known as the method of lines, allows to easily switch from one ...
gnzlbg's user avatar
  • 1,075
15 votes
3 answers
6k views

Finding which triangles points are in

Suppose I have a 2D mesh consisting of nonoverlapping triangles $\{T_k\}_{k=1}^N$, and a set of points $\{p_i\}_{i=1}^M \subset \cup_{k=1}^N T_K$. What is the best way to determine which triangle each ...
Nick Alger's user avatar
  • 3,143
3 votes
1 answer
2k views

Rigid Body Elements

I am currently developing structural FEM solver in FORTRAN. My question is about Rigid Body Elements (Multi Point Constraints). In NASTRAN there is RBE2 element defined by one independent and one or ...
Vladimir Smajevic's user avatar
2 votes
0 answers
3k views

COMSOL - Implementing Perfectly Matched Layers (3D)

If there are any COMSOL gurus here, I want to simulate a 3D RF system infinite in 1 direction. I realize I need to use PML's, but am unsure how to configure them. it seems that simply scaling them by ...
Hopeful's user avatar
  • 21
2 votes
0 answers
78 views

Complexity of direct solvers? [duplicate]

Possible Duplicate: How to reorder variables to produce a banded matrix of minimum bandwidth? What is the time and space complexity of direct sparse solvers (e.g., UMFPACK, SUPERLU, PARDISO, etc.)...
Thomas W.'s user avatar
  • 810
30 votes
9 answers
2k views

What is a good way to run parameter studies in C++

The problem I'm currently working on a Finite Element Navier Stokes simulation and I would like to investigate the effects of a variety of parameters. Some parameters are specified in an input file ...
Matija Kecman's user avatar
8 votes
5 answers
745 views

Recommended Route for Mastering Inverse PDE Problems

I would like to master Inverse PDE Problems particularly with the use of Finite Elements. My problem is I don't know where to start. Should I begin by reading a book on Inverse Problems or on PDE-...
Jason's user avatar
  • 81
9 votes
3 answers
2k views

Construction of $C^1$/$H^2$-conforming finite element basis for triangular or tetrahedral mesh

In the paper Hierarchical Conforming Finite Element Methods for the Biharmonic Equation, P. Oswald claimed Clough-Tocher type elements has $C^1$-continuity while being a cubic polynomial on each ...
Shuhao Cao's user avatar
  • 2,552
7 votes
1 answer
302 views

Forced viscous damping in elastodynamics

I have an 2D elastodynamics problem, that is a problem which is driven by the Cauchy equation: $$\rho\ddot u-\mathrm{div}\sigma=\rho f$$ where $u$ is the displacement, $\sigma$ the Cauchy stress ...
Benjamin's user avatar
  • 265
11 votes
2 answers
7k views

FEM: singularity of the stiffness matrix

I'm solving the differential equation $$ \left( \sigma^{2}(x) u ''(x) \right)'' = f(x), \;\;\; 0 \leqslant x \leqslant 1 $$ with initial conditions $u(0) = u(1) = 0$, $u''(...
Appliqué's user avatar
  • 445
5 votes
2 answers
561 views

What does symmetrize mean? (imposing multifreedom constraints to stiffness matrix)

I have a small FEM implementation program. And I want to add imposing multifreedom constraints (MFC) feature to it. The theory of master-slave method is given here (page 10 for general case). ...
danny_23's user avatar
  • 501
10 votes
2 answers
372 views

What about this simple error estimate for linear PDE?

Let $\Omega$ be a convex polygonally bounded Lipschitz domain in $\mathbb R^2$, let $f \in L^2(\Omega)$. Then the solution of the Dirichlet problem $\Delta u = f$ in $\Omega$, $\operatorname{trace} u ...
shuhalo's user avatar
  • 3,620
8 votes
2 answers
150 views

Initial guesses for perturbed linear systems

Suppose you solve a linear system $Au = f$ by an iterative method, e.g. conjugate gradients or Richardson iteration. Then you try to solve a linear system that is slightly perturbed in the matrix and ...
shuhalo's user avatar
  • 3,620
9 votes
3 answers
2k views

Standard format for finite element meshes

Does there exist a standard format for finite element meshes which is widely used in the industry? Thanks!
Benjamin's user avatar
  • 265
33 votes
9 answers
6k views

Modern resources for learning FEM

I need to get started using Finite Element Methods. I am about to start reading Numerical solutions of partial differential equations by the finite element method by Claes Johnson, but it's dated 1987....
a06e's user avatar
  • 1,729
25 votes
4 answers
4k views

How to incorporate the boundary conditions with the Galerkin method?

I've been reading some resources on the web about Galerkin methods to solve PDEs, but I'm not clear about something. The following is my own account of what I have understood. Consider the following ...
a06e's user avatar
  • 1,729
4 votes
0 answers
355 views

How does one handle the source term in the Shallow Water Equations when using the discontinuous galerkin method? [closed]

I use the discontinuous galerkin method to solve the steady flow 1D shallow water equations with a bump at the bottom. This flow is frictionless. I use the runge-kutta method to approximate the time ...
Thida's user avatar
  • 57
7 votes
1 answer
184 views

Conjugate Gradient with Hierarchical Basis Functions: How can the hierarchical base be decomposed?

I'm trying to implement a Conjugate Gradient solver using Hierarchical Basis Functions, following this paper. In section 3 the paper says that the hierarchical basis matrix $S$ can be decomposed into ...
rsp1984's user avatar
  • 435
8 votes
2 answers
2k views

Shape regularity in higher dimensions

In Finite Element theory, and other methods in scientific computing for PDEs, one uses meshes which fulfill several regularity criteria, many of them being equivalent. It is of interest to have ...
shuhalo's user avatar
  • 3,620
7 votes
2 answers
2k views

Reference implementation of Nédélec-Elements

Does anybody know of an implementation of Nédélec elements that does not come along with a huge bulk of additional software? Is there a small library written in a language like Python, Matlab, or ...
shuhalo's user avatar
  • 3,620
6 votes
2 answers
2k views

Which libraries have good implementations of Basis splines?

I'm looking to use the finite element method with B-splines as my function basis. Which C/C++ libraries have good B-spline support? Specifically, I'm looking for an implementation of a stable ...
Dan's user avatar
  • 3,355
19 votes
1 answer
6k views

How to derive the Weak Formulation of a Partial Differential Equation for Finite Element Method?

I have taken a basic introduction to Finite Element Method, which did not emphasize a sophisticated understanding of a 'weak formulation'. I understand that with the galerkin method, we multiply both ...
Paul's user avatar
  • 12k
11 votes
1 answer
337 views

What are the possible numerical schemes for a diffusion equation with a nonlinear reaction term?

For some simple convex domain $\Omega$ in 2D, we have some $u(x)$ satisfying the following equation: $$ -\mathrm{div}(A\nabla u)+cu^n = f $$ with certain Dirichlet and/or Neumann boundary conditions....
Shuhao Cao's user avatar
  • 2,552
10 votes
1 answer
254 views

Quadrature rules, methodologies, and references

There is at least one quite comprehensive encyclopaedia of quadrature rules that doesn't seem to have been updated in quite a while and has restricted access. This source refers to several classical ...
6 votes
0 answers
580 views

How to do FEM in sector elements?

Suppose we have an element in the polar coordinate as $0\leq r\leq 1$, $0\leq\theta\leq \alpha$. What are the correct basis functions in this element? For this kind of mesh, there are a lot of ...
Hui Zhang's user avatar
  • 1,319
10 votes
2 answers
2k views

Are 8 Gauss points required for second order hexahedral finite elements?

Is it possible to get second order accuracy for hexahedral finite elements with fewer than 8 Gauss points without introducing unphysical modes? A single central Gauss point introduces an unphysical ...
Geoffrey Irving's user avatar
5 votes
2 answers
2k views

Solution oscillations with a small timestep in backward Euler

I am using backward Euler in a FEM scheme for a convection-diffusion problem. On a given mesh, I can take arbitrarily large time steps, as expected. But if I decrease time step, at some point it will ...
Dominik Lark's user avatar
3 votes
2 answers
1k views

Looking for a library or algorithms to perfom clipping 3D unstructured meshes by a set of surfaces

We have a 3D (volume) unstructured, possibly hybrid, degenerative irregular mesh data structure that we are capable of generating (mostly composed of hexahedra and general polyhedra, using a mix of ...
Marc-Olivier Titeux's user avatar
9 votes
1 answer
713 views

What numerical quadrature to choose to integrate a function with singularities?

For example, I would like to numerically compute the $L^2$-norm of $\displaystyle u = \frac{1}{(x^2+y^2+z^2)^{1/3}}$ in some domain that includes zero, I tried Gauss quadrature and it fails, it is ...
Shuhao Cao's user avatar
  • 2,552
8 votes
2 answers
242 views

Is there one general approach to build a projection methods for different problems?

My question is probably going to be too general to answer it with a couple words. Could you please suggest a good reading in that case. Projection methods are used to reduce size of the solution space ...
danny_23's user avatar
  • 501
53 votes
5 answers
11k views

What are criteria to choose between finite-differences and finite-elements

I am used to thinking of finite-differences as a special case of finite-elements, on a very constrained grid. So what are the conditions on how to choose between Finite Difference Method (FDM) and ...
shuhalo's user avatar
  • 3,620
5 votes
5 answers
1k views

Introductions to hp-FEM

do you know good introductions into or surveys $hp$-adaptive finite elements? In particular I do not know how the heuristics for choosing spatial refinement or increased polynomial degree are ...
shuhalo's user avatar
  • 3,620
13 votes
1 answer
2k views

What are possible methods to solve compressible Euler equations

I would like to write my own solver for compressible Euler equations, and most importantly I want it to work robustly in all situations. I would like it to be FE based (DG is ok). What are the ...
Ondřej Čertík's user avatar
9 votes
2 answers
778 views

Coupling FEM DG methods to Riemann solvers

Are there any good papers and or codes that couple discontinuous galerkin finite element solvers with Riemann solvers? I need to explore coupling elliptic and hyperbolic problems but most splitting ...
aterrel's user avatar
  • 3,624
14 votes
1 answer
460 views

What spatial discretizations work for incompressible flow with anisotropic boundary meshes?

High Reynolds number flows produce very thin boundary layers. If wall resolution is used in Large Eddy Simulation, the aspect ratio may be on the order of $10^6$. Many methods become unstable in this ...
Jed Brown's user avatar
  • 25.6k

1
21 22 23 24
25