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

9
votes
3answers
121 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 ...
3
votes
1answer
87 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
111 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
80 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
142 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
113 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
116 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
72 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
155 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
153 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
129 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
67 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
43 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
39 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
89 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
49 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
86 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 ...
3
votes
1answer
100 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
45 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
159 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( ...
1
vote
0answers
143 views

MATLAB: incorrect computation with pdepe for parabolic-elliptic system

I have found that MATLAB solves 1D parabolic-elliptic system incorrectly by using pdepe function. Here is a system: $$ u_t = u_{xx} + 2, $$ $$ 0 = v_{xx} + u. $$ Boundary conditions: $$ ...
5
votes
1answer
135 views

Numerical method for solving PDE with non-linear boundary conditions

I have the following problem: \begin{align} \frac{\partial w}{\partial t} = \frac{1}{r^2} \frac{\partial}{\partial r}\left( r^2 D_1 \frac{\partial w}{\partial r} \right) \\ \\ \frac{dr_d}{dt} = ...
0
votes
0answers
72 views

Eigenvalue analysis of preconditioned partial differential operator

today, I encountered a confused problem by accident, but I have no ideas to deal with it fully. The question can be described as follows: for example, when we need to use FDM/FEM to discrete the ...
2
votes
2answers
74 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
0answers
66 views

4th order tensor [closed]

I'm new with FEniCS and Python and I'm stuck with this issue: is there a way to write a 4th order tensor in an easy way to implement? I have to compute the following stiffnes tensor: $A_{ijkl}= ...
1
vote
1answer
56 views

What are the good testing problems for hyperbolic equation?

I read the whole list of this question: Where can one obtain good data sets/test problems for testing algorithms/routines? But the answers are in different areas and I want to ask a specific area. I ...
3
votes
2answers
73 views

Adaptive mesh refinement algorithms and the difference between AMR and moving mesh

I'm working on my thesis and a part of it has to do with adaptive mesh refinement. As a computer science major, I'm not too familiar with this field. The best way I can put my knowledge of AMR is: I ...
0
votes
0answers
44 views

Solvers for nonlinear parabolic PDEs

Could you please advise some programs or libraries for solving parabolic PDEs (or its systems) in 1D, 2D and 3D, for example, with the method of lines? The system of parabolic PDEs can be nonlinear in ...
2
votes
0answers
59 views

Solvers for stiff initial value ODEs with sparse Jacobian

What ODE solvers are optimized for solving stiff systems with sparse Jacobian? Such systems appear, for instance, when a parabolic PDE is discretized in space using typical finite difference or finite ...
1
vote
2answers
210 views

How can I solve this 1D nonlinear, variable-coefficient hyperbolic PDE?

I need to solve the following hyperbolic equation in x and phi co-ordinates $$\frac{\partial \left ( -s/f \right )}{\partial \varphi }+\frac{\partial \left ( 1/f \right )}{\partial x}=0$$ $$\varphi ...
4
votes
2answers
358 views

Crank-Nicolson method for solving nonlinear parabolic PDEs

Is the Crank-Nicolson method appropriate for solving a system of nonlinear parabolic PDEs like $\partial u/\partial t - a\Delta u + u^4 = 0$ ? I tried to apply this method for solving such system but ...
-1
votes
1answer
51 views

Libraries with the method of lines for parabolic PDEs [closed]

Could you please advise some programs or libraries for solving parabolic PDEs (or its systems) in 1D, 2D and 3D, for example, with the method of lines? The system of parabolic PDEs can be nonlinear in ...
3
votes
2answers
152 views

How to set up a shock tube problem such that the solution includes a shock with a specified Mach number

One of the famous and convenient test cases for shock wave modeling is the 1D Sod's shock tube. This is a Riemann problem for the compressible Euler equations of gas dynamics. The initial set up has ...
0
votes
2answers
171 views

Numerical method of lines for solving PDEs

Could you please advise some literature about the numerical method of lines (MOL) for parabolic PDEs? It is a method of solving PDEs with discretizing only by space but not by time. A system of ODEs ...
2
votes
1answer
186 views

Solving Advection (Convection) - Diffusion - Reaction Partial Differential Equation in Python

I am looking for library written in Python which will enable me to solve the coupled nonlinear equations which looks like: I need the library which will enable me to couple this solver to other ...
2
votes
1answer
119 views

What is the exact formulation of compressible Euler equation of gas dynamics in polar coordinates with artificial diffusion in 2D?

The interested equation is advection-diffusion equation. One of the canonical example is Navier-Stokes equations. However, I would like to let the coefficient of diffusion constant goes to zero, ...
1
vote
2answers
65 views

Looking for some literature about numerical solution of _coupled_ PDEs

There's a vast amount of literature about the numerical solution of partial differential equations. But in none of my books or lectures did I find anything about coupled systems. Internet and library ...
3
votes
1answer
111 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 + ...
5
votes
2answers
179 views

Implementation of nonlinear term in FEM

Although there are similar questions, I am also struggling with the implementation of the following term in "my own code" by Finite Element Method, namely, $\nabla \phi \cdot \nabla \phi$. $\phi$ is ...
3
votes
1answer
151 views

Numerically solve a PDE in Python with a term calculated by coarse-graining

I'm trying to solve a PDE in Python of the form, $\dfrac{\partial c(\mathbf{x}, t)}{\partial t} = \mathrm{D} \nabla^2 c(\mathbf{x}, t) -\gamma \rho(\mathbf{x}, t) c(\mathbf{x}, t)$ where $c$ ...
3
votes
1answer
159 views

Partial differential equations with octave [closed]

I need to find a numerical solution for $-\Delta U = f$, on the $\Omega = [0,1]^2$, with $ U|_{\partial \Omega} = 0$. I found a method: POISSONFD in ...
5
votes
1answer
181 views

Is there a jump condition for this PDE? ( Brinkman model , piecewise constant permeability)

The Brinkman equations for steady flow of an incompressible fluid through rigid porous solid are: $-\dfrac{\mu_0}{k}\mathbf{v} - \mathrm{grad}p + \mu_0 \Delta \mathbf{v} =0$ and ...
1
vote
1answer
119 views

How to solve coupled steady laminar diffusion flame jet problem? [closed]

I am trying to solve governing equations of laminar diffusion flame jet for steady state case. In the next step, I will solve for unsteady case. I have non-dimensional continuity, axial momentum, ...
2
votes
2answers
642 views

Visualising Maxwell's equations using MATLAB

I asked this question to help me understand what is going on in one of Maxwell's equations. I am happy with following through the maths on paper now, but would like to use MATLAB to take it one step ...
0
votes
0answers
168 views

Finite differences and Neumann boundary conditions

I am dealing with a highly nonlinear system of two PDEs. I already have a code to solve the system in case of Dirichlet boundary conditions. The explicit system is: $$ \begin{eqnarray*} \partial_{t}u ...
6
votes
2answers
281 views

Understanding the cost of adjoint method for pde-constrained optimization

I'm trying to understand how the adjoint-based optimization method works for a PDE constrained optimization. Particularly, I'm trying to understand why the adjoint method is more efficient for ...
1
vote
2answers
285 views

use Matlab's PDE toolbox to solve PDE with variable coefficients [closed]

I'm new to the PDE toolbox in Matlab. From the PDE specification window of the toolbox, it looks like one can only solve PDE with constant coefficients. How can ...
2
votes
3answers
246 views

Books on mathematical foundation of finite element methods

After reading three books about finite element method, with two of them covering also finite volume and grid generation, I found myself lost when I have to discuss these topics with library developers ...
3
votes
1answer
163 views

Finite Elements Weak Formulation generalization

I am struggling with an equation that represents the Weak form of Galerkin method: $ \phi^{T}F(\textbf{u})\sim \int_{\Omega}^{ } \phi.f_{0}(\mathit{u},\nabla \mathit{u}) + ...