CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

Dirichlet boundary condition for additional variable on the wall

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   July 23, 2012, 11:05
Default Dirichlet boundary condition for additional variable on the wall
  #1
Member
 
Ftab
Join Date: Sep 2011
Posts: 53
Rep Power: 5
ftab is on a distinguished road
Hi,

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?

Thanks!
ftab is offline   Reply With Quote

Old   July 23, 2012, 13:42
Default
  #2
Member
 
Felggv's Avatar
 
Felipe Gobbi
Join Date: Apr 2012
Location: Brazil
Posts: 76
Rep Power: 5
Felggv is on a distinguished road
Have you tried creating a Source Point on the walls?

Don't know if it will serve you... just a hint.
Felggv is offline   Reply With Quote

Old   July 23, 2012, 18:39
Default
  #3
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,803
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Felipe is correct, you can do this using a source term on the boundary surface. That will work.
ghorrocks is offline   Reply With Quote

Old   July 24, 2012, 07:12
Default
  #4
Member
 
Ftab
Join Date: Sep 2011
Posts: 53
Rep Power: 5
ftab is on a distinguished road
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
ftab is offline   Reply With Quote

Old   July 24, 2012, 10:37
Default
  #5
Member
 
Felggv's Avatar
 
Felipe Gobbi
Join Date: Apr 2012
Location: Brazil
Posts: 76
Rep Power: 5
Felggv is on a distinguished road
Hello,

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?

Bye
Felggv is offline   Reply With Quote

Old   July 24, 2012, 11:00
Default
  #6
Member
 
Ftab
Join Date: Sep 2011
Posts: 53
Rep Power: 5
ftab is on a distinguished road
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:
http://www.vizworld.com/2010/06/mit-...create-stents/

or this one:
http://maths.dur.ac.uk/~dma0mpj/summ...lo_Zunino.html

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.
ftab is offline   Reply With Quote

Old   July 24, 2012, 13:45
Default
  #7
Member
 
Felggv's Avatar
 
Felipe Gobbi
Join Date: Apr 2012
Location: Brazil
Posts: 76
Rep Power: 5
Felggv is on a distinguished road
Well,

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.

Bye
Felggv is offline   Reply With Quote

Old   July 24, 2012, 17:25
Default
  #8
Member
 
Ftab
Join Date: Sep 2011
Posts: 53
Rep Power: 5
ftab is on a distinguished road
Felipe,
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.
ftab is offline   Reply With Quote

Old   July 24, 2012, 19:26
Default
  #9
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,803
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
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?
ghorrocks is offline   Reply With Quote

Old   July 25, 2012, 06:11
Default
  #10
Member
 
Ftab
Join Date: Sep 2011
Posts: 53
Rep Power: 5
ftab is on a distinguished road
Glenn,
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.
ftab is offline   Reply With Quote

Old   July 25, 2012, 12:14
Default
  #11
Member
 
Felggv's Avatar
 
Felipe Gobbi
Join Date: Apr 2012
Location: Brazil
Posts: 76
Rep Power: 5
Felggv is on a distinguished road
ftab,

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
Felggv is offline   Reply With Quote

Old   July 25, 2012, 18:31
Default
  #12
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,803
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
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.
ghorrocks is offline   Reply With Quote

Old   January 27, 2013, 13:23
Default
  #13
Member
 
bitak
Join Date: Jul 2009
Posts: 56
Rep Power: 7
bitak is on a distinguished road
Hello.
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!!!!!!!!!!!
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 is offline   Reply With Quote

Old   January 27, 2013, 13:24
Default
  #14
Member
 
bitak
Join Date: Jul 2009
Posts: 56
Rep Power: 7
bitak is on a distinguished road
How I can source term?
bitak is offline   Reply With Quote

Reply

Tags
additional variable, wall boundary condition

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
vorticity boundary condition bearcharge Main CFD Forum 0 May 14, 2010 11:32
Variable temperature wall boundary condition Virginie FLUENT 2 October 2, 2007 09:23
UDF hook on the WALL boundary condition Luke FLUENT 0 June 7, 2006 11:54


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