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

Is there anyone who has experience using dsmcFoam solver?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 23, 2011, 09:23
Default Is there anyone who has experience using dsmcFoam solver?
  #1
spk
Member
 
Join Date: Aug 2009
Posts: 67
Rep Power: 8
spk is on a distinguished road
Hi guys,
I'm new in OpenFOAM and i would like to solve a 2D orifice geometry using the dsmcFoam solver. Has anyone used this solver before?
I will honestly appreciate any advice or help!!!
Thanks
spk is offline   Reply With Quote

Old   October 24, 2011, 05:13
Default
  #2
Senior Member
 
Join Date: Mar 2011
Posts: 174
Rep Power: 6
anon_a is on a distinguished road
Well, you need to provide some more information.
What are your boundary conditions?
The pressure BC is not implemented by default and you have to program it
(modify src/lagrangian/dsmc/submodels/InflowBoundaryModel/FreeStream/FreeStream.C to read n from a dictionary for each incoming patch).

Also, are you sure you need this solver? What is your Kn number upstream?
The dsmcFoam solver can not deal with 2D AXISYMMETRIC flows,
which means that you have to turn to 3D, increasing the computational cost.
Just trying to help you avoid expensive calculations.

If you need some specific information, let me know
anon_a is offline   Reply With Quote

Old   October 24, 2011, 05:37
Default
  #3
spk
Member
 
Join Date: Aug 2009
Posts: 67
Rep Power: 8
spk is on a distinguished road
Hi anon_a,
thank you for your answer!
The boundary conditions are:
pressure inlet: absolute pressure 1000 Pa
pressure outlet: vacuum
orifice diameter 0.5 mm
The gas density and the orifice size correspond to the entire regime from free molecule to continuum flow.
If you want and have time i could send you via email an isodensity contours of gas through the orifice solving this problem with Prof. Bird's DS2V program.
I just want to know can i solve the same problem with openFoam?
and if yes, is it easy? because i am new user in OpenFoam!!
Thanks a lot!
spk
spk is offline   Reply With Quote

Old   October 25, 2011, 08:10
Default
  #4
Senior Member
 
Join Date: Mar 2011
Posts: 174
Rep Power: 6
anon_a is on a distinguished road
Yes, you can solve the problem you mention with dsmcFoam.
It needs a little bit of modification but it can be done relatively easy.

Two quick ways of simulating open boundaries with specified pressure
were given at the same time I was replying to you in this thread:

http://www.cfd-online.com/Forums/ope...am-vacuum.html

The first method is only appropriate for vacuum BCs.
The second way is better (but requires the compilation of the code each time you change the boundary number density).

Personally, I would advise to study the code and try to change it in a general manner, using dictionaries. Deeper knowledge really helps to fine tune your runs.
anon_a is offline   Reply With Quote

Old   December 1, 2011, 19:37
Default Boundary conditions dsmcFoam!
  #5
spk
Member
 
Join Date: Aug 2009
Posts: 67
Rep Power: 8
spk is on a distinguished road
I would like help about boundary conditions!

First about inlet. I have pressure 1000 Pa and velocity 0 m/s.
So in folder 0, in file boundaryU i wrote for inlet (in order to define velocity):
type fixed value;
value uniform (0 0 0);
In folder 0, in file rhoN i wrote for inlet (in order to define pressure):
type fixed value;
value uniform <particle number density corresponding to the pressure value>;

Then about outlet. I have vacuum!
in file boundaryU i wrote for outlet:
type calculated;
value uniform (0 0 0);

Are these right?? How can i define the pressure in vacuum?
spk is offline   Reply With Quote

Old   December 2, 2011, 02:58
Default
  #6
Senior Member
 
Join Date: Mar 2011
Posts: 174
Rep Power: 6
anon_a is on a distinguished road
Did you read the link I gave you?
Did you have trouble applying it?
anon_a is offline   Reply With Quote

Old   December 2, 2011, 05:45
Default
  #7
spk
Member
 
Join Date: Aug 2009
Posts: 67
Rep Power: 8
spk is on a distinguished road
Hi anon_a!

I think i cant apply the second method because i dont know the desired outlet number density.
The first method is something that i proposed with the previous post?
spk is offline   Reply With Quote

Old   December 2, 2011, 08:26
Default
  #8
Senior Member
 
Join Date: Mar 2011
Posts: 174
Rep Power: 6
anon_a is on a distinguished road
The number density that you provide to DSMC in the boundaries is not the number density that you expect to have there. Instead, it is the number density that characterizes the incoming stream. Therefore, a vacuum boundary condition means that there are no INCOMING molecules from that surface (and not "no molecules at all").

Therefore, for a vacuum BC, you just need to specify that no molecules are coming in the domain, which can be done with the methods of that post.

EDIT: and this is not done by the rhoN, rhoM, dsmcRhoN files.
anon_a is offline   Reply With Quote

Old   December 2, 2011, 08:58
Default
  #9
spk
Member
 
Join Date: Aug 2009
Posts: 67
Rep Power: 8
spk is on a distinguished road
In the file FreeStream.C what value i have to enter for the desired inlet number density and outlet?
The only value known is the input particle number density for inlet and vacuum for outlet.
Please help i'm new OpenFOAM user!
spk is offline   Reply With Quote

Old   February 20, 2013, 10:36
Default how to intialise mean free path and knudsen number
  #10
New Member
 
murali
Join Date: Feb 2013
Location: pondicherry
Posts: 10
Rep Power: 4
murali is on a distinguished road
hi foamers


im new here to OpenFoam this is my knusden number 0.113 and mean free path at .0068m this is my dimension .1x.06x.06m .and then where i going mention this value in dsmcFoam.(0,constant,system)
murali is offline   Reply With Quote

Old   February 21, 2013, 10:02
Smile hypersonic flow corner
  #11
New Member
 
murali
Join Date: Feb 2013
Location: pondicherry
Posts: 10
Rep Power: 4
murali is on a distinguished road
hi foamers

i am doing my project related to hypersonic flow, therefore anyone know the hypersonic flow corner dsmc code or have?
murali is offline   Reply With Quote

Old   June 26, 2013, 09:53
Default
  #12
New Member
 
Carl
Join Date: Dec 2012
Posts: 14
Rep Power: 4
jiaoxiaolei is on a distinguished road
Hi, murali, there is a case in dsmcFoam, it is the same as yours, you can have a look
jiaoxiaolei is offline   Reply With Quote

Old   July 2, 2013, 04:48
Default hypersonic flow corner
  #13
New Member
 
murali
Join Date: Feb 2013
Location: pondicherry
Posts: 10
Rep Power: 4
murali is on a distinguished road
hi carl


thanks to reply which case is similiar to my cases.
murali is offline   Reply With Quote

Reply

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
Is it possible to accelerate the ode solver of reactingFoam? pajofego OpenFOAM Running, Solving & CFD 2 August 10, 2014 05:38
Quarter Burner mesh with periosic condition SamCanuck FLUENT 2 August 31, 2011 11:34
Could you comare StarCD with CFX 5?Help, please... Suteh CFX 54 November 7, 2001 21:12
CFX 5.5 Roued CFX 1 October 2, 2001 16:49
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 16:25


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