CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

overset running blow up

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 26, 2021, 04:35
Exclamation overset running blow up
  #1
New Member
 
Join Date: Nov 2020
Posts: 15
Rep Power: 5
OREOMax is on a distinguished road
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.
__________________
My email adress is
825543420@qq.com
or
gwf-1000rd@outlook.com
OREOMax is offline   Reply With Quote

Old   June 26, 2021, 04:45
Default
  #2
New Member
 
Join Date: Nov 2020
Posts: 15
Rep Power: 5
OREOMax is on a distinguished road
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;
}
__________________
My email adress is
825543420@qq.com
or
gwf-1000rd@outlook.com
OREOMax is offline   Reply With Quote

Old   June 26, 2021, 04:49
Default
  #3
New Member
 
Join Date: Nov 2020
Posts: 15
Rep Power: 5
OREOMax is on a distinguished road
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
__________________
My email adress is
825543420@qq.com
or
gwf-1000rd@outlook.com
OREOMax is offline   Reply With Quote

Old   June 26, 2021, 04:55
Exclamation
  #4
New Member
 
Join Date: Nov 2020
Posts: 15
Rep Power: 5
OREOMax is on a distinguished road
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)
OREOMax is offline   Reply With Quote

Reply

Tags
overset overinterdymfoam


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 07:34.