Questions tagged [discretization]
The process of representing a continuum space with a finite set of points/elements
-1
votes
0answers
48 views
Finite difference method for conservative form of equations
My question is about how do we discretize the equations in the conservative form using finite difference method.
I'm trying to solve Euler equations in conservative form.
$$
\frac{\partial u}{\...
-1
votes
0answers
42 views
How to discretize continuity equation with velocity calculated using Darcy's law?
$$
\partial_t(\epsilon_g\rho_g)+\partial_x\cdot(\epsilon_g\rho_g\mathbf{v}_g)=\Pi
$$
I want to program normal continuity equation and Darcy's law to calculate velocity.
$$
\mathbf{v}_g=-\frac{1}{\...
0
votes
1answer
46 views
Changing the domain of a 3D Finite Difference code from cube to sphere
I have an explicit FD (Finite Difference) code for diffusion/heat on a PDE in a cuboid domain, and it works fine. I would like to update the discretized equations and change the code so as to solve ...
4
votes
1answer
96 views
Matrix Representation of a Discretization for a Partial Differential Equation
I want to discretize the following problem
\begin{cases}
\mu \nabla^2u+(\lambda+\mu)\nabla \nabla\cdot u = \rho \frac{\partial^2u }{\partial t^2 } + \beta \frac{\partial u}{\partial t}\\
u(...
6
votes
1answer
106 views
How do I integrate a function defined over an arbitrary area?
Let's say, I have a compact area $S$ (for example a circle, a square or some arbitrary polygon) and a function $f: S \rightarrow \mathbb{R}$. I want to numerically calculate the Integral
$$ \int_S f(\...
2
votes
1answer
75 views
Stability Analysis
The partial differential equation,
\begin{align}
\dfrac{\partial f}{\partial x} + a(x)\dfrac{\partial f}{\partial y} = 0
\qquad & f(0,y) = f(L_1,y) = c_0e^{-y} \\
& f(x,0) = c_0 \;,\; f(x,L_2) ...
5
votes
1answer
132 views
Finite Differencing schemes for Convection-Diffusion equation
I'm using the Convection(/advection)-Diffusion(-Reaction) equation to calculate the temperature over time in different hydraulic parts like a pipe or a heat exchanger.
The flow/convection is always 1D,...
2
votes
1answer
60 views
Is there Von Neumann stability analysis for 9-point laplacian like we have for the 5-point Laplacian?
For spatial accuracy in 2-D Laplace equation, a 9-point stencil is better than a 5-point one.
$$\partial_tq= r\left(\partial^2_x q + \partial^2_y q\right)$$
for FTCS (forward-time, central-space) ...
1
vote
1answer
49 views
How to isolate and test time discretization order of accuracy
I have a code that uses both spatial and time discretization/integration. For convergence analysis, I am wondering how one would test the order of accuracy of their ${time}$ integration scheme? I ...
0
votes
0answers
75 views
Strange behavior when modeling deflection (e.g., cantilever) problems using linear elasticiy
I am observing some strange behavior in my code as I model deflection such as cantilever problems. I am using the linear elasticity equations with finite volume methods.
I am attaching a figure at ...
0
votes
1answer
22 views
Calculating volume of a discretised diffuse interface object
Suppose I have a spherical object projected onto a discrete square mesh. The dicretised circle can be represented by filling a logical matrix such that voxels in the interior of the sphere are filled ...
2
votes
1answer
55 views
Discretisation of logarithmic derivative: Deriving the formula
I'm reading a paper where they use a discrete approximation of a logarithmic mass growth rate as follows:
$$ \frac{d \log M}{d \log t} \approx \frac{(t_B + t_A)(M_B - M_A)}{(t_B - t_A)(M_B + M_A)}$$
...
1
vote
0answers
45 views
Decrease in slope during convergence analysis
I am using the method of manufactured solutions to perform the order of accuracy testing. I am using a cube for the testing. The cube is size 1m on all sides.
I used 5 refinements:
$dx = dy = dz = ...
1
vote
0answers
36 views
Limit to volume change in a discretized mathematical model?
I have set up a mathematical model describing the diffusion of ozone out of a gas bubble. The bubble is surrounded by a thin gas film. So actually, the model describes the diffusion of ozone through ...
2
votes
1answer
71 views
FEM-Laplace with Dirichlet in only a few points: Nonsingular operator?
Let's consider the FEM discretization of the Laplace operator without boundary conditions, i.e.,
$$
a(u,v) = \int_\Omega \nabla u \cdot \nabla v - \int_\Gamma (n\cdot \nabla u) v.
$$
For one-...
5
votes
2answers
126 views
Residual norm of PDE discretization: correspondence in the continuous problem?
Solving a linear PDE like
$$
\Delta u = f \quad\text{on } \Omega,\\
n\cdot \nabla u = 0 \quad\text{on } \Gamma,
$$
with Finite Elements usually goes like this:
Create the discretization $Au=b$ via
$$
...
1
vote
0answers
56 views
Global truncation error behavior at fixed time step
I am trying to solve the following diffusion equation problem:
$\frac{\partial f}{\partial t}=\frac{\partial (D\frac{\partial f}{\partial x})}{\partial x}+S$
$D=1+x^{2}+\sin(x)$
$f(x,0)=1 , f(0,t)...
1
vote
3answers
76 views
Discretization Error amplification instead of stagnation to machine precision
I wrote a code on Python 2.7.5 to solve numerically the following differential equation.
$\frac{\partial^2f}{\partial x^2}=-S$
$S=\pi^{2}\sin(\pi x)$
S is chosen that way in order to have $f= \sin(\...
1
vote
1answer
48 views
Problems with deriving an equation for a finite-difference scheme given in the journal paper
I'm reading this paper and trying to follow everything that the author has done.
A Comparative Study of Finite Volume Method and Finite Difference Method for Convection-Diffusion Problem
But there ...
2
votes
0answers
79 views
Monotone, monotonicity preserving, LED, TVD, NVD, bounded, stable and stability preserving discretization schemes [closed]
When it comes to discretization schemes for finite volume method, the following terms can be found in literature:
monotone schemes
monotonicity preserving schemes
local extremum diminishing schemes
...
1
vote
5answers
203 views
Don't we care about the numerical diffusion in the diffusion term?
In the context of the solution of advection-diffusion equations by finite volume method, many numerical schemes, papers and book chapters are dedicated to address the numerical diffusion and/or ...
4
votes
0answers
97 views
How to account for the interface between two different phases in a discretized diffusion model?
I have tried to set up a model for the diffusion of a gas into a liquid. The two media are next to each other and the geometry is spherical because the system should simulate the diffusion out of a ...
3
votes
1answer
95 views
Discontinuous Galerkin FEM : Control points are mid-points of edges instead of nodes
I am thinking to use discontinuous galerkin FEM (DGFEM) method to estimate discontinuous displacement field $u: \Omega \rightarrow \mathbb{R}^2$ at the crack surface of a material.
The domain is ...
0
votes
1answer
38 views
Implementation of stochastic cellular automata
In my problem, I have a lattice with a stochastic cellular automaton. In order to simplify a bit, let's say it is 1D. In my system, each node can be type A, B or C. A way to represent the system and ...
1
vote
1answer
111 views
Can this nonlinear advection-diffusion equation be discretized as to only have to solve SPD systems?
Consider a nonlinear advection-diffusion equation of the form
$$
\frac{\partial u}{\partial t} = \nabla \cdot (a(u) \nabla b(u) - \vec{c}(u)u) \tag{1}
$$
on a rectangular domain with Dirichlet ...
1
vote
0answers
178 views
Implement Robin boundary condition (finite volume)
I have a PDE equation with Robin Boundary condition in an annulus system and I should solve it by finite volume method:
\begin{align}
\frac{\partial T_f}{\partial t} - k \left(\frac{\...
1
vote
0answers
78 views
Discrete operator textbooks
This will be a vague question.
When I was writing a finite element matrix assembly routine, a colleague noticed that I had a bug in my code because the sparsity pattern of the one of the blocks didn'...
8
votes
1answer
359 views
Discrete wave simulation - absorbing boundaries?
I wrote a simple 2D wave simulation using the following equations:
$$\frac{\partial^2 u}{\partial t^2}=c^2\nabla^2u$$
Where $\nabla^2$ is the discrete laplace operator using a Von Neumann neighborhood ...
1
vote
1answer
210 views
Finite Volume Polar Discretization: Lengths
Given a uniform polar grid, as in the figure below:
and a FV discretization of a gradient for example:
$\frac{\partial p}{\partial \varphi} = 0$
$\Delta r \frac{p_e - p_w}{\Delta \varphi} = 0$
My ...
2
votes
0answers
75 views
Literature on numerical solving based on multiple meshes?
Consider solving a differential equation system for 1D, 2D, or 3D.
It involves various input and output "field" variables, which, more often than not, correspond to various physical quantities ...
5
votes
0answers
255 views
Order of accuracy of FVM discretization
I've recently got interested in CFD and started a small project by solving the radial Reynoldsequation. Why the Reynoldsequation? I recently encountered it through my studies and somehow got stuck :)
...
1
vote
0answers
240 views
Rhie and Chow Pressure Velocity Coupling
In a collocated grid, which velocity is used in the convective term in the momentum equation? Is the Rhie and Chow constructed face velocity or an average of the adjacent cell center values(in a ...
1
vote
0answers
71 views
Three steps of pde numerical solution and nonlinear equation
I'm very new here.
I'm trying to solve nonlinear elliptic equation
$$
(n(u)u')' = f(u)
$$
and face with crucial misunderstanding.
As I suppose, the procedure of solving some nonlinear equation ...
-4
votes
1answer
148 views
Use Finite Difference Discretization to find approximate solution to the Poisson's equation
I've just been introduced to the Poisson's equation. I've never had the need to dealt with PDE, so I'm a bit lost.
Apparently we can compute an approximate solution of the Poisson's equation
$$\frac{...
2
votes
2answers
113 views
Combining trapezoidal rule with upwind scheme
I want to discretize and numerically solve the equation:
\begin{equation}
v(k)\dfrac{\partial f}{\partial z} + F(z)\dfrac{\partial f}{\partial k} = \alpha f(z,k) \,\,.
\end{equation}
Discretization ...
1
vote
1answer
248 views
Mass Lumping in case of Dirichlet boundary conditions
I'm currently trying to implement a FEM to solve a type of wave equation with homogeneous Dirichlet boundary conditions by using standard $\mathcal{P}_1$ triangle elements and an explicit scheme for ...
1
vote
0answers
34 views
Inner Products vs Discretizations of Functions
Let $f:(0,1)\rightarrow [0,1]$ be a $\mathcal{C}^{\infty}$ function.
We say that a function $D_r^f:\{1,...,r\}\rightarrow [0,1]$ is an $r$-discretization of $f$ if $$D_r^f(j) = \frac{1}{|a(j)-b(j)|}\...
4
votes
1answer
122 views
Space-time Galerkin of Burgers changes the convection speed
tldr: Can space-time Galerkin schemes applied to convection-diffusion problems lead to effects on the convection velocity?
For time $t\in (0,1)$ and the spatial variable $\xi \in (0,1)$, I am ...
4
votes
1answer
150 views
Motivation behind Collocation Method
In previous question "Motivation behind Galerkin method", Paul gives a good and easy-understanding explanation indicating that the Galerkin method is a kind of projection method. Can anyone explain ...
1
vote
0answers
42 views
Second and Higher Order Order Corrector in Spectral Deferred Correction
I am trying to work out a second order or higher order correction for the method of Spectral deferred Correction (SDC). Specifically using as a corrector a second order or third order multi-step.
In ...
4
votes
1answer
134 views
High order unconditionally stable discretization for a scalar hyperbolic PDE
In order to numerically solve the following differential equation:
\begin{equation}
\text{Fr}\{f\} := v(k)\dfrac{\partial f(z,k)}{\partial z} - F(z) \dfrac{\partial f(z,k)}{\partial k} = -\dfrac{f-...
3
votes
1answer
350 views
Second order interpolation scheme
On a grid I am having the values of a physical quantity say for example Temperature, at the E,W,N,S and P node all of them being calculated using a second order discretization scheme. I want a second ...
2
votes
1answer
355 views
How to define residual in multigrid approach?
I wish to solve the two-dimensional Navier Stokes equations using multigrid method on a collocated grid using the Predictor-Corrector method mentioned below. But first, let me elaborate on what I had ...
0
votes
2answers
48 views
How to choose a good distribution for visualizing phase changes in the nature of the roots of a quadratic equation
I'm not sure if this SE site is the best one for this question, so let me know where it should be moved to if you think it doesn't belong here.
After learning about the quadratic formula, I'm ...
1
vote
0answers
44 views
Does scaling factor affect discretization?
Suppose I want to solve the below equation numerically.
$$
\frac{dy}{dx}=y
$$
I'd like to normalize the space discretization by choosing
$$
a\bar{x}=x
$$
where I assume $\bar{x}$ is unity.
Then the ...
3
votes
1answer
95 views
Apply second order finite difference discretization for mixed boundary condition
I want to solve the problem below
\begin{equation}
\begin{aligned}
\eta u-\Delta u &=f, &\text{in $\Omega$}\\
\end{aligned}
\end{equation}
where $\Omega=(0,1)\times (...
1
vote
1answer
369 views
Discretization method for advection equation without numerical diffusion
Given the advection equation for an incompressible flow field
$$\frac{\partial c}{\partial t} + \mathrm{Pe} \frac{\partial c}{\partial x} = 0$$
what would the best method be for discretizing this ...
4
votes
1answer
105 views
Boltzmann equation and equilibrium distribution function
The free-streaming term in stationary 1D Boltzmann equation satisfies the equilibrium distribution function, that is:
\begin{equation}
\text{Fr}\{f\} := v(k)\dfrac{\partial f(x,k)}{\partial x}
- \...
1
vote
0answers
36 views
Stability in Discretization of 1D Stationary Boltzmann equation
I want to discretize and numerically solve the following PDE:
\begin{equation}
v(k)\dfrac{\partial f}{\partial x} + E(x)\dfrac{\partial f}{\partial k} = S\{f\}
\end{equation}
using finite volume (box ...
8
votes
1answer
484 views
Finite difference recursion and higher order
This may be a trivial question, but I've always wondered...
For the classical, central finite difference schemes, if I'm interested in determining the second derivative, does applying the first ...
