周期性边界条件cyclicAMI的问题?
-
各位大佬,我设置了Mesh的边界为cyclicAMI
inlet { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 1200; startFace 254800; matchTolerance 0.0001; neighbourPatch outlet; transformType translational; separation (-7 0 0); //method faceAreaWeightAMI; } outlet { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 1200; startFace 256000; matchTolerance 0.0001; neighbourPatch inlet; transformType translational; separation (7 0 0); }
那么在0文件中U、k、epsion、p该怎么设置呢,我把inlet/oulet都设置为了cyclicAMI,
inlet { type cyclicAMI; } outlet { type cyclicAMI; }
但是一运行就报错,我想问一下,原来给定inlet的fixvalue,outlet的0梯度,在设置为cyclicAMI type后,这些进口速度什么的去哪给定呢?是在internalField uniform 0.015;这里吗???
报错如下:一改cyclicAMI就会报这个bounding k epsion问题
Time = 13 smoothSolver: Solving for Ux, Initial residual = 2.42791e-08, Final residual = 2.42791e-08, No Iterations 0 smoothSolver: Solving for Uy, Initial residual = 0.000811094, Final residual = 3.38517e-05, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.737597, Final residual = 0.0203503, No Iterations 2 GAMG: Solving for p, Initial residual = 0.890921, Final residual = 0.0244256, No Iterations 2 time step continuity errors : sum local = 6.70555e+80, global = 3.10688e+64, cumulative = 3.10688e+64 smoothSolver: Solving for epsilon, Initial residual = 2.70667e-12, Final residual = 2.70667e-12, No Iterations 0 smoothSolver: Solving for k, Initial residual = 0.0213526, Final residual = 0.000592874, No Iterations 2 bounding k, min: -1.30004e+125 max: 2.00379e+126 average: 1.2466e+122 ExecutionTime = 2.03 s ClockTime = 2 s Time = 14 smoothSolver: Solving for Ux, Initial residual = 3.62239e-09, Final residual = 3.62239e-09, No Iterations 0 smoothSolver: Solving for Uy, Initial residual = 1.07394e-09, Final residual = 1.07394e-09, No Iterations 0 smoothSolver: Solving for Uz, Initial residual = 3.38562e-09, Final residual = 3.38562e-09, No Iterations 0 [stack trace] ============= #1 Foam::sigFpe::sigHandler(int) in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so #2 ? in /lib/x86_64-linux-gnu/libc.so.6 #3 Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so #4 Foam::PCG::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so #5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so #6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so #7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so #8 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so #9 Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so #10 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so #11 Foam::fvMatrix<double>::solve() in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/bin/simpleFoam #12 ? in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/bin/simpleFoam #13 ? in /lib/x86_64-linux-gnu/libc.so.6 #14 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #15 ? in ~/OpenFOAM/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/bin/simpleFoam ============= Floating point exception