CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Instantaneous point release / advection-diffusion (

april0100 February 20, 2013 15:11

Instantaneous point release / advection-diffusion
1 Attachment(s)
Hi All,

I am simulating what should be a very simple problem: laminar flow in a 2-d channel with scalar transport. I've amended icoFoam to include the scalar transport equation:


fvScalarMatrix TEqn
      + fvm::div(phi, T)
      - fvm::laplacian(DT, T) == f

where 'f' is a matrix that contains 0 everywhere at all times except at the time and location of the instantaneous release.

My problem is that the scalar 'T' is not properly being advected from the domain. That is, after some time the concentration of my scalar levels off and does not decrease further as shown in the attached plot.

I am thinking I may have a boundary condition problem. Currently I'm using zeroGradient boundaries on all sides of the channel for T. The inlet velocity is parabolic with a maximum centerline value of 2 m/s. If anyone has any ideas where I might be erring I would greatly appreciate it. Thanks,


april0100 March 30, 2013 12:59

In case anyone else has this problem, the solution is to set the inlet boundary condition to a fixed value of 0.

All times are GMT -4. The time now is 19:49.