CFD Online Discussion Forums

CFD Online Discussion Forums (
-   CFX (
-   -   Dirichlet boundary condition for additional variable on the wall (

ftab July 23, 2012 11:05

Dirichlet boundary condition for additional variable on the wall

In a simulation, I have been asked to put the unity Dirichlet Boundary condition on the wall for concentration.

I have defined volumetric additional variable(AV) with dimension of kg m^-3 and set up the case. The problem is how to set the boundary condition of unit concentration on the wall? CFX asks for constant value of AV (kgm^-3), but obviously the wall is a surface not a volume. Hence, I am puzzled with physical meaning of it and how to apply a boundary condition which would be equivalent to Dirichlet BC of unity on the wall.
In fluent I remember, we simply put the constant value of UDS equal to 1. What would be the setting to get similar results in CFX?


Felggv July 23, 2012 13:42

Have you tried creating a Source Point on the walls?

Don't know if it will serve you... just a hint.

ghorrocks July 23, 2012 18:39

Felipe is correct, you can do this using a source term on the boundary surface. That will work.

ftab July 24, 2012 07:12

Thanks Felipe and Glenn for reply.
If I need to use the source for this boundary, what should be the setting in Boundary details then? Value (kg m^-3) or Flux equal to zero?
Setting the source in dimension of kg m^-2 s^-1 in a steady solution will introduce huge amount of additional variable (concentration) inside the domain and having seudo transient solution approach of CFX for steady cases, the concentration goes up to 10^8 or so!!
Thanks for your help

Felggv July 24, 2012 10:37


In my opinion, Steady State simulations are way too harder to model than transient ones.

Would it be valid to do a transient run and come up with an expression that would cease the source of your contaminant in time just to study this case in a simpler model before running it in steady state?

Just for us to understand, when you set your wall with boundary sources, which options do you have? Would you attach a picture?

What is you material?


ftab July 24, 2012 11:00

Hi Felipe,
It is a simple cylindrical artery with stent inside, which is eluting drug both inside the wall and lumen of the artery.
Check this image for instance:

or this one:

To be honest, I have not set the boundary source, but the flux in boundary details, which has exactly the same dimension and got the stupid results as mentioned.

Felggv July 24, 2012 13:45


What are you trying to study? How it mixes with blood flow? The only information you have is the concentration of the drug and it only goes to the blood flow by concentration difference?

If you have more information you can model it some other way.

Tell us more about your problem, now that we see the picture we can understand better what you want to do.


ftab July 24, 2012 17:25

The purpose is simply to model the transport of drug(concentration AV) from the stent into blood flow and wall both with convection and diffusion.
And due to the concentration difference (1 on the stent and 0 in wall and lumen) the drug is transported. The equations which are being solved are Mass and momentum together with AV transport.
The flow is simply laminar, incomp and steady for now.

ghorrocks July 24, 2012 19:26

Are you sure the Dirichlet BC is the appropriate one? Would a constant flux be better? Or what about a convective type BC, but based on AV concentration, not temperature?

ftab July 25, 2012 06:11

you are right! If it was transient solution it would fit. Even for the steady one, but as I said I am trying to replicate a reference paper who has used this BC.

Felggv July 25, 2012 12:14


As Glenn always recommends, begin with a simpler model, that's why I asked you if you only had the concentration, because if you had the flux it would be easier to replicate the situation!

Even if that's not what you want, other BCs should do the job.

Good luck

ghorrocks July 25, 2012 18:31

If you are sure you need a dirichlet condition on AV concentration then here's how to do it :)

* Set a source term for the AV equation on the boundary to impose the condition on.
* Set the source term equation to -1.0e3*(1.0-(AV concentration variable) and the source term coefficient to 1.0e3. Please note I am recalling this equation from memory so might have got some signs wrong or other little error, but that is the basic idea.

Where 1.0e3 is a "big" number relative to your AV concentration, 1.0 is the value you want for the dirichlet boundary condition and "AV concentration variable" is the AV concentration variable however you have defined it.

This uses a source term to pull the local value to the set value (1.0), and the source term coefficient helps convergence.

bitak January 27, 2013 13:23

I have a problem like you.:)
I want to reproduce an article but I don't know how Should put a boundary condition for wall!!!!!!!!!!!:confused:
my geometry is a tube( d=7mm and L=15.4 cm ).

Boundary conditions for concentration of LDL at the inlet, outlet and wall, sequentially, are as follows:
C(0, y) = C0 ( inlet concentration)

@ outlet flux = 0

on the wal : CwVw (wall concentration * infilteration velocity)
I really be happy and will be thanks if you answer to me.

bitak January 27, 2013 13:24

How I can source term?

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