[Solved] How can I change the density of fluid phase in "cfdemSolverIB"??

Submitted by bhwang on Tue, 02/19/2019 - 04:15

Hi, I'm the beginner for using "cfdemSolverIB", and I want to know changing 'density' of fluid phase.
(I'm using the source code version (3.5.0), not the workbench version.)
In case of 'non-resolved CFD-DEM' (such as cfdemSolverPIso, ...), the 'rho' file is included in the '0' folder in CFD.
However, there is no data for density of fluid phase in the tutorial 'twospheresGlowinski'.

How can I change the density of fluid phase in "cfdemSolverIB"??
Thank you

Hwang

alice's picture

alice | Mon, 02/25/2019 - 14:07

Hello Hwang,

like in pressure based OpenFOAM solvers the density is not used as an explicit quantity, i.e. the pressure in the p-file is actually pressure/density and the viscosity is chosen accordingly in the transportProperties file.

Cheers,

Alice

bhwang | Mon, 03/04/2019 - 19:13

and, what I really know was that how to change the density of fluid.
And in my case (CFDEM 3.5.1, not the workbench version), the solver 'cfdemSolverIB' has a "fixed density" in the 'createFields.H'.
=====================================================
Info<< "\nCreating dummy density field rho = 1\n" << endl;
volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("0", dimensionSet(1, -3, 0, 0, 0), 1.0)
);

==============================================

So, I can change the density of fluid part in this file (createFields.H).
Thanks again, and I changed the title by adding "[solved]".

achuth1992 | Sat, 03/09/2019 - 18:29

Insert a rho file in the 0 folder of the case. This will do what is required.

Bests
Achuth

bhwang | Tue, 03/12/2019 - 08:43

if then, the last line (with comma) in the "createFields.H" shoul be commented such as,

=====================================================
Info<< "\nCreating dummy density field rho = 1\n" << endl;
volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
mesh //,
//dimensionedScalar("0", dimensionSet(1, -3, 0, 0, 0), 1.0)
);

==============================================
Thanks.

Hwang