lagrangian cloud::SU() returns different
Hi all,I noticed the following when coupling lagrangian cloud with eulerian continuum. In the momentum source calculator function, MomentumCloud::SU() (for OF8+, otherwise it would probably be KinematicCloud::SU()), we see the following code:
Code:
if (solution_.semiImplicit("U")) Where UCoeff and UTrans are the self-explanatory fields, derived by adding up the interaction from all member particles of the cloud. The returned object is an fvVectorMatrix, with the first and last are explicit components, and the middle term, the coefficient term. Now, if there is no "semi-implicit" flag, ie the source is purely an explicit term, the code returns: Code:
else Which is quite analogous to the explicit part of semi-implicit source above, but it does not have a division by cell volume. It appears to me that in general, we supply fvVectorMatrix terms on a per unit volume basis. What is then the reason behind difference in the divisor in the two cases. Is the second formulation correct? If it is, can someone help me understand why it is so? I would expect that the two should not even run with the same solver/case, based on dimensional considerations, but they do. However, I cannot find any other point in lagrangian libraries, where the semi-implicit flag is used to change the way in which UTrans is calculated differently depending of the flag. I would appreciate any help in enlightening me about this. Thanks |
In the explicit formulation, it doesn't construct fvMatrix but assigns the field to the source of the fvMatrix, hence no division by volume.
|
Hey,
Thank you for your answer! :) You say that Quote:
but, in the explicit case's code above doesn't Code:
tmp<fvVectorMatrix> tfvm(new fvVectorMatrix(U, dimForce)); Thanks again. |
Hi,
There is the section (3.4.9, P-43) about the source terms in Programmer's Guide. https://sourceforge.net/projects/openfoam/files/v2112/ With the implicit method, the term will be integrated over a control volume. This explanation is about a flow field, not about a lagrangian cloud. However, there may be some relations with your question. |
All times are GMT -4. The time now is 08:26. |