CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Drag Model for twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 23, 2020, 04:21
Default Drag Model for twoPhaseEulerFoam
  #1
New Member
 
Patrizio
Join Date: Feb 2019
Location: Paris
Posts: 5
Rep Power: 7
lollico is on a distinguished road
Good morning everyone,
I would like to know if it is possibile to use the "PlessisMasliyahDrag particleForce" in an eulerian-eulerian solver like twophaseeulerfoam. The goal is to find the best drag model to match experimental data of metal particle injected in a convergent divergent supersonic nozzle (the continuum media is gas).
Otherwise are there other suggestions about other drag model?
Thank you.
Patrizio
lollico is offline   Reply With Quote

Old   September 29, 2020, 06:19
Default
  #2
Senior Member
 
Join Date: Dec 2019
Location: Cologne, Germany
Posts: 355
Rep Power: 8
geth03 is on a distinguished road
if you go to
solvers -> multiphase -> twoPhaseEulerFoam -> interfacialModels ->
dragModels,
you will see all the available drag models for the twoPhaseEulerFoam solver.

The Plessis Masliyah Drag Model is not available, but you could create a copy of the Ergun drag model and then replace the constants (150 and 1.75) with the mathematical expressions of the Plessis-Masliyah Model.
geth03 is offline   Reply With Quote

Old   September 29, 2020, 06:56
Default
  #3
New Member
 
Patrizio
Join Date: Feb 2019
Location: Paris
Posts: 5
Rep Power: 7
lollico is on a distinguished road
Quote:
Originally Posted by geth03 View Post
if you go to
solvers -> multiphase -> twoPhaseEulerFoam -> interfacialModels ->
dragModels,
you will see all the available drag models for the twoPhaseEulerFoam solver.

The Plessis Masliyah Drag Model is not available, but you could create a copy of the Ergun drag model and then replace the constants (150 and 1.75) with the mathematical expressions of the Plessis-Masliyah Model.

Thank you very much, I will try this!
lollico is offline   Reply With Quote

Old   October 1, 2020, 09:58
Default Drag Model for twoPhaseEulerFoam
  #4
New Member
 
Patrizio
Join Date: Feb 2019
Location: Paris
Posts: 5
Rep Power: 7
lollico is on a distinguished road
Hi,
Following your suggestion (I am not an expert in C++) I changed Ergun.c from :
________________________________

(4.0/3.0)
*(
150
*max
(
scalar(1) - pair_.continuous(),
pair_.continuous().residualAlpha()
)/max(pair_.continuous(), pair_.continuous().residualAlpha())
+ 1.75
*pair_.Re()
);
_____________________________________
to:

return
(4.0/3.0)
*(
( 26.8*pow3(pair_.continuous())
/(
sqr(cbrt(scalar(1) - pair_.continuous()))
*(1.0 - cbrt(scalar(1) - pair_.continuous()))
*sqr(1.0 - sqr(cbrt(scalar(1) - pair_.continuous())))))

*max
(
scalar(1) - pair_.continuous(),
pair_.continuous().residualAlpha()
)/max(pair_.continuous(), pair_.continuous().residualAlpha())
+ (sqr(pair_.continuous())
/sqr(1.0 - sqr(cbrt(scalar(1) - pair_.continuous()))))

*pair_.Re()
);



where A=
(26.8*pow3(pair_.continuous())
/(
sqr(cbrt(scalar(1) - pair_.continuous()))
*(1.0 - cbrt(scalar(1) - pair_.continuous()))
*sqr(1.0 - sqr(cbrt(scalar(1) - pair_.continuous())))
))

And B=

sqr(pair_.continuous())
/sqr(1.0 - sqr(cbrt(scalar(1) - pair_.continuous())));


__________________________________________________ ____


I don't know if I well understood the variables, but it should work.
There is something I should change in the .H files too ?
Many thanks.
Cordially,
Patrizio
lollico is offline   Reply With Quote

Old   October 2, 2020, 04:51
Default
  #5
Senior Member
 
Join Date: Dec 2019
Location: Cologne, Germany
Posts: 355
Rep Power: 8
geth03 is on a distinguished road
you can't just change the equation and expect the code to work.

to make it clear i will write down the steps you need to follow:

1. create a copy of the solver you want to modify and make changes only there.
if you don't know how to do it, check out the 3 weeks series, there is a tutorial on how to do that.

2. if you want to add a drag model, copy one and change
i) the name of the directory, the names in the .C and .H - files.
ii) change the equation

3. add your new model to the files-file in the Make-directory. look how it is done for other models and add your model.

4. compile your modified stuff whether it is the solver or the library.
geth03 is offline   Reply With Quote

Old   November 23, 2020, 05:57
Default
  #6
New Member
 
Patrizio
Join Date: Feb 2019
Location: Paris
Posts: 5
Rep Power: 7
lollico is on a distinguished road
Thank you, following your instructions I succeeded to extend the lagrangian model to the twophaseeulerfoam euler-euler model.
Have a nice day!





Quote:
Originally Posted by geth03 View Post
you can't just change the equation and expect the code to work.

to make it clear i will write down the steps you need to follow:

1. create a copy of the solver you want to modify and make changes only there.
if you don't know how to do it, check out the 3 weeks series, there is a tutorial on how to do that.

2. if you want to add a drag model, copy one and change
i) the name of the directory, the names in the .C and .H - files.
ii) change the equation

3. add your new model to the files-file in the Make-directory. look how it is done for other models and add your model.

4. compile your modified stuff whether it is the solver or the library.
lollico is offline   Reply With Quote

Reply

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
Finding Drag Force from Skin Friction Rob Wilk Main CFD Forum 0 May 8, 2020 07:04
How to compile a new drag force in Lagrangian model lilinmin OpenFOAM Programming & Development 4 June 9, 2016 05:15
Drag model implementation vbchris OpenFOAM Running, Solving & CFD 1 March 31, 2015 09:28
about Subgrid-scale model impecca OpenFOAM Running, Solving & CFD 4 December 20, 2013 11:36
multiphase model and drag law Yasmail AKARIOUH FLUENT 0 April 29, 2008 08:44


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