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

pimpleDyMFoam - rho not found

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

Like Tree1Likes
  • 1 Post By kerim

Reply
 
LinkBack Thread Tools Display Modes
Old   April 5, 2017, 06:42
Default pimpleDyMFoam - rho not found
  #1
New Member
 
M. Amirul
Join Date: Mar 2017
Location: Malaysia
Posts: 20
Rep Power: 2
maasyraf3 is on a distinguished road
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
Attached Files
File Type: zip viv_test17.zip (182.8 KB, 2 views)
__________________
Best Regards,
maasyraf3
maasyraf3 is offline   Reply With Quote

Old   April 5, 2017, 07:22
Default
  #2
Senior Member
 
Join Date: Mar 2015
Posts: 244
Rep Power: 5
KateEisenhower is on a distinguished road
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
KateEisenhower is offline   Reply With Quote

Old   April 5, 2017, 07:43
Default
  #3
New Member
 
M. Amirul
Join Date: Mar 2017
Location: Malaysia
Posts: 20
Rep Power: 2
maasyraf3 is on a distinguished road
Quote:
Originally Posted by KateEisenhower View Post
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;
}
);
}

// ************************************************** *********************** //
__________________
Best Regards,
maasyraf3

Last edited by maasyraf3; April 5, 2017 at 17:37.
maasyraf3 is offline   Reply With Quote

Old   April 5, 2017, 07:50
Default
  #4
New Member
 
M. Amirul
Join Date: Mar 2017
Location: Malaysia
Posts: 20
Rep Power: 2
maasyraf3 is on a distinguished road
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);
}





// ************************************************** *********************** //
__________________
Best Regards,
maasyraf3

Last edited by maasyraf3; April 5, 2017 at 17:39.
maasyraf3 is offline   Reply With Quote

Old   April 5, 2017, 17:49
Default My case file -simulate two cylinders in tandem in pimpleDyMFoam
  #5
New Member
 
M. Amirul
Join Date: Mar 2017
Location: Malaysia
Posts: 20
Rep Power: 2
maasyraf3 is on a distinguished road
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;
}
);
}

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




__________________
Best Regards,
maasyraf3
maasyraf3 is offline   Reply With Quote

Old   April 20, 2017, 19:11
Default My case file -simulate two cylinders in tandem in pimpleDyMFoam
  #6
New Member
 
abdikerim kurbanaliev
Join Date: Jun 2010
Location: Kyrgyzstan, Osh
Posts: 19
Rep Power: 9
kerim is on a distinguished road
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
Attached Images
File Type: jpg Screenshot from 2017-04-21 05-00-12.jpg (70.6 KB, 15 views)
File Type: jpg Screenshot from 2017-04-21 05-06-35.jpg (116.0 KB, 14 views)
kerim is offline   Reply With Quote

Old   April 20, 2017, 20:32
Default
  #7
New Member
 
M. Amirul
Join Date: Mar 2017
Location: Malaysia
Posts: 20
Rep Power: 2
maasyraf3 is on a distinguished road
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
__________________
Best Regards,
maasyraf3
maasyraf3 is offline   Reply With Quote

Old   May 7, 2017, 06:12
Default
  #8
New Member
 
M. Amirul
Join Date: Mar 2017
Location: Malaysia
Posts: 20
Rep Power: 2
maasyraf3 is on a distinguished road
Quote:
Originally Posted by kerim View Post
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 View Post
The p




U





pointDisplacement





blockMeshDict







dynamicMeshDict






controlDict










Sent from my M631Y using CFD Online Forum mobile app
__________________
Best Regards,
maasyraf3
maasyraf3 is offline   Reply With Quote

Old   May 7, 2017, 23:49
Default
  #9
New Member
 
abdikerim kurbanaliev
Join Date: Jun 2010
Location: Kyrgyzstan, Osh
Posts: 19
Rep Power: 9
kerim is on a distinguished road
Hi,
Please see the attached file. There is all case file.
Best,
Kerim
Attached Files
File Type: gz viv_test17.tar.gz (176.9 KB, 5 views)
maasyraf3 likes this.
kerim is offline   Reply With Quote

Old   May 11, 2017, 09:21
Default
  #10
New Member
 
Anders Utnes
Join Date: May 2017
Location: Norway
Posts: 6
Rep Power: 2
Alasir is on a distinguished road
Try changing "rhoName" to "rho". The code for this was updated at some point.
Alasir is offline   Reply With Quote

Old   June 5, 2017, 13:05
Default
  #11
New Member
 
Join Date: May 2017
Posts: 8
Rep Power: 2
Eefoe is on a distinguished road
Quote:
Originally Posted by Alasir View Post
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.
Eefoe is offline   Reply With Quote

Reply

Tags
cylinder, pimpledymfoam, viv

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Gmsh installation on terminal help spitfire Main CFD Forum 3 August 23, 2016 10:43
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 04:34
Problems in compiling paraview in Suse 10.3 platform chiven OpenFOAM Installation 3 December 1, 2009 08:21
OpenFOAM15 paraFoam bug koen OpenFOAM Bugs 19 June 30, 2009 10:46
Installation OF1.5-dev ttdtud OpenFOAM Installation 46 May 5, 2009 02:32


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