How to evaluate the force of particles to the fluid in Lagrangian method

marialhm · November 8, 2018, 04:09

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, 10: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, 11: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, 21:48

Quote:

Originally Posted by NewKid

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, 00: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, 21:46

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

NewKid · November 13, 2018, 22:28

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

November 8, 2018, 04:09	How to evaluate the force of particles to the fluid in Lagrangian method	#1
marialhm Member Maria Join Date: Jul 2013 Posts: 84 Rep Power: 12	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

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Physics of an incompressible fluid	Jonas Holdeman	Main CFD Forum	87	April 22, 2019 10:40
fluid flow but temperature raises from nowhere?	mxcfd	STAR-CCM+	5	September 16, 2014 05:46
which solver for particle tracking through a fluid for lagrangian based method?	ranasa	OpenFOAM Running, Solving & CFD	0	August 9, 2014 03:09
Volume of Fluid method for multiple fluid flows	j.patel	Main CFD Forum	6	March 7, 2014 22:13
Terrible Mistake In Fluid Dynamics History	Abhi	Main CFD Forum	12	July 8, 2002 09:11

November 8, 2018, 10:41		#2
NewKid Senior Member LT Join Date: Dec 2010 Posts: 104 Rep Power: 15	Hi, maria, regarding your thread, seems what I do is similar to your issue. How is your case going now?

November 8, 2018, 11:04		#3
NewKid Senior Member LT Join Date: Dec 2010 Posts: 104 Rep Power: 15	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?

November 9, 2018, 00:50		#5
NewKid Senior Member LT Join Date: Dec 2010 Posts: 104 Rep Power: 15	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.

November 13, 2018, 21:46		#6
NewKid Senior Member LT Join Date: Dec 2010 Posts: 104 Rep Power: 15	Hi, Maria, do you understand dimensionedField and 'dimForce' in KinematicCloudI.H?

November 13, 2018, 22:28		#7
NewKid Senior Member LT Join Date: Dec 2010 Posts: 104 Rep Power: 15	fvOptions seems another option to hold the source? https://openfoam.org/release/2-2-0/fv-options/