Linked Questions

1 vote
0 answers
232 views

Closed (Robin) boundaries in advection-diffusion equation with FDM

I am solving the equation $$ \frac{\partial \phi}{\partial t} = \frac{\partial}{\partial x} \left( D \frac{\partial \phi}{\partial x} + v\phi \right) $$ using finite differences. I want to include ...
rugermini's user avatar
3 votes
1 answer
1k views

Implementing Robin Boundary condition (finite difference)

I'm interested in applying Robin boundary condition to a convection-diffusion problem in 1D. In the following system, $$\frac{\partial C}{\partial t} = D\frac{\partial ^2 C}{\partial x^2} - v\frac{\...
Natasha's user avatar
  • 421
2 votes
1 answer
494 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 (...
Gabriele's user avatar
  • 173
3 votes
1 answer
1k views

Implementation of gradient zero boundary conditon in advection-diffusion equation

My question is about Finite Element Method. I want to know how to implement "gradient zero" conditions to advection-diffusion equations in conservative form like, $\frac{\partial \rho}{\partial t} + ...
Phoenix Kyoma's user avatar
5 votes
3 answers
4k views

Conservation of Mass in 1D Advection-Diffusion Equation

My long-term goal is to numerically solve the 1D advection-diffusion equation of the form: $$\frac{\partial u}{\partial t}=\frac{\partial }{\partial x}\left( v(x,t) u+D\frac{\partial u}{\partial x}\...
tom's user avatar
  • 231
8 votes
1 answer
879 views

Conservative finite-difference expression for the advection equation

Following on from the earlier question I am trying to derive a finite-difference scheme for the advection equation which is conservative. It was suggested that for advection equation with variable ...
boyfarrell's user avatar
  • 5,409
14 votes
1 answer
2k views

Can the advection equation with variable velocity be conservative?

I am trying to understand the advection equation with variable velocity coefficient a bit better. In particular I don't understand how the equation can be conservative. The advection equation, $$ \...
boyfarrell's user avatar
  • 5,409
6 votes
1 answer
3k views

Trouble implementing Neumann boundary conditions because the ghost points cannot be eliminated

Neumann boundary conditions are implemented by introducing ghost points outside the domain and then using the boundary conditions to eliminate the ghost points. For example, see this question. I ...
boyfarrell's user avatar
  • 5,409
5 votes
3 answers
2k views

No flux boundaries for mixed hyperbolic parabolic PDE

I read this post, "Conservation of a physical quantity when using Neumann boundary conditions applied to the advection-diffusion equation" and although it is the same type of equation it does not fit ...
Moritz's user avatar
  • 343
8 votes
3 answers
4k views

Open boundary conditions with the advection-diffusion equation

Following on from my previous equation I'm would like to apply open boundary condition to the advection-diffusion equation (with reaction term), $$ \frac{\partial \phi}{\partial t} = \frac{\partial}{\...
boyfarrell's user avatar
  • 5,409
34 votes
2 answers
6k 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 ...
boyfarrell's user avatar
  • 5,409