|
[Sponsors] |
March 1, 2016, 06:30 |
Explicit treatment of drag term
|
#1 |
New Member
Akash
Join Date: Jun 2014
Location: Oslo
Posts: 29
Rep Power: 11 |
Hi All
I am modelling solidification process and so far constructed a solver with phase change with fixed solid. Now I am interested in solid movement (solid grain motion ) and need to model drag term. The source term is written as follows - Mld = (3*gs*rhom*Cd)/(4*dg) *mag(Ul-Us)* (Ul-Us) gs - solid fraction rhom - mixture density of solid and liquid Cd - drag coefficient dg - diameter of grain Ul - liquid velocity Us - solid velocity I am pasting the partial code for UlEqn.H pertaining to the drag term volScalarField drag("drag", ( (3.*gs)/(4.*dg) )*rhom*Cd); and in fvVectorMatrix UlEqn + fvm::SuSp(drag*mag(Ul-Us), Ul) - fvc::SuSp(drag*mag(Ul-Us), Us) since I am placing mag(Ul-Us) in the fvm part, the calculation is being done implicitly leading to prohibitively slower calculation. The trick would be to treat the term mag(Ul-Us) explicitly. This could be done using fvc::, but I am not sure of how to do this. Can anyone suggest me how to proceed with this? Cheers Akash |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
wrong SU2 calculation for lift and drag coefficient for NAC4421 | mechy | SU2 | 7 | January 9, 2017 05:18 |
Drag coefficient_cylinder | B.Stas | FloEFD, FloWorks & FloTHERM | 1 | August 26, 2014 07:27 |
Pressure drag calculation | lc05 | Main CFD Forum | 2 | November 1, 2010 07:50 |
ATTENTION! Reliability problems in CFX 5.7 | Joseph | CFX | 14 | April 20, 2010 15:45 |
Drag coefficient for parcels in dieselFoam | sebastian_vogl | OpenFOAM Running, Solving & CFD | 5 | December 31, 2008 12:19 |