CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   pimpleDyMFoam - rho not found (https://www.cfd-online.com/Forums/openfoam-solving/185904-pimpledymfoam-rho-not-found.html)

maasyraf3 April 5, 2017 07:42

pimpleDyMFoam - rho not found
 
1 Attachment(s)
Hi all - very good day.

I'd like to get your help here guys to look into this..

After run it, appear "cannot find rho".

below is the comments

theuser@oscae03:~/OpenFOAM/theuser-4.1/run/tutorials/incompressible/pimpleDyMFoam/final_test17$ pimpleDyMFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.1
Exec : pimpleDyMFoam
Date : Apr 06 2017
Time : 03:05:57
Host : "oscae03"
PID : 19972
Case : /home/theuser/OpenFOAM/theuser-4.1/run/tutorials/incompressible/pimpleDyMFoam/final_test17
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: displacementLaplacian
Selecting motion diffusion: inverseDistance

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
No MRF models present

Reading/calculating face velocity Uf

No finite volume options present

Courant Number mean: 3.94318e-05 max: 0.000470908
--> FOAM Warning :
From function void Foam::timeControl::read(const Foam::dictionary&)
in file db/functionObjects/timeControl/timeControl.C at line 89
Reading "/home/theuser/OpenFOAM/theuser-4.1/run/tutorials/incompressible/pimpleDyMFoam/final_test17/system/controlDict.functions.probes" from line 56 to line 70
Using deprecated 'outputControl'
Please use 'writeControl' with 'writeInterval'
--> FOAM Warning :
From function void Foam::timeControl::read(const Foam::dictionary&)
in file db/functionObjects/timeControl/timeControl.C at line 89
Reading "/home/theuser/OpenFOAM/theuser-4.1/run/tutorials/incompressible/pimpleDyMFoam/final_test17/system/controlDict.functions.forces" from line 75 to line 96
Using deprecated 'outputControl'
Please use 'writeControl' with 'writeInterval'
forces forces:
Not including porosity effects
forceCoeffs forces:
Not including porosity effects
--> FOAM Warning :
From function void Foam::timeControl::read(const Foam::dictionary&)
in file db/functionObjects/timeControl/timeControl.C at line 89
Reading "/home/theuser/OpenFOAM/theuser-4.1/run/tutorials/incompressible/pimpleDyMFoam/final_test17/system/controlDict.functions.fieldAverage1" from line 101 to line 120
Using deprecated 'outputControl'
Please use 'writeControl' with 'writeInterval'
fieldAverage fieldAverage1:
Starting averaging at time 0


Starting time loop

Courant Number mean: 3.94298e-05 max: 0.000470884
deltaT = 5.99952e-06
Time = 5.99952e-06

GAMG: Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 0

GAMG: Solving for cellDisplacementy, Initial residual = 1, Final residual = 6.39574e-06, No Iterations 5

GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0114705, No Iterations 7
time step continuity errors : sum local = 2.53103e-09, global = -3.6467e-18, cumulative = -3.6467e-18

PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 5.61624e-08, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 6.19727e-08, No Iterations 1

GAMG: Solving for p, Initial residual = 1, Final residual = 0.00913275, No Iterations 8
time step continuity errors : sum local = 2.09605e-09, global = -3.64661e-18, cumulative = -7.29331e-18

PIMPLE: iteration 2
smoothSolver: Solving for Ux, Initial residual = 1.62711e-05, Final residual = 4.687e-13, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1.97001e-05, Final residual = 4.24567e-13, No Iterations 1

GAMG: Solving for p, Initial residual = 0.419213, Final residual = 7.03197e-07, No Iterations 28
time step continuity errors : sum local = 5.63715e-13, global = -3.64672e-18, cumulative = -1.094e-17
ExecutionTime = 0.11 s ClockTime = 1 s



--> FOAM FATAL ERROR:
Could not find rho

From function void Foam::functionObjects::forces::initialise()
in file forces/forces.C at line 196.

FOAM exiting

KateEisenhower April 5, 2017 08:22

Hi,

I think this error is related to a wrong entry in the functions section in your controlDict. What does it look like?

Rho should be defined in the following manner:

Code:

rho rhoInf;
rhoInf 1.199;

Best regards,

Kate

maasyraf3 April 5, 2017 08:43

Quote:

Originally Posted by KateEisenhower (Post 643773)
Hi,

I think this error is related to a wrong entry in the functions section in your controlDict. What does it look like?

Rho should be defined in the following manner:

Code:

rho rhoInf;
rhoInf 1.199;

Best regards,

Kate


Hi Kate- I already did the correct way but still happen the error.

The file controlDict

Quote:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application pimpleDyMFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 40;

deltaT 5e-06;

writeControl adjustableRunTime;

writeInterval 0.1;

purgeWrite 0;

writeFormat binary;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep yes;

maxCo 0.5;

functions
{
probes
{
type probes;
functionObjectLibs ("libsampling.so");
enabled true;
outputControl timeStep;
outputInterval 1;
probeLocations
(
( 0.0 0.0 0.0 )
( 10.0 0.0 0.0 )

);

fields
(
p
);
}

forces
{
type forceCoeffs;
functionObjectLibs ( "libforces.so" );
outputControl timeStep;
outputInterval 1;
patches
(
cylinder
);
directForceDensity no;

pName p;
UName U;
rhoName rhoInf;
//log true;
rhoInf 1;//1000; //994.5;
CofR (0 0 0); //( 0.080415 0 0 );
liftDir ( 0 1 0 );
dragDir ( 1 0 0 );
pitchAxis ( 0 0 1 );
magUInf 1.0; //0.1;
lRef 1.0; //0.01787;
Aref 1.0; //0.0003193369;
}

fieldAverage1
{
type fieldAverage;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;
fields
(
U
{
mean on;
prime2Mean on;
base time;
}

p
{
mean on;
prime2Mean on;
base time;
}
);
}

// ************************************************** *********************** //

maasyraf3 April 5, 2017 08:50

Hi kate , is it something to do with my dynamic meshing?

Quote:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solver displacementLaplacian;


displacementLaplacianCoeffs
{
diffusivity inverseDistance 1(cylinder);
}





// ************************************************** *********************** //

maasyraf3 April 5, 2017 18:49

My case file -simulate two cylinders in tandem in pimpleDyMFoam
 
The p

Quote:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{

topAndBottom
{
type zeroGradient;
}

outlet
{
type zeroGradient;
}

inlet
{
type zeroGradient;
}

frontAndBack
{
type empty;
}

cylinder
{
type zeroGradient;
}

}

// ************************************************** *********************** //

U

Quote:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (2 0 0);

boundaryField
{
topAndBottom
{
type slip;
}

inlet
{
type fixedValue;
value uniform (2 0 0);
}

outlet
{
type zeroGradient;
}

frontAndBack
{
type empty;
}


cylinder
{
type movingWallVelocity;
value uniform (0 0 0);
}

}

// ************************************************** *********************** //


pointDisplacement


Quote:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class pointVectorField;
location "0.01";
object pointDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 0 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
topAndBottom
{
type uniformFixedValue;
uniformValue (0 0 0);
}

inlet
{
type uniformFixedValue;
uniformValue (0 0 0);
}

outlet
{
type uniformFixedValue;
uniformValue (0 0 0);
}

cylinder
{
type oscillatingDisplacement;
omega 1.256;
amplitude (0 0.5 0);
value uniform (0 0 0);
}

frontAndBack
{
type empty;
}
}


blockMeshDict




Quote:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.1;

vertices
(
//z=0
//inner circle1
(5 0 0) //0
(3 3 0) //1
(0 5 0) //2
(-3 3 0) //3
(-5 0 0) //4
(-3 -3 0) //5
(0 -5 0) //6
(3 -3 0) //7

//outer circle1
(15 0 0)//8
(10 10 0)//9
(0 15 0)//10
(-10 10 0)//11
(-15 0 0)//12
(-10 -10 0)//13
(0 -15 0)//14
(10 -10 0)//15

//rectangle1
(75 0 0) //16
(75 10 0) //17
(75 50 0) //18
(10 50 0) //19
(0 50 0) //20
(-10 50 0) //21
(-25 50 0) //22
(-25 10 0) //23
(-25 0 0) //24
(-25 -10 0) //25
(-25 -50 0) //26
(-10 -50 0) //27
(0 -50 0) //28
(10 -50 0) //29
(75 -50 0) //30
(75 -10 0) //31


//z=1
//inner circle1
(5 0 1) //32
(3 3 1) //33
(0 5 1) //34
(-3 3 1) //35
(-5 0 1) //36
(-3 -3 1) //37
(0 -5 1) //38
(3 -3 1) //39

//outer circle1
(15 0 1) //40
(10 10 1) //41
(0 15 1) //42
(-10 10 1) //43
(-15 0 1) //44
(-10 -10 1) //45
(0 -15 1) //46
(10 -10 1) //47

//rectangle1
(75 0 1) //48
(75 10 1) //49
(75 50 1) //50
(10 50 1) //51
(0 50 1) //52
(-10 50 1) //53
(-25 50 1) //54
(-25 10 1) //55
(-25 0 1) //56
(-25 -10 1) //57
(-25 -50 1) //58
(-10 -50 1) //59
(0 -50 1) //60
(10 -50 1) //61
(75 -50 1) //62
(75 -10 1) //63

//z=0
//inner circle2
(105 0 0)//64
(103 3 0)//65
(100 5 0)//66
(97 3 0)//67
(95 0 0)//68
(97 -3 0)//69
(100 -5 0)//70
(103 -3 0)//71

//outer circle2
(115 0 0)//72
(110 10 0)//73
(100 15 0)//74
(90 10 0)//75
(85 0 0)//76
(90 -10 0)//77
(100 -15 0)//78
(110 -10 0)//79

//rectangle2
(200 0 0)//80
(200 10 0)//81
(200 50 0)//82
(110 50 0)//83
(100 50 0)//84
(90 50 0)//85
(75 50 0)//18
(75 10 0)//17
(75 0 0)//16
(75 -10 0)//31
(75 -50 0)//30
(90 -50 0)//91
(100 -50 0)//92
(110 -50 0)//93
(200 -50 0)//94
(200 -10 0)//95

//inner circle2
(105 0 1)//96
(103 3 1)//97
(100 5 1)//98
(97 3 1)//99
(95 0 1)//100
(97 -3 1)//101
(100 -5 1)//102
(103 -3 1)//103

//outer circle2
(115 0 1)//104
(110 10 1)//105
(100 15 1)//106
(90 10 1)//107
(85 0 1)//108
(90 -10 1)//109
(100 -15 1)//110
(110 -10 1)//111

//rectangle2
(200 0 1)//112
(200 10 1)//113
(200 50 1)//114
(110 50 1)//115
(100 50 1)//116
(90 50 1)//117
(75 50 1)//50
(75 10 1)//49
(75 0 1)//48
(75 -10 1)//63
(75 -50 1)//62
(90 -50 1)//123
(100 -50 1)//124
(110 -50 1)//125
(200 -50 1)//126
(200 -10 1)//127
);

blocks
(
hex (1 0 8 9 33 32 40 41) (10 10 1) simpleGrading (1 8 1) //Block 0
hex (2 1 9 10 34 33 41 42) (10 10 1) simpleGrading (1 8 1) //Block 1
hex (3 2 10 11 35 34 42 43) (10 10 1) simpleGrading (1 8 1) //Block 2
hex (4 3 11 12 36 35 43 44) (10 10 1) simpleGrading (1 8 1) //Block 3
hex (5 4 12 13 37 36 44 45) (10 10 1) simpleGrading (1 8 1) //Block 4
hex (6 5 13 14 38 37 45 46) (10 10 1) simpleGrading (1 8 1) //Block 5
hex (7 6 14 15 39 38 46 47) (10 10 1) simpleGrading (1 8 1) //Block 6
hex (0 7 15 8 32 39 47 40) (10 10 1) simpleGrading (1 8 1) //Block 7
hex (8 16 17 9 40 48 49 41) (10 10 1) simpleGrading (2 1 1) //Block 8
hex (9 17 18 19 41 49 50 51) (10 10 1) simpleGrading (2 2 1) //Block 9
hex (10 9 19 20 42 41 51 52) (10 10 1) simpleGrading (1 2 1) //Block 10
hex (11 10 20 21 43 42 52 53) (10 10 1) simpleGrading (1 2 1) //Block 11
hex (23 11 21 22 55 43 53 54) (10 10 1) simpleGrading (1 2 1) //Block 12
hex (24 12 11 23 56 44 43 55) (10 10 1) simpleGrading (1 1 1) //Block 13
hex (25 13 12 24 57 45 44 56) (10 10 1) simpleGrading (1 1 1) //Block 14
hex (26 27 13 25 58 59 45 57) (10 10 1) simpleGrading (1 0.5 1) //Block 15
hex (27 28 14 13 59 60 46 45) (10 10 1) simpleGrading (1 0.5 1) //Block 16
hex (28 29 15 14 60 61 47 46) (10 10 1) simpleGrading (1 0.5 1) //Block 17
hex (29 30 31 15 61 62 63 47) (10 10 1) simpleGrading (2 0.5 1) //Block 18
hex (15 31 16 8 47 63 48 40) (10 10 1) simpleGrading (2 1 1) //Block 19
hex (65 64 72 73 97 96 104 105) (10 10 1) simpleGrading (1 8 1) //Block 20
hex (66 65 73 74 98 97 105 106) (10 10 1) simpleGrading (1 8 1) //Block 21
hex (67 66 74 75 99 98 106 107) (10 10 1) simpleGrading (1 8 1) //Block 22
hex (68 67 75 76 100 99 107 108) (10 10 1) simpleGrading (1 8 1) //Block 23
hex (69 68 76 77 101 100 108 109) (10 10 1) simpleGrading (1 8 1) //Block 24
hex (70 69 77 78 102 101 109 110) (10 10 1) simpleGrading (1 8 1) //Block 25
hex (71 70 78 79 103 102 110 111) (10 10 1) simpleGrading (1 8 1) //Block 26
hex (64 71 79 72 96 103 111 104) (10 10 1) simpleGrading (1 8 1) //Block 27
hex (72 80 81 73 104 112 113 105) (10 10 1) simpleGrading (2 1 1) //Block 28
hex (73 81 82 83 105 113 114 115) (10 10 1) simpleGrading (2 2 1) //Block 29
hex (74 73 83 84 106 105 115 116) (10 10 1) simpleGrading (1 2 1) //Block 30
hex (75 74 84 85 107 106 116 117) (10 10 1) simpleGrading (1 2 1) //Block 31
hex (17 75 85 18 49 107 117 50) (10 10 1) simpleGrading (1 2 1) //Block 32
hex (16 76 75 17 48 108 107 49) (10 10 1) simpleGrading (1 1 1) //Block 33
hex (31 77 76 16 63 109 108 48) (10 10 1) simpleGrading (1 1 1) //Block 34
hex (30 91 77 31 62 123 109 63) (10 10 1) simpleGrading (1 0.5 1) //Block 35
hex (91 92 78 77 123 124 110 109) (10 10 1) simpleGrading (1 0.5 1) //Block 36
hex (92 93 79 78 124 125 111 110) (10 10 1) simpleGrading (1 0.5 1) //Block 37
hex (93 94 95 79 125 126 127 111) (10 10 1) simpleGrading (2 0.5 1) //Block 38
hex (79 95 80 72 111 127 112 104) (10 10 1) simpleGrading (2 1 1) //Block 39


);

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(22 23 55 54)
(23 24 56 55)
(24 25 57 56)
(25 26 58 57)
);
}

