|
[Sponsors] |
June 26, 2021, 04:35 |
overset running blow up
|
#1 |
New Member
Join Date: Nov 2020
Posts: 15
Rep Power: 6 |
Hi,guys
I am struggling in a overset problem.I've check my mesh several times and I think its ok.celltypes,zoneID,overset almost anywhere are exactly the same just like the tutorials. But my running is still blow up,first I try overInterdymfoam,copied fvschemes and solutions for my own case,also the 0/initial alpha.water,p_rgh,U and so on. |
|
June 26, 2021, 04:45 |
|
#2 |
New Member
Join Date: Nov 2020
Posts: 15
Rep Power: 6 |
object alpha.water;
} // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 0 0]; internalField uniform 1; boundaryField { #includeEtc "caseDicts/setConstraintTypes" up { type overset; } inlet { type fixedValue; value uniform 1; } outlet { type inletOutlet; inletValue uniform 1; } fixedwalls { type zeroGradient; } down { patchType wall; type zeroGradient; } back0 { type wedge; } back { type wedge; } front0 { type wedge; } front { type wedge; } axis { type empty; } object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0]; internalField uniform 0; boundaryField { #includeEtc "caseDicts/setConstraintTypes" oversetPatch { type overset; } up { //patchType overset; type overset; value $internalField; } inlet { type zeroGradient; } outlet { type fixedValue; value $internalField; } fixedwalls { type zeroGradient; } down { type zeroGradient; } back0 { type wedge; } back { type wedge; } front0 { type wedge; } front { type wedge; } axis { type empty; } object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0]; internalField uniform (1 0 0); boundaryField { #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedValue; value uniform (1 0 0); } outlet { type pressureInletOutletVelocity; value uniform (1 0 0); //meaningless } fixedwalls { type zeroGradient; } up//@@ { type overset; } down { patchType wall; type movingWallVelocity; value uniform (0 0 0); } back0 { type wedge; } back { type wedge; } front0 { type wedge; } front { type wedge; } axis { type empty; } object zoneID; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 0 0 0 0]; internalField uniform 0; boundaryField { #includeEtc "caseDicts/setConstraintTypes" up { type overset; } front0 { type wedge; } back0 { type wedge; } inlet { type zeroGradient; } outlet { type zeroGradient; } fixedwalls { type zeroGradient; } front { type wedge; } back { type wedge; } down { patchType wall; type zeroGradient; } axis { type empty; } |
|
June 26, 2021, 04:49 |
|
#3 |
New Member
Join Date: Nov 2020
Posts: 15
Rep Power: 6 |
And the exception message goes like
Execution time for mesh.update() = 4.66 s smoothSolver: Solving for alpha.water, Initial residual = 0.0211346821908, Final residual = 0.0135779401161, No Iterations 1000 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.0183158744867, Final residual = 0.0132939251548, No Iterations 1000 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.0145174915134, Final residual = 0.000141800548143, No Iterations 145 time step continuity errors : sum local = 2.36996719524e-07, global = -1.36752080751e-07, cumulative = 4.2624379384e-05 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.00168907824048, Final residual = 3.02691821844e-10, No Iterations 349 time step continuity errors : sum local = 1.15558499179e-11, global = -1.034834113e-11, cumulative = 4.26243690356e-05 PIMPLE: iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 0.0357528511759, Final residual = 0.0135027766473, No Iterations 1000 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.0183358150928, Final residual = 0.0134021256367, No Iterations 1000 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000792210553907, Final residual = 7.793206267e-06, No Iterations 202 time step continuity errors : sum local = 2.97630519529e-06, global = -1.02830716398e-06, cumulative = 4.15960618717e-05 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000733350123913, Final residual = 4.04437451455e-10, No Iterations 342 time step continuity errors : sum local = 9.02037418296e-10, global = -8.18616151815e-10, cumulative = 4.15952432555e-05 ExecutionTime = 1667.02 s ClockTime = 1668 s Courant Number mean: 0.0860046423749 max: 25.8662257931 Interface Courant Number mean: 0 max: 0 deltaT = 1.73812108246e-155 Time = 1.20785726e-07 PIMPLE: iteration 1 cellVolumeWeight : detected 2 mesh regions zone:0 nCells:87400 zone:1 nCells:2425 Marking patch-cells on zone 0 Marking patch-cells on zone 1 After patch analysis : nCells : 89825 other : 87946 patch : 1382 overset: 497 Creating mesh-to-mesh addressing for region0 and region0 regions using cellVolumeWeight Overlap volume: 5.00255075476e-10 cellVolumeWeight : detected 3 mesh regions after overset Front : fraction:1 size:1602 Overset analysis : nCells : 89825 calculated : 59166 interpolated : 1443 hole : 29216 Execution time for mesh.update() = 4.6 s smoothSolver: Solving for alpha.water, Initial residual = 0.0204861849052, Final residual = 0.0127746202415, No Iterations 1000 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for alpha.water, Initial residual = 0.0178901469289, Final residual = 0.012371690307, No Iterations 1000 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0 #3 Foam:ILUPreconditioner:recondition(Foam::Field <double>&, Foam::Field<double> const&, unsigned char) const at ??:? #4 Foam::PBiCGStab::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #5 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:? #7 Foam::fvMatrix<double>::solveSegregatedOrCoupled(F oam::dictionary const&) at ??:? #8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:? #9 Foam::SolverPerformance<double> Foam::dynamicOversetFvMesh::solve<double>(Foam::fv Matrix<double>&, Foam::dictionary const&) const at ??:? #10 Foam::dynamicOversetFvMesh::solve(Foam::fvMatrix<d ouble>&, Foam::dictionary const&) const at ??:? #11 ? in ~/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/bin/overInterDyMFoam #12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #13 ? in ~/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/bin/overInterDyMFoam |
|
June 26, 2021, 04:55 |
|
#4 |
New Member
Join Date: Nov 2020
Posts: 15
Rep Power: 6 |
Then I try overPimpleDyMFoam and remove alpha.water,all settings follow the wolf dynamics overset tutorials --flapping aerofoil.
log↓ Create mesh for time = 0 Selecting dynamicFvMesh dynamicOversetFvMesh Selecting motion solver: multiSolidBodyMotionSolver Applying solid body motion to entire mesh Selecting solid-body motion function linearMotion Applying solid body motion linearMotion to 5820 points of cellZone move PIMPLE: no residual control data found. Calculations will employ 2 corrector loops Reading field p Reading field U Reading/calculating face flux field phi Creating cellMask field to block out hole cells Creating interpolatedCells field Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes Reading/calculating face velocity Uf No MRF models present No finite volume options present Courant Number mean: 0.000232462285904 max: 0.00115752728397 Starting time loop Courant Number mean: 0.000232462285904 max: 0.00115752728397 deltaT = 1.2e-08 Time = 1.2e-08 inverseDistance : detected 2 mesh regions zone:0 nCells:87400 voxels:(299 299 1) bb:(-2.57139161993e-08 -2.57139161993e-08 -0.000261742038109) (0.0250000257139 0.00599431504341 0.000261742038109) zone:1 nCells:2425 voxels:(299 299 1) bb:(0.00459218180196 9.98880241227e-05 -0.00019630404118) (0.020800136798 0.00449573379516 0.00019630404118) Overset analysis : nCells : 89825 calculated : 2284 interpolated : 2077 (interpolated from local:2076 mixed local/remote:0 remote:0) hole : 85464 PIMPLE: iteration 1 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0 #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) at ??:? #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) at ??:? #5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) at ??:? #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) at ??:? #7 Foam::PBiCGStab::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #8 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #9 Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:? #10 Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:? #11 Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<do uble> >&, Foam::dictionary const&) const at ??:? #12 Foam::SolverPerformance<Foam::Vector<double> > Foam::dynamicOversetFvMesh::solve<Foam::Vector<dou ble> >(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at ??:? #13 Foam::dynamicOversetFvMesh::solve(Foam::fvMatrix<F oam::Vector<double> >&, Foam::dictionary const&) const at ??:? #14 ? in ~/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/bin/overPimpleDyMFoam #15 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #16 ? in ~/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/bin/overPimpleDyMFoam Floating point exception (core dumped) |
|
Tags |
overset overinterdymfoam |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Overset mesh interpolation extremely time consuming | JohnMartinGodo | OpenFOAM Running, Solving & CFD | 12 | February 24, 2023 04:01 |
How to set up an overset mesh ? (Moving square inside a rectangular box) | jordan.c | STAR-CCM+ | 5 | September 6, 2019 09:53 |
OpenFOAM v1706 – aeroFoil tutorial | Fieschi | OpenFOAM Post-Processing | 12 | February 23, 2018 10:18 |
Body Coupling in overset region | sjtuwjh | STAR-CCM+ | 1 | July 28, 2014 05:34 |
[Star-CCM+ v9] Dynamic Overset Mesh basic problem | johncoffey | STAR-CCM+ | 5 | June 17, 2014 15:14 |