Hydrofoil in Openfoam

silver.blaze · July 16, 2022, 14:02

Dear Altruist,

I am quite new in OpenFOAM. I am doing turbulence flow pattern of NACA0012 hydrofoil. I have done it in Ansys Fluent. Now I am tring to solve it in OpenFOAM.com to validate the answer. I wanto share the procedures:

1. I have save the ansys fluent mesh in ASCII and used a command fluentMeshtoFoam and covert the mesh file.

2. I have updated the boundary file, p and U file. ( My Re = 2.98X10^6, nu = 0.001002991, rho = 997 kg/m^3, Speed = 3 m/s)

3. Now, simpleFoam command used and execute the case with this error (functionObjects::div contErr writing field: div(phi)). Then foam file generated.

4. Open the foam file in paraview. Extraxt Cp for lower and upper surface. Plot the Cp vs x/c. but didnot get the same result like Ansys, the difference is huge.

Now my questions are:

1. Whats wrong I have done?

2. should I update other files according my Re, V, nu value?

3. How can I solve this?

Kindly help me

My files :

P:

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            freestreamPressure;
        freestreamValue $internalField;
    }

    outlet
    {
        type            freestreamPressure;
        freestreamValue $internalField;
    }

    hydrofoil-wall
    {
        type            zeroGradient;
    }

    frontAndBack
    {
        type            empty;
    }
}

U

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (3 0 0);

boundaryField
{
    inlet
    {
        type            freestreamVelocity;
        freestreamValue $internalField;
    }

    outlet
    {
        type            freestreamVelocity;
        freestreamValue $internalField;
    }

    hydrofoil-wall
    {
        type            fixedValue;
		value			uniform (0 0 0); 
    }

    frontAndBack
    {
        type            empty;
    }
}

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

Transport Properties

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

transportModel  Newtonian;

rho             997;

nu              0.001002991;

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

nut:

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0.14;

boundaryField
{
    inlet
    {
        type            freestream;
        freestreamValue uniform 0.14;
    }

    outlet
    {
        type            freestream;
        freestreamValue uniform 0.14;
    }

    hydrofoil-wall
    {
        type            nutUSpaldingWallFunction;
        value           uniform 0;
    }

    frontAndBack
    {
        type            empty;
    }
}

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

nuTilda

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      nuTilda;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0.14;

boundaryField
{
    inlet
    {
        type            freestream;
        freestreamValue uniform 0.14;
    }

    outlet
    {
        type            freestream;
        freestreamValue uniform 0.14;
    }

    hydrofoil-wall
    {
        type            fixedValue;
        value           uniform 0;
    }

    frontAndBack
    {
        type            empty;
    }
}

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

Fouch · July 17, 2022, 12:04

Hi,

Could you upload the whole case ?

BR

silver.blaze · July 21, 2022, 16:29

Quote:

Originally Posted by Fouch

Hi,

Could you upload the whole case ?

BR

Sorry for my late reply. Can I send it to you in the inbox?

geth03 · July 21, 2022, 17:13

Can you provide more information?
1. Checkmesh result
2. Residuals
3. FvSchemes file
4. Fvsolution file

kandelabr · July 21, 2022, 18:33

looks your 'nu' is the dynamic viscosity but simpleFoam takes kinematic viscosity and doesn't need density...

silver.blaze · July 22, 2022, 00:56

Quote:

Originally Posted by geth03

Can you provide more information?
1. Checkmesh result
2. Residuals
3. FvSchemes file
4. Fvsolution file

Thanks for your response. Here is the result.

1. checkMesh Result:

Code:

Mesh stats
    points:           60700
    internal points:  0
    faces:            120350
    internal faces:   59650
    cells:            30000
    faces per cell:   6
    boundary patches: 4
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     30000
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology
    hydrofoil-wall      200      400      ok (non-closed singly connected)
    inlet               300      602      ok (non-closed singly connected)
    outlet              200      402      ok (non-closed singly connected)
    frontAndBackPlanes  60000    60700    ok (non-closed singly connected)

Checking faceZone topology for multiply connected surfaces...
    No faceZones found.

