Partial differential equations (PDEs) are equations that relate the partial derivatives of a function of more than one variable. This tag is intended for questions on modeling phenomena with PDEs, solving PDEs, and other related aspects.

learn more… | top users | synonyms

1
vote
0answers
29 views

scipy.integrate.ode ignores boundary conditions

I am trying to solve the 1-dimensional diffusion problem numerically using method of lines: $$ \frac{\partial c}{\partial t} =D \frac{\partial^2 c}{\partial z^2},$$ where the right hand side is ...
2
votes
1answer
88 views

Computational methods for finding the energy eigenvalues of the time-independent Schrodinger equation with arbitrary potential

I have seen in some papers that the energy levels in some arbitrary potential are specified. How can one find the energy levels in such arbitrary potentials. For example, $V(x)=\sin^2(x/2)$ with ...
10
votes
1answer
104 views

Can an approximated Jacobian with finite differences cause instability in the Newton method?

I have implemented a backward-Euler solver in python 3 (using numpy). For my own convenience and as an exercise, I also wrote a small function that computes a finite difference approximation of the ...
2
votes
0answers
50 views

Numerical integration when solving PDE: Simpsons rule and high frequency noise

I am solving a PDE, and one of the intermediate steps is to numerically integrate a function over a compact interval. The function is represented on a linearly spaced grid. I am using Simpsons rule ...
0
votes
1answer
56 views

How can PDE matrices be identified?

I need to include experimental results for lots of PDE (partial differential equation) matrices in my research work. How can I identify PDE matrices? For example, matrices in the UFL Sparse Matrix ...
0
votes
0answers
50 views

The centered difference operator for fractional function

Recently, I come to a question about the 2nd order centered finite difference approximation of a fractional function, more precisely, we set $\delta^{2}_x u(x,t) = u(x_{i-1},t)-2u(x_i,t) + ...
4
votes
1answer
61 views

How can I compare errors in PDE solvers with non-uniform grids?

Is there a standard approach to testing codes with refined regions? Specifically, I am interested in testing whether the refinement is working correctly. For the sake of simplicity, let's consider a ...
0
votes
1answer
87 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: $$ ...
1
vote
0answers
87 views

Tips on improving stability in numerical scheme for non-linear PDE

I am solving a non-linear second order system of PDEs in two variables. The equations are too complicated to write out here, but an essential feature is that there is a propagating wave which then ...
6
votes
3answers
89 views

Stability criterion for waves in anisotropic solids

The equations of motion for an elastic solid are given by $$\begin{align} &\nabla \cdot \boldsymbol{\sigma} + \mathbf{f} = \rho \ddot{\mathbf{u}}\\ &\boldsymbol{\sigma} = ...
0
votes
0answers
31 views

Matlab PDE toolbox

I'm new at matlab and I can't figure out how to make external temperature not constant but as a function. I wanted to use pde toolbox, I choose parabolic type. My external temperature should be ...
1
vote
1answer
66 views

Instability of pdepe in Matlab… boundary conditions?

here is a Matlab beginner banging his head on the wall... I am trying to solve a system of partial differential equations in Matlab, with both derivatives in time and space domains. I am using the ...
3
votes
1answer
133 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 ...
0
votes
1answer
81 views

Solving a PDE using Matlab (with varying initial conditions)

I want to solve a 1-D heat conduction PDE using Matlab which looks like $$ \rho c_p \dfrac{\partial T}{\partial t} = \dfrac{\partial}{\partial z}\left( \lambda \dfrac{\partial T}{\partial z} \right), ...
0
votes
2answers
84 views

Explain this multivariate differential identity

$$ \frac{\partial|\nabla\phi|^2}{\partial\phi}=-2\nabla\cdot\nabla\phi$$ I would very appreciate that you help me . Please do it in detail, I am quite not good at such problems. There is something ...
4
votes
3answers
146 views

Solving Laplace's equation on a domain with moving boundary

