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

Please help me...

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 16, 2022, 23:09
Default Please help me...
  #1
New Member
 
DongCheon
Join Date: Aug 2022
Posts: 21
Rep Power: 3
Lee97 is on a distinguished road
I try to solve naca 0015 3D wing. Using pimpleFoam in steady-state.

But I can't run simulation until end.

Time step error goes huge and calculation is stopped.

Please help what is the problems....

I show my settings and grids...

Please help..

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

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

internalField uniform (45.8879 3.2088 0);

boundaryField
{
Airfoil
{
type noSlip;
}
FarField
{
type freestream;
freestreamValue uniform (45.8879 3.2088 0);
value uniform (45.8879 3.2088 0);
}
Inlet
{
type freestream;
freestreamValue uniform (45.8879 3.2088 0);
value uniform (45.8879 3.2088 0);
}
Outlet
{
type zeroGradient;
}
}


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

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

internalField uniform 0;

boundaryField
{
Airfoil
{
type zeroGradient;
}
FarField
{
type freestreamPressure;
value uniform 0;
}
Inlet
{
type freestreamPressure;
value uniform 0;
}
Outlet
{
type zeroGradient;
}
}


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

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

internalField uniform 6.3788e-05;

boundaryField
{
Airfoil
{
type fixedValue;
value uniform 0;
}
FarField
{
type freestream;
freestreamValue uniform 6.3788e-05;
value uniform 6.3788e-05;
}
Inlet
{
type freestream;
freestreamValue uniform 6.3788e-05;
value uniform 6.3788e-05;
}
Outlet
{
type freestream;
freestreamValue uniform 6.3788e-05;
value uniform 6.3788e-05;
}
}


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

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

internalField uniform 1.5947e-05;

boundaryField
{
Airfoil
{
type nutUSpaldingWallFunction;
value uniform 0;
}
FarField
{
type freestream;
freestreamValue uniform 1.5947e-05;
value uniform 1.5947e-05;
}
Inlet
{
type freestream;
freestreamValue uniform 1.5947e-05;
value uniform 1.5947e-05;
}
Outlet
{
type freestream;
freestreamValue uniform 1.5947e-05;
value uniform 1.5947e-05;
}
}


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

simulationType RAS;

RAS
{
RASModel SpalartAllmaras;

turbulence on;

printCoeffs on;
}

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

transportModel Newtonian;

nu [0 2 -1 0 0 0 0] 1.5947e-05;

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

application pimpleFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 5;

deltaT 1e-8;

writeControl adjustableRunTime;

writeInterval 0.1;

purgeWrite 0;

writeFormat ascii;

writePrecision 10;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

adjustTimeStep no;





///////////

functions
{
#includeFunc residuals
}

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

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwind grad(U);
div(phi,nuTilda) bounded Gauss linearUpwind grad(nuTilda);
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

wallDist
{
method meshWave;
}


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

solvers
{
p
{
solver GAMG;
tolerance 1e-06;
relTol 0.01;
smoother GaussSeidel;
}
pFinal
{
$p;
relTol 0;
}

U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0.01;
}
UFinal
{
$U;
relTol 0;
}

nuTilda
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0.01;
}
nuTildaFinal
{
$nuTilda;
relTol 0;
}

}

PIMPLE
{
nNonOrthogonalCorrectors 0;
nCorrectors 1;
nOuterCorrectors 1;
pRefCell 0;
pRefValue 0;
}

relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
nuTilda 0.7;
}
}


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

Quote:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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 : /opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/checkMesh -allGeometry -allTopology
Date : Sep 17 2022
Time : 11:51:33
Host : "anode0"
PID : 1480142
Case : /home/users/DongCheon/3D_Finite_wing/NACA0015_Wing_SA
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 polyMesh for time = 0

Enabling all (cell, face, edge, point) topology checks.

Enabling all geometry checks.

Time = 0

Mesh stats
points: 868532
faces: 4639727
internal faces: 4607299
cells: 1977729
faces per cell: 4.675577898
boundary patches: 4
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 578266
prisms: 112402
wedges: 0
pyramids: 67176
tet wedges: 0
tetrahedra: 1219885
polyhedra: 0

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Topological cell zip-up check OK.
Face-face connectivity OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology Bounding box
Airfoil 23302 21460 ok (closed singly connected) (0 -0.07500506158 0) (1 0.07500506158 6.6)
FarField 6084 6240 ok (non-closed singly connected) (-15 -25 -23) (45 25 27)
Inlet 1521 1600 ok (non-closed singly connected) (-15 -25 -23) (-15 25 27)
Outlet 1521 1600 ok (non-closed singly connected) (45 -25 -23) (45 25 27)

Checking geometry...
Overall domain bounding box (-15 -25 -23) (45 25 27)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (1.491714331e-17 6.618226196e-16 9.501032657e-18) OK.
***High aspect ratio cells found, Max aspect ratio: 68694.60481, number of cells 185822
<<Writing 185822 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 4.360592817e-09. Maximum face area = 1.972386602. Face area magnitudes OK.
Min volume = 2.39262343e-11. Max volume = 0.7156414278. Total volume = 149999.3262. Cell volumes OK.
Mesh non-orthogonality Max: 89.99789344 average: 21.09421162
*Number of severely non-orthogonal (> 70 degrees) faces: 83474.
Non-orthogonality check OK.
<<Writing 83474 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
***Max skewness = 13.75897281, 2317 highly skew faces detected which may impair the quality of the results
<<Writing 2317 skew faces to set skewFaces
Coupled point location match (average 0) OK.
***Error in face tets: 153 faces with low quality or negative volume decomposition tets.
<<Writing 153 faces with low quality or negative volume decomposition tets to set lowQualityTetFaces
*Edges too small, min/max edge length = 7.999979852e-06 2.72060009, number too small: 322015
<<Writing 322015 points on short edges to set shortEdges
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : min = 0.7607182977 average = 0.9998524182
*There are 8 faces with ratio between projected and actual area < 0.8
Minimum ratio (minimum flatness, maximum warpage) = 0.7607182977
<<Writing 8 warped faces to set warpedFaces
Cell determinant (wellposedness) : minimum: 2.881311859e-14 average: 0.8050691697
***Cells with small determinant (< 0.001) found, number of cells: 512017
<<Writing 512017 under-determined cells to set underdeterminedCells
Concave cell check OK.
Face interpolation weight : minimum: 0.005552204927 average: 0.4435016131
***Faces with small interpolation weight (< 0.05) found, number of faces: 12184
<<Writing 12184 faces with low interpolation weights to set lowWeightFaces
Face volume ratio : minimum: 0.005568238408 average: 0.8150794223
***Faces with small volume ratio (< 0.01) found, number of faces: 370
<<Writing 370 faces with low volume ratio cells to set lowVolRatioFaces

Failed 6 mesh checks.

End





Lee97 is offline   Reply With Quote

Reply

« Previous Thread | Next Thread »

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


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

Contact Us - CFD Online - Privacy Statement - Top