Checking basic cellZone addressing...
    No cellZones found.

Checking geometry...
    Overall domain bounding box (-4.851962 -5 -0.1790477) (10 5 0.1790477)
    Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
    Mesh has 2 solution (non-empty) directions (1 1 0)
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (-4.998397e-18 -7.594004e-19 8.795399e-19) OK.
    Max cell openness = 2.390399e-16 OK.
    Max aspect ratio = 178.1052 OK.
    Minimum face area = 9.417987e-06. Maximum face area = 0.2519532.  Face area magnitudes OK.
    Min volume = 3.372538e-06. Max volume = 0.04976618.  Total volume = 49.31023.  Cell volumes OK.
    Mesh non-orthogonality Max: 42.75253 average: 6.294525
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.8074126 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

2. Residuals:
I am sorry, I am not sure, where it is.

3. fvSchemes File:

Code:

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

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

ddtSchemes
{
	default	steadyState;
}
gradSchemes
{
	default	Gauss linear;
}
divSchemes
{
	default	            bounded Gauss upwind;
	div(phi,U)	        bounded Gauss upwind;
	div(phi,k)	        bounded Gauss upwind;
	
	//div(phi,epsilon)	bounded Gauss upwind;
    div(phi,omega)	bounded Gauss upwind;
	
	div(phi,R)          bounded Gauss upwind;	
	div(R)	    		Gauss linear;
	
	div(phi,nuTilda)   	bounded Gauss upwind;
	div((nuEff*dev2(T(grad(U)))))	Gauss linear;
	
}
laplacianSchemes
{
	default	Gauss linear corrected;
}
interpolationSchemes
{
	default	linear;
}
snGradSchemes
{
	default	corrected;
}
fluxRequired
{
	default	no;
	p	;
}
wallDist
{
	method	meshWave;
	nRequired	true;
}

4. fvSolution File:

Code:

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

solvers
{
	U
	{
		solver	PBiCGStab;
		tolerance	1E-8;
		relTol	0.1;
		preconditioner	DILU;
		maxIter	100;
	}

	p
	{
		solver	GAMG;
		smoother	GaussSeidel;
		tolerance	1E-8;
		relTol	0.1;
		cacheAgglomeration	true;
		nFinestSweeps	2;
		maxIter	100;
		nPreSweeps	0;
		nPostSweeps	1;
		agglomerator	faceAreaPair;
		nCellsInCoarsestLevel	10;
		mergeLevels	1;
	}

	"(k|epsilon|nuTilda)"

	{
		relTol	0.1;
		preconditioner	DILU;
		tolerance	1E-8;
		maxIter	100;
		solver	PBiCGStab;
	}
	
	Phi
	{
		solver	GAMG;
		smoother	GaussSeidel;
		tolerance	1E-10;
		relTol	01;
		nPreSweeps	0;
		nPostSweeps	1;
		nFinestSweeps	1;
		nCellsInCoarsestLevel	10;
		mergeLevels	2;
		agglomerator	faceAreaPair;
		cacheAgglomeration	true;
		maxIter	100;
	}
}

SIMPLE
{
	nNonOrthogonalCorrectors	1;
	consistent	true;
	pRefCell	0;
	pRefValue	0;
	residualControl
	{
		p	1E-6;
		epsilon	1E-6;
		k	1E-6;
		U	1E-6;
	}
}

relaxationFactors
{
	p	0.75;
	epsilon	0.5;
	k	0.5;
	U	0.75;
}

potentialFlow
{
	nNonOrthogonalCorrectors	10;
	PhiRefCell	0;
	PhiRefValue	0;
}

silver.blaze · July 22, 2022, 00:58

Quote:

Originally Posted by kandelabr

looks your 'nu' is the dynamic viscosity but simpleFoam takes kinematic viscosity and doesn't need density...

Thanks, I have updated it already, but there is no major change

Code:

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

transportModel  Newtonian;

//rho             998.2002;

nu              1e-6;

Kindly help me.

geth03 · July 22, 2022, 03:22