Consider a function $X(\xi,\nu)$, $2\pi$ periodic in $\xi$ satisfying $$\nabla^2 X = 0$$ in a domain $D$ with $\nabla = (\partial_{\xi},\partial_{\nu})$. If I know the values of $X$ on the boundary ...
1
vote
2answers
177 views

Python, numpy and complex functions (PDE's)

Update 4 I have almost given up on getting this right. This is the solution to the time-independent Schrodinger's equation, so the analytical solution is: $\psi(x,t) = \psi(x,0)e^{\frac{-iE ...
0
votes
1answer
77 views

Solving PDE with state and time dependent boundary conditions

I am interested in solving the following PDE (heat equation): $$\frac{\partial u}{\partial t} = \kappa \frac{\partial ^2 u}{\partial x^2}$$ In order to solve it, I discretize space uniformly into $N$ ...
0
votes
0answers
55 views

Solution to PDE with differential boundary conditions

I have the following equations $$ a_t(x,t)=1-a(x,t)b(x,t)^\gamma+D_1a_{xx}(x,t) $$ and $$ b_t(x,t)=\alpha(a(x,t)b(x,t)^\gamma -b(x,t))+D_2b(x,t)_xx $$ where $a,b:]0;4\pi[\times \mathrm{R}_+ ...
0
votes
0answers
76 views

parallel computatioan of a PDE in MATLAB

I want to solve a 1-D PDE $(\partial_{tt} + \alpha\partial_t)u(x,t)=\partial_{xx}u(x,t)-\sin(u(x,t))+f$, using method of lines and for this I defined a spatial grid of about n~1000 points. Since my ...
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, ...
8
votes
1answer
134 views

Can the method of lines be used to discretize all PDEs?

I have found that the method of lines is a very natural way to think about the discretization of PDE's. Therefore I always default to that mindset when presented with a new set of equations. I have ...
2
votes
0answers
92 views

1 D Diffusion equation FDM with different layers

I'm trying to solve this particular equation $\frac{\partial u}{\partial t} = \frac{\partial}{\partial x} \big[D_{i}(x)\frac{\partial u}{\partial x} \big] + S(x,t)$ where the $i$ index denotes ...
1
vote
1answer
50 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$ ...
0
votes
0answers
33 views

Isolated solving of coupled system of PDE [duplicate]

I would like to solve 3 differential equations for 3 unknowns. So I wrote a MATLAB code, which solves (using the '\' operator) these equation using a linear system of equations (in which the 3 ...
4
votes
2answers
129 views

Usability of upwind finite difference schemes

NOTE: I asked this on Mathematics Stack Exchange and there were no answers. So, I thought I might try here. Upwind schemes like the classic "upstream" scheme, can be used to solve, for example, the ...
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 ...
4
votes
0answers
89 views

numerical analysis of a partial integro-differential equation

I have to numerically solve a nonlinear partial integro-differential equation. This is my equation, $$\frac{\partial y(x,t)}{\partial t}=\int_{-\infty}^\infty K_0(|x-u|) \frac{\partial^2 ...
9
votes
3answers
164 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
166 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
2answers
133 views

Well-posedness of a linear elasticity problem and Navier-Cauchy equation

I read a master thesis on a topic I'm interested too. This work concern the solution of the displacement equation of motion for a homogeneous, elastic, isotropic material: $$\rho \ddot{\mathbf{u}} - ...
1
vote
1answer
89 views

What does “strongly conservative” mean in the context of numerical methods?

I have a homework problem that asks me to show that 1st order unwinding or central differencing can give a strongly conservative, consistent scheme for the 1-D Burger's Equation using a finite volume ...
3
votes
3answers
146 views

Algorithms for radiation treatment planning

I have a medical physics problem - I want to maximise the dose absorbed by a brain tumour whilst minimising the dose in the rest of the brain, especially certain organs, such as the pituitary gland, ...
4
votes
3answers
126 views

Variable viscosity Stokes equation

One very efficient way to solve Stokes equation with periodic boundary conditions \begin{equation} -\eta \nabla^{2} \bf{v} + \nabla p = f \\ \nabla \cdot \bf{v} = 0 \end{equation} is using the ...
1
vote
0answers
154 views

Stationary 2D/3D Navier-Stokes source code

Trying to solve stationary Navier-Stokes problem for incompressible laminar Newtonian fluid. I've found a couple solutions for instationary Navier-Stokes equations (like FeniCS examples or CFD ...
1
vote
1answer
89 views

Minimization of The Blind Deconvolution Functional

I want to minimize the functional of teh Blind Deconvolution model as given in: Total Variation Blind Deconvolution by Chan and Wong. Their model is given by: $$ z = h \ast u + \eta $$ Where $ \ast ...
3
votes
2answers
304 views

The rate of convergence for finite difference methods for Poisson's equation with piecewise constant data

I am solving the following PDE; $$ \nabla^2 u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = \rho, $$ where $\rho(0.5,0.5) = 2$ (zero elsewhere), $0\leq x,y\leq1$ and the ...
2
votes
2answers
169 views

FEM for a nonlinear parabolic PDE

I'm looking to numerically compute the solution to $$ k(x,u) \partial_t u - \Delta u = f \quad\quad\text{ in } \Omega \times [0,T]$$ where $k$ is a continuous but nonlinear (in $u$) real-valued ...
3
votes
2answers
232 views

Neumann Boundary Condition at r=0 in Polar Coordinates (Numerical BCs)

I have asked a question in this regard earlier. I am trying to solve the following equation in Polar Co-ordinates: $$ u_t - (u_{rr} + \frac{1}{r} * u_r + \frac{1}{\theta} * u_{\theta\theta} + bu) = ...
0
votes
0answers
86 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
44 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
vote
1answer
47 views

Application of CLAWPACK to Richards' equation

I'm looking to solve the Richards' equation. This models water flow in porous media and is a nonlinear, possibly degenerative, parabolic differential equation that takes the form $\partial_t ...
-1
votes
1answer
42 views

Can you give some information for rothe method [closed]

I want to learn a numerical method for PDEs other than finite difference method. After some research on internet i have found Rothe method and it looks interesting to me. Unfortunately, i couldn't ...
0
votes
2answers
104 views

Numerical Solution of non-linear diffusion equation using Finite Differencing

I'm trying to solve the following non-linear diffusion equation: $$ \frac{\partial}{\partial t} u(x,t)= \frac{\partial^{2}}{\partial x^{2}}u(x,t)^{3}$$ $$ -1\leq x \leq1, t \geq 0 $$ with the boundary ...
0
votes
1answer
50 views

How can we compute the exact value of some rotational symmetric n-by-n convolution kernel

The commonly used $3\times 3$ Laplacian convolution kernel $\frac{1}{6}\begin{bmatrix}1 & 4 & 1\\4 & -20 & 4\\1 & 4 & 1\end{bmatrix}$ is only an approximation of rotational ...
2
votes
2answers
97 views

Can the conservative form of the advection equation be re-written by replacing the velocity term with an integral over all other points in space?

Suppose I have a 1D advection equation in conservation (divergence) form $\partial_t u(x,t) = -\partial_x [v(x)u(x,t)],$ where $u$ is a conserved quantity in space, and $v$ gives the velocity of the ...
4
votes
1answer
129 views

Inverse advection-diffusion problem, solving for a drift coefficient with experimental data?

I am investigating a physical process where I believe the 1-D advection-diffusion equation: \begin{equation} \frac{\partial u}{\partial t} = -\frac{\partial}{\partial x}[\mu(x,t) u(x,t)] + ...
1
vote
0answers
64 views

Frozen coefficient method (von Neumann stability analysis)

Earlier it was considered that frozen coefficients method for Neumann stability analysis for finite difference scheme is more heuristic than rigorous. But I have read some information in a book by ...
1
vote
1answer
164 views

What kind of test cases are convenient to use for testing the code for Euler equations of gas dynamics in polar coordinates?

Consider Euler equation of gas dynamics in polar coordinates as $$ \left( \begin{array}{ccc} \rho \\ \rho u_{r} \\ \rho u_{\theta} \\ E \end{array} \right)_{t} + \frac{1}{r} \left( ...