View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0000464 | OpenFOAM | Bug | public | 2012-03-16 20:12 | 2012-08-13 10:43 |
| Reporter | Assigned To | ||||
| Priority | normal | Severity | feature | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| Platform | Linux | OS | OpenSUSE | OS Version | 11.4 |
| Summary | 0000464: implementation of mapping for clouds (and possibly parcels/particles) | ||||
| Description | I've been trying to work with dynamic meshes and lagrangrian sprays and I keep getting the following error that crashes my simulation: --> FOAM FATAL ERROR: Not implemented From function cloud::autoMap(const mapPolyMesh&) in file fields/cloud/cloud.C at line 65. FOAM aborting It would be really helpful to include a mapper for clouds as this would open up a lot of functionality. | ||||
| Steps To Reproduce | 1. Compile the attached solver sprayFoamDynamicMesh (included in the tarball) 2. Run sprayFoamDynamicMesh on the attached case aachenBombDynamicMesh (based on the regular tutorial, also included in the tarball) | ||||
| Additional Information | The attached solver and case are based on the following discussion on the forum: http://www.cfd-online.com/Forums/openfoam-solving/94493-adaptive-local-mesh-refinement.html See also the following report: http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2010/anneKoesters/anneKoestersReport.pdf I'm also including the stack trace from the debug version when the refinement begins: #0 Foam::error::printStack(Foam::Ostream&) at ~/OpenFOAM/OpenFOAM-2.1.x/src/OSspecific/POSIX/printStack.C:201 #1 Foam::error::abort() at ~/OpenFOAM/OpenFOAM-2.1.x/src/OpenFOAM/lnInclude/error.C:249 #2 Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) at ~/OpenFOAM/OpenFOAM-2.1.x/src/OpenFOAM/lnInclude/errorManip.H:85 #3 Foam::cloud::autoMap(Foam::mapPolyMesh const&) at ~/OpenFOAM/OpenFOAM-2.1.x/src/OpenFOAM/fields/cloud/cloud.C:66 #4 Foam::mapClouds(Foam::objectRegistry const&, Foam::mapPolyMesh const&) at ~/OpenFOAM/OpenFOAM-2.1.x/src/OpenFOAM/lnInclude/mapClouds.H:52 #5 Foam::fvMesh::mapFields(Foam::mapPolyMesh const&) at ~/OpenFOAM/OpenFOAM-2.1.x/src/finiteVolume/fvMesh/fvMesh.C:489 #6 Foam::fvMesh::updateMesh(Foam::mapPolyMesh const&) at ~/OpenFOAM/OpenFOAM-2.1.x/src/finiteVolume/fvMesh/fvMesh.C:674 #7 Foam::dynamicRefineFvMesh::refine(Foam::List<int> const&) at ~/OpenFOAM/OpenFOAM-2.1.x/src/dynamicFvMesh/dynamicRefineFvMesh/dynamicRefineFvMesh.C:264 #8 Foam::dynamicRefineFvMesh::update() at ~/OpenFOAM/OpenFOAM-2.1.x/src/dynamicFvMesh/dynamicRefineFvMesh/dynamicRefineFvMesh.C:1186 #9 at ~/OpenFOAM/OpenFOAM-2.1.x/applications/solvers/lagrangian/sprayFoamDynamicMesh/sprayFoam.C:80 #10 __libc_start_main in "/lib/libc.so.6" #11 at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/i386/elf/start.S:122 Aborted | ||||
| Tags | No tags attached. | ||||
|
2012-03-16 20:12
|
|
|
|
It would be nice to get feedback on this bug, but I realize I have to give before I can get ;) Looking at the lagrangian libraries, there is an autoMap function in $FOAM_SRC/lagrangian/basic/Cloud that should be called when autoMap is called from here $FOAM_SRC/OpenFOAM/fields/cloud. Or at least it was in the older versions of OpenFOAM when Foam::cloud::autoMap was a pure virtual function. I tried changing Foam::cloud::autoMap to a pure virtual but that didn't work. So there is likely more to the linking/templates than I can see. |
|
|
The cloud::autoMap is the fall-back version. Your cloud class should define its own virtual autoMap function which then can call the Cloud::autoMap(trackData&, const mapPolyMesh&) to do the particles and any additional stuff. The key is that the Cloud::autoMap needs trackData argument so you have to supply that somehow. |
|
|
This is already implemented in our internal development line - however, there are additional changes still required that relate to the improved tracking algorithm introduced in version 2.0. |
|
|
Great news on the impending release. Are there major issues with just making the substitution (similar to your development line)? I'd like to test the change myself, as this is some functionality I'd like to try out, so that when the real version hits I can be ready to use it right away. Or if you need testers I would also be happy to assist :) Because of the amount of templatization involved, an example would be instructive. Keep up the great work! |
|
|
I've pushed a couple of commits 0ffb024 and 9c8ce25 that add the autoMap functionality to the clouds and clear the tets used by the tracking algorithm when the mesh is updated - hopefully these should help |
|
|
Awesome! I'll pull them now and give them a test. Should I report anything under here or file a new issue? |
|
|
I'll leave the thread open for now for anything directly related to the recent commits |
|
|
The mapping looks like its working well, there are some issues with the tets after a refine operation (using dynamicRefineFvMesh). liquidSpray0 is a basicSprayCloud and a gasSpray0 is a coalCloud. The error I get is: Courant Number mean: 0.000213884 max: 0.0845929 deltaT = 1.36242e-06 Time = 9.86376e-05 s Selected 24 cells for refinement out of 168100. Refined from 168100 to 168268 cells. Cloud<ParticleType>::autoMap(TrackData&, const mapPolyMesh&) for lagrangian cloud gasSpray0 Cloud<ParticleType>::autoMap(TrackData&, const mapPolyMesh&) for lagrangian cloud liquidSpray0 --> FOAM FATAL ERROR: cell, tetFace and tetPt search failure at position (0.000108393 0.0550254 0.000119704) for requested cell 93295 From function void Foam::particle::initCellFacePt() in file /home/gandalf/OpenFOAM/OpenFOAM-2.1.x/src/lagrangian/basic/lnInclude/particleI.H at line 758. My guess is the problem is with the liquid spray, as at this point there are only parcels from the liquidSpray0. I have a debug version running that I will upload the log for if you want it. |
|
|
Thanks - the log file would be useful |
|
2012-08-03 17:15
|
foam.log (58,048 bytes)
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.x-69fb1b62fd1e
Exec : sprayDyMFoam
Date : Aug 03 2012
Time : 08:46:43
Host : "tolkien01"
PID : 10316
Case : /home/gandalf/OpenFOAM/OpenFOAM-2.1.x/tutorials/lagrangian/sprayFoam/aachenBombRefine
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Selecting dynamicFvMesh dynamicRefineFvMesh
Reading g
Creating combustion model
Selecting combustion model PaSR<psiChemistryCombustionModel>
Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>>
Selecting chemistryReader chemkinReader
ODEChemistryModel: Number of species = 5 and reactions = 1
Selecting ODE solver SIBS
Creating component thermo properties:
multi-component carrier - 5 species
liquids - 1 components
solids - 0 components
Reading field U
Reading/calculating face flux field phi
Creating turbulence model
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
Prt 1;
}
Creating field dpdt
Creating field kinetic energy K
Constructing sources
No field sources present
Constructing reacting cloud
Constructing particle forces
Selecting particle force sphereDrag
Constructing cloud functions
none
Selecting dispersion model none
Selecting injection model coneNozzleInjection
Constructing 3-D injection
Selecting distribution model RosinRammler
Selecting patch interaction model standardWallInteraction
Selecting surface film model none
Selecting U integration scheme Euler
Selecting heat transfer model RanzMarshall
Selecting T integration scheme analytical
Selecting composition model singlePhaseMixture
Selecting phase change model liquidEvaporationBoil
Participating liquid species:
C7H16
Selecting AtomizationModel none
Selecting BreakupModel ReitzDiwakar
Selecting StochasticCollisionModel none
Average parcel mass: 2.4e-10
Selecting radiationModel none
Courant Number mean: 0 max: 0
PIMPLE: Operating solver in PISO mode
Starting time loop
Courant Number mean: 0 max: 0
deltaT = 2.94118e-06
Time = 2.94118e-06
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 58 new parcels
Cloud: sprayCloud
Current number of parcels = 58
Current mass in system = 6.86832e-10
Linear momentum = (2.48554e-11 -9.24816e-09 7.76549e-11)
|Linear momentum| = 9.24852e-09
Linear kinetic energy = 7.11574e-08
Rotational kinetic energy = 0
Total number of parcels added = 58
Total mass introduced = 6.86839e-10
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 8.6095e-15
D10, D32, Dmax (mu) = 69.0532, 93.5336, 144.356
Liquid penetration 95% mass (m) = 6.42925e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (0 -9.38861e-11 0), (8.47866e-13 0 1.22648e-12)
UCoeff min/max = 0, 6.66618e-12
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 9.42829e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 3.83676e-08, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 9.42842e-08, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 1, Final residual = 9.43118e-08, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.809069, Final residual = 7.63048e-08, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -1.40149e-06, 0
hsCoeff min/max = 0, 2.92499e-09
DILUPBiCG: Solving for hs, Initial residual = 0.998334, Final residual = 9.44958e-08, No Iterations 2
T gas min/max = 799.765, 800
DICPCG: Solving for p, Initial residual = 0.999999, Final residual = 0.0965782, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.71251e-10, global = -1.21366e-10, cumulative = -1.21366e-10
DICPCG: Solving for p, Initial residual = 0.0428956, Final residual = 4.93705e-07, No Iterations 17
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.53569e-15, global = -5.70368e-16, cumulative = -1.21367e-10
DILUPBiCG: Solving for epsilon, Initial residual = 0.00241979, Final residual = 2.64192e-07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 4.63677e-08, No Iterations 2
ExecutionTime = 171.62 s ClockTime = 172 s
Courant Number mean: 9.34524e-08 max: 0.000139867
deltaT = 3.36134e-06
Time = 6.30252e-06
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 67 new parcels
Cloud: sprayCloud
Current number of parcels = 125
Current mass in system = 2.52497e-09
Linear momentum = (-3.30581e-10 -6.64294e-08 7.40131e-11)
|Linear momentum| = 6.64302e-08
Linear kinetic energy = 9.9954e-07
Rotational kinetic energy = 0
Total number of parcels added = 125
Total mass introduced = 2.52501e-09
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 5.59051e-14
D10, D32, Dmax (mu) = 67.4864, 91.1281, 146.808
Liquid penetration 95% mass (m) = 0.000148969
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.16057e-11 -1.6836e-09 -3.34498e-12), (0 0 0)
UCoeff min/max = 0, 5.60121e-11
DILUPBiCG: Solving for Ux, Initial residual = 0.367917, Final residual = 1.08466e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.91455, Final residual = 1.06953e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.365244, Final residual = 1.08551e-08, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.57861, Final residual = 8.32285e-09, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.564599, Final residual = 8.64363e-09, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -7.9049e-06, 0
hsCoeff min/max = 0, 1.65429e-08
DILUPBiCG: Solving for hs, Initial residual = 0.615201, Final residual = 1.27401e-08, No Iterations 2
T gas min/max = 798.379, 800
DICPCG: Solving for p, Initial residual = 0.64303, Final residual = 0.0497098, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.73343e-09, global = -6.78791e-10, cumulative = -8.00158e-10
DICPCG: Solving for p, Initial residual = 0.0262455, Final residual = 5.32277e-07, No Iterations 18
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.3937e-14, global = -5.48745e-15, cumulative = -8.00163e-10
DILUPBiCG: Solving for epsilon, Initial residual = 0.00128521, Final residual = 1.6628e-07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.432937, Final residual = 1.66135e-08, No Iterations 2
ExecutionTime = 188.66 s ClockTime = 189 s
Courant Number mean: 2.96943e-07 max: 0.00120253
deltaT = 3.9725e-06
Time = 1.0275e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 79 new parcels
Cloud: sprayCloud
Current number of parcels = 204
Current mass in system = 6.14558e-09
Linear momentum = (-6.45333e-10 -2.49111e-07 -1.4119e-09)
|Linear momentum| = 2.49116e-07
Linear kinetic energy = 5.7658e-06
Rotational kinetic energy = 0
Total number of parcels added = 204
Total mass introduced = 6.14577e-09
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 2.26456e-13
D10, D32, Dmax (mu) = 64.0571, 89.4653, 144.392
Liquid penetration 95% mass (m) = 0.000320529
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-9.54018e-12 -1.11601e-08 -1.3036e-10), (0 0 0)
UCoeff min/max = 0, 2.454e-10
DILUPBiCG: Solving for Ux, Initial residual = 0.338231, Final residual = 2.81292e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.687614, Final residual = 1.73343e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.344486, Final residual = 2.83836e-08, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.432171, Final residual = 1.40554e-08, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.429964, Final residual = 1.42596e-08, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -2.81043e-05, 0
hsCoeff min/max = 0, 5.93312e-08
DILUPBiCG: Solving for hs, Initial residual = 0.456155, Final residual = 2.35866e-08, No Iterations 2
T gas min/max = 793.481, 800
DICPCG: Solving for p, Initial residual = 0.549988, Final residual = 0.0549441, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.00624e-08, global = -4.51785e-09, cumulative = -5.31802e-09
DICPCG: Solving for p, Initial residual = 0.034133, Final residual = 6.8206e-07, No Iterations 21
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.37269e-13, global = -4.09025e-14, cumulative = -5.31806e-09
DILUPBiCG: Solving for epsilon, Initial residual = 0.00141658, Final residual = 2.43786e-07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.257031, Final residual = 1.6415e-08, No Iterations 2
ExecutionTime = 205.78 s ClockTime = 206 s
Courant Number mean: 9.28205e-07 max: 0.00666354
deltaT = 4.41389e-06
Time = 1.46889e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 88 new parcels
Cloud: sprayCloud
Current number of parcels = 292
Current mass in system = 1.20085e-08
Linear momentum = (4.127e-09 -6.61355e-07 1.70989e-09)
|Linear momentum| = 6.6137e-07
Linear kinetic energy = 2.07704e-05
Rotational kinetic energy = 0
Total number of parcels added = 292
Total mass introduced = 1.20091e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 7.411e-13
D10, D32, Dmax (mu) = 41.0838, 79.6885, 144.408
Liquid penetration 95% mass (m) = 0.000528636
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (0 -4.59268e-08 0), (5.96149e-10 0 1.94492e-10)
UCoeff min/max = 0, 7.54989e-10
DILUPBiCG: Solving for Ux, Initial residual = 0.309687, Final residual = 5.80298e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.49843, Final residual = 1.39279e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.304605, Final residual = 5.97195e-08, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.360428, Final residual = 4.31513e-08, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.359916, Final residual = 4.27498e-08, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -7.86123e-05, 0
hsCoeff min/max = 0, 1.69867e-07
DILUPBiCG: Solving for hs, Initial residual = 0.367344, Final residual = 5.30469e-08, No Iterations 2
T gas min/max = 780, 800
DICPCG: Solving for p, Initial residual = 0.486128, Final residual = 0.0354199, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.60414e-08, global = -1.19429e-08, cumulative = -1.72609e-08
DICPCG: Solving for p, Initial residual = 0.0273464, Final residual = 6.35212e-07, No Iterations 23
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.01073e-13, global = -1.30345e-13, cumulative = -1.7261e-08
DILUPBiCG: Solving for epsilon, Initial residual = 0.00146717, Final residual = 2.99171e-07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.181318, Final residual = 2.08313e-08, No Iterations 2
ExecutionTime = 222.96 s ClockTime = 223 s
Courant Number mean: 2.83293e-06 max: 0.0275381
deltaT = 5.04444e-06
Time = 1.97333e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 100 new parcels
Cloud: sprayCloud
Current number of parcels = 392
Current mass in system = 2.10799e-08
Linear momentum = (9.56914e-09 -1.4751e-06 1.56034e-08)
|Linear momentum| = 1.47522e-06
Linear kinetic energy = 5.94005e-05
Rotational kinetic energy = 0
Total number of parcels added = 392
Total mass introduced = 2.1082e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 2.33825e-12
D10, D32, Dmax (mu) = 27.9901, 72.5879, 144.431
Liquid penetration 95% mass (m) = 0.000837753
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (0 -8.71829e-08 0), (9.16166e-10 0 1.4076e-09)
UCoeff min/max = 0, 1.05381e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.291732, Final residual = 9.92258e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.420606, Final residual = 3.851e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.294936, Final residual = 9.6699e-08, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.345213, Final residual = 3.91466e-07, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.345061, Final residual = 3.86469e-07, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000102813, 0
hsCoeff min/max = 0, 2.23981e-07
DILUPBiCG: Solving for hs, Initial residual = 0.317572, Final residual = 3.37885e-07, No Iterations 2
T gas min/max = 764.653, 800
DICPCG: Solving for p, Initial residual = 0.363512, Final residual = 0.0313712, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.10689e-08, global = -2.93565e-08, cumulative = -4.66175e-08
DICPCG: Solving for p, Initial residual = 0.0277295, Final residual = 7.58246e-07, No Iterations 26
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.6323e-12, global = -3.87776e-13, cumulative = -4.66179e-08
DILUPBiCG: Solving for epsilon, Initial residual = 0.00157133, Final residual = 4.08499e-07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.1492, Final residual = 5.50521e-08, No Iterations 2
ExecutionTime = 240.27 s ClockTime = 240 s
Courant Number mean: 8.53039e-06 max: 0.0888915
deltaT = 5.04444e-06
Time = 2.47778e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 100 new parcels
Cloud: sprayCloud
Current number of parcels = 492
Current mass in system = 3.26793e-08
Linear momentum = (2.21837e-08 -2.76592e-06 -8.63988e-09)
|Linear momentum| = 2.76602e-06
Linear kinetic energy = 0.000136179
Rotational kinetic energy = 0
Total number of parcels added = 492
Total mass introduced = 3.2685e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 6.20352e-12
D10, D32, Dmax (mu) = 23.5977, 67.0428, 142.619
Liquid penetration 95% mass (m) = 0.00120868
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.17232e-12 -1.63964e-07 -1.90347e-09), (1.08475e-09 0 1.05003e-10)
UCoeff min/max = 0, 2.39823e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.244788, Final residual = 5.97451e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.314307, Final residual = 6.70634e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.247943, Final residual = 1.24733e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.286766, Final residual = 7.97117e-07, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.286722, Final residual = 7.89413e-07, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000252942, 0
hsCoeff min/max = 0, 5.99369e-07
DILUPBiCG: Solving for hs, Initial residual = 0.229877, Final residual = 5.83031e-07, No Iterations 2
T gas min/max = 744.108, 800.006
DICPCG: Solving for p, Initial residual = 0.323419, Final residual = 0.0239176, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.18555e-07, global = -4.82875e-08, cumulative = -9.49054e-08
DICPCG: Solving for p, Initial residual = 0.0242483, Final residual = 7.02696e-07, No Iterations 26
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.02498e-12, global = -7.32253e-13, cumulative = -9.49061e-08
DILUPBiCG: Solving for epsilon, Initial residual = 0.00152889, Final residual = 7.06084e-07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.127782, Final residual = 7.2883e-08, No Iterations 2
ExecutionTime = 257.64 s ClockTime = 258 s
Courant Number mean: 1.89261e-05 max: 0.188386
deltaT = 2.52222e-06
Time = 2.73e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 50 new parcels
Cloud: sprayCloud
Current number of parcels = 542
Current mass in system = 3.94262e-08
Linear momentum = (3.47807e-08 -3.65409e-06 3.07719e-09)
|Linear momentum| = 3.65426e-06
Linear kinetic energy = 0.000199232
Rotational kinetic energy = 0
Total number of parcels added = 542
Total mass introduced = 3.94354e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 9.64344e-12
D10, D32, Dmax (mu) = 22.8845, 66.0858, 141.31
Liquid penetration 95% mass (m) = 0.00141654
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.66483e-10 -1.03091e-07 -3.18872e-10), (1.06944e-09 0 6.63136e-10)
UCoeff min/max = 0, 1.47933e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.126808, Final residual = 3.28206e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.127639, Final residual = 3.06469e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.127783, Final residual = 2.79724e-08, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.117677, Final residual = 2.1584e-07, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.117669, Final residual = 2.16308e-07, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000147714, 0
hsCoeff min/max = 0, 3.8061e-07
DILUPBiCG: Solving for hs, Initial residual = 0.080231, Final residual = 1.58141e-07, No Iterations 2
T gas min/max = 730.158, 800.013
DICPCG: Solving for p, Initial residual = 0.196334, Final residual = 0.0125171, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.80612e-08, global = -4.90818e-08, cumulative = -1.43988e-07
DICPCG: Solving for p, Initial residual = 0.0106665, Final residual = 7.22987e-07, No Iterations 12
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.90372e-12, global = -1.25463e-12, cumulative = -1.43989e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.000946986, Final residual = 5.67837e-07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.0858123, Final residual = 4.20183e-08, No Iterations 2
ExecutionTime = 274.46 s ClockTime = 274 s
Courant Number mean: 1.25406e-05 max: 0.114338
deltaT = 2.06364e-06
Time = 2.93636e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 41 new parcels
Cloud: sprayCloud
Current number of parcels = 583
Current mass in system = 4.54156e-08
Linear momentum = (2.77707e-08 -4.48242e-06 -1.04128e-08)
|Linear momentum| = 4.48252e-06
Linear kinetic energy = 0.000262078
Rotational kinetic energy = 0
Total number of parcels added = 583
Total mass introduced = 4.54289e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 1.37777e-11
D10, D32, Dmax (mu) = 22.4595, 64.0942, 139.989
Liquid penetration 95% mass (m) = 0.00159155
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-4.38526e-10 -1.0887e-07 -9.57572e-10), (1.06281e-09 0 2.43636e-10)
UCoeff min/max = 0, 1.42029e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.102998, Final residual = 1.13931e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.102827, Final residual = 9.02004e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.103683, Final residual = 1.10033e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.0957183, Final residual = 7.86697e-07, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.0957138, Final residual = 8.21449e-07, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.0001344, 0
hsCoeff min/max = 0, 3.53093e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0642986, Final residual = 5.25292e-07, No Iterations 2
T gas min/max = 726.877, 800.016
DICPCG: Solving for p, Initial residual = 0.11215, Final residual = 0.00505941, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.40848e-08, global = -2.44978e-08, cumulative = -1.68487e-07
DICPCG: Solving for p, Initial residual = 0.00537725, Final residual = 9.85439e-07, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.66704e-12, global = -2.24726e-12, cumulative = -1.68489e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.00140632, Final residual = 2.04466e-06, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.120519, Final residual = 3.87412e-07, No Iterations 2
ExecutionTime = 291.21 s ClockTime = 291 s
Courant Number mean: 1.3028e-05 max: 0.106824
deltaT = 1.87603e-06
Time = 3.12397e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 37 new parcels
Cloud: sprayCloud
Current number of parcels = 620
Current mass in system = 5.12258e-08
Linear momentum = (1.69548e-08 -5.31244e-06 1.14537e-08)
|Linear momentum| = 5.31247e-06
Linear kinetic energy = 0.000328393
Rotational kinetic energy = 0
Total number of parcels added = 620
Total mass introduced = 5.12449e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 1.97628e-11
D10, D32, Dmax (mu) = 21.8154, 62.1078, 138.616
Liquid penetration 95% mass (m) = 0.00175031
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-3.1184e-10 -1.0379e-07 -2.42166e-10), (5.41135e-10 0 1.59562e-09)
UCoeff min/max = 0, 1.21747e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0974742, Final residual = 6.10677e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.106244, Final residual = 4.59366e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.098691, Final residual = 6.39474e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.0985098, Final residual = 3.83525e-06, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.0985065, Final residual = 3.52769e-06, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000119989, 0
hsCoeff min/max = 0, 3.17301e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0962426, Final residual = 3.56018e-06, No Iterations 2
T gas min/max = 733.921, 800.025
DICPCG: Solving for p, Initial residual = 0.128375, Final residual = 0.00907655, No Iterations 3
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.26624e-08, global = -5.12554e-08, cumulative = -2.19745e-07
DICPCG: Solving for p, Initial residual = 0.00988771, Final residual = 6.67112e-07, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.02712e-12, global = -2.03123e-12, cumulative = -2.19747e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.00229067, Final residual = 7.69595e-06, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.146113, Final residual = 2.77459e-06, No Iterations 2
ExecutionTime = 307.96 s ClockTime = 308 s
Courant Number mean: 1.45637e-05 max: 0.10161
deltaT = 1.70548e-06
Time = 3.29452e-05
Selected 0 cells for refinement out of 168100.
Selected 0 split points out of a possible 0.
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 34 new parcels
Cloud: sprayCloud
Current number of parcels = 654
Current mass in system = 5.68082e-08
Linear momentum = (4.02628e-08 -6.13554e-06 -1.91919e-08)
|Linear momentum| = 6.1357e-06
Linear kinetic energy = 0.000396941
Rotational kinetic energy = 0
Total number of parcels added = 654
Total mass introduced = 5.6836e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 2.83631e-11
D10, D32, Dmax (mu) = 19.4152, 59.3464, 143.326
Liquid penetration 95% mass (m) = 0.0018948
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-4.62395e-10 -1.04345e-07 -2.73912e-09), (1.67291e-09 0 4.62437e-10)
UCoeff min/max = 0, 1.0747e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0868041, Final residual = 1.15625e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.105333, Final residual = 7.18322e-08, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.0876434, Final residual = 1.13178e-06, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.111547, Final residual = 8.75853e-06, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.111544, Final residual = 1.00797e-07, No Iterations 3
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -9.60544e-05, 0
hsCoeff min/max = 0, 2.57001e-07
DILUPBiCG: Solving for hs, Initial residual = 0.101794, Final residual = 8.52745e-06, No Iterations 2
T gas min/max = 744.184, 800.037
DICPCG: Solving for p, Initial residual = 0.0967298, Final residual = 0.00567746, No Iterations 3
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.35679e-08, global = -3.74624e-08, cumulative = -2.57209e-07
DICPCG: Solving for p, Initial residual = 0.00608439, Final residual = 5.84036e-07, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.20856e-12, global = -2.01934e-12, cumulative = -2.57211e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.00285522, Final residual = 5.92026e-08, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.13297, Final residual = 6.74214e-06, No Iterations 2
ExecutionTime = 326.28 s ClockTime = 326 s
Courant Number mean: 1.59187e-05 max: 0.09728
deltaT = 1.70548e-06
Time = 3.46506e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 34 new parcels
Cloud: sprayCloud
Current number of parcels = 688
Current mass in system = 6.26754e-08
Linear momentum = (1.902e-08 -7.03297e-06 -3.38234e-08)
|Linear momentum| = 7.03308e-06
Linear kinetic energy = 0.000475468
Rotational kinetic energy = 0
Total number of parcels added = 688
Total mass introduced = 6.27162e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 4.14912e-11
D10, D32, Dmax (mu) = 16.5973, 57.55, 141.334
Liquid penetration 95% mass (m) = 0.00203765
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-6.32474e-10 -1.35158e-07 -1.61225e-09), (8.47988e-10 0 6.34746e-10)
UCoeff min/max = 0, 1.29996e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0820898, Final residual = 1.21541e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.086671, Final residual = 9.11667e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0827013, Final residual = 1.27178e-06, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.10564, Final residual = 9.12078e-06, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.105637, Final residual = 1.86113e-07, No Iterations 3
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000111468, 0
hsCoeff min/max = 0, 2.97527e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0782431, Final residual = 9.65057e-06, No Iterations 2
T gas min/max = 748.091, 800.063
DICPCG: Solving for p, Initial residual = 0.0775124, Final residual = 0.00450934, No Iterations 3
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.06439e-08, global = -3.53593e-08, cumulative = -2.9257e-07
DICPCG: Solving for p, Initial residual = 0.00470496, Final residual = 5.38059e-07, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.70416e-12, global = -2.50004e-12, cumulative = -2.92573e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.0036728, Final residual = 1.09868e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.102366, Final residual = 5.20097e-06, No Iterations 2
ExecutionTime = 343.09 s ClockTime = 343 s
Courant Number mean: 1.91274e-05 max: 0.0987231
deltaT = 1.70548e-06
Time = 3.63561e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 34 new parcels
Cloud: sprayCloud
Current number of parcels = 722
Current mass in system = 6.88248e-08
Linear momentum = (-8.67466e-09 -8e-06 -4.69766e-09)
|Linear momentum| = 8e-06
Linear kinetic energy = 0.000563199
Rotational kinetic energy = 0
Total number of parcels added = 722
Total mass introduced = 6.88856e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 6.16011e-11
D10, D32, Dmax (mu) = 14.5895, 55.6291, 141.429
Liquid penetration 95% mass (m) = 0.00220233
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-2.93565e-09 -1.55501e-07 -1.89032e-09), (1.05282e-09 0 1.42049e-09)
UCoeff min/max = 0, 1.44379e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0849134, Final residual = 2.11049e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0828639, Final residual = 4.49448e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0846775, Final residual = 1.67493e-06, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.109398, Final residual = 1.0439e-07, No Iterations 3
DILUPBiCG: Solving for O2, Initial residual = 0.109397, Final residual = 1.11474e-07, No Iterations 3
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000129119, 0
hsCoeff min/max = 0, 3.52261e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0601803, Final residual = 7.51948e-06, No Iterations 2
T gas min/max = 747.997, 800.077
DICPCG: Solving for p, Initial residual = 0.0730554, Final residual = 0.00395071, No Iterations 3
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.24738e-08, global = -3.62694e-08, cumulative = -3.28842e-07
DICPCG: Solving for p, Initial residual = 0.00395339, Final residual = 4.65841e-07, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.91023e-12, global = -2.7463e-12, cumulative = -3.28845e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.0048954, Final residual = 1.39163e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.0872726, Final residual = 4.38874e-06, No Iterations 2
ExecutionTime = 359.96 s ClockTime = 360 s
Courant Number mean: 2.30643e-05 max: 0.101556
deltaT = 1.51599e-06
Time = 3.78721e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 30 new parcels
Cloud: sprayCloud
Current number of parcels = 752
Current mass in system = 7.45243e-08
Linear momentum = (6.27935e-09 -8.92363e-06 -1.87639e-08)
|Linear momentum| = 8.92365e-06
Linear kinetic energy = 0.000650954
Rotational kinetic energy = 0
Total number of parcels added = 752
Total mass introduced = 7.46124e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 8.87846e-11
D10, D32, Dmax (mu) = 12.9467, 54.144, 141.986
Liquid penetration 95% mass (m) = 0.00235394
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.28828e-09 -1.55784e-07 -1.41838e-09), (1.20849e-09 0 1.16865e-09)
UCoeff min/max = 0, 1.37915e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.077595, Final residual = 1.25677e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0731772, Final residual = 5.9657e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0768393, Final residual = 1.11595e-06, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.109527, Final residual = 1.37845e-07, No Iterations 3
DILUPBiCG: Solving for O2, Initial residual = 0.109527, Final residual = 1.64348e-07, No Iterations 3
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.00012303, 0
hsCoeff min/max = 0, 3.51194e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0481022, Final residual = 2.46643e-06, No Iterations 2
T gas min/max = 748.578, 800.082
DICPCG: Solving for p, Initial residual = 0.059477, Final residual = 0.00586695, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.01876e-08, global = -6.23104e-08, cumulative = -3.91155e-07
DICPCG: Solving for p, Initial residual = 0.00545318, Final residual = 6.03011e-07, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.01855e-11, global = -3.99012e-12, cumulative = -3.91159e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.00597507, Final residual = 1.27867e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.0713387, Final residual = 3.35765e-06, No Iterations 2
ExecutionTime = 376.81 s ClockTime = 377 s
Courant Number mean: 2.41581e-05 max: 0.0940003
deltaT = 1.51599e-06
Time = 3.93881e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 30 new parcels
Cloud: sprayCloud
Current number of parcels = 782
Current mass in system = 8.04409e-08
Linear momentum = (5.98232e-09 -9.90764e-06 -1.51737e-08)
|Linear momentum| = 9.90765e-06
Linear kinetic energy = 0.000747276
Rotational kinetic energy = 0
Total number of parcels added = 782
Total mass introduced = 8.05677e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 1.27535e-10
D10, D32, Dmax (mu) = 11.4735, 52.7039, 142.684
Liquid penetration 95% mass (m) = 0.00251162
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.23604e-09 -1.59596e-07 -1.85843e-09), (1.60894e-09 0 1.38335e-09)
UCoeff min/max = 0, 1.33714e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0750896, Final residual = 7.07386e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0719852, Final residual = 5.33043e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.075036, Final residual = 7.33536e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.102729, Final residual = 1.08384e-07, No Iterations 3
DILUPBiCG: Solving for O2, Initial residual = 0.102728, Final residual = 3.4246e-07, No Iterations 3
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000140034, 0
hsCoeff min/max = 0, 3.93686e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0459169, Final residual = 2.92925e-06, No Iterations 2
T gas min/max = 742.183, 800.085
DICPCG: Solving for p, Initial residual = 0.0563951, Final residual = 0.00562846, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.00401e-07, global = -7.04533e-08, cumulative = -4.61613e-07
DICPCG: Solving for p, Initial residual = 0.00528073, Final residual = 5.41269e-07, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.05504e-11, global = -4.13642e-12, cumulative = -4.61617e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.00800659, Final residual = 1.89887e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.0674135, Final residual = 3.52194e-06, No Iterations 2
ExecutionTime = 393.62 s ClockTime = 394 s
Courant Number mean: 2.8311e-05 max: 0.0997731
deltaT = 1.51599e-06
Time = 4.09041e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 30 new parcels
Cloud: sprayCloud
Current number of parcels = 812
Current mass in system = 8.65607e-08
Linear momentum = (5.88892e-12 -1.09502e-05 5.07803e-09)
|Linear momentum| = 1.09502e-05
Linear kinetic energy = 0.000853145
Rotational kinetic energy = 0
Total number of parcels added = 812
Total mass introduced = 8.67514e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 1.9149e-10
D10, D32, Dmax (mu) = 10.3552, 51.4864, 138.287
Liquid penetration 95% mass (m) = 0.00263581
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.50647e-09 -1.68639e-07 -2.09207e-09), (1.70284e-09 0 1.77933e-09)
UCoeff min/max = 0, 1.35379e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0720314, Final residual = 6.13108e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.06978, Final residual = 5.22537e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0721308, Final residual = 6.079e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.112133, Final residual = 1.84136e-07, No Iterations 3
DILUPBiCG: Solving for O2, Initial residual = 0.112132, Final residual = 7.41345e-07, No Iterations 3
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000152436, 0
hsCoeff min/max = 0, 4.39888e-07
DILUPBiCG: Solving for hs, Initial residual = 0.043013, Final residual = 2.52253e-06, No Iterations 2
T gas min/max = 735.844, 800.098
DICPCG: Solving for p, Initial residual = 0.0535091, Final residual = 0.00505205, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.03714e-07, global = -7.07083e-08, cumulative = -5.32325e-07
DICPCG: Solving for p, Initial residual = 0.00478496, Final residual = 4.86865e-07, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.08904e-11, global = -4.24368e-12, cumulative = -5.32329e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.0106358, Final residual = 2.66404e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.0642187, Final residual = 3.48501e-06, No Iterations 2
ExecutionTime = 410.44 s ClockTime = 411 s
Courant Number mean: 3.29324e-05 max: 0.105941
deltaT = 1.29942e-06
Time = 4.22035e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 25 new parcels
Cloud: sprayCloud
Current number of parcels = 837
Current mass in system = 9.19316e-08
Linear momentum = (1.69657e-08 -1.18961e-05 -7.56753e-09)
|Linear momentum| = 1.18961e-05
Linear kinetic energy = 0.000952955
Rotational kinetic energy = 0
Total number of parcels added = 837
Total mass introduced = 9.22273e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 2.96487e-10
D10, D32, Dmax (mu) = 9.15515, 51.5598, 144.192
Liquid penetration 95% mass (m) = 0.00274436
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.61118e-09 -1.51188e-07 -1.89482e-09), (1.46205e-09 0 1.67868e-09)
UCoeff min/max = 0, 1.29011e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0592101, Final residual = 2.72358e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0570737, Final residual = 3.85805e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0592149, Final residual = 3.0256e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.121565, Final residual = 1.17685e-07, No Iterations 3
DILUPBiCG: Solving for O2, Initial residual = 0.121564, Final residual = 2.00152e-07, No Iterations 3
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000144666, 0
hsCoeff min/max = 0, 4.24544e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0371872, Final residual = 1.99989e-06, No Iterations 2
T gas min/max = 730.969, 800.1
DICPCG: Solving for p, Initial residual = 0.0443314, Final residual = 0.00330739, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.57399e-08, global = -5.69659e-08, cumulative = -5.89295e-07
DICPCG: Solving for p, Initial residual = 0.00329679, Final residual = 3.14236e-07, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.73475e-12, global = -3.20372e-12, cumulative = -5.89298e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.0117704, Final residual = 1.9136e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.0523725, Final residual = 1.94378e-06, No Iterations 2
ExecutionTime = 427.25 s ClockTime = 427 s
Courant Number mean: 3.18656e-05 max: 0.0955914
deltaT = 1.29942e-06
Time = 4.35029e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 25 new parcels
Cloud: sprayCloud
Current number of parcels = 861
Current mass in system = 9.74444e-08
Linear momentum = (4.74701e-08 -1.28535e-05 -3.34553e-08)
|Linear momentum| = 1.28536e-05
Linear kinetic energy = 0.00105365
Rotational kinetic energy = 0
Total number of parcels added = 862
Total mass introduced = 9.78083e-08
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 3.64733e-10
D10, D32, Dmax (mu) = 17.1539, 51.1778, 142.341
Liquid penetration 95% mass (m) = 0.00285642
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.60956e-09 -1.38499e-07 -1.9771e-09), (4.27006e-09 0 1.79104e-09)
UCoeff min/max = 0, 1.39478e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0563238, Final residual = 4.08265e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0556011, Final residual = 2.12413e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0557207, Final residual = 3.90029e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.0930332, Final residual = 9.30786e-06, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.093033, Final residual = 1.33326e-07, No Iterations 3
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000153484, 0
hsCoeff min/max = 0, 4.589e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0353239, Final residual = 1.95405e-06, No Iterations 2
T gas min/max = 727.252, 800.096
DICPCG: Solving for p, Initial residual = 0.0402545, Final residual = 0.00308799, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.87884e-08, global = -6.16298e-08, cumulative = -6.50928e-07
DICPCG: Solving for p, Initial residual = 0.00314768, Final residual = 3.46942e-07, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.43494e-12, global = -4.24494e-12, cumulative = -6.50932e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.0143832, Final residual = 2.58414e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.0497904, Final residual = 1.98766e-06, No Iterations 2
ExecutionTime = 443.98 s ClockTime = 444 s
Courant Number mean: 3.5687e-05 max: 0.101388
deltaT = 1.08285e-06
Time = 4.45858e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 21 new parcels
Cloud: sprayCloud
Current number of parcels = 882
Current mass in system = 1.02135e-07
Linear momentum = (7.46445e-08 -1.36486e-05 -8.74937e-09)
|Linear momentum| = 1.36488e-05
Linear kinetic energy = 0.00113631
Rotational kinetic energy = 0
Total number of parcels added = 883
Total mass introduced = 1.0253e-07
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 3.9568e-10
D10, D32, Dmax (mu) = 16.519, 49.8796, 144.81
Liquid penetration 95% mass (m) = 0.00294498
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-1.89907e-09 -1.28808e-07 -2.3018e-09), (3.4545e-09 0 1.6802e-09)
UCoeff min/max = 0, 1.17197e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0443761, Final residual = 2.85356e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0457828, Final residual = 1.2962e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0438123, Final residual = 1.26766e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.0520051, Final residual = 2.81375e-06, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.0520064, Final residual = 4.14832e-06, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000129977, 0
hsCoeff min/max = 0, 4.03009e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0274896, Final residual = 6.40496e-07, No Iterations 2
T gas min/max = 724.82, 800.111
DICPCG: Solving for p, Initial residual = 0.0316172, Final residual = 0.0015172, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.17936e-08, global = -3.28635e-08, cumulative = -6.83796e-07
DICPCG: Solving for p, Initial residual = 0.00157688, Final residual = 9.61899e-07, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.79949e-11, global = -1.39185e-11, cumulative = -6.8381e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.0145192, Final residual = 1.68639e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.0399482, Final residual = 9.66014e-07, No Iterations 2
ExecutionTime = 460.66 s ClockTime = 461 s
Courant Number mean: 3.25453e-05 max: 0.0869848
deltaT = 1.08285e-06
Time = 4.56686e-05
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Solving cloud sprayCloud
--> Cloud: sprayCloud
Added 21 new parcels
Cloud: sprayCloud
Current number of parcels = 903
Current mass in system = 1.06884e-07
Linear momentum = (6.76286e-08 -1.44473e-05 1.70602e-08)
|Linear momentum| = 1.44474e-05
Linear kinetic energy = 0.00121938
Rotational kinetic energy = 0
Total number of parcels added = 904
Total mass introduced = 1.07315e-07
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Mass transfer phase change = 4.31838e-10
D10, D32, Dmax (mu) = 15.7124, 48.7406, 142.293
Liquid penetration 95% mass (m) = 0.00302726
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UTrans min/max = (-2.07749e-09 -1.29391e-07 -2.19472e-09), (2.67548e-09 0 2.28963e-09)
UCoeff min/max = 0, 1.19883e-09
DILUPBiCG: Solving for Ux, Initial residual = 0.0428498, Final residual = 1.72397e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0444204, Final residual = 1.37431e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0429048, Final residual = 1.34647e-07, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0.0416281, Final residual = 2.05073e-06, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 0.0416289, Final residual = 2.27837e-06, No Iterations 2
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
hsTrans min/max = -0.000130387, 0
hsCoeff min/max = 0, 4.15105e-07
DILUPBiCG: Solving for hs, Initial residual = 0.0264585, Final residual = 5.1751e-07, No Iterations 2
T gas min/max = 722.747, 800.124
DICPCG: Solving for p, Initial residual = 0.0311986, Final residual = 0.00159239, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.7641e-08, global = -3.74995e-08, cumulative = -7.21309e-07
DICPCG: Solving for p, Initial residual = 0.00168385, Final residual = 2.02216e-07, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.39762e-12, global = -2.96175e-12, cumulative = -7.21312e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.0168877, Final residual = 2.26053e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.038646, Final residual = 9.9899e-07, No Iterations 2
ExecutionTime = 477.41 s ClockTime = 478 s
Courant Number mean: 3.55657e-05 max: 0.0896867
deltaT = 1.08285e-06
Time = 4.67515e-05
Selected 63 cells for refinement out of 168100.
Refined from 168100 to 168541 cells.
Cloud<ParticleType>::autoMap(TrackData&, const mapPolyMesh&) for lagrangian cloud sprayCloud
|
|
2012-08-03 17:15
|
|
|
|
I've uploaded the logfile as well as the sprayCloud_positions.obj file that was spat out. Hope this is useful. |
|
|
Have you recompiled the OpenFOAM library after pulling commit 9c8ce25? |
|
|
I did an Allwmake, but I just noticed that the OpenFOAM library didn't fully recompile because of the polynomial commit. Did a full clean and compile and it appears that helped. The mapping appears to work, though there is still a subsequent crash that I'm not sure if its linked to the lagrangian or not. I'll run some more tests and get back to you. |
|
|
I'm going to leave the cases running over the weekend, as lagrangian stuff in debug always takes forever to run. This is what the optimized version gave: Courant Number mean: 3.54942e-05 max: 0.0897122 deltaT = 1.08285e-06 Time = 4.67515e-05 Selected 66 cells for refinement out of 168100. Refined from 168100 to 168562 cells. Selected 0 split points out of a possible 66. DICPCG: Solving for pcorr, Initial residual = 0.987988, Final residual = 5.1753e-07, No Iterations 11 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 Solving cloud sprayCloud #0 Foam::error::printStack(Foam::Ostream&) in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigSegv::sigHandler(int) in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #4 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #5 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #6 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #7 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #8 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #9 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #10 __libc_start_main in "/lib/libc.so.6" #11 at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/i386/elf/start.S:122 Segmentation fault |
|
|
It's likely that this relates to the injection model as they generally cache the injection cell/tet/tet-face - which injection model are you using? this has been updated in our internal line, but not pushed to 2.1.x (as yet) |
|
|
For this simulation I was using the coneNozzle. The simulation hangs when solving for the sprayCloud when I'm running the code in debug mode, so unfortunately I haven't got a debug stack trace for you yet. Any advice for this, I know there is something to do with the array bounds checking that does odd things to lagrangian solving... |
|
|
You could try using the 'disc' injectionMethod (if not already doing so) as this will re-calculate the injection cell/tet addressing information |
|
|
Under coneNozzleCoeffs, the injection method is listed as disc. Is there another setting I should change, or try another injectionModel? EDIT: I tried the injectionMethod = point, and I didn't crash but the code hangs at evolving the spray. |
|
|
I increased the refinement frequency to each iteration, and the cloud evolution can at least move to after parcel injection (I get similar results by changing the distribution model to uniform). Selected 12 cells for refinement out of 168100. Refined from 168100 to 168184 cells. Selected 0 split points out of a possible 12. Solving cloud sprayCloud --> Cloud: sprayCloud Added 30 new parcels #0 Foam::error::printStack(Foam::Ostream&) in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigSegv::sigHandler(int) in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #4 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #5 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #6 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #7 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #8 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #9 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #10 in "/home/gandalf/OpenFOAM/OpenFOAM-2.1.x/platforms/linuxGccDPOpt/bin/sprayDyMFoam" #11 __libc_start_main in "/lib/libc.so.6" #12 at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/i386/elf/start.S:122 Segmentation fault Still trying to get some better debug info. Will keep posting interesting results. |
|
|
The final piece appears to be the mapping of the cloud source fields - this is now working in our development line, and I'll look to push similar changes into 2.1.x soon |
|
|
I've pushed the field mapping under commit 8ee747a - there may be some issues to do with the cell unrefinement, but this may need to be put on hold as an item for the next release. |
|
|
Thanks for the commits! Good news is that the code stops crashing, but its hangs right after parcels are added. Not sure if this is an issue with setup, will do my best to narrow the factors. Do you have an internal tutorial case that you have been testing against that works? If so could I try it out? |
|
2012-08-10 09:47
|
|
|
|
I've been running a slightly modified version of your original attachments (sprayFoamDynamic_2.tgz) - all should work well based on commit bfe8c6e. |
|
|
Very nice, the new version is working on my machine as well. Great job! One question I have is if this will work for moving meshes as well? In solvers like pimpleDyMFoam, there is a flux correction step that occurs (usually involving calls to meshPhi). But in compressible solvers I haven't seen one that does this. This may not be the best place to continue discussion, so please feel free to email me: mturcios (at) westport (dot) com. I don't know if this issue counts as fixed/resolved. I'll start testing this on some real world cases and give some feedback, if unrefinement is an issue I can run some tests for you if you'd like. Once again, thanks for being so responsive to these reports! |
|
|
No problem - the mapping feature was lost in some of the upgrades to version 2.0 and it's good to get it up and running again. These updates should work for moving mesh cases; however, there is still the limitation that not all injection models have been updated to accommodate the mapping, which will be resolved in the future. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2012-03-16 20:12 |
|
New Issue | |
| 2012-03-16 20:12 |
|
File Added: sprayFoamDynamic.tar.gz | |
| 2012-07-26 00:01 |
|
Note Added: 0001504 | |
| 2012-07-26 10:26 |
|
Note Added: 0001506 | |
| 2012-07-26 12:27 |
|
Note Added: 0001507 | |
| 2012-07-26 18:02 |
|
Note Added: 0001509 | |
| 2012-07-27 00:31 |
|
Note Edited: 0001509 | |
| 2012-08-02 12:25 |
|
Note Added: 0001537 | |
| 2012-08-02 17:16 |
|
Note Added: 0001542 | |
| 2012-08-02 18:05 |
|
Note Added: 0001543 | |
| 2012-08-02 21:56 |
|
Note Added: 0001544 | |
| 2012-08-02 22:10 |
|
Note Edited: 0001544 | |
| 2012-08-03 15:37 |
|
Note Added: 0001551 | |
| 2012-08-03 17:15 |
|
File Added: foam.log | |
| 2012-08-03 17:15 |
|
File Added: sprayCloud_positions.obj | |
| 2012-08-03 17:16 |
|
Note Added: 0001552 | |
| 2012-08-03 17:34 |
|
Note Added: 0001553 | |
| 2012-08-03 18:48 |
|
Note Added: 0001554 | |
| 2012-08-03 21:25 |
|
Note Added: 0001555 | |
| 2012-08-06 09:47 |
|
Note Added: 0001556 | |
| 2012-08-06 09:48 |
|
Note Edited: 0001556 | |
| 2012-08-07 18:42 |
|
Note Added: 0001569 | |
| 2012-08-08 12:47 |
|
Note Added: 0001576 | |
| 2012-08-08 17:43 |
|
Note Added: 0001577 | |
| 2012-08-08 17:48 |
|
Note Edited: 0001577 | |
| 2012-08-08 18:20 |
|
Note Edited: 0001577 | |
| 2012-08-08 22:21 |
|
Note Added: 0001578 | |
| 2012-08-09 10:37 |
|
Note Added: 0001579 | |
| 2012-08-09 15:49 |
|
Note Added: 0001580 | |
| 2012-08-09 20:42 |
|
Note Added: 0001582 | |
| 2012-08-10 09:47 |
|
File Added: sprayFoamDynamic_2.tgz | |
| 2012-08-10 09:58 |
|
Note Added: 0001583 | |
| 2012-08-10 18:20 |
|
Note Added: 0001584 | |
| 2012-08-13 10:42 |
|
Note Added: 0001585 | |
| 2012-08-13 10:43 |
|
Status | new => resolved |
| 2012-08-13 10:43 |
|
Fixed in Version | => 2.1.x |
| 2012-08-13 10:43 |
|
Resolution | open => fixed |
| 2012-08-13 10:43 |
|
Assigned To | => user2 |
