The tag has no usage guidance.

learn more… | top users | synonyms

5
votes
1answer
67 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 ...
1
vote
1answer
53 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
0answers
39 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
1answer
44 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 $(...
6
votes
2answers
70 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{...
4
votes
1answer
41 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 $...
3
votes
1answer
64 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 ...
1
vote
0answers
20 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 ...
2
votes
1answer
43 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\...
3
votes
1answer
81 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 governign PDEs are: \begin{equation} \left. \begin{aligned} C_1\frac{\...
0
votes
1answer
109 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 ...
1
vote
0answers
68 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 ...
3
votes
1answer
182 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 ...
3
votes
1answer
91 views

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

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. Dirichlet boundaries ...
1
vote
2answers
84 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 ...
4
votes
1answer
86 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)$. ...
0
votes
0answers
58 views

Injection Vs Full Restriction in Dirichlet-Neumann 3-D Multigrid

I have implemented the Multigrid method for a Mixed Dirichlet-Neumann boundary value problem where $\nabla^{2}{u}=0$, $u = 1+x+y+z$ for Dirichlet and $\frac{\partial e}{\partial n} = 1$ for Neumann ...
1
vote
0answers
46 views

Monotonic convergence of Newton's method for boundary value problems

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
1answer
76 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 ...
9
votes
1answer
815 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 ...
1
vote
2answers
105 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 ...
2
votes
1answer
58 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 ...
0
votes
1answer
56 views

Equivalence of linear elasticity and biharmonic equations: variational formulation

Wikipedia tells me that the equations for linear elasticity and biharmonic equations have the same solution for Dirichlet boundary condition. How do you show the equivalence in the variational ...
0
votes
1answer
82 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
169 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
116 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
0answers
46 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 (...
2
votes
2answers
79 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 ...
3
votes
1answer
85 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
175 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)...
0
votes
1answer
117 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 ...
2
votes
1answer
192 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 ...
4
votes
0answers
79 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 $-\...
3
votes
2answers
185 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 ...
3
votes
1answer
59 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
322 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: ...
5
votes
2answers
312 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 ...
7
votes
2answers
507 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 ...
17
votes
3answers
2k 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 (...
2
votes
1answer
119 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 \...
5
votes
2answers
1k 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 ...
4
votes
3answers
803 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 ...
8
votes
1answer
212 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)...
3
votes
0answers
266 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}$ (...
3
votes
1answer
583 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 ...
7
votes
2answers
404 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 ...