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

How the momentum source term due to droplet evaporation is calculated?

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

LinkBack Thread Tools Search this Thread Display Modes
Old   July 3, 2020, 10:57
Default How the momentum source term due to droplet evaporation is calculated?
New Member
Join Date: Jul 2016
Posts: 11
Rep Power: 8
pappo1890 is on a distinguished road
Hi Foamers.

The momentum source terms due to the dispersed phase can be identified as those due to the mass transfer and the forces acting on the particle. In literature, this term is expressed as in the figure in the attachment (Experiments and Numerical Simulations of Diluted Spray Combustion, B. Merci, D. Roekaerts, A. Saidiki., p.17)

I would like to have feedback about this expression and the one is used in OpenFOAM.

In OpenFOAM the momentum source term is written in the UEqn.H as parcels.SU(U).
SU(volVectorField& u) const is defined in line 563 of ~/OpenFOAM2.4.0/src/lagrangianba/intermediate/clouds/Templates/KinematicCloud/KinematicCloudI.H:

tmp<fvVectorMatrix> tfvm(new fvVectorMatrix(U, dimForce));
fvVectorMatrix& fvm = tfvm();

fvm.source() = -UTrans()/(this->db().time().deltaT());

return tfvm;
UTrans() is updated twice when solving the lagrangian libraries. First in the KinematicCloud, for the momentum spray source term due to the force acting on the particle and later in parcels/Templates/ReactingParcel/ReactingParcel.C (line 689):

Code:[cellI] += dm*U0;
U0 is (line 579):
const vector& U0 = this->U_;
and at line 731 the this-> U_ pointer is calculated as the particle velocity.

The expression does not completely look like that found in literature, where the difference between the carrier and the particle velocity is considered. Could, please someone double-check if my reading of the following lines of code is correct?

Thanks in advance
Attached Images
File Type: jpg momsource.JPG (13.2 KB, 27 views)

Last edited by pappo1890; July 3, 2020 at 19:01.
pappo1890 is offline   Reply With Quote


evaporation, openfoam, spray source term, su(u), utrans

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
[] swak4foam compiling issues on a cluster saj216 OpenFOAM Installation 5 January 17, 2023 17:05
[Other] How to use finite area method in official OpenFOAM 2.2.0? Detian Liu OpenFOAM Meshing & Mesh Conversion 4 November 3, 2015 04:04
[] problem when installing foam-extend-1.6 Thomas pan OpenFOAM Installation 7 September 9, 2015 22:53
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 20:08
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 18:51

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