outlet
{
type patch;
faces
(
(94 95 127 126)
(95 80 112 127)
(80 81 113 112)
(81 82 114 113)
);
}

topAndBottom
{
type patch;
faces
(
(82 83 115 114)
(83 84 116 115)
(84 85 117 116)
(85 18 50 117)
(18 19 51 50)
(19 20 52 51)
(20 21 53 52)
(21 22 54 53)

(26 27 59 58)
(27 28 60 59)
(28 29 61 60)
(29 30 62 61)
(30 91 123 62)
(91 92 124 123)
(92 93 125 124)
(93 94 126 125)
);
}


cylinder
{
type wall;
faces
(
(1 0 32 33)
(2 1 33 34)
(3 2 34 35)
(4 3 35 36)
(5 4 36 37)
(6 5 37 38)
(7 6 38 39)
(0 7 39 32)

(65 64 96 97)
(66 65 97 98)
(67 66 98 99)
(68 67 99 100)
(69 68 100 101)
(70 69 101 102)
(71 70 102 103)
(64 71 103 96)

);
}

frontAndBack
{
type empty;
faces
(
//circular faces
(0 1 9 8)
(32 40 41 33)

(1 2 10 9)
(33 41 42 34)

(2 3 11 10)
(34 42 43 35)

(3 4 12 11)
(35 43 44 36)

(4 5 13 12)
(36 44 45 37)

(5 6 14 13)
(37 45 46 38)

(6 7 15 14)
(38 46 47 39)

(7 0 8 15)
(39 47 40 32)

(64 65 73 72)
(96 104 105 97)

(65 66 74 73)
(97 105 106 98)

(66 67 75 74)
(98 106 107 99)

(67 68 76 75)
(99 107 108 100)

(68 69 77 76)
(100 108 109 101)

(69 70 78 77)
(101 109 110 102)

(70 71 79 78)
(102 110 111 103)

(71 64 72 79)
(103 111 104 96)

//rectangular faces
(8 9 17 16)
(40 48 49 41)

(9 19 18 17)
(41 49 50 51)

(9 10 20 19)
(41 51 52 42)

(10 11 21 20)
(42 52 53 43)

(11 23 22 21)
(43 53 54 55)

(11 12 24 23)
(43 55 56 44)

(12 13 25 24)
(44 56 57 45)

(13 27 26 25)
(45 57 58 59)

(13 14 28 27)
(45 59 60 46)

(14 15 29 28)
(46 60 61 47)

(15 31 30 29)
(47 61 62 63)

(15 8 16 31)
(47 63 48 40)

(72 73 81 80)
(104 112 113 105)

(73 83 82 81)
(105 113 114 115)

(73 74 84 83)
(105 115 116 106)

(74 75 85 84)
(106 116 117 107)

(75 17 18 85)
(107 117 50 49)

(75 76 16 17)
(107 49 48 108)

(76 77 31 16)
(108 48 63 109)

(77 91 30 31)
(109 63 62 123)

(77 78 92 91)
(109 123 124 110)

(78 79 93 92)
(111 110 124 125)

(79 95 94 93)
(111 125 126 127)

(79 72 80 95)
(111 127 112 104)
);
}
);

