|June 6, 2012, 06:25||
Inlet Boundary conditions for pressure driven flow in DSMC
Jesús González Serrano
Join Date: Nov 2011
Posts: 8Rep Power: 6
I have been lately trying to program inlet boundary conditions in order to solve a test case.
The boundary conditions implemented are described in the following papers:
I used as a base the Freestream.H and Freestream.C files and I modified them in order to implement the stated inflow conditions.
As well I modified the cloud in order to calculate velocity, temperature (only translational because I run the test case for Neon) and pressure at each time step because it is needed by the inflow model.
The following problems appeared when I compiled and run the test case:
1. In the inlet, after a couple of steps, I obtain negative velocities. Even forcing in the initialisation that all molecules have positive velocities, after few steps velocities in the inlet turn negative. It happens always and I haven't found any reasonable explanation.
2. When using the velocity distribution for entering molecules from the papers, wrong pressure profiles are obtained. But when leaving the original velocity treatment promising results are obtained.
As well an approach in order to calculate self-diffusion coefficient (I am not really sure about this one) and viscosity in each step is done.
CASE 1: Leaving the original velocity distribution for entering molecules
CASE 2: Using the distribution velocity for entering from the papers
I hope anyone has any idea, I think it is really interesting if we could obtain proper boundary conditions to simulate pressure driven flows with OpenFoam.
Last edited by Jesgonse; June 6, 2012 at 09:48.
|June 6, 2012, 10:29||
Join Date: Mar 2011
Posts: 174Rep Power: 7
I am generally interested and have been working on such flows for some time.
Currently, I don't have enough time to check the paper and discuss the BCs but I would suggest that you check your time step.
If the time step is too large, molecules enter in pulses, while they are supposed to enter in a continuous fashion in time.
If you inject them all together in sparsely distributed times, the number of injecting molecules will be unrealistically high and a large part of them will exit the boundaries very fast. This may be why you witness this effect (at least this is what I have seen in the original dsmcFoam code).
Let me know if there are any changes with a smaller time step.
|February 1, 2014, 03:32||
Join Date: Aug 2013
Posts: 8Rep Power: 4
Could you solve your problem?
If so, would you please guide how I can apply these BCs in dsmcFoam?
How should I change FreeStream?
By the way I'm a newbie in C++. Please help me.
Thanks in advance
|October 3, 2014, 09:52||
is the problem solved??
Join Date: May 2014
Posts: 1Rep Power: 0
I too am facing a similar problem. Did you find any solution to this problem?
Please do enlighten me with the solution.
|boundary condition, dsmc openfoam dsmcfoam, dsmcfoam, openfoam|
|Thread||Thread Starter||Forum||Replies||Last Post|
|Proper Pressure Boundary Conditions for Buoyant Flow||mchurchf||OpenFOAM||0||March 25, 2010 13:16|
|steam flow in a pipe driven by a pressure gradient between inlet and outlet||SalvoCalvo||Main CFD Forum||0||March 11, 2010 07:52|
|Domain Reference Pressure and mass flow inlet boundary||AdidaKK||CFX||12||February 11, 2010 21:28|
|Pulsating flow with non-reflective inlet boundary||cc1000||CFX||6||April 27, 2009 08:10|
|New topic on same subject - Flow around race car||Tudor Miron||CFX||15||April 2, 2004 06:18|