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.
895
questions
3
votes
1
answer
176
views
How can we calculate mixed derivatives numerically using the Chebyshev derivative matrix?
Using the Chebyshev derivative matrix $D$, we can numerically approximate the first and second derivative of a function by doing matrix multiplication:
$${df(x) \over dx} = Df(x) \tag 1$$
$${d^2f(x) \...
1
vote
0
answers
103
views
Method of lines for a mixed PDE
I am trying to solve the following PDE using the method of lines to discretize space, and then solve it as system of ODEs at each point in space using ODE15s:
subject to
and initial condition $w(z,t=...
1
vote
2
answers
196
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 ...
0
votes
1
answer
80
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, ...
2
votes
1
answer
145
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 ...
2
votes
1
answer
99
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 ...
0
votes
1
answer
56
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
votes
2
answers
279
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 ...
1
vote
0
answers
278
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 ...
2
votes
0
answers
108
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}...
0
votes
0
answers
53
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 ...
0
votes
0
answers
81
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....
0
votes
0
answers
43
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}}...
11
votes
1
answer
1k
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 ...
2
votes
2
answers
1k
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 ...