If you use simplefoam you should track the residual behaviour. Look in the tutorials for simplefoam, you should find the files in one. You can also search in google how to do that, there is at least one video.

Gauss upwind is only first order accurate, so if your first run converges with upwind you can use this as starting point for second order schemes like gauss linearUpwindV grad(U) for div(phi,U) for example.

May i ask which turbulence model you use? If k-epsilon you did not upload the boundary files.

SimpleC algorithm: use 0.9 or 0.95 for all equation URF and optionally 0.9 for the fields.

silver.blaze · July 22, 2022, 03:29

Quote:

Originally Posted by geth03

If you use simplefoam you should track the residual behaviour. Look in the tutorials for simplefoam, you should find the files in one. You can also search in google how to do that, there is at least one video.

Gauss upwind is only first order accurate, so if your first run converges with upwind you can use this as starting point for second order schemes like gauss linearUpwindV grad(U) for div(phi,U) for example.

May i ask which turbulence model you use? If k-epsilon you did not upload the boundary files.

SimpleC algorithm: use 0.9 or 0.95 for all equation URF and optionally 0.9 for the fields.

Thanks for your quick response.

Yes, I am using k-epsilon.
Boundary file:

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
    hydrofoil-wall
    {
        type            wall;
        inGroups        1(wall);
        nFaces          200;
        startFace       59650;
    }
    inlet
    {
        type            patch;
        nFaces          300;
        startFace       59850;
    }
    outlet
    {
        type            patch;
        nFaces          200;
        startFace       60150;
    }
    frontAndBackPlanes
    {
        type            empty;
        inGroups        1(empty);
        nFaces          60000;
        startFace       60350;
    }
)

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

silver.blaze · July 22, 2022, 03:30

Quote:

Originally Posted by geth03

If you use simplefoam you should track the residual behaviour. Look in the tutorials for simplefoam, you should find the files in one. You can also search in google how to do that, there is at least one video.

Gauss upwind is only first order accurate, so if your first run converges with upwind you can use this as starting point for second order schemes like gauss linearUpwindV grad(U) for div(phi,U) for example.

May i ask which turbulence model you use? If k-epsilon you did not upload the boundary files.

SimpleC algorithm: use 0.9 or 0.95 for all equation URF and optionally 0.9 for the fields.

I can send you my case file.

geth03 · July 22, 2022, 03:35

Quote:

Originally Posted by silver.blaze

Thanks for your quick response.

Yes, I am using k-epsilon.
Boundary file:

Code:

FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
    hydrofoil-wall
    {
        type            wall;
        inGroups        1(wall);
        nFaces          200;
        startFace       59650;
    }
    inlet
    {
        type            patch;
        nFaces          300;
        startFace       59850;
    }
    outlet
    {
        type            patch;
        nFaces          200;
        startFace       60150;
    }
    frontAndBackPlanes
    {
        type            empty;
        inGroups        1(empty);
        nFaces          60000;
        startFace       60350;
    }
)

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

Sorry, i meant the boundary conditions for k and epsilon.

silver.blaze · July 22, 2022, 03:40

Quote:

Originally Posted by geth03

Sorry, i meant the boundary conditions for k and epsilon.

Sorry,
here it is:

k

Code:

FoamFile
{
	version             2.0;
	format              ascii;
	class               volScalarField;
	location            "0";
	object              k;
}


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

internalField	uniform 0.03375;

boundaryField
{
	hydrofoil-wall
	{
		type	kqRWallFunction;
		value	$internalField;
	}

	inlet
	{
		type	turbulentIntensityKineticEnergyInlet;
		value	$internalField;
		intensity	0.03;
	}
	outlet
	{
		type    zeroGradient;
		//type	inletOutlet;
		//inletValue	uniform 1.0;
	}
}

Epsilon

Code:

FoamFile
{
	version             2.0;
	format              ascii;
	class               volScalarField;
	location            "0";
	object              epsilon;
}


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

internalField	uniform 0.000102
;

boundaryField
{
	hydrofoil-wall
	{
		type	epsilonWallFunction;
		value	$internalField;
	}
	inlet
	{
		type	turbulentMixingLengthDissipationRateInlet;
		value	$internalField;
		mixingLength	0.3;
	}
	outlet
	{
		type	zeroGradient;
	}
}

