
[Sponsors] 
January 21, 2011, 06:12 
boundary condition at the wall

#1 
Member
Join Date: Apr 2009
Posts: 34
Rep Power: 9 
Hello,
I would like to model mass transport using convectiondiffusion equation. I want to apply a boundary condition of this type: − D ∇c n + u c n = Js, where c is the concentration, u the velocity vector and n the normal vector of the domain. Any help? Thanks in advance 

January 21, 2011, 07:35 

#2 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,397
Rep Power: 96 
Are you proposing solving for c (or the flux of c) at the boundary (or something like that)? Sounds tricky. I think you will probably need to write a fortran subroutine for the boundary condition to do this.


January 21, 2011, 09:10 

#3 
Member
Join Date: Apr 2009
Posts: 34
Rep Power: 9 
I want to solve for the concentration of a solute, but Js is a solute flux across the boundary (it is described by a simple CEL expression).
Can you give more details about what to program with fortran code, because i am not very familiar. Thanks anyway 

January 22, 2011, 06:02 

#4 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,397
Rep Power: 96 
What is Js a function of?
What is D? What is known and what is not known in your equation? 

January 22, 2011, 06:14 

#5 
Member
Join Date: Apr 2009
Posts: 34
Rep Power: 9 
thanks again for your answer
Js is a function of wall shear stress (known from solving NavierStokes) D is the diffusivity (known) u is the fluid velocity given by the NavierStokes (unknown but it is solved) c is the concentration (unknown and it is solved using the convectiondiffusion equation) I want to calculate c at the wall which is described by the above equation. If you see − D ∇c n + u c n = Js is like convectiondiffusion equation (transport equation in CFX: ∇(− D ∇c + u c) = 0) 

January 22, 2011, 06:29 

#6 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,397
Rep Power: 96 
If you can recast it into a form like the transport equation then you can simply use an additional variable with a flux boundary = Js (or some function like it). But if that is not possible you have a much harder problem on your hand and you will need to write a solver in fortran to do it. I don't think example fortran comes with CFX any more, you would have to talk to the support people to get one I guess.


January 22, 2011, 07:17 

#7 
Member
Join Date: Apr 2009
Posts: 34
Rep Power: 9 
Since Js describe the flux across the wall, could i use the flux in boundary condition?
But what is the mathematical formulation of this boundary condition? So that i could recast it in my form? P.S. the positive flux in describes direction into the domain. If i want it to move out of the domain, should i use Js? 

January 23, 2011, 05:49 

#8 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,397
Rep Power: 96 
For available boundary condition options see the documentation.
For what positive and negative does, give it a try and see what happens. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
isothermal wall boundary condition  Neil  Main CFD Forum  3  November 9, 2015 03:34 
How to set rotating wall as boundary condition  babu  FLUENT  7  December 14, 2013 05:48 
wall boundary condition in wall fired boiler  bharath  CFX  0  August 9, 2010 21:40 
Change wall boundary condition in multiphase?  Jacques  FLUENT  0  November 27, 2008 07:26 
Variable temperature wall boundary condition  Virginie  FLUENT  2  October 2, 2007 09:23 