mergePatchPairs
(
);

// ************************************************** *********************** //

dynamicMeshDict



Quote:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solver displacementLaplacian;


displacementLaplacianCoeffs
{
diffusivity inverseDistance 1(cylinder);
}





// ************************************************** *********************** //

controlDict



Quote:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application pimpleDyMFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 40;

deltaT 5e-06;

writeControl adjustableRunTime;

writeInterval 0.1;

purgeWrite 0;

writeFormat binary;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep yes;

maxCo 0.5;

functions
{
probes
{
type probes;
functionObjectLibs ("libsampling.so");
enabled true;
outputControl timeStep;
outputInterval 1;
probeLocations
(
( 0.0 0.0 0.0 )
( 10.0 0.0 0.0 )

);

fields
(
p
);
}

forces
{
type forceCoeffs;
functionObjectLibs ( "libforces.so" );
outputControl timeStep;
outputInterval 1;
patches
(
cylinder
);
directForceDensity no;

pName p;
UName U;
rhoName rhoInf;
//log true;
rhoInf 1;//1000; //994.5;
CofR (0 0 0); //( 0.080415 0 0 );
liftDir ( 0 1 0 );
dragDir ( 1 0 0 );
pitchAxis ( 0 0 1 );
magUInf 1.0; //0.1;
lRef 1.0; //0.01787;
Aref 1.0; //0.0003193369;
}

fieldAverage1
{
type fieldAverage;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;
fields
(
U
{
mean on;
prime2Mean on;
base time;
}

p
{
mean on;
prime2Mean on;
base time;
}
);
}

// ************************************************** *********************** //




:)

