|
[Sponsors] |
BrowninanMotionForce of particles in laminar flow |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 28, 2016, 05:08 |
BrowninanMotionForce of particles in laminar flow
|
#1 |
Member
Bijan Darbari
Join Date: Nov 2015
Posts: 48
Rep Power: 10 |
Hello
In DPM model, Brownian force is one of the major forces, which act on a particles. Look at the following link ( line 194-197 ): https://github.com/OpenFOAM/OpenFOAM...nMotionForce.C For laminar regime, the Brownian force is defined by : Code:
const scalar rhoRatio = p.rho()/p.rhoc(); const scalar s0 = 216*muc*sigma*Tc/(sqr(mathematical::pi)*pow5(dp)*(rhoRatio)*cc); f = eta*sqrt(mathematical::pi*s0/dt); For laminar regime, the correlation for this force is presented in the following : note:where subscripts "f" and "p" denote the fluid and particle, respectively. Moreover, rho is for continuous phase. As shown in this figure, the square of rho_p (particle density) shoud appears at the equation of S0 but I had supposed that the power one of rho_p hed been appeared in the (S0) of the code. - Is this a bug in the code ?? ( If i don't mistake. ) - what is muc?? ( dynamic or kinematic viscosity of continuous phase ) Please kindly guide me. Last edited by bijan darbari; March 30, 2016 at 10:30. |
|
March 28, 2016, 10:29 |
muc
|
#2 |
New Member
Rahand
Join Date: Mar 2016
Posts: 10
Rep Power: 10 |
Just to your second question:
mu is usually the dynamic viscosity and muc in this cas the dynamic viscosity of the continous phase. |
|
March 30, 2016, 10:31 |
|
#3 |
Member
Bijan Darbari
Join Date: Nov 2015
Posts: 48
Rep Power: 10 |
Please kindly guide me.
|
|
April 17, 2016, 01:26 |
|
#4 |
Member
Bijan Darbari
Join Date: Nov 2015
Posts: 48
Rep Power: 10 |
for more information :
http://www.openfoam.org/mantisbt/view.php?id=2036 |
|
June 7, 2016, 07:37 |
|
#5 |
Member
Evangelos
Join Date: Sep 2011
Posts: 87
Rep Power: 14 |
I think that there are some mistakes in the code of BrownianMotionForce.C
a) dt in the code is not the time step (as the reference paper) of the simulation. b) the Boltzmann constant is in eVK^-1 but the reference paper uses the constant in J/K !!! The results from the simulation (with the values above) seems to be logical |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
CFX Treatment of Laminar and Turbulent Flows | Jade M | CFX | 18 | September 15, 2022 07:08 |
Laminar Isothermal Flow in a duct | HectorRedal | Main CFD Forum | 29 | June 2, 2012 07:04 |
Optimize for laminar flow, assume it valid for turbulent flow? | Chander | CFX | 15 | November 6, 2011 05:06 |
Optimize for laminar flow, assume it valid for turbulent flow? | Chander | Main CFD Forum | 2 | October 24, 2011 07:43 |
Validation for laminar Disperse phase flow | shashwat | Main CFD Forum | 0 | April 4, 2008 02:20 |