CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Verification & Validation

Validating Brownian motion in OpenFOAM

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

Like Tree1Likes
  • 1 Post By raptorishungry

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 22, 2019, 20:05
Default Validating Brownian motion in OpenFOAM
  #1
New Member
 
Adithya Gurumurthy
Join Date: Jun 2017
Posts: 17
Rep Power: 4
raptorishungry is on a distinguished road
Hello Foamers!

I am trying to validate OpenFOAM's Brownian motion force. For this purpose, I am following the process highlighted in this link .
HTML Code:
https://bugs.openfoam.org/view.php?id=2153
I understand that the source code is updated by henry. But I am having trouble validating it, I tried all possible combinations (changing particle density, carrier fluid, mean free path, temperature of the fluid/particle) but unfortunately not able to pull it off. Also I do not understand what henry means in the last comment of this post
HTML Code:
https://bugs.openfoam.org/view.php?id=2153
where he says "RMS of the magnitude of the displacement in the spherical distribution.", as of now I think solving that ambiguity might help me get closer to validating it.

Could someone please shed some light on this?

Thanks,
adithya
raptorishungry is offline   Reply With Quote

Old   March 6, 2019, 13:51
Default
  #2
New Member
 
Adithya Gurumurthy
Join Date: Jun 2017
Posts: 17
Rep Power: 4
raptorishungry is on a distinguished road
Hey all, I found out the solution to this after performing quite a lot of computational runs and varying constants here and there to narrow down on the discrepancy. The results will not match if you follow this thread
HTML Code:
https://bugs.openfoam.org/view.php?id=2153
This is because the sphereDrag in OpenFOAM needs to be corrected by the Cunningham slip correction factor. I modified OpenFOAM's sphereDrag to a new force following the formulation in ANSYS for Stokes-Cunningham drag
HTML Code:
https://www.sharcnet.ca/Software/Ansys/16.2.3/en-us/help/flu_th/x1-62400016.7.3.html
I also modified the BrownianMotionFoce in OpenFOAM according to Omid et al.
HTML Code:
https://www.tandfonline.com/doi/pdf/10.1080/02786820802587888
Making these changes validated Brownian motion force in OpenFOAM.
Daniel_Khazaei likes this.
raptorishungry is offline   Reply With Quote

Old   September 4, 2019, 10:25
Default j have problem in validating brownian force openfoam
  #3
New Member
 
Join Date: Jun 2019
Posts: 1
Rep Power: 0
jafarbakhshi1374 is on a distinguished road
hello dear
Adithya Gurumurthy
i have the same problem and my results have about 50 precent deviation by exact answere (sqrt(2Dt)) so would you please send the changes y have done to validate this force?
jafarbakhshi1374 is offline   Reply With Quote

Old   October 18, 2019, 14:45
Default
  #4
Member
 
Join Date: Jun 2016
Posts: 35
Rep Power: 5
xuegy is on a distinguished road
Hello Adithya

Could you please share your modification? I also have problem validating the Brownian motion force. Seems like the code from OpenFOAM is always smaller by sqrt(2).
xuegy is offline   Reply With Quote

Old   October 18, 2019, 18:00
Default
  #5
New Member
 
Adithya Gurumurthy
Join Date: Jun 2017
Posts: 17
Rep Power: 4
raptorishungry is on a distinguished road
Hello jafarbakshi and xuegy, sorry for the late reply. I have attached my files below.
Attached Files
File Type: c cDrag.C (2.7 KB, 3 views)
File Type: h cDrag.H (3.3 KB, 0 views)
File Type: c myBMF.C (7.3 KB, 2 views)
File Type: h myBMF.H (4.3 KB, 0 views)
File Type: h myBMFI.H (1.5 KB, 0 views)
raptorishungry is offline   Reply With Quote

Old   October 21, 2019, 13:09
Default
  #6
Member
 
Join Date: Jun 2016
Posts: 35
Rep Power: 5
xuegy is on a distinguished road
Quote:
Originally Posted by raptorishungry View Post
Hello jafarbakshi and xuegy, sorry for the late reply. I have attached my files below.
I just tried your code. The mathematics equation are basically the same as OpenFOAM cubic distribution(comment in the source code). It has the same problem verifying the Brownian motion. It's always diffusing sqrt(2) times slower than the continuum model.
xuegy is offline   Reply With Quote

Old   October 21, 2019, 17:27
Default
  #7
New Member
 
Adithya Gurumurthy
Join Date: Jun 2017
Posts: 17
Rep Power: 4
raptorishungry is on a distinguished road
It worked for me. What are you validating your model with? could you post your case set-up/plots?
raptorishungry is offline   Reply With Quote

Old   October 22, 2019, 14:28
Default
  #8
Member
 
Join Date: Jun 2016
Posts: 35
Rep Power: 5
xuegy is on a distinguished road
Quote:
Originally Posted by raptorishungry View Post
It worked for me. What are you validating your model with? could you post your case set-up/plots?
I'm trying the easiest 1D case. A 2um*2um*20um box(1*1*100 mesh) is filled with water. Half of the box has 100nm spherical particles initially (dimensionless concentration "2") while another half is empty (0). The continuum model solves the diffusion equation while the discrete model simulates Brownian motion and converts the Lagrangian field to the void fraction.

My results shows that the discrete model is always slower by sqrt(2), e.g. the result of discrete model at t=11.3 s converges to the result of continuum model at t=16 s.

0.0.png

Please correct me if I make any mistake!
xuegy is offline   Reply With Quote

Reply

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
External management of solid motion (using dynamicmesh) maxou1993 Main CFD Forum 0 July 28, 2015 12:37
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
How to incorporate facezone motion using dynamicmeshmotion in openfoam? Aj Nair OpenFOAM Programming & Development 1 March 16, 2014 00:12
Automatic Mesh Motion solver michele OpenFOAM Running, Solving & CFD 10 September 26, 2005 09:21
Brownian motion Sank Main CFD Forum 7 September 11, 2004 17:36


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