|
[Sponsors] |
Running rhoSimpleFoam for Pressure Driven Flow |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 21, 2018, 14:53 |
Running rhoSimpleFoam for Pressure Driven Flow
|
#1 |
Member
Yijiu Jiang
Join Date: Jan 2013
Location: Michigan, US
Posts: 49
Rep Power: 13 |
Hello everyone,
I am running rhoSimpleFoam for a pressure driven flow model. The model is a tank with fixed wall temperature and cooler air with room temperature is blown into the tank. The inlet and outlet pressure are 9 and 6 bar. I use kEpsilon model and wallfunction for this model. The mesh was generated by snappyHexMesh. But I got the following error in 11th iteration. Does anyone have idea on what this message refer to? The BCs and fvScheme and fvSolution are shown below. Thank you in advance for your help! Time = 11 GAMG: Solving for Ux, Initial residual = 0.04648656, Final residual = 7.246233e-05, No Iterations 1 GAMG: Solving for Uy, Initial residual = 0.04807173, Final residual = 5.169749e-05, No Iterations 1 GAMG: Solving for Uz, Initial residual = 0.02602888, Final residual = 1.604139e-05, No Iterations 1 GAMG: Solving for e, Initial residual = 0.02056888, Final residual = 0.0001085386, No Iterations 1 #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:? #4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:? #5 ? in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam" #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 ? in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam" ~ Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object alphat; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { defaultFaces { type empty; } tankfull_inlet { type calculated; value uniform 0; } tankfull_outlet { type calculated; value uniform 0; } tankfull_needle { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } tankfull_tank { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object epsilon; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -3 0 0 0 0]; internalField uniform 65000; boundaryField { defaultFaces { type empty; } tankfull_inlet { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.005; value uniform 0.0028; } tankfull_outlet { type zeroGradient; } tankfull_needle { type epsilonWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 65000; } tankfull_tank { type epsilonWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 65000; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object k; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 85; boundaryField { defaultFaces { type empty; } tankfull_inlet { type turbulentIntensityKineticEnergyInlet; intensity 0.05; value uniform 0.0028; } tankfull_outlet { type zeroGradient; } tankfull_needle { type kqRWallFunction; value uniform 85; } tankfull_tank { type kqRWallFunction; value uniform 2; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object nut; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { defaultFaces { type empty; } tankfull_inlet { type calculated; value uniform 0; } tankfull_outlet { type calculated; value uniform 0; } tankfull_needle { type nutkWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 0; } tankfull_tank { type nutkWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 0; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 800000; boundaryField { defaultFaces { type empty; } tankfull_inlet { type fixedValue; value uniform 900000; } tankfull_outlet { type fixedValue; value uniform 600000; } tankfull_needle { type zeroGradient; } tankfull_tank { type zeroGradient; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 383.0; boundaryField { defaultFaces { type empty; } tankfull_inlet { type fixedValue; value uniform 293.0; } tankfull_outlet { type inletOutlet; value uniform 293.0; inletValue uniform 293.0; } tankfull_needle { type zeroGradient; } tankfull_tank { type fixedValue; value uniform 473.0; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // //Uinlet (10 0 0); dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { defaultFaces { type empty; } tankfull_inlet { type pressureInletOutletVelocity; value uniform (0 0 0); } tankfull_outlet { type pressureInletOutletVelocity; value uniform (0 0 0); } tankfull_needle { type fixedValue; value uniform (0 0 0); } tankfull_tank { type fixedValue; value uniform (0 0 0); } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / 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 upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; div(phi,e) bounded Gauss upwind; div(phi,epsilon) bounded Gauss upwind; div(phi,k) bounded Gauss upwind; div(phid,p) Gauss upwind; div(phi,Ekp) bounded Gauss upwind; div((phi|interpolate(rho)),p) Gauss upwind; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / 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-08; relTol 0.1; smoother GaussSeidel; nCellsInCoarsestLevel 20; } "(U|e|k|epsilon)" { solver GAMG; tolerance 1e-08; relTol 0.1; smoother GaussSeidel; nCellsInCoarsestLevel 20; } } SIMPLE { nNonOrthogonalCorrectors 0; pMinFactor 0.1; pMaxFactor 2; transonic yes; consistent yes; residualControl { p 1e-3; U 1e-3; e 1e-3; // possibly check turbulence fields "(k|epsilon|omega)" 1e-3; } } relaxationFactors { fields { p 0.7; } equations { p 0.7; U 0.3; e 0.6; k 0.6; epsilon 0.6; } } // ************************************************************************* // Last edited by y_jiang; August 24, 2018 at 11:03. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Issues on the simulation of high-speed compressible flow within turbomachinery | dowlee | OpenFOAM Running, Solving & CFD | 11 | August 6, 2021 06:40 |
problem about pressure driven flow | yhaomin2007 | OpenFOAM Running, Solving & CFD | 8 | June 4, 2021 07:51 |
sonicFoam - pressure driven pipe: flow continuity violation and waveTransmissive BC | Endel | OpenFOAM Running, Solving & CFD | 3 | September 11, 2014 16:29 |
Running a feedback control simulation of the flow around a cylinder | gasagna | OpenFOAM Running, Solving & CFD | 0 | July 16, 2014 10:03 |
LTSInterFoam and pressure driven flow | Andrea_85 | OpenFOAM | 3 | August 21, 2012 08:44 |