大涡模拟哪种差值格式收敛性最好啊？ Ｇamma ? vanLee?
M
mark
@mark
帖子
-
差值格式 -
calculated 边界条件怎么用？有谁知道啊？这样设置好像不行啊@cfd-china
-
calculated 边界条件怎么用？有谁知道啊？PIMPLE: iteration 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensioned<double> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(double const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#7 Foam::phaseChangeTwoPhaseMixture::vDotAlphal() const at ??:?
#8 ? at ??:?
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? at ??:? -
calculated 边界条件怎么用？有谁知道啊？FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //dimensions [1 -1 -2 0 0 0 0];
internalField uniform 3e+06;
boundaryField
{
INLET
{
type fixedValue;
value uniform 10e+06;
}WALL { type zeroGradient; } OUTLET { type waveTransmissive; value uniform 3e+06; field p; gamma 1.17; phi phi; rho rho; psi psi; lInf 0.05; fieldInf 3e+06; }}
-
calculated 边界条件怎么用？有谁知道啊？FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //dimensions [1 -1 -2 0 0 0 0];
internalField uniform 3e+06;
boundaryField
{
INLET
{
type calculated;
}WALL { type zeroGradient; } OUTLET { type calculated; }}
// ************************************************************************* //
-
calculated 边界条件怎么用？有谁知道啊？This boundary condition
-
射流算例速度发散了，有谁知道怎么解决吗速度和压力设置不对？
-
射流算例速度发散了，有谁知道怎么解决吗/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object alpha.fuel; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 0 0 0 0]; internalField uniform 0; boundaryField { INLET { type fixedValue; value uniform 1; } WALL { type zeroGradient; } OUTLET { type zeroGradient; } } // ************************************************************************* // -
射流算例速度发散了，有谁知道怎么解决吗/*------------none--------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 3e+06; boundaryField { INLET { type fixedValue; value uniform 100e+06; } WALL { type zeroGradient; } OUTLET { type zeroGradient; } } // ************************************************************************* // -
射流算例速度发散了，有谁知道怎么解决吗/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / 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 3e+06; boundaryField { INLET { type fixedValue; value uniform 100e+6; } WALL { type zeroGradient; } OUTLET { type zeroGradient; } } // ************************************************************************* // -
射流算例速度发散了，有谁知道怎么解决吗速度
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.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 (0 0 0); boundaryField { INLET { type zeroGradient; } WALL { type zeroGradient; } OUTLET { type zeroGradient; } } // ************************************************************************* // -
射流算例速度发散了，有谁知道怎么解决吗
大口进，小口出 -
射流算例速度发散了，有谁知道怎么解决吗/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { alpha.fuel { nAlphaCorr 1; nAlphaSubCycles 1; cAlpha 1; } ".*(rho|rhoFinal)" { solver diagonal; } pcorr { solver GAMG; preconditioner { preconditioner GAMG; tolerance 1e-05; relTol 0; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } tolerance 1e-05; relTol 0; maxIter 100; } p_rgh { solver GAMG; tolerance 1e-07; relTol 0.01; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } p_rghFinal { solver GAMG; preconditioner { preconditioner GAMG; tolerance 1e-07; relTol 0; nVcycles 2; smoother GaussSeidel; nPreSweeps 2; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } tolerance 1e-07; relTol 0; maxIter 20; } "(U|T|k|B|nuTilda).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-08; relTol 0; } "e.*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-8; relTol 0; minIter 1; } } PIMPLE { momentumPredictor yes; transonic yes; nOuterCorrectors 5; nCorrectors 1; nNonOrthogonalCorrectors 2; residualControl { U { tolerance 1e-07; relTol 0; } p_rgh { tolerance 1e-07; relTol 0; } T { tolerance 1e-07; relTol 0; } k { tolerance 1e-07; relTol 0; } } relaxtionFactors { fields { p_rgh 0.4; p_rghFinal 1; } equations { "(U|T|k|B|nuTilda)" 0.4; "(U|T|k|B|nuTilda)Final" 1; } } // ************************************************************************* // -
射流算例速度发散了，有谁知道怎么解决吗*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { default Gauss upwind; div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss linear; div(rhoPhi,U) Gauss upwind; div(phi,thermo:rho.fuel) Gauss upwind; div(phi,thermo:rho.air) Gauss upwind; div(rhoPhi,T) Gauss upwind; div(rhoPhi,K) Gauss upwind; div(phi,p) Gauss upwind; div(phi,k) Gauss upwind; div((muEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear uncorrected; } interpolationSchemes { default linear; } snGradSchemes { default uncorrected; } fluxRequired { default no; p_rgh; pcorr; } // ************************************************************************* // -
射流算例速度发散了，有谁知道怎么解决吗FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application myFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 1.8e-04; deltaT 1e-09; writeControl adjustableRunTime; writeInterval 1e-07; purgeWrite 0; writeFormat binary; writePrecision 8; writeCompression uncompressed; timeFormat general; timePrecision 6; runTimeModifiable yes; adjustTimeStep yes; maxCo 0.5; maxAlphaCo 0.5; maxDeltaT 1e-09; // ************************************************************************* // -
源项与压力耦合的时候，sp su怎样设置？每个网格内的压力不一样，怎样使源项在网格中不同，而不是用整体的源项呢？
-
MULES方法求解alpha方程出错，有谁知道怎样修改吗？2 sp su 怎样赋值给每个网格，而不是整体的源项？
-
MULES方法求解alpha方程出错，有谁知道怎样修改吗？1 , sp su 该怎样修改？
-
MULES方法求解alpha方程出错，有谁知道怎样修改吗？
-
MULES方法求解alpha方程出错，有谁知道怎样修改吗？