CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Explicit treatment of drag term

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

LinkBack Thread Tools Display Modes
Old   March 1, 2016, 07:30
Default Explicit treatment of drag term
New Member
Join Date: Jun 2014
Location: Oslo
Posts: 29
Rep Power: 5
pakanatiakash is on a distinguished road
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?

pakanatiakash is offline   Reply With Quote


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
wrong SU2 calculation for lift and drag coefficient for NAC4421 mechy SU2 7 January 9, 2017 06: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 08: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 13:19

All times are GMT -4. The time now is 04:55.