kandelabr · July 22, 2022, 13:46

Just out of curiosity, you can try using fvSchemes as recommended by wolfDynamics tutorials/courses (change omega to epsilon):

Code:

gradSchemes
{
    default cellLimited leastSquares 1.0;
    grad(U) cellLimited leastSquares 1.0;
}

divSchemes
{
    default         none;

    div(phi,U) bounded Gauss linearUpwindV grad(U);
    div(phi,k) bounded Gauss upwind;
    div(phi,omega) bounded Gauss upwind;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited 0.5;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.5;
}

wallDist
{
    method meshWave;
}

But my guess is that there's something wrong with your ParaView business.
If you need lift and drag, you can use the forceCoeffs functionObject but if you insist on using PV, please share that pipeline as I had similar enterprise there in the past and was stabbed in the back multiple times. I'm lucky I survived at all

silver.blaze · July 22, 2022, 14:15

Quote:

Originally Posted by kandelabr

Just out of curiosity, you can try using fvSchemes as recommended by wolfDynamics tutorials/courses (change omega to epsilon):

Code:

gradSchemes
{
    default cellLimited leastSquares 1.0;
    grad(U) cellLimited leastSquares 1.0;
}

divSchemes
{
    default         none;

    div(phi,U) bounded Gauss linearUpwindV grad(U);
    div(phi,k) bounded Gauss upwind;
    div(phi,omega) bounded Gauss upwind;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited 0.5;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.5;
}

wallDist
{
    method meshWave;
}

But my guess is that there's something wrong with your ParaView business.
If you need lift and drag, you can use the forceCoeffs functionObject but if you insist on using PV, please share that pipeline as I had similar enterprise there in the past and was stabbed in the back multiple times. I'm lucky I survived at all

Thanks for your valuable reply. I can send you the case file.

Fouch · July 22, 2022, 17:01

Hi,

I ran your case with the following configuration file :

fvSchemes :

Code:

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

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

ddtSchemes
{
    default    steadyState;
}
gradSchemes
{
    default    Gauss linear;
}
divSchemes
{
    default                bounded Gauss upwind;
//    div(phi,U)            bounded Gauss upwind;
    div(phi,U)            bounded Gauss linearUpwind grad(U);
    div(phi,k)            bounded Gauss upwind;
    
    //div(phi,epsilon)    bounded Gauss upwind;
    div(phi,omega)    bounded Gauss upwind;
    
    div(phi,R)          bounded Gauss upwind;    
    div(R)                Gauss linear;
    
    div(phi,nuTilda)       bounded Gauss upwind;
    div((nuEff*dev2(T(grad(U)))))    Gauss linear;
    
}
laplacianSchemes
{
    default    Gauss linear corrected;
}
interpolationSchemes
{
    default    linear;
}
snGradSchemes
{
    default    corrected;
}
fluxRequired
{
    default    no;
    p    ;
}
wallDist
{
    method    meshWave;
    nRequired    true;
}

fvSolution :

Code:

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

