All Questions
Tagged with pde finite-volume
43
questions
0
votes
0
answers
56
views
What is the best finite volume method for the following equation?
I'm trying to create a partial differential equation that approximates 1-D climate in a rocky planet's atmosphere, which accounts for energy transport via radiation and convection. I am only ...
0
votes
1
answer
51
views
numerical schemes for 1D PDE: for smaller grid size there is an increased roundoff error, larger size more truncation, so sweet spot in between?
I had a discussion with a colleague today. He claimed that usually for a general numerical scheme for solving a general 1D PDE, for smaller grid size there is an increased roundoff error because of ...
4
votes
1
answer
148
views
Burger's equation (PDE) does not work with downwind difference?
I'm working on implementing the discretised Burger's equation. I am quite confused as to why it does not work when using a step-function and downwind difference formula. When using a step-function and ...
4
votes
0
answers
126
views
FVM for non-regular domain with triangular mesh
Setup
The 1D convection-diffusion equation is given by:
\begin{equation}\tag{1}
\frac{\partial u}{\partial t} + v \frac{\partial u}{\partial x} - \mu \frac{\partial^2 u}{\partial x^2} = 0,
\end{...
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
0
answers
64
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 ...
2
votes
0
answers
72
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 ...
0
votes
0
answers
58
views
Discretizing Multi-species Ion Exchange Equations by Finite Volume Method
I'm solving a system of multispecies ion exchange equations (diffusion+drift fluxes) in 1-d spherical domain using finite volume method to obtain the ion concentrations at the next time step. After ...
1
vote
0
answers
162
views
Two-dimensional modelling of a flate-plate reactor
I am trying to simulate the unsteady pollutant concentration along the reactor by solving the 2nd order PDE below with the stated BC and IC. which method is appropriate to solve convection-diffusion-...
0
votes
0
answers
76
views
The relation between PDE order and discretization order
In Jasak's Ph.D. thesis (2000), a notion is given about discretization of a transport equation:
For good accuracy, it is necessary for the order of the discretization to
be equal to or higher than the ...
1
vote
1
answer
365
views
How to compute turbulent energy cascade
I need to compute the kinetic energy cascade using a finite volume solution in an equally spaced grid. I wonder if it is more correct to first compute the kinetic energy in the space (or time) domain, ...
2
votes
1
answer
60
views
Discretization with non-constant matrix containg entries form unknown vector
Consider a system of PDEs
$$
\begin{cases}
u_t = \nabla \cdot (D(u)\nabla u) + \frac{c}{K_U+c}u-ku\\
c_t = d_c\Delta c -\frac{\nu_U c}{K_U + c}u
\end{cases}
$$
with some boundary conditions. Here, $D(...
1
vote
0
answers
66
views
Discretizing a parabolic PDE with finite volume method
I want to discretize the following parabolic PDE:
$$u_t = \nabla\cdot(\alpha(x)\nabla u)- \beta u\\
x\in\Omega \subset \mathbb{R}^2\\
\partial_n u = 0\\
u(t,0) = u_0(x)\ge 0, \alpha(x)>0$$
Given ...
1
vote
0
answers
106
views
Determine truncation error of PDE discretization
The equation is $$\frac{\partial}{\partial x}\left(u\frac{\partial u}{\partial x}\right)=f(x)\\ 0<x<1, u(0)=u(1)=0$$
I'm discretizing this PDE using FVM as follows:
$0=x_0=x_{1/2}<x_1<x_{...
0
votes
1
answer
348
views
Elliptic PDE finite volume method with Dirichlet boundary condition
I want to discretize the following equation using a Finite Volume Method
$$\nabla \cdot (a(x)\nabla u)=f(x)\\x\in \Omega \subset \mathbb{R}^2
\\u_{|\partial\Omega}=g$$
I'm using Voronoi cells here: ...