|
[Sponsors] |
September 15, 2022, 10:22 |
modified stochasticDispersionRAS
|
#1 |
New Member
Charles Guaquiere
Join Date: Sep 2022
Posts: 14
Rep Power: 3 |
Hello,
I try to modify stochastic dispersion model on Open FOAM to improved my results concerning particles deposition. I found a tutorial from Chalmers University: http://www.tfd.chalmers.se/~hani/kur..._Xu_Report.pdf I changed the stochasticDispersionRAS.C and the stochasticDispersionRAS.H like in the tutorial but it used Open FOAM 4 and I use Open FOAM 9. My problem is on the following part of the code (from line 96 to 103): const word turbName= IOobject::groupName ( turbulenceModel:: propertiesName, this->owner().U().group() ); const turbulenceModel& turbModel = obr.lookupObject<turbulenceModel>(turbName); I changed all "turbulenceModel" by "momentumTransportModel" concidering the modifications between Open FOAM 4 et 9 but "propertiesName" is no longer defined in momentumTransportModel.C and momentumTransportModel.H like it was on Open FOAM 4. I don't know how to modify this part of the code to make it works on Open FOAM 9. I don't have enough experience developing code on OpenFoam yet, so please can someone help me ? Kind regards, |
|
September 23, 2022, 11:20 |
|
#2 |
Senior Member
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 147
Rep Power: 17 |
Any idea where to find some helpful pointers on this particular point ?
|
|
September 23, 2022, 22:05 |
|
#3 | |
New Member
Chen Xiaoxiao
Join Date: Jun 2018
Location: China
Posts: 6
Rep Power: 7 |
Quote:
so you can use the following simplified form: Code:
const word turbName= IOobject::groupName ( "momentumTransport.<phase name>" ); |
||
September 26, 2022, 05:32 |
|
#4 |
New Member
Charles Guaquiere
Join Date: Sep 2022
Posts: 14
Rep Power: 3 |
Thanks for your reply,
Unfortunatly it does not work, are you sure it is momentumTransport.<phase name> and not momentumTransportModel.<phase name> ? Because OpenFOAM suggested me to use "momentumTransportModel_H" instead of "momentumTransport" when I compile the code. Also, how can I find my phase name to replace in your simplified form because "phase name" is not declared. Thank you in advance for your answer |
|
September 26, 2022, 06:56 |
|
#5 | |
New Member
Chen Xiaoxiao
Join Date: Jun 2018
Location: China
Posts: 6
Rep Power: 7 |
Quote:
|
||
September 26, 2022, 08:26 |
|
#6 |
New Member
Charles Guaquiere
Join Date: Sep 2022
Posts: 14
Rep Power: 3 |
Sorry it is my first time on CFD Online.
I think I found the solution by compararing files like shearStress.C or yPlus.C, so I changed my previous code into: Code:
const word turbName= IOobject::groupName ( momentumTransportModel::typeName, phaseName_ ); Code:
yP_ ( IOobject ( "yP", this->owner_.mesh().time().typeName(), this->owner_.mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), this->owner_.mesh(), dimensionedScalar("yP", dimless, 0.0) ) error_message.jpg i'm not sure what that means. |
|
September 26, 2022, 08:33 |
|
#7 | |
New Member
Chen Xiaoxiao
Join Date: Jun 2018
Location: China
Posts: 6
Rep Power: 7 |
Quote:
|
||
September 26, 2022, 09:45 |
|
#8 |
New Member
Charles Guaquiere
Join Date: Sep 2022
Posts: 14
Rep Power: 3 |
Do you have any idea what i can modify to make it works on OpenFOAM 9 ? I tried modifications but nothing works.
thank you in advance |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[General] Extracting ParaView Data into Python Arrays | Jeffzda | ParaView | 30 | November 6, 2023 21:00 |
Problem while compiling modified simpleFoam solver | mattiafc | OpenFOAM Programming & Development | 0 | May 26, 2020 20:01 |
why in solve(UEqn == -fvc::grad(p)) UEqn.source() is not modified and UEqn.H() chang | mAlletto | OpenFOAM | 1 | January 1, 2019 11:59 |
Modified rhoCentralFoam slip boundary fails in parallel | ChrisA | OpenFOAM Programming & Development | 0 | June 25, 2014 20:28 |
Modified Equation for CFX algorithm | Craig Johansen | CFX | 0 | August 27, 2004 23:02 |