solvers
{
    U
    {
        solver    PBiCGStab;
        tolerance    1E-8;
        relTol    0.1;
        preconditioner    DILU;
        maxIter    100;
        minIter    2;
    }

    p
    {
        solver    GAMG;
        smoother    GaussSeidel;
        tolerance    1E-8;
        relTol    0.1;
        cacheAgglomeration    true;
        nFinestSweeps    2;
        maxIter    100;
        nPreSweeps    0;
        nPostSweeps    1;
        agglomerator    faceAreaPair;
        nCellsInCoarsestLevel    10;
        mergeLevels    1;
    }

    "(k|epsilon|nuTilda)"
    //"(k|omega|nuTilda)"
    {
        relTol    0.1;
        preconditioner    DILU;
        tolerance    1E-8;
        maxIter    100;
        solver    PBiCGStab;
    }
    
    Phi
    {
        solver    GAMG;
        smoother    GaussSeidel;
        tolerance    1E-10;
        relTol    01;
        nPreSweeps    0;
        nPostSweeps    1;
        nFinestSweeps    1;
        nCellsInCoarsestLevel    10;
        mergeLevels    2;
        agglomerator    faceAreaPair;
        cacheAgglomeration    true;
        maxIter    100;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors    1;
    consistent    true;
    pRefCell    0;
    pRefValue    0;
    residualControl
    {
        p    1E-6;
        epsilon    1E-6;
        k    1E-6;
        U    1E-6;
    }
}

relaxationFactors
{
//    p    0.75;
    epsilon    0.5;
    k    0.5;
    U    0.75;
}

potentialFlow
{
    nNonOrthogonalCorrectors    10;
    PhiRefCell    0;
    PhiRefValue    0;
}

I did not compared the result with yours but I used a 2nd order schemes that should be more accurate.

Test it and tell me if the result is better.

Cheers.

silver.blaze · July 27, 2022, 02:05

Quote:

Originally Posted by Fouch

Hi,

I ran your case with the following configuration file :

fvSchemes :

Code:

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

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

ddtSchemes
{
    default    steadyState;
}
gradSchemes
{
    default    Gauss linear;
}
divSchemes
{
    default                bounded Gauss upwind;
//    div(phi,U)            bounded Gauss upwind;
    div(phi,U)            bounded Gauss linearUpwind grad(U);
    div(phi,k)            bounded Gauss upwind;
    
    //div(phi,epsilon)    bounded Gauss upwind;
    div(phi,omega)    bounded Gauss upwind;
    
    div(phi,R)          bounded Gauss upwind;    
    div(R)                Gauss linear;
    
    div(phi,nuTilda)       bounded Gauss upwind;
    div((nuEff*dev2(T(grad(U)))))    Gauss linear;
    
}
laplacianSchemes
{
    default    Gauss linear corrected;
}
interpolationSchemes
{
    default    linear;
}
snGradSchemes
{
    default    corrected;
}
fluxRequired
{
    default    no;
    p    ;
}
wallDist
{
    method    meshWave;
    nRequired    true;
}

fvSolution :

Code:

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

solvers
{
    U
    {
        solver    PBiCGStab;
        tolerance    1E-8;
        relTol    0.1;
        preconditioner    DILU;
        maxIter    100;
        minIter    2;
    }

    p
    {
        solver    GAMG;
        smoother    GaussSeidel;
        tolerance    1E-8;
        relTol    0.1;
        cacheAgglomeration    true;
        nFinestSweeps    2;
        maxIter    100;
        nPreSweeps    0;
        nPostSweeps    1;
        agglomerator    faceAreaPair;
        nCellsInCoarsestLevel    10;
        mergeLevels    1;
    }

    "(k|epsilon|nuTilda)"
    //"(k|omega|nuTilda)"
    {
        relTol    0.1;
        preconditioner    DILU;
        tolerance    1E-8;
        maxIter    100;
        solver    PBiCGStab;
    }
    
    Phi
    {
        solver    GAMG;
        smoother    GaussSeidel;
        tolerance    1E-10;
        relTol    01;
        nPreSweeps    0;
        nPostSweeps    1;
        nFinestSweeps    1;
        nCellsInCoarsestLevel    10;
        mergeLevels    2;
        agglomerator    faceAreaPair;
        cacheAgglomeration    true;
        maxIter    100;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors    1;
    consistent    true;
    pRefCell    0;
    pRefValue    0;
    residualControl
    {
        p    1E-6;
        epsilon    1E-6;
        k    1E-6;
        U    1E-6;
    }
}

relaxationFactors
{
//    p    0.75;
    epsilon    0.5;
    k    0.5;
    U    0.75;
}

potentialFlow
{
    nNonOrthogonalCorrectors    10;
    PhiRefCell    0;
    PhiRefValue    0;
}

I did not compared the result with yours but I used a 2nd order schemes that should be more accurate.

Test it and tell me if the result is better.

Cheers.

Sorry!

Could you please send me the case file?

Fouch · July 27, 2022, 15:06

Could you explain what's wrong ?

Download link for the case : https://nextcloud.cm-en-transition.fr/s/kjnYJkF2RrHWbKR

silver.blaze · July 27, 2022, 17:16

Quote:

Originally Posted by Fouch

Could you explain what's wrong ?

Download link for the case : https://nextcloud.cm-en-transition.fr/s/kjnYJkF2RrHWbKR

Cp vs X/c graph's trend or pattern is okay, but the main problem is value. the value is not correct of Cp.
Kindly help me

July 17, 2022, 12:04		#2
Fouch Member Join Date: Feb 2020 Posts: 79 Rep Power: 6	Hi, Could you upload the whole case ? BR silver.blaze likes this.

July 21, 2022, 18:33	Hydrofoil in Openfoam	#5
kandelabr Senior Member Nejc Join Date: Feb 2017 Location: Slovenia Posts: 196 Rep Power: 9	looks your 'nu' is the dynamic viscosity but simpleFoam takes kinematic viscosity and doesn't need density... __________________ www.damogranlabs.com

July 22, 2022, 13:46		#13
kandelabr Senior Member Nejc Join Date: Feb 2017 Location: Slovenia Posts: 196 Rep Power: 9	Just out of curiosity, you can try using fvSchemes as recommended by wolfDynamics tutorials/courses (change omega to epsilon): Code: gradSchemes { default cellLimited leastSquares 1.0; grad(U) cellLimited leastSquares 1.0; } divSchemes { default none; div(phi,U) bounded Gauss linearUpwindV grad(U); div(phi,k) bounded Gauss upwind; div(phi,omega) bounded Gauss upwind; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited 0.5; } interpolationSchemes { default linear; } snGradSchemes { default limited 0.5; } wallDist { method meshWave; } But my guess is that there's something wrong with your ParaView business. If you need lift and drag, you can use the forceCoeffs functionObject but if you insist on using PV, please share that pipeline as I had similar enterprise there in the past and was stabbed in the back multiple times. I'm lucky I survived at all silver.blaze likes this. __________________ www.damogranlabs.com

July 27, 2022, 15:06		#17
Fouch Member Join Date: Feb 2020 Posts: 79 Rep Power: 6	Could you explain what's wrong ? Download link for the case : https://nextcloud.cm-en-transition.fr/s/kjnYJkF2RrHWbKR silver.blaze likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
OpenFOAM 5.0 Released	CFDFoundation	OpenFOAM Announcements from OpenFOAM Foundation	11	June 5, 2018 23:48
OpenFOAM Training Jan-Jul 2017, Virtual, London, Houston, Berlin	CFDFoundation	OpenFOAM Announcements from Other Sources	0	January 4, 2017 06:15
OpenFOAM Training, London, Chicago, Munich, Sep-Oct 2015	cfd.direct	OpenFOAM Announcements from Other Sources	2	August 31, 2015 13:36
OpenFOAM Foundation releases OpenFOAM 2.2.2	opencfd	OpenFOAM Announcements from ESI-OpenCFD	0	October 14, 2013 07:18
Adventure of fisrst openfoam installation on Ubuntu 710	jussi	OpenFOAM Installation	0	April 24, 2008 14:25

July 21, 2022, 17:13		#4
geth03 Senior Member Join Date: Dec 2019 Location: Cologne, Germany Posts: 363 Rep Power: 8	Can you provide more information? 1. Checkmesh result 2. Residuals 3. FvSchemes file 4. Fvsolution file

July 22, 2022, 03:22		#8
geth03 Senior Member Join Date: Dec 2019 Location: Cologne, Germany Posts: 363 Rep Power: 8	If you use simplefoam you should track the residual behaviour. Look in the tutorials for simplefoam, you should find the files in one. You can also search in google how to do that, there is at least one video. Gauss upwind is only first order accurate, so if your first run converges with upwind you can use this as starting point for second order schemes like gauss linearUpwindV grad(U) for div(phi,U) for example. May i ask which turbulence model you use? If k-epsilon you did not upload the boundary files. SimpleC algorithm: use 0.9 or 0.95 for all equation URF and optionally 0.9 for the fields.