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

3
votes
2answers
123 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
54 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
78 views

Partial differential equations with octave [on hold]

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
82 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
71 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
133 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
111 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
200 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
73 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
190 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
131 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}) + ...
7
votes
1answer
120 views

Finite element convergence rates for mixed problems

I've coded up a Stokes Flow problem using finite elements and am in the process of verifying that it works. I'm just not sure what convergence rate I should be expecting as I globally refine the mesh. ...
4
votes
3answers
149 views

vector PDEs on manifolds

What are the subtleties involved in solving vector PDEs on manifolds? Can someone suggest a reference summarizing the problems involved? Specifically I want to solve a vector Helmholtz equation with ...
1
vote
2answers
149 views

wavelet for numerical partial differential equations

Is there a good introduction into wavelet Galerkin schemes for numerical partial (and ordinary) differential equations?
1
vote
2answers
127 views

How to compute $\Delta u$ on the boundary of the biharmonic equation?

Let $u$ be the answer of a PDE.Is there any relationship between $u,\frac{\partial u} {\partial n}$ and $\Delta u$. I have the values of $u$ and $\frac{\partial u} {\partial n}$ on $\partial \Omega$ ...
2
votes
1answer
96 views

Time Integration of a nonlinear reaction-diffusion system

I want to solve the following system of nonlinear reaction-diffusion equations (Schnakenberg Turing) using FEM methods (such as deal.ii): $$ \partial_{t} u = \Delta u + \gamma\left(a-u+u²v\right)$$ ...
0
votes
0answers
28 views

Second order PDE with minimum condition

I need to solve 2nd order PDE of $F=F(X,Y)$, where X and Y are functions of t. I have got to point where I transformed my PDE to canonical with logarithmical transformation for X and Y to get constant ...
2
votes
0answers
33 views

How should I choose the knot sequence when using B-splines as a basis for solving a PDE?

I'm looking to solve the Schrödinger equation with a basis made of a tensor product of basis splines. A number of papers describe calculations made with a program designed this way, but they never ...
0
votes
1answer
47 views

Guidelines for choosing manufactured solutions for numerical PDE schemes

When testing a numerical method for a PDE, I know that it's often useful to compare it to a known analytical solution. If none is available, one can always 'manufacture' a solution, substitute it ...
9
votes
1answer
420 views

What are the relative benefits of using Adams-Moulton over Adams-Bashforth algorithm?

I am solving a system of two coupled PDE's in two spatial dimensions and in time computationally. Since the function evaluations are expensive, I would like to use a multistep method (initialised ...
3
votes
2answers
125 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 ...
1
vote
2answers
61 views

Performance metrics to compare initial-boundary value problem solutions

I am comparing the performance several finite difference methods of solving an initial-boundary value problem. There are several dimensions to this comparison: Number of cells Number of timesteps ...
2
votes
0answers
105 views

Integration of nonlinear PIDE via spectral methods

At the mean-field level, the dynamics of a polariton condensate can be described by a type of nonlinear Schrodinger equation (Gross-Pitaevskii-type), for a classical (complex-number) wavefunction ...
4
votes
1answer
255 views

Coupled nonlinear PDEs with time dependence on the RHS

I would like to numerically solve the following system of 2 coupled partial differential equations for the unknown functions $\psi_X(x,y,t)$ and $\psi_C(x,y,t)$: $\partial_t \psi_X = -i\psi_C - ...
5
votes
1answer
97 views

Solving the quadratic in the Fast Marching Method

The Fast Marching Method is a way of solving the Eikonal Equation on a discrete grid, essentially just computing a wavefront speading out from initial points, IE: The idea is that we want to ...
3
votes
1answer
107 views

Solving a system of nonlinear PDEs by minimization

I have two coupled nonlinear partial differential equations of the form: $ \begin{align} \dot{u} -f(u,u',u'',v,v',v'')=0 \\ \dot{v} -g(u,u',u'',v,v',v'')=0 \end{align} $ The boundary conditions are ...
3
votes
0answers
140 views

A Question About Weak Forms in Fenics

