Questions tagged [elliptic-pde]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
1
vote
2answers
102 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
2answers
92 views

Simulating the heat equation with insulating material

My plan is to solve the heat equation in the right half portion of the domain, while having the left half completely isolated with constant temperature. To do so, I model the left half with a very low ...
1
vote
1answer
88 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 \...
1
vote
1answer
237 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
2answers
169 views

Specific questions for 2-D Multigrid

I am simulating $\nabla^{2}u=0$ with mixed Dirichlet-Neumann boundary conditions on 2-D using 2-Grid method. Dirichlet ...
1
vote
1answer
238 views

An example of mixed elliptic problem using lowest-order Raviart Thomas element

I try to solve the following mixed second order elliptic PDE in the domain $D=[0, 1]^2$ \begin{eqnarray*} v+\nabla p=&0 \quad &\text{in} \quad D,\\ \text{div}(v)=&1/2 \quad &\...
1
vote
1answer
116 views

A doubt in Multigrid V-cycle

Assume I have 3 levels of grids. Finest Grid = level 2, Coarser Grid = level 1, Coarsest Grid = level 0. Relax $u$ on $Au = b$ at level 2 for 3 times. Find residual $r2$ at level 2, then restrict to ...
1
vote
0answers
45 views

When to discretize nonlinear Poisson Equation

I am trying to solve a nonlinear poisson equation of the form: $u_{xx} + f(u_y)u_{yy} = 0$. In trying to get a handle on this problem, it seems like there are two approaches. I could either ...
1
vote
0answers
73 views

Determine truncation error of PDE discretization

The equation is $$\frac{\partial}{\partial x}\left(u\frac{\partial u}{\partial x}\right)=f(x)\\ 0<x<1, u(0)=u(1)=0$$ I'm discretizing this PDE using FVM as follows: $0=x_0=x_{1/2}<x_1<x_{...
1
vote
0answers
102 views

PDE discretization on triangular domain

