Questions related to solving the advection-diffusion equation using numerical methods, including derivation and implementation of boundary conditions.

94 questions
Filter by
Sorted by
Tagged with
28 views

### Does the time-dependent 1D advection-diffusion with point sources have an analytical solution?

I am looking for the analytical solution of 1-dimensional advection-diffusion equation with several point sources, Q, along the axial length of a cylinder through which the fluid flow occurs. Neumann ...
218 views
+50

### How to implement point source or volume source in finite element implementations

I'm trying to do a simple implementation to study the advection-diffusion-reaction dynamics in a straight pipe. I have points positioned along the length of the pipe (blue dots in the image above). I ...
88 views

### Implementing Dirichlet BC for the Advection-Diffusion equation using a second-order Upwind Scheme finite difference discretization

i am implementing a Matlab code to solve the following equation numerically : $$(\frac{\partial c}{\partial t} =-D_{e} \frac{\partial^2 c}{\partial z^2} +U_{z}\frac{\partial c}{\partial z})$$ with ...
46 views

29 views

### Semi-analytical/empirical modelling of wall boundary conditions in advection-diffusion-reaction equation with distributed source

Let's suppose I need to numerically solve a 3D steady-state transport equation of the form $$\nabla \cdot (\mathbf{u} c) = \nabla \cdot (D \nabla c) - \lambda c + S$$ where $c$ is the transported ...
289 views

528 views

3k views

### BDF vs implicit Runge Kutta time stepping

Are there any reasons for why one should choose high order implicit Runge Kutta (IMRK) over BDF time stepping? BDF seems much easier to me as $q$ stage IMRK needs $q$ linear solves per time step. ...
40 views

143 views

106 views

263 views

### Cauchy problem with a change of variables

Consider an advection-diffusion problem $$u_t+au_{xx}+bu_{x}=0,\quad x>0.$$ Now I want to remove the drift and rewrite the problem for the new domain. I use the change of variables $y=x-bt$, and ...
758 views

### Impose Neumann Boundary Condition in advection-diffusion equation 1D

when solving the advection equation in 1D that is: $$\frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0$$ with $u'(t,0) = 0$ and $u(t,L) = 0$ , $u(0,x) = u_{0}$ one numerical ...
412 views

### How write a integration loop in fortran, leapfrog scheme to solvind PDE (advection)?

I want to resolve numerically this equation using of difference finite method with Leapfrog Scheme $$\frac{\partial{u}}{\partial t}+ v \frac{\partial{u}}{\partial x}= 0$$ I'm trying to write a code ...
334 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 ...
427 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,...
54 views

### Finite difference Neumann boundary conditions: uneven weighting of edge nodes?

Originally asked this on math.stackexchange, but I figure it's also appropriate here. I'm reading through some finite difference code for a diffusion equation and came across something odd for the ...
803 views

### Methods of solving non-linear advection-diffusion systems beyond Newton-Raphson?

I'm working on a project where I have two adv-diff coupled domains through their respective source terms (one domain adds mass, the other subtracts mass). For brevity, I'm modeling them in steady ...
53 views

### Jacobian Elements for Coupled Drift-Diffusion System using Vertex-Centered Finite Volume

I'm trying to solve the fully coupled drift-diffusion system using Newton's Method. Although I eventually plan to potentially use a Jacobian-Free Newton-Krylov approach, this is still something that I ...
237 views

### Closed boundary conditions in finite difference method for diffusive-advective equation

I am implementing a finite difference method in solving the diffusive-advective equation: $$u_t + v \cdot u_x = D\cdot u_{xx}$$ (v, D are constants). Planning to use the operator splitting method (...
215 views

### Finite Element Stabilization for Drift-Diffusion/Advection-Diffusion Equations

I've tried my best to look through the relevant suggested similar questions when posting this, and hopefully this contains enough new material to not be considered a duplicate. I'm currently trying ...
104 views

### finite differences on a slanted grid — advection diffusion equation

I used pretty much all my expertise with finite differences to solve an advection-diffusion equation with space-dependent coefficient with a grid in the $x$,$z$ domain with regular spacings. Something ...
1k views

### How do you debug numerical code, what could be source of this oscillatory error?

Quiet a lot of insight can be gained form experience, I was just wondering if anybody has seen something similar to this before. The plot shows the initial condition (green) for the advection-...
873 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 ...
294 views

### Should particles in Smoothed Particle Hydrodynamics (SPH) always move during a simulation?

Or can they just be used as an interpolation points and use some other "transported property" which are just evolved and propagated from boundary conditions like for eg. heat conduction through a ...
413 views

### 1D k-epsilon turbulence model in a turbidity current

I am trying to implement a 1D k-$\varepsilon$ turbulent model for a turbidity current, hence the conservation equation for $c$. I'm solving for the variables $u,c,k$ and $\varepsilon$. The remaining ...
1k views

### CFL condition in Stokes equation

Does the CFL condition play any role in a pure Stokes flow, i.e. convective term is neglibile, or vanishing? If not, what is the "equivalent" condition for stability? I have read something about the ...