Questions tagged [pde]
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.
866
questions
1
vote
2
answers
150
views
Modeling contamination diffusion in a draining container, part 2
Part 1, but I'll repeat here. This time we'll move the top boundary.
I have a container of water with initial volume $V_0$ and height $L$, open at the top which receives a constant mass flux $\phi$. A ...
- 197
0
votes
1
answer
69
views
Determining the importance of different parameters in a simulation
Suppose that I have a function of, say, three parameters, $f(p_1,p_2,p_3)$ whose output is a field(s) (e.g. velocity field) and is dependent on some real-valued parameters (e.g. viscosity, density, ...
- 623
2
votes
1
answer
135
views
Modeling contamination diffusion in a draining container
I asked this on the Physics site, but it fits much better here, and I now can ask a more specific question.
For the scope of the problem, I have a 250-mL bottle filled with pure water and a sampling ...
- 197
2
votes
0
answers
54
views
Symmetrization of Laplacian Matrix Operator (finite volumes)
The aim is to construct symmetric Laplacian matrix $L$ for 2D square domain.
I have the following discretization of second derivatives on non-uniform grid (will skip the steps of derivation, any ...
- 364
0
votes
1
answer
53
views
Time discretisation after splitting a 4th order equation
Suppose we have a fourth-order parabolic PDE
$$\partial_t u + a(t)\Delta( \Delta u) - div( b(x,t)\nabla u) =0$$
We split the equation into two second-order equations by introducing $w = \Delta u$ thus
...
- 3
3
votes
2
answers
130
views
Solving systems of advection-diffusion-reaction equations with finite element methods
I have been doing a lot of self-study on numerically solving PDEs so that I can solve system of linear and nonlinear Advection-Diffusion-Reaction (ADR) systems on complex meshes.
I have been watching ...
- 287
1
vote
0
answers
67
views
Using solve_ivp for a PDE: how to handle multiple time-dependent variables?
I am trying to build a Python code that solves a set of coupled differential equations which will be spatially discretized by the method of lines advancing in time. I am planning to use ...
- 11
2
votes
0
answers
73
views
How can we symbolically working out $\phi^4$ theory green's function/propagator and consequences in python?
I am having some difficulty calculating Green's function symbolically in Python for $\phi^4$ theory.
The specific rendition of the $\phi^4$ theory I have in mind can be written as follows.
$\mathcal{L}...
- 121
0
votes
0
answers
44
views
PETSC: Solving a simpler PDE results in error while solving the original equation works in snes/tutorials/ex13.c
In snes/tutorials/ex13.c,
there is a function SetupPrimalProblem(),
which sets up the $f_0$ and $f_1$ in ...
- 101
0
votes
0
answers
64
views
How to get damping matrix for structural model in FE analysis
I need to implement in C a method of obtaining transient solution of damped FE models based on modal results for a structural model (imported CAD geometry) defined with hysteretic (structural) damping....
- 1
0
votes
0
answers
35
views
Verification of a Function Definition in Python
I want to write a function $f$ and it is defined as $f = - \nabla \cdot(|\nabla u|^{p-2} \nabla u) $ and I exact solution $u(x) = \tilde{u}(r) = 1 - \frac{p-1}{p-2} \left( s^{p/{p-1}} - (1-s)^{p/{p-1}}...
- 143
11
votes
1
answer
980
views
Is using iterative methods to solve a linear system always superior to inversing the matrix?
I have a silly question. Is it always more computationally efficient to use iterative methods to solve for some matrix $A$, $Ax=b$, where $x$ and $b$ change but $A$ stays constant, compared to ...
- 121
2
votes
2
answers
347
views
Problems solving 2D heat equation using physics-informed neural networks
I am trying to solve 2D heat equation using the physics-informed neural networks approach. The training loss is decreasing, but my final network outputs make no sense. I am using Python/Pytorch.
2D ...
2
votes
1
answer
111
views
How to numerically solve differential equations involving sines, cosines and inverses of the unknown function?
I'm very new to finite difference method and I am just introduced to methods of solving differential equation using finite difference method via sparse matrix method.
I find that the main idea is to ...
- 35
0
votes
0
answers
41
views
A question from boundedness property of LMT
I have the following PDE
$$
u_t = k \Delta u + \alpha u H(u-c)
$$
I am trying to show the boundedness property to apply LMT. I get confused with the estimates of third term.
The space I have
$$V = \{...
- 143
3
votes
1
answer
156
views
Numerical scheme for the level set equation that solves inverse mean curvature flow problems
I am considering the problem of simulating the evolution of an interface given as a curve in 2D (or surface in 3D) that evolves according to a velocity specified at the interface of the form:
$$\vec{v}...
- 33
3
votes
0
answers
129
views
Discontinuous Galerkin for transport equation with non-constant advection
This question is mainly an inquiry about the usefulness of Discontinuous Galerkin (DG) for the time-independent transport equation of the form
$$\sigma u+\beta\cdot\nabla u =f,\;\;\;\text{on }\Omega\...
- 131
2
votes
1
answer
77
views
Differential Equation with Forced Behavior
I'm attempting to solve a strange differential equation problem. My goal is to know if there are kinds of ODE solver packages to solve this kind of problem.
I'm solving a 1D Partial Differential ...
- 297
0
votes
0
answers
53
views
Solving Laplace for Velocity Potential in Constricted Channel
I am trying to solve the 2D Laplace equation numerically to give the velocity potential of a fluid flowing in a channel with a constriction:
$u_{xx}+u_{yy}=0$
There is a constriction in the channel at ...
3
votes
0
answers
135
views
Form of nonlinear diffusion equation
Consider the following nonlinear diffusion problem,
$$
\frac{\partial u}{\partial t} = x^{-2}\frac{\partial}{\partial x}\left(x^2 u^4 \frac{\partial u}{\partial x}\right), \quad 0 < x < 1
$$
We ...
- 454
0
votes
0
answers
63
views
Analytical Equation of the gaussian 1D wave equation with periodic Boundary condition
I am trying to validate the 1D analytical wave equation with a numerical solution with periodic boundary conditions. I have implemented the periodic boundary condition for the numerically calculated ...
- 121
0
votes
1
answer
114
views
How do I use pdepe for a first order parabolic PDE with only one boundary condition?
I am trying to use Matlab's pdepe.m to solve the first order parabolic PDE
$$\frac{\partial u}{\partial x}+\frac{\partial u}{\partial x}=x$$
I have not had trouble coding the argument of pdepe @pdefun:...
6
votes
2
answers
273
views
Continuous vs discontinuous space-time FEM
What are some reasons for approximating a (e.g. parabolic) PDE using the space-time method, with continuous finite elements in time, vs discontinuous finite elements in time?
Are there e.g. ...
- 127
1
vote
0
answers
84
views
Accuracy of the Crank-Nicolson method for non-linear, inhomogeneous heat equation
I am currently coding a solution to the following PDE:
$\frac{\partial T }{\partial t} =\frac{\partial}{\partial \theta}(A(\theta ,\phi )\frac{\partial T }{\partial \theta}) +\frac{\partial }{\partial ...
- 11
0
votes
0
answers
105
views
2D wave equation is numerically unstable using Finite Difference Method
I'm working with simulating both the heat and wave equation in 2D in a Python code. When simulating the heat equation, I learned about the CFL which I used to get a numerical stable solution.
I found ...
- 101
0
votes
0
answers
91
views
Solving 2D Heat Equation with Input using central finite difference method
So I have implemented central finite difference method for solving the 2D heat equation. When I leave all initial and boundary conditions as 0s, but apply an input uniform across the entire space or ...
1
vote
1
answer
86
views
Simulating First Order Hyperbolic PDE with Finite Difference Scheme
I am trying to simulate a hyperbolic PDE with some control given by the following:
$$u_t(x, t) = u_x(x, t) + \theta(x) u(0, t)$$
with boundary conditions:
$$u(1, t) = U(t) = \int_0^1 k(1-y) u(y, t)dy$$...
- 31
2
votes
2
answers
216
views
Errors imposing boundary conditions weakly with DG
I am using interior penalty discontinuous Galerkin to solve a simple Laplace problem:
\begin{align*}
\nabla u=0
\end{align*}
with prescribed 0 and 1 Dirichlet boundary conditions on opposite edges of ...
- 395
1
vote
0
answers
133
views
Solving PDE on a non-uniform grid with Crank-Nicolson scheme
I am solving a 1D diffusion-type equation with the finite-difference Crank-Nicolson (CN) scheme, and I need to densify the spatial grid around the central point. One could change the spatial variable ...
- 21
0
votes
0
answers
80
views
Best approach to solve this system of equations?
I have the following 1D (in space, that is) system of equations I would like to solve:
\begin{equation}
\rho_{fs}\frac{\partial x_{fs}}{\partial t} = h_m\left(W_a - W_{fs}\right) - D_{eff}\left(\...
- 9
2
votes
1
answer
363
views
Meaning of Degree of Freedom in FEM
Assume we want to solve the Poisson eq. with the FEM on some Domain $\Omega$, i.e.
$$\begin{cases} -\Delta u = f, \; \Omega\\
u = 0, \; \partial \Omega \end{cases}$$
For the sake of the discussion let ...
- 123
1
vote
0
answers
75
views
Linear PDE solution with constraints
Consider the following linear PDE:
$$\nabla_q V(q) - M_d(q)M^{-1}(q)\nabla_q V_d(q) = 0,$$
where $V(q)$ and $M(q)$ are known and $M_d(q)$ is a grey box function (e.x., $M_d(q)$ is fitted using a ...
- 11
1
vote
0
answers
66
views
Discrete laplacian 9 point
I am trying to write a code for 9 point discrete laplacian. I would like to write a matrix and solve the linear system $AU=F$ with gradient conjugate method.
I wrote the matrix this way
...
- 41
3
votes
1
answer
89
views
Upwind scheme with periodic conditions
I am struggling with this assignment.
I have to write an upwind scheme for the following PDE:
$$u_t+a Du=0 \quad\mathrm{on}\;(-1,3)$$
$a$ is said to be positive, the initial condition is $\sin(2\pi x)$...
- 41
0
votes
0
answers
17
views
Kolmogorov $\epsilon$-entropy, $n$-width, and $\epsilon$-capacity and applications
What is the relationship between Kolmogorov $\epsilon$-entropy, Kolmogorov $n$-width, and Kolmogorov $\epsilon$-capacity of a set $M$ in a metric space $X$? (The $\epsilon$-capacity here is the ...
- 1
1
vote
0
answers
160
views
Solving second order coupled differential equations in python
as I have to design a reactor and therefore have to get its length x, I have to solve the following differential equations:
$$D_{eg}\tfrac{d^2A_g}{dx^2}-u_g\tfrac{dA_g}{dx} = k_la_b\left(\tfrac{A_g}{...
- 11
0
votes
0
answers
59
views
Can I use MOL to solve 2D steady state PDE in terms of r and z spatial coordinates?
recently I need to solve a 2D steady state PDE equation.
It’s not time dependent, and the only two independent variables are z and r direction.
So far for this solution, I was thinking using Method ...
- 1
1
vote
1
answer
30
views
Applying spectral method to a damped, driven 2D bending-mode wave equation on an irregular domain with heterogeneous boundary conditions
I'm trying to model some two-dimensional waves, and am unsure how to combine my boundary conditions with spectral methods. The PDE I'd like to explore resembles the equations for damped, driven ...
- 153
0
votes
1
answer
155
views
Can successive over-relaxation (SOR) method deal with ill-posed PDE BVP?
Recently, I've been struggling to understand the limitations and capabilities of the successive over-relaxation (SOR) method for boundary value problems which are ill-posed, such as, for instance, ...
- 53
2
votes
2
answers
170
views
Poisson equation with discontinuous variable coefficient
Trying to solve numerically the 2D Poisson equation with variable diffusion coefficient $K$, that in general can be discontinuous.
$$ \frac{\partial}{\partial x} ( K(x,y)\frac{\partial C}{\partial x}) ...
- 435
1
vote
0
answers
85
views
Crank-Nicolson vs Spectral Methods for the TDSE
The time-dependent Schroedinger equation (TDSE) depends linearly on the system's initial state $\vert \psi(0) \rangle$, such that the solution can be generally written as
$$ \vert \psi(t) \rangle = \...
- 111
1
vote
1
answer
199
views
Convergence stall when solving 2D Poisson PDE with pure Neumann boundaries (finite differences)
I recently started coding a small library of 2D PDE solvers (time dependent and time independent), and my first attempt was a 2D Poisson equation of the form:
$$\nabla(\epsilon\nabla\varphi)=\nabla\...
- 53
4
votes
0
answers
220
views
Spurious oscillations in solving diffusion problems using finite elements
I've been struggling with this problem for a while so I hope someone can help me here.
I'm trying to solve the McNabb-Foster equations for hydrogen diffusion in metals using a simple 1D finite element ...
- 41
1
vote
1
answer
99
views
One dimensional $C^1$ finite elements
I tried to solve a one dimensional biharmonic equation with finite elements. I wanted to use a conforming approach (as I simply do not know a lot about other approaches) and therefore was looking for ...
user43841
0
votes
1
answer
61
views
Numerical solution of PDE with uniform initial condition
I have a PDE like this
$$
\frac{\partial h}{\partial t} = \bigg(\frac{\dot{L}}{L}\bigg)x\frac{\partial h}{\partial x} - \alpha\bigg[h^3\frac{\partial^3 h}{\partial x^3}\bigg]
$$
With boundary and ...
- 63
2
votes
1
answer
945
views
How do you handle the singularity in polar or cylindrical coordinates?
Governing equations in polar or cylindrical coordinates often have terms with $\frac{1}{r}$ involved. At $r = 0$, such terms blow up to become a "singularity." The Cartesian version of such ...
- 29
3
votes
1
answer
146
views
orthogonal basis functions on arbitrary domains and boundary conditions
I'm interested in solving an inverse coefficient problem for a PDE.
Let's say the field to be estimated is $\theta$.
The conventional approach would be to use a finite element discretization for $\...
- 9,813
4
votes
2
answers
150
views
Numerical estimation of eigenfunctions of Laplacian
Consider the Laplace equation,
$$
\nabla^2 f(r,\theta,\phi) = 0
$$
in spherical coordinates. We know that the solution to this equation can be derived analytically, and is given by,
$$
f(r,\theta,\phi)...
- 1,038
2
votes
0
answers
52
views
Centered finite volume scheme for an advective term on an unstructured/irregular/non-uniform grid
Consider the continuity equation
$$\frac{\partial u}{\partial t} + \frac{\partial \Phi}{\partial x} = 0$$
$$\Phi = au + b\frac{\partial u}{\partial x}$$
Suppose I want to solve the above using ...
- 297
0
votes
1
answer
117
views
Solving PDE with a non-linear constraint in MATLAB
I am trying to solve a DAE with a non-linear constraint. The governing equations have the following form.
The second equation is a constraint and it must be satisfied everywhere. Is there a way to ...
