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

maxwellianThermal wall interaction

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 2, 2011, 17:17
Default maxwellianThermal wall interaction
  #1
Member
 
Daniel
Join Date: Jul 2010
Location: California
Posts: 39
Rep Power: 7
hyperion is on a distinguished road
could someone explain the final expression for the reflected velocity vector in the diffuse wall interaction model in dsmcFOAM. It looks like it's supposed to be the most probable molecular thermal speed times a random vector, but the most probable molecular thermal speed is (2kT/m)^1/2. Looks like they left out the 2??? Thanks.

00109 U = 00110 sqrt(CloudType::kb*T/mass) 00111 *( 00112 rndGen.GaussNormal()*tw1 00113 + rndGen.GaussNormal()*tw2 00114 - sqrt(-2.0*log(max(1 - rndGen.scalar01(), VSMALL)))*nw 00115 );
hyperion is offline   Reply With Quote

Old   June 12, 2011, 10:25
Default
  #2
New Member
 
Fabian
Join Date: Jan 2010
Posts: 27
Rep Power: 7
gautcho is on a distinguished road
/E: Problem 2 is solved. Multiplication with the inverse of the factor scales a gaussNormal-dist exaktly to the Maxwellian-dist.



I was also thinking about this these days and I am also not so sure about it, but I have some thoughts, we can discuss over:

1. The first expression equals sqrt(RT), which is the inverse of the factor of the Maxwellian distribution, which I think the name of this model comes from.
2. I THINK there is a try to implement a randomly generated velocity based on the Maxwellian dist. But there are some problems:
2a. The first two expressions based on gaussNormal is a N(0,1)-dist with a factor sqrt(RT), but it should be a N(0, sqrt(RT))-dist in my opinion.
2b. The last expression is in fact also the gaussNormal-dist. This is the formula you can get gaussNormal numbers, BUT there is an pi missing and also the same problem as in 2a (my prof helped me with that).
2c. The same problems took place at the initialise-method in dsmcCloud.C. Perhaps you can have a look. The only difference is, that there is no thing like the 3rd expression...
3. I have found these formulas in Shen's excellente book "Rarefied gas dynamics: fundamentals, simulations and micro flows" from 2005 beginning on page 135. But to be honest: I didnt get all the points, I think...


So perhaps, we could figure out something together...

Last edited by gautcho; June 14, 2011 at 09:18.
gautcho is offline   Reply With Quote

Reply

Tags
diffuse, dsmcfoam, maxwellianthermal, wall interaction

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
tmerge utility creates unwanted interface/walls comes in the final mesh dinesh2n@gmail.com OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 11 January 20, 2012 07:23
Patches for OpenFOAM 1.7 on MacOS X gschaider OpenFOAM Installation on Windows, Mac and other Unsupported Platforms 101 September 21, 2011 05:37
UDF for wall slipping HFLUENT Fluent UDF and Scheme Programming 0 April 27, 2011 12:03
Question about bcdefw.f for wall temperature bc. Jimmy CD-adapco 10 March 18, 2008 16:28
X-Y plot of Yplus in Fluent 5.3 Luo FLUENT 24 April 11, 2000 06:07


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