Given the 2D Poisson equation $$\Delta u = f\\ u(x,0) = g_1(x), 0<x<1\\u(0,y) = g_2(y), 0<y<1\\ \partial_n u (x, 1-x) =0, 0<x<1$$ defined on the domain $\Omega := \{(x,y) \in \...
1
vote
0answers
61 views

scaling in discretized PDE system

I want to solve the following system via Matlab $\Omega=(0,1)^2$ $$\Delta y=\frac{1}{\alpha} p$$ $$ -\Delta p= y -1 $$ $$p|_{\partial \Omega}=0,~y|_{\partial \Omega}=0$$ using ...
1
vote
0answers
179 views

FiPy with derivative source terms

I have a coupled nonlinear PDE system in 1 spatial dimension, in which I want to solve using FiPy. The dependent variables are $n$ and $T$: \begin{align} \frac{\partial n}{\partial t} \,&=\, D\,\...
1
vote
0answers
99 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
0answers
53 views

Decreasing - increasing - stabilising $l_{2}$ norm

Let $\bar{x}$ denote the analytical solution of a PDE. Let $x^{(k)}$ be the solution at the $k^{th}$ iteration. The initial guess for the solution is $x^{(0)} = 0$. Let $r_{0} = ||\bar{x}-x^{(0)}||_{2}...
1
vote
1answer
94 views

3 questions on FEM to solve elliptic PDE with homogeneous and mixed boundary conditions

Assume we use FEM with piecewise linear finite elements to discretize the BVP over $\omega = (0,1)$: $-u''+ bu' + u = 2x$, $u(0) = u(1) = 0$ for parameter $b\in R$. Given a mesh $T = \left\{x_i\right\}...
1
vote
0answers
73 views

Manufacturing a solution for non-smooth coefficients in elliptic problems

This question is a continuation of this answer (I can't comment) If we were going to manufacture a solution for a problem with discontinuous coefficients, I understand that the solution should have ...
1
vote
0answers
32 views

Gradient convergence on a checkerboard domain in finite element

I've got an elliptic problem in a domain with two materials in a checkerboard arrangement. Both materials are highly dissimilar, being one much more compliant than the other. If I use classic linear ...
1
vote
0answers
118 views

understanding interior penalty jump of basis function

If cardinal basis functions are used (i.e. $\psi_{ij}=1$ iff $i=j$, and 0 otherwise) in interior penalty methods for elliptic equations such as SIPG & NIPG, shoudn't the jump in basis functions ...
1
vote
0answers
88 views

Monotonic convergence of Newton's method for boundary value problems [closed]

I’m interested in solving nonlinear elliptic boundary value problems of the type $$ -a\Delta u + f(u) = 0, $$ $$ u|_\Gamma = u_0 $$ by Newton’s method when its convergence is global and monotonic. ...
1
vote
0answers
507 views

Finite Difference in Polar Co-ordinates

Background Please note that I am duplicating the question on scicomp. I have already asked this in math. I am trying to come up with a scheme in Polar Co-ordinates for the following PDE: PDE I am ...
0
votes
2answers
548 views

Correctly setting boundary condition for periodic linear elasticity problem

From an old, wise engineering book Peterson's Stress Concentration Factors (http://eu.wiley.com/WileyCDA/WileyTitle/productCd-0470048247.html page 324) I've got the following problem: There is 2D ...
0
votes
1answer
130 views

Calculation of error

I have written a code in which I find the approximation of the solution of this elliptic problem. I calculated the error using the following part of code: http://pastebin.com/7b5mmuRW but I get the ...
0
votes
1answer
123 views

Imposing zero mean condition in FEM

I wanted to solve a periodic elliptic equation of the form $$-\nabla\cdot(A\nabla u)=-\nabla\cdot F$$ on $Y=[0,2\pi)^d$ using FreeFem++, where $A$ and $F$ are $Y$-periodic. The space of solutions is $...
0
votes
1answer
202 views

jump conditions for Poisson/Darcy equation in primal form versus mixed form

Consider the Darcy equation, $$\mathbf{v} + \dfrac{k}{\mu_0}\nabla p = \mathbf{f} \\ \mathrm{div}\; \mathbf{v} = 0$$ If the coefficient $k$ is piecewise constant across an interface $\Gamma$ in the ...
0
votes
1answer
100 views

Elliptic PDE: Proving that a 2nd order accurate discretization of the 2nd derivative of the unknown is 2nd order accurate for the unknown itself

For an ODE: $\frac{dy}{dt}=f(y(t),t)$ The Euler Explicit scheme reads: $y_{n+1}=y_{n}+\Delta tf_n$ and it can be easily shown with a Taylor expension that: $y_{n+1}=y_{n}+\Delta t \frac{dy}{dt}|...
0
votes
1answer
251 views

Discontinuous Galerkin - Inhomogeneous Dirichlet B.C. for 1D Poisson Equation

I am trying to get some code working for the 1D Poisson equation using the textbook: Nodal Discontinuous Galerkin Methods Algorithms, Analysis, and Applications. I use the following formulation (for ...
0
votes
1answer
191 views

Bifurcation of linear PDE

I have a linear elliptic PDE (unfortunately not allowed to be shown here) with a constant parameter $\epsilon$ giving the stable solutions qualitatively shown by the functions below. As we smoothly ...
0
votes
1answer
158 views

Simulation of Laplace Equation in 3-D with mixed BC of Dirichlet-Neumann

I am simulating a Mixed-Boundary value (Dirichlet-Neumann) problem using Finite Difference Method on a unit 3-D cube such that the left, lower, and front plane have $u=u(x,y,z)=1$ (Dirichlet) and ...
0
votes
1answer
666 views

How to impose Neumann boundary conditions in interior penalty DG method

Consider the following two point BVP: $$ -u''(x)=f(x),~~~u(0)=u(1)=0. $$ An interior penalty DG method for this BVP that weakly imposes homogeneous Dirichlet boundary conditions is of the form: $$ \...
0
votes
0answers
57 views

Difficulty to prove the coercivity of $a(u, v)=\int_{\Omega} \nabla u \cdot \nabla v$

I'm stuck at proving the coercivity of the bilinear functional in the variational formulation of the problem: \begin{array}{c} -\nabla^{2} u=f \quad \text { in } \Omega \\ u=g_{D} \text { on } \...
0
votes
1answer
139 views

Singularity in gradient caused by Dirichlet boundary condition

I am looking for a mathematical explanation for the singularity caused by a Dirichlet boundary condition partially imposed at a boundary. For instance $$ \nabla^2u=0 ~ \text{in}~\Omega $$ where $\...
0
votes
0answers
48 views

Existence and uniquness of solution of FVM for Poisson equation

I'm discretizing the following Poisson equation using FVM where the domain $\Omega$ of the solution is a regular hexagon of side $1$ centered about the origin. $$\Delta u =k,\text{ $k$ constant}\\ \...
0
votes
1answer
199 views

Elliptic PDE finite volume method with Dirichlet boundary condition

I want to discretize the following equation using a Finite Volume Method $$\nabla \cdot (a(x)\nabla u)=f(x)\\x\in \Omega \subset \mathbb{R}^2 \\u_{|\partial\Omega}=g$$ I'm using Voronoi cells here: ...
0
votes
0answers
80 views

Finite element method for Surface integrals using polar coordinates

I am trying to solve a 2D elliptic PDE (see complete electrode model for electrical impedance tomography) using the finite element method (FEM) over a circular region $\Omega$. I have discretized the ...
0
votes
0answers
224 views

How to form the stiffness matrix for the Poisson equation using a spectral method

This is a follow-up question to How do I form the Chebyshev differentiation matrix in MATLAB? The goal of the following code is to solve the Poisson problem: ...
0
votes
0answers
52 views

What is a reasonable manufactured solution to test finite difference method? [duplicate]

What is a reasonable manufactured solution to test the following equation against its finite difference approximation? I want it to look like a Cosine function about $0$, rotated about $Z$ axis (...

1
2