CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   How to evaluate the force of particles to the fluid in Lagrangian method (https://www.cfd-online.com/Forums/openfoam-programming-development/210946-how-evaluate-force-particles-fluid-lagrangian-method.html)

marialhm November 8, 2018 05:09

How to evaluate the force of particles to the fluid in Lagrangian method
 
Dear all,

I'm using a Eulerian-Lagrangian solver to simulate gas-solid flow. And in the UEqn.H file of the solver, there is a term of parcels.SU(U). I know that this term is used to expand the momentum equation for two-way coupling. But how to evaluate this term?

Thanks in advance.

Maria

NewKid November 8, 2018 11:41

Hi, maria, regarding your thread, seems what I do is similar to your issue.
How is your case going now?

NewKid November 8, 2018 12:04

Hi, I just read the source code of kinematicCloudI.H which has the implementation of SU(volVectorField& U).
In the implementation of SU(), a switch named coupled determines the two-way coupling works indeed.
What does Vdt() mean in this function?

marialhm November 8, 2018 22:48

1 Attachment(s)
Quote:

Originally Posted by NewKid (Post 714624)
Hi, I just read the source code of kinematicCloudI.H which has the implementation of SU(volVectorField& U).
In the implementation of SU(), a switch named coupled determines the two-way coupling works indeed.
What does Vdt() mean in this function?

I also have the same findings, the equation for SU is shown in the attached figure. But I don't know the meaning of each term, do you have any clue?

About Vdt, I think it is related to the time step, but don't know the specific meaning either...

NewKid November 9, 2018 01:50

Hi, maria, nice to see your response.
I have no clue to the terms in the implementation of SU() in kinematicCloudI.H, in fact I don't know how the coupling works.
But I have a model to describe the interactions between the gas and the liquid, what I need to do is to write a new function/routine in kinematicCloudI.H to implement my model, and replace the parcels.SU(U) with my new function.
I do hope it's a successful attempt.

NewKid November 13, 2018 22:46

Hi, Maria, do you understand dimensionedField and 'dimForce' in KinematicCloudI.H?

NewKid November 13, 2018 23:28

fvOptions seems another option to hold the source?
https://openfoam.org/release/2-2-0/fv-options/


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