CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

dsmcFoam - micro-hole limiting the number of dsmc particles

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

LinkBack Thread Tools Search this Thread Display Modes
Old   June 25, 2015, 06:50
Default dsmcFoam - micro-hole limiting the number of dsmc particles
New Member
AdriÓ Salvador
Join Date: Jun 2015
Posts: 4
Rep Power: 11
Araist is on a distinguished road
Dear all,

I'm running one hell of a case where I have to simulate mollecules of Helium passing through a skimmer (micro-hole).

I use an axi-symmetric geometry (using two symmetry planes at a small angle) which is composed of two regions: the inlet region, where the number density is bigger and the outlet region, where is much smaller.

The reason for the drop of the number density is that between the two regions there is a wall with a micrometer-range hole centered at the axis of symmetry. The particles can only pass from the inlet region to the outlet through this small hole.

The problem is that, within a reasonable computational time, i find it practically impossible to have enough dsmc particles in the outlet region (I'm aiming for 2-3 particles per cell at least)

I have thought about the following strategies:

1. Try to find a much more efficient parallelization procedure, I will try if scotch increases my efficiency enough. However I believe I would need at least a x20 increase of speed to be able to compute enough particles.

2. Time-averaging: I have not enough particles per cell, but if I let the program run enough time I do have enough particles that have been in that cell under a stable regime: is then time averaging (under convergence) matematically equivalent to averaging over the particles?

3. Divide the computation in two parts: first compute the fields up to the hole (skimmer), and use the results as boundary conditions for the part where I have a much smaller number of particles. Then change the nequivalentparticles in the second region.

4. The fields after the inlet seem reasonable and symmetric around the axis. However, due to the area-reduction of the inlet section near the line corresponding to radius=0, the number of inserted dsmc particles tends to 0 as well, reducing even further the number of dsmc particles that actually go across my hole. I'm uncertain if OpenFOAM is doing it in the right way or maybe the particle reduction close to the axis is too big.

The attached images show the inlet dsmc particles at different times of simulations, (em is the negative exponent ie 3_7_em6 = 3.7e-6)

5. Another option would be, considering that the mean free path after the hole is very small due to the minimum number of dsmc particles and the narrow velocity distribution, increase the size of the cells in the transversal direction. I have calculated that if the cell would be increased about 10 times, the number of dsmc particles should be enough. However I'm not sure if that cell would then not suffice the condition L<0.5*mean free path.

Have any of you worked in a similar problem? what do you do when you have regions with a much smaller number of dsmc particles? What strategy would you follow?

Any help will be very welcomed!

Attached Images
File Type: jpg screenshot_t_5_em11.jpg (42.9 KB, 30 views)
File Type: jpg screenshot_t_3_5_em9.jpg (52.7 KB, 25 views)
File Type: jpg screenshot_t_3_7_em6.jpg (67.8 KB, 24 views)

Last edited by Araist; June 25, 2015 at 07:57.
Araist is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Cluster ID's not contiguous in compute-nodes domain. ??? Shogan FLUENT 1 May 28, 2014 15:03
decomposePar pointfield flying OpenFOAM Running, Solving & CFD 28 December 30, 2013 15:05
No of simulated particles and No of real molecules in DSmc Simulation archeoptyrx OpenFOAM Running, Solving & CFD 0 June 19, 2013 15:55
AMI interDyMFoam for mixer danny123 OpenFOAM Running, Solving & CFD 4 June 19, 2013 04:49
Coal combustion. Particles number at inlet. Georg CFX 2 November 29, 2008 08:33

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