
[Sponsors] 
Transport a concentration field with convection equation 

LinkBack  Thread Tools  Search this Thread  Display Modes 
December 19, 2020, 13:07 
Transport a concentration field with convection equation

#1 
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 3 
Assume that we solve a simple convection equation for an incompressible flow, by which a concentration scalar field (between 0 and 1) is transported with the flow field.
The geometry is a simple 2d channel with the height of H connected to a bigger channel with Y=2H, and flow enters the domain from one side of the bigger channel and exits from the smaller channel. I set the inlet velocity boundary condition with a constant velocity and i set the inlet concentration boundary condition to 1. My question is about the conservation of the concentration field. Intuitively, when concentration field enters the smaller channel, its value should increase because of the incompressibility. But we know that the concentration field cannot exceed 1. Can anybody please clarify a bit about this confusion? Thanks. 

December 19, 2020, 13:28 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,146
Rep Power: 66 
Quote:
Why should it increase? You are considering a scalar tracer that cannot diffuse and the velocity field is divergencefree. Thus, the equation is dC/dt + div (vC) = 0 > DC/Dt=0 Hence, the concentration C cannot exceed the initial values. Conservation is guarateed if your numerical method is conservative but you can have undershoot and overshoot using high order method without flux limiters. 

December 19, 2020, 17:39 

#3  
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 3 
Quote:
Another short question: If we impose the noslip boundary condition on the walls for the velocity, should we expect a boundary layer for concentration field, similar to the one we have for velocity? or the concentration field 'fills' the whole domain after reaching the steady situation? 

December 20, 2020, 04:30 

#4  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,146
Rep Power: 66 
Quote:
The advection equation says the answer, the concentration is transported along the trajectory. If the value of the concentration is uniform (for example 1) and starts from the steady inlet BC you will see a transient where the concentration profile is somehow similar to a BL but not in the values that remains uniform. After the transient is ended, the whole BL region will be filled by the uniform concentration. 

December 20, 2020, 15:37 

#5  
New Member
Shabi
Join Date: Dec 2020
Posts: 16
Rep Power: 3 
Quote:
Because if the concentration is gonna fill the whole domain with the value of 1 (where the velocity field is not uniform), I think we might get the above conclusion that the number of particles increase until infinity. 

December 20, 2020, 15:50 

#6  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,146
Rep Power: 66 
Quote:
Last edited by FMDenaro; December 21, 2020 at 15:32. 

Tags 
concentration, convection, transport equation 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Domain Reference Pressure and mass flow inlet boundary  AdidaKK  CFX  75  August 20, 2018 05:37 
Solving Userdefined Scalar Transport Equation properly  meepokman  Fluent UDF and Scheme Programming  0  July 10, 2018 07:57 
Adding diffusion term to interFoam transport equation  Gearb0x  OpenFOAM Programming & Development  2  August 30, 2017 07:51 
Boundary conditions for a newly added transport equation  Divyaprakash  OpenFOAM Running, Solving & CFD  0  May 2, 2016 05:32 
transport equation for SGS scalar variance  chinthakindi  Main CFD Forum  0  April 27, 2004 09:26 