Questions tagged [advection]

To move in some direction (as a fluid does in a pipe). Often contrasted with diffusion, which is a spreading out without necessarily having any movement of the field as a whole.

Filter by
Sorted by
Tagged with
0
votes
1answer
72 views

Implementing structured grid boundary conditions using NumPy arrays?

I am making a toy code in Python to solve the advection equation $$u_t + cu_x = 0$$ with, for example, periodic boundary conditions. Background information The numerical grid is specified like this: ...
2
votes
3answers
505 views

Finite volume piecewise linear 2D advection develops instability

I'm developing a finite volume solver for the simple twodimensional advection equation with constant velocities $u, v$ and constant mesh spaces $\Delta x$: $$ \frac{\partial \rho}{\partial t} + u \...
0
votes
0answers
36 views

Is semi-Lagrangian 1D advection identical to upwind Euler, when $|u|<\Delta x/\Delta t$?

This looked true to me, I worked through the algebra for 1D and confirmed it. But no one seems to mention it, so maybe I'm missing something... I'm using Bridson's SIGGRAPH 2007 course notes. [5MB ...
1
vote
1answer
47 views

How to calculate the analytical solution of linear advection equation with Dirichlet's boundary conditions?

I am trying to find the solution of linear advection equation of the form: $\frac {\partial c}{\partial t}+u\frac {\partial c}{\partial x}=0$ $c(x,0)=0$ $c(0,t)=\{c_0 \ \text{for}\ t \leq t_1 \text{...
0
votes
0answers
46 views

Computationally obtaining the convergence rate of upwind scheme for Advection equation

The Advection equation (with velocity = 1) is $${\partial u \over \partial x} + {\partial u \over \partial t} = 0$$ I am trying to solve the equation with periodic BC. One of the ways to numerically ...
9
votes
5answers
643 views

How can I derive a bound on the spurious oscillations in the numerical solution of the 1D advection equation?

Suppose I had the following periodic 1D advection problem: $\frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0$ in $\Omega=[0,1]$ $u(0,t)=u(1,t)$ $u(x,0)=g(x)$ where $g(x)$ has a ...
2
votes
1answer
112 views

Why the numerical solution of advection-dominant problem is challenging

In many CFD text books, usually there is a dedicated chapter for advection term discretization. Why discretization of such term in advection-dominated problems and near the discontinuities is ...
32
votes
2answers
4k views

Strange oscillation when solving the advection equation by finite-difference with fully closed Neumann boundary conditions (reflection at boundaries)

I am trying to solving the advection equation but have a strange oscillation appearing in the solution when the wave reflects from the boundaries. If anybody has seen this artefact before I would be ...
2
votes
0answers
67 views

Efects from the boundary in advection equation [duplicate]

I am implementing the advection equation $u_x+(1/c)u_t=0$ following a Crank-Nicholson finite difference scheme. The equation for this is \begin{eqnarray*} -\frac{\gamma}{4} w_{n-3 j+1} + w_{n-2 j+1} ...
1
vote
0answers
140 views

Inflow and outflow boundary conditions for advection-diffusion equation

I'm trying to solve this advection-diffusion equation (ADE): $$\frac{\partial \phi}{\partial t} + \nabla \cdot (-D \nabla \phi + \mathbf{u} \phi) = 0$$ In fact, this ADE framework is coupled to a ...
26
votes
2answers
5k views

Is Crank-Nicolson a stable discretization scheme for Reaction-Diffusion-Advection (convection) equation?

I am not very familiar with the common discretization schemes for PDEs. I know that Crank-Nicolson is popular scheme for discretizing the diffusion equation. Is also a good choice for the advection ...
14
votes
4answers
6k views

Boundary conditions for the advection equation discretized by a finite difference method

I am trying to find some resources to help explain how to choose boundary conditions when using finite difference methods to solve PDEs. The books and notes which I currently have access to all say ...
8
votes
1answer
11k views

How to discretize the advection equation using the Crank-Nicolson method?

The advection equation needs to be discretized in order to be used for the Crank-Nicolson method. Can someone show me how to do that?
1
vote
1answer
273 views

Advection equation in 2D using finite differences - the scheme works, but the pulse loses “energy”

I am trying to solve the following equation $ \partial_t g(x,y,t) = - v\left(\partial_x + \partial_y\right) g(x,y,t)$ using finite differences (here $v>0$). The equation is also solvable ...
1
vote
1answer
366 views

Is there a general analytic solution to 1D advection of velocity, $u_t=-uu_x$?

This is to help me relate continuous and discrete, predict what my scheme should be doing, and move toward using the method of manufactured solutions. There's a solution for constant velocity $c$ ...
7
votes
1answer
258 views

Solving the Advection Equation with Forcing using the Discontinuous Galerkin Method

I've been learning about the Discontinous Galkerin Method by reading the book by Hesthaven and Warburton and have ran into a problem with the advection equation with forcing $u_t + u_x = g(x,t)$ ...
3
votes
4answers
1k views

Advection equation with finite difference: importance of forward, backward or centered difference formula for the first derivative

I am a newbie in finite difference methods, so I apologize in advance if the question is trivial. I am trying to solve the advenction equation, i.e. $\frac{\partial \phi(x,t)}{\partial t} + v \frac{\...
0
votes
3answers
1k views

Neumann boundary conditions for the upwind scheme applied to the advection equation (Python)

I'm trying to solve the linear advection equation $$u_{t} = cu_{x}, \\ x \in [x_{0}, x_{e}], \quad t \in (0, T], \quad c \in \mathbb{R} \\ u(x,0) = f(x)$$ Note that for $c > 0$, the solution is ...
1
vote
0answers
47 views

Extrapolating to non-fluid cells (for Shallow Water Equations), for a shore/beach?

The water height $h$ and 2d velocity field $(u,w)$ are "extrapolated to non-fluid-cells, i.e., setting $h$ equal to the value in the nearest fluid cell." [Bridson] I'm using finite differences. ...
2
votes
1answer
150 views

Finite difference method not working for advection PDE with negative coefficient

I'm trying to solve a very simple advection PDE $\frac{\partial u}{\partial t}+c\frac{\partial u}{\partial x}=0$ where $c<0$. I have been able to implement a simple Modelica code to solve the ...
1
vote
0answers
225 views

Crank-Nicolson scheme in space for advection equation

Consider the equation $$\frac{\partial}{\partial t}v(t,x)=\frac{\partial}{\partial x}v(t,x)$$, for $t,x\in\mathbb{R}$. I'd like to solve this equation forward in space and backward in time, ...
1
vote
0answers
179 views

Why can I not solve the negative advection equation (backwards in time)?

Suppose we have the negative, inhomogeneous advection equation: $$\left(\frac{\partial}{\partial x}-\frac{1}{c}\frac{\partial}{\partial t}\right)v(t,x)=u(t,x)\qquad(t\in\mathbb{R}_{+},x\in\mathbb{R})$$...
1
vote
1answer
433 views

Crank-Nicolson method for inhomogeneous advection equation

Suppose we have the inhomogeneous advection equation $$\left(\frac{\partial}{\partial x}+\frac{1}{c}\frac{\partial}{\partial t}\right)u(t,x)=v(t,x)$$ for $u,v:\mathbb{R}\times\mathbb{R}\to\mathbb{R}$ (...
25
votes
1answer
5k views

Conservation of a physical quantity when using Neumann boundary conditions applied to the advection-diffusion equation

I don't understand the different behaviour of the advection-diffusion equation when I apply different boundary conditions. My motivation is the simulation of a real physical quantity (particle density)...
1
vote
0answers
91 views

Discontinuos Galerkin Method - inhomogeneous advection problem

I'm currently trying to get into this topic. I've learned that the basic scheme for the advection problem ($D_{x}u+a*D_{t}u=0$) can be solved in a scheme like $$ M^{k}\frac{d}{dt}u^{k}_{h}-(S^{k})^{T}...
2
votes
1answer
87 views

Deposition model in laminar flow

I have a chamber full with a fluid flowing horizontally in laminar regime from one side to the other. It carries a suspension with concentration $c$. This suspension also falls to the bottom of the ...
2
votes
0answers
109 views

Numerical diffusion in during advection of a free surface in an FE context

I am currently working on a project where a two-phase flow is considered. The phases are described using a level set approach and a signed distance function from the interface between the phases where ...
1
vote
3answers
177 views

When is it safe to ignore the diffusion term in an advection-diffusion equation?

Given the one dimensional equation: $\epsilon\frac{\partial^2u}{\partial x^2} +\frac{\partial u}{\partial x} = 0 $ with $0\le\epsilon \ll1$ with boundary conditions $u(0) = 0$ and $u(1) = 2$, we ...
2
votes
2answers
110 views

How to support or contradict a hypothesis on unconditional stability using numerical optimization

The main motivation behind my next question is that I think I derived a higher order numerical scheme for linear advection equation that is unconditionally stable using Von Neumann stability analysis. ...
5
votes
2answers
825 views

Advice on numerical solution for 2D hyperbolic PDE with zero flux boundary conditions

I would like to numerically solve a hyperbolic PDE of the form $\frac{\partial\theta_t}{\partial t}(x,y)+\frac{\partial\left[\theta_t \gamma_t^x\right]}{\partial x}(x,y)+\frac{\partial\left[\theta_t \...
1
vote
0answers
212 views

Simple steady-state advection problem: do I need FVM with upwind scheme?

I have a 2D (x,y) scalar advection problem that describes net blowing snow ($q$) transport at a point. This takes the form $$q = A - F*\nabla\cdot(q {\bf \hat u}),$$ where A ($kg\cdot m^{-2}\cdot s^{...
1
vote
1answer
281 views

Upwind difference for velocity in staggered grid

I am reading the paper, http://math.mit.edu/~gs/cse/codes/mit18086_navierstokes.pdf In the paper, the nonlinear term is treated as mix of central central difference and upwind difference using a ...
1
vote
0answers
233 views

Numerical solution of non-linear advection equation other than inviscid burgers

I am solving a non-linear advection equation of the form $u_t + f(u)_x = 0$ where $f(u)$ is a complicated function of $u$. I am solving this equation using a first order fully implicit scheme (...
1
vote
0answers
179 views

Spherical Advection Discretization (boundary nodes)

Consider the spherical advection problem: describing the conservation of a property $u$ in a closed spherical domain. $$ \frac{\partial u}{\partial t}+\frac{1}{r^2}\frac{\partial }{\partial r}\left(r^...
1
vote
0answers
60 views

Direction-splitting for SSP-RK schemes

What are the implications of applying a direction-splitting within each stage of an SSP-RK scheme? For instance, given a standard advective transport type equation: $$ \partial_{t}Q + \operatorname{...
1
vote
1answer
400 views

More Smearing with decreasing timestep in advection problems

I find it kind of counter intuitive, that the result of an advection gets more smeared out at the borders when decreasing the timestep (which should make it more accurate). Let there be a equally ...
3
votes
0answers
251 views

Corner Transport Upwind for Linear Advection in Arbitrary Velocity Field

I need to implement a 3D version of the Corner Transport Upwind (CTU) finite volume method (in python); and so I've been reading Leveque, "Finite Volume Methods for Hyperbolic Problems" which I think ...
2
votes
1answer
676 views

Numerically computing the advection equation

I am trying to write a program to compute the advection equation. $$u_t +u_x = 0$$ I use the spectral method for the spatial derivative $u_x$ and the leapfrog method for the time derivative $u_t$. ...
1
vote
1answer
438 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 ...
1
vote
1answer
155 views

CFD: Doubt with time convergence in advection fully implicit upwind scheme

I'm trying to solve an advection - convection problem using an implicit upwind scheme - you can see here the finite difference discretization used. I start the model (built from scratch on Scilab) ...
0
votes
1answer
286 views

boundary conditions of linear advection problem

I am solving the 1D advection problem given by: $$\frac{\partial u}{\partial t}=-c\frac{\partial u}{\partial x}$$ where c is the wave speed, and u is the unknown field variable, and x and t are time ...
0
votes
1answer
283 views

Find cfl condition

We have the advection equation $u_t+a u_x=0, a>0, 0<t<T_f, x \in \mathbb{R}$ with initial condition $u(0,x)=u_0(x)$. Suppose that we have the following sheme: I want to find the CFL ...
4
votes
2answers
909 views

Advection equation using the finite element method

I want to solve this simple advection equation using the finite element method. $$\frac{dc}{dt}+v\cdot\nabla c = f$$ What's the best FEM discretization for this? I have tried using the standard/...
3
votes
2answers
150 views

what do zero real parts of eigenvalues mean? Any good references?

I am solving a 1D advection problem of the the form $$dQ/dt=[A]Q$$ where {Q} is the vector of unknowns and [A] is the matrix of coefficients of spatial discretisation. I have worked out the ...
7
votes
2answers
292 views

Stabilization of convection-dominated flow and turbulence modeling

Are stabilization techniques for convection-dominated flows like SUPG+PSPG, interior penalty methods, etc. able to handle turbulent flows without tubulence model being employed, at least up to some ...
2
votes
0answers
136 views

what do positive real parts of eigenvalues mean?

I am solving a 1D advection problem of the the form $$d{Q}/dt = [A]{Q}$$ where {Q} is the vector of unknowns and [A] is the matrix of coefficients of spatial discretisation. I have worked out the ...
2
votes
0answers
109 views

Calculating theoretical order of accuracy of least squares fit advection scheme

I'm familiar with finding the order of accuracy using von Neumann analysis for finite difference schemes formulated using Taylor series expansions. But is there a similar technique for finding the ...
2
votes
1answer
2k views

Implementation of 1D Advection in Python using WENO and ENO schemes [closed]

I'm trying to implement 1D advection solver using WENO and ENO schemes. \begin{equation} \frac{\partial u}{\partial t} + \frac{\partial f(u)}{\partial x} =0 \end{equation} where: \begin{...
7
votes
1answer
242 views

Implementation of convection scheme given by normalized variable diagram

In finite difference and finite volume methods, convection schemes (upwind, central, quick, ...) are usually shown in a normalized variable diagram. The diagram gives the normalized face variable as a ...
4
votes
2answers
449 views

Usability of upwind finite difference schemes

NOTE: I asked this on Mathematics Stack Exchange and there were no answers. So, I thought I might try here. Upwind schemes like the classic "upstream" scheme, can be used to solve, for example, the ...