Is it possible to use test and trial functions from two different function spaces (defined over two different meshes) in a single weak form? Under what conditions can I do this (eg., each term in the ...
1
vote
3answers
157 views

Surface normals integration

I am trying to reconstruct a 3D surface given the normals of the unknown surface. Reading through this paper on section 4 they say [...] denote the surface by $z(x,y)$. The directions of the ...
1
vote
1answer
79 views

Numerically solving the polar poisson equation [closed]

I want to solve the Poisson equation for a 2D polar system: $$\Delta_r f(r) = u(r)$$ with the Laplace operator: $ \Delta_r f(r) = \frac{1}{r}\frac{\partial}{\partial r} \left[r ...
2
votes
1answer
61 views

Imposing symmetry plane boundary condtition

I want to impose symmetry plane boundary condition for a solid mechanics problem. I googled around and found out that in many places people say to "forbid displacemnts out of symmentry plane and ...
2
votes
0answers
94 views

Segregated solving of a tightly coupled system of PDEs

To compute the evolution of a free surface between two incompressible, immiscible liquids, two tightly coupled equations have to be solved, the volume fraction advection and the Navier-Stokes ...
5
votes
0answers
96 views

Strategy for solving a non-trivial differential equation

I would like to numerically solve an equation of a type as shown below. Does anyone of you have an idea how to approach such a problem? Any links to literature or for further reading would be greatly ...
1
vote
1answer
226 views

Parallel 2d heat equation (implicit timestepping) using MPI

I am trying to solve the time dependent heat equation with backward euler timestepping and second order space finite differences. This results in a Poisson system needing to be inverted. In serial ...
2
votes
1answer
116 views

Von Neumann stability analysis in 3d

I need to get a stability criterion for the numerical scheme for equation $$\frac{\partial u}{\partial t}-\frac{\partial^2 u}{\partial x^2}-\frac{\partial^2 u}{\partial y^2}-\frac{\partial^2 ...
2
votes
1answer
116 views

Courant Friedrichs Lewy condition - how to get it?

I am interested, how can we get CFL condition for every type of PDE? It's known that for 1st order linear equation $$\frac{\partial u}{\partial t}+a\frac{\partial u}{\partial x}=0 $$ CFL is get from ...
1
vote
0answers
127 views

Derivation of a Higher Order Compact Alternating Direction Implicit Method

I dont understand how this Higher Order Compact ADI scheme, which is fourth order in time and space, for the wave equation is derived: I go through the following Using Tylor's expansion $u(t+h,x,y)$ ...
1
vote
0answers
61 views

Locally conservative method for differential generalized eigenvalue problem

I have to approximate the smallest eigenvalue of the following generalized eigenvalue problem $$ - \nabla \cdot D(x) \nabla p(x) + \alpha(x) p(x) = \lambda \beta(x) p(x) $$ over a domain like ...
4
votes
2answers
109 views

Discretization method for a reaction dominated elliptic PDE

I'm working with an elliptic reaction diffusion PDE of the form $$-k\nabla^2u+cu=0$$ I've noticed that when the reaction term dominates over the diffusion (i.e. $c>>k$), the true (exact) ...
7
votes
2answers
154 views

The effect of decoupling a coupled system of PDEs

I asked a somewhat similar question previously but perhaps it might have been too specific for anyone to really answer. Here is a bit more general of a question that I am struggling with. Consider the ...
4
votes
1answer
147 views

User Friendly GUI based 3d PDE solving software

I'm searching for software which can solve 3d IBVP like PDEtool in Matlab (user-friendly GUI, easy to learn how to use). What could you advise to me? (ANSYS,COMSOL,OpenFOAM are too difficult for me)
6
votes
2answers
160 views

Spectral Methods in time

I was reading up on Spectral Methods for PDEs. In all the descriptions I read, while the position component is approximated via a Fourier series or other methods, the time component is still ...
6
votes
0answers
113 views

Sequential approach to solving coupled PDEs

I'm dealing with a coupled system of three transient, non-linear convection-diffusion equations. Let's just say to simplify the problem that they take the following form: $$ ...
2
votes
1answer
113 views

One finite difference scheme

There is PDE: $$\frac{\partial u(r,\varphi,\psi,t)}{\partial t}=\operatorname{div}A(r,\varphi,\psi)\nabla u +f(r,\varphi,\psi,t) $$ We solve numerically IBVP for the ball $B_{1}(0)\subset ...
1
vote
0answers
57 views

Stability question (finite difference): dealing with corner nodes

Consider one initial boundary value problem for sphere. $$\frac{\partial u}{\partial t}=\operatorname{div}A\nabla u +f$$ Here is explicit numerical scheme (we consider that it is stable): ...
2
votes
1answer
202 views

Implementation of gradient zero boundary conditon in advection-diffusion equation

My question is about Finite Element Method. I want to know how to implement "gradient zero" conditions to advection-diffusion equations in conservative form like, $\frac{\partial \rho}{\partial t} + ...
1
vote
0answers
41 views

how to find outward normal for robin codition [duplicate]

I wrote the code to fix this problem; now I should validate it using a test function to see the error in my code. I can't figure out how starting from u can derive the function g on the edge, in ...
-1
votes
1answer
162 views

Pde problem with robin boundary condition

I have my pde 2D problem with robin condition (form: du/dn +ku=g) to solve with matlab. i have the exact function u and I want to find the function g in robin condition. How can i do it? thanks for ...
2
votes
1answer
84 views

variational formulation

I would like to minimize: $$J = \int_{\Omega} \|\nabla u - \nabla g\|^2 + \lambda \|\frac{\partial u}{\partial t} + \nabla u.v||^2 ~\text{dx dy dt}$$ where $u(x,y,t)$ is the unknown function, the ...
0
votes
1answer
127 views

Calculating Divergence in COMSOL

Is it computationally safe and accurate to use the following equation in COMSOL to compute the divergence of the vector quantity J (instead of using its general built-in equations that have $\nabla$ ...
9
votes
3answers
372 views

How to deal with curved boundary condition when using finite difference method?

I'm trying to learn about numerically solving PDE by myself. I've been beginning with finite difference method(FDM) for some time because I heard that FDM is the fundament of numerous numerical ...