Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [elliptic-pde]

The tag has no usage guidance.

22
votes
3answers
3k views

What is the purpose of using integration by parts in deriving a weak form for FEM discretization?

When going from the strong form of a PDE to the FEM form it seems one should always do this by first stating the variational form. To do this you multiply the strong form by an element in some (...
16
votes
1answer
5k views

What is the general idea of Nitsche's method in numerical analysis?

I know that the Nitsche's method is a very attractive methods since it allows to take into account Dirichlet type boundary conditions or contact with friction boundary conditions in a weak way without ...
10
votes
1answer
319 views

$L^2$-convergence of finite element method when right hand side is only in $H^{-1}$ (Poisson eqn)

I know that the piecewise linear finite element approximation $u_h$ of $$ \Delta u(x)=f(x)\quad\text{in }U\\ u(x)=0\quad\text{on }\partial U $$ satisfies $$ \|u-u_h\|_{H^1_0(U)}\leq Ch\|f\|_{L^2(U)} $...
8
votes
2answers
534 views

Helmholtz and Biharmonic equation examples with exact solution

I'm looking for examples of Helmholtz and Biharmonic equations in Cartesian co-ordinates with exact solutions, in order to compare my numerical solutions with it. I was able to find quite a few ...
8
votes
1answer
274 views

Asymptotic convergence of the solution to a parabolic pde to the solution of an elliptic pde

Suppose I have the parabolic system $$u_t=\nabla\cdot(k(x)\nabla u)+f,\quad (x,t)\in\Omega\times I$$ with Dirichlet boundary conditions $$u=g, \quad x\in\partial\Omega$$ and initial condition $$u(x,t)...
8
votes
0answers
424 views

Numerical implementation of the Dirichlet-to-Neumann map

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
318 views

Increasing V-cycles for constant Coarsest Grid Size and increasing Fine Grid size

Problem statement I implemented geometric multigrid for $-\nabla^{2}=f$ where $f=\frac{3\pi^{2}}{4}sin \frac{\pi x}{2} sin \frac{\pi y}{2} sin \frac{\pi z}{2}$ on $\Omega \in [0,1]$ on a unit cube. ...
7
votes
3answers
3k views

Role of boundary conditions (e.g. periodic) in Poisson equation

Given 3D Poisson equation $$ \nabla^2 \phi(x, y, z) = f(x, y, z) $$ and the right hand side and the domain, am I free to impose any boundary conditions (BC) on the function $\phi$, or do they have to ...
7
votes
2answers
708 views

Is there a good tutorial or textbook-like source on implementing ENO/WENO with limiters in one (and more than one) dimension?

I've inherited a finite volume code that does a second-order discretization of flux terms for a set of mixed parabolic-elliptic equations with discontinuous diffusion coefficients. The impression I ...
6
votes
2answers
117 views

Stokes Equation in “two-fold saddle point” form?

Are there papers that deal with the (nondimensionalized) Stokes equation for incompressible fluid flow in a "doubly mixed" form like the following? \begin{align*} 0&=\underline{\epsilon} + \frac{...
5
votes
4answers
778 views

biharmonic equation

I want to solve the biharmonic equation numerically, that is: $$\Delta^2 u=f~~in~~\Omega$$ $$u=g_1~~on ~~\partial \Omega$$ $$\frac {\partial u}{\partial n}=g_2~~on ~~\partial \Omega$$ Using Green's ...
5
votes
1answer
88 views

Finite-volume method applied to a particular advection equation

I'm trying to apply the finite-volume method (FVM), with which I'm not so familiar, so a simple 1D PDE equation. I already asked this question in the math stackexchange, but was told that it could be ...
5
votes
2answers
355 views

Finite element discretization of Reaction-diffusion problem with Dirac source term

I'm writing a code using continuous piecewise linear finite elements on triangular grids to solve the diffusion-reaction problem. the source function f is a Dirac mass at the center. How can i compute ...
5
votes
1answer
66 views

Stabilization parameter for an elliptic equation

I simply want to solve the elliptic equation: $$ -\kappa \nabla^2 u + u = f $$ where $f\in [0,1]$. When using continuous Galerkin with Lagrange elements, I have noticed that $\kappa$ has to be greater ...
5
votes
3answers
253 views

How to impose boundary condition with mixed derivatives?

I have the biharmonic equation on a 2D rectangular domain $\Omega$ with the following boundary conditions: $\Delta^2 u = f$ on $\Omega$ $\nabla u \bullet \mathbf{n}=0$ on $\partial \Omega$ (1) $u_{...
5
votes
0answers
128 views

Elliptical problem with Robin BC

Working in finite differences, I am using a transformation on the temperature variable $\Theta = \int_{T0}^T \kappa(T)dT$ to linearize the steady-state heat equation into a Poisson equation $-\...
4
votes
3answers
1k views

Online Poisson Solver

I'm wondering if anyone can point to a browser-based FEM (or other) 2D PDE solver for simple elliptic problems. It seem like there ought to be a javacript implementation that would allow for the ...
4
votes
1answer
1k views

variational formulation of linear elasticity

First I'm not 100% sure I'm on the good stack for asking my question. I would like to get a bilinear form for linear elasticity that separate a rotational part from a pure divergence part, so starting ...
4
votes
1answer
163 views

Finite element error for second order ODE at nodes equal to zero

I coded a finite element method with linear basis elements for the problem $$-u'' = f(x), x\in[0,1], u(0) = u(1) = 0$$ The nodes are uniformly spaced and I will denote them as $x_i$. I initially ...
4
votes
1answer
128 views

Comments needed on the doubts of PDEs in moving boundary problems

We know that in classical two-phase Stefan problems, let's say in the temperature distribution of ice-water problem here, the governing PDEs are: \begin{equation} \left. \begin{aligned} C_1\frac{\...
4
votes
1answer
246 views

Points on the interface

We consider the problem $\left\{\begin{matrix} k(x)\Delta u(x)=f(x) & \text{ in } \Omega\\ u=0 & \text{ in } \Gamma \end{matrix}\right.$ where $\Omega \subset \mathbb{R}^2$ open and ...
4
votes
1answer
83 views

In mixed elliptic formulation, what are the weakest requirements to ensure the flux is in $H^1$?

In the book Mixed Finite Element Methods and Applications by Boffi, Brezzi, and Fortin there is a pretty long discussion about why the Raviart-Thomas (RT) projection is only defined for functions in $...
4
votes
1answer
112 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 ...
4
votes
1answer
123 views

Computing Fourier representation of space dependent advection operator via FFT

Consider the following equation on the circle: $$\dfrac{\partial p(x,t)}{\partial t} = a(x)\dfrac{\partial p(x,t)}{\partial x} \equiv L(p) \enspace ,$$ where $L$ is the operator acting on $p(x,t)$. ...
4
votes
0answers
273 views

Neumann-Neumann boundary intersection

I am trying to solve a Vertex Centered, Finite Difference, Poisson equation $-\nabla^{2}u=f$ with Dirichlet boundaries on the left and bottom and Neumann boundaries on the top and right using ...
3
votes
1answer
206 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 ...
3
votes
1answer
256 views

How do I avoid divide-by-zero when solving the Poisson equation with Fourier transforms?

I wanted to try to implement part of the method in the following article using Fourier transforms. http://www.shodor.org/media/content/jocse/student_submissions/nocito2010/nocito2010_pdf Right now I ...
3
votes
1answer
407 views

Jump condition for elliptic equation in standard finite element method

Consider the elliptic PDE $ -\mathrm{div}(k(x)\nabla u) = f(x) \in \Omega$ and $u = 0$ on $\partial \Omega$ If $k(x)$ is piecewise constant across an interface $\Gamma$ in the domain, then we ...
3
votes
1answer
100 views

Fast methods to solve an elliptic PDE if high accuracy is needed only in part of the domain

Does someone know a method to get cheap approximation of harmonic problems (and possibly local approximations)? Let me explain: I need to compute the solution of an harmonic problem \begin{equation} ...
3
votes
1answer
254 views

steady state solution from parabolic problem vs solution of elliptic problem

My question is related, but not a duplicate of Asymptotic convergence of the solution to a parabolic pde to the solution of an elliptic pde Suppose I solve the parabolic PDE: $u_t = \Delta u + f(x,t)...
3
votes
1answer
74 views

Library to compute eigenvalues of the Laplace operator in a polyhedral domain

What library can one use to compute efficiently the lowest eigenvalues of the Laplace operator in a polyhedral domain in $R^3$? For the application I have in mind one has to consider very acute ...
3
votes
0answers
88 views

$L^2$ norm error estimates of conforming FEM about Poisson’s equation with mixed boundary conditions

Consider Poisson’s equation $$- \Delta u = f{\rm\qquad{ in }}\;\Omega $$ with following mixed boundary cconditions $$u = g{\rm\qquad{ on }}\;\Gamma \subset \partial \Omega $$ $$\frac{{\...
3
votes
0answers
92 views

Unstable convergence of a Poisson equation

What could be the reason that the solution of a Poisson equation is smooth when obtained by an iterative solver, only if the maximum residual is set to a high value (e.g. 0.1)? When the maximum ...
3
votes
0answers
58 views

Solving numerically a linearized system of elliptic (?) Navier-Stokes equation (Shallow Water Derived)

For my PhD Thesis, my advisor asked me to build a solver inspired from the article "Optimal Control Theory Applied to an Objective Analysis of a Tidal Current Mapping by HF Radar, J-L Devenon, 1989". ...
3
votes
1answer
298 views

General case Kutta condition

I'm working on a 2D inviscid fluid simulation using a "panel method", with Potential being used to enforce the no-through boundary condition. I'm trying to incorporate the Kutta condition, which says ...
3
votes
0answers
319 views

Solving PDE or eigenvalue problems without FEM

Do you know any methods/solvers for PDE or eigenvalue problems like $\begin{cases} \Delta u= 0\ (\text{ or }\lambda u) & \text{ in }\Omega \\ u =0 & \text{ on }\partial \Omega \end{cases}$ (...
2
votes
2answers
215 views

$O(h^2)$ convergence for Elliptic PDE

I am trying to solve an elliptic PDE in 2-D: $$-\nabla^{2} u = 20tanh(10x-5)(10-10tanh(10x-5)^2) = f$$ I know that the solution is $u = tanh(10x-5)$ but I am unable to get $O(h^2)$ solution with a ...
2
votes
2answers
97 views

Infinite Function Value on Dirichlet Boundary

I have been working on a multigrid solution to a non-homogeneous Dirichlet boundary value problem. However, the function goes to infinity on the boundary. This causes numerical overflow errors to be ...
2
votes
2answers
81 views

Elliptic equation with finite volume and unstructured high order geometry

I have found that in unstructured mesh, discretizing the laplacian operator with finite volumes requires special care, as given in An Introduction to Computational Fluid Dynamics: The Finite Volume ...
2
votes
1answer
98 views

Boundary elements method — calculation of solid angle

I am developing a BEM code based on a deal.ii tutorial, see https://www.dealii.org/8.3.0/doxygen/deal.II/step_34.html . Consider the Poisson equation $$ \Delta u=-f\,, $$ and its Green's function $G\...
2
votes
1answer
121 views

Finding the frequencies of vibration of a circular and square drum

I want to find the frequencies of vibration of a circular and square drum. To do this, I need to solve a 2-dimensional wave equation (PDE) with boundary conditions. Every method that I have researched ...
2
votes
1answer
144 views

Solving the elliptic eigenproblem with periodic boundary conditions

Given an energy level $\mu$, I'm looking to calculate the eigenvectors corresponding to the time-independent Schrodinger operator on the torus (that is, periodic boundary conditions)-- $H = -h^2 \...
2
votes
0answers
68 views

Are there known accuracy issues between 2D axisymmetric and 3D solutions?

In my full 3D solutions I am solving for the potential throughout a $100\times 200\times 200$ grid. Inside is a ring electrode set to -5V via a Dirichlet boundary condition, and surrounded on all ...
2
votes
0answers
95 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 ...
2
votes
1answer
333 views

FreeFem user-defined function [closed]

I'm trying to solve the steady-state heat equation with a space-dependent thermal diffusivity in FreeFem. I.e. $$\nabla\cdot(\kappa(x,y)\nabla T) = 0$$ I'm hoping to read in from a file locations $(...
2
votes
0answers
388 views

assembly matrices in finite element method

I'm trying to construct the right–hand side of my 2D Poisson's equation in Matlab. I used the vertex rule in order to approximate the integral: ...
1
vote
1answer
179 views

How to correctly define the flux in a finite volume method for Poisson's equation with a piecewise constant material

How do we correctly define the flux in a finite volume method applied to Poisson's equation where we have a piecewise constant material? Specifically, say we have the equation \begin{align*} -\nabla\...
1
vote
2answers
94 views

Building minimization optimization problem for 2nd-order elliptic PDE

I am solving elliptic PDE problem, for which, Euler scheme looks as following: $$ \nabla [\gamma ( |\nabla u|^2) \nabla u] = 0,$$ where $$\gamma(|\nabla u|^2) = (1 + |\nabla u|^2)^{-1/2}. $$ I am ...
1
vote
1answer
169 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 &= ...
1
vote
1answer
75 views

Does a mixed method solve this elliptic pde exactly if the source function is piecewise polynomial?

Let $\Omega\subset \mathbb{R}^d$, $d\in \{2,3\}$ be an open bounded polygonal/polyhedral set. Suppose I want to solve the following pde \begin{align*} \vec{q}+\vec{\nabla}u &=0\,&x&\in \...