kerim April 20, 2017 20:11

My case file -simulate two cylinders in tandem in pimpleDyMFoam
 
Greetings from Kyrgyzstan!
Dear maasyraf3,

Please see attached files. I have made very little changes in your controlDict, others are OK. Test case works fine.

Regards,
Kerim

maasyraf3 April 20, 2017 21:32

Hi Kerim! Tq vm for th modification

One more thing, the "rho" is equal to 1, means the rho of cylinder is 1kg/m3 or else?



Sent from my M631Y using CFD Online Forum mobile app

maasyraf3 May 7, 2017 07:12

Quote:

Originally Posted by kerim (Post 645734)
Greetings from Kyrgyzstan!
Dear maasyraf3,

Please see attached files. I have made very little changes in your controlDict, others are OK. Test case works fine.

Regards,
Kerim

Hi Kerim- what do you change in this?
Because i follow yours but still appear the what is rho.
My OF is 4.0
Quote:

Originally Posted by maasyraf3 (Post 643846)
The p




U





pointDisplacement





blockMeshDict







dynamicMeshDict






controlDict









:)


Sent from my M631Y using CFD Online Forum mobile app

kerim May 8, 2017 00:49

1 Attachment(s)
Hi,
Please see the attached file. There is all case file.
Best,
Kerim

Alasir May 11, 2017 10:21

Try changing "rhoName" to "rho". The code for this was updated at some point.

Eefoe June 5, 2017 14:05

Quote:

Originally Posted by Alasir (Post 648474)
Try changing "rhoName" to "rho". The code for this was updated at some point.

That worked for me. I was going through a tutorial for OF 2.x with OF 4.1, and changing rhoName to rho fixed this error for me. One of the many things that needed to change to make the tutorial work properly.


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