|
[Sponsors] |
March 27, 2014, 05:42 |
manualInjection model in sprayFoam
|
#1 |
New Member
Join Date: Dec 2012
Posts: 12
Rep Power: 13 |
Hi,
i'm trying to use the manualInjection model with the sprayFoam Solver. Unfortunately the simulation crashes prior the first iteration. This is the corresponding error message: Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.2.x-0f3a889dd39b Exec : sprayFoam Date : Mar 27 2014 Time : 11:35:58 Host : "philipp-Home" PID : 5779 Case : /home/philipp/cfd/OF_22x/EtF1/full_Spray_EtF1_8 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 Reading g Creating combustion model Selecting combustion model PaSR<psiChemistryCombustion> Selecting chemistry type { chemistrySolver ode; chemistryThermo psi; } Selecting thermodynamics package { type hePsiThermo; mixture reactingMixture; transport sutherland; thermo janaf; energy sensibleEnthalpy; equationOfState perfectGas; specie specie; } Selecting chemistryReader chemkinReader Reading CHEMKIN thermo data in new file format chemistryModel: 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 kOmegaSST kOmegaSSTCoeffs { alphaK1 0.85034; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.85616; Prt 1; gamma1 0.5532; gamma2 0.4403; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; } Creating field dpdt Creating field kinetic energy K No finite volume options present Constructing reacting cloud Constructing particle forces Selecting particle force sphereDrag Constructing cloud functions none Constructing particle injection models Creating injector: model1 Selecting injection model manualInjection Constructing 3-D injection Selecting distribution model RosinRammler Selecting dispersion model none 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: C2H5OH Selecting AtomizationModel none Selecting BreakupModel ReitzDiwakar Selecting StochasticCollisionModel none #0 Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/sprayFoam" #4 in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/sprayFoam" #5 in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/sprayFoam" #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/sprayFoam" Gleitkomma-Ausnahme (Speicherabzug geschrieben) Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format binary; class dictionary; location "constant"; object SprayCloudProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solution { active true; coupled true; transient yes; cellValueSourceCorrection on; sourceTerms { schemes { rho explicit 1; U explicit 1; Yi explicit 1; h explicit 1; radiation explicit 1; } } interpolationSchemes { rho cell; U cellPoint; thermo:mu cell; T cell; Cp cell; p cell; kappa cell; } integrationSchemes { U Euler; T analytical; } } constantProperties { T0 294; // place holders for rho0 and Cp0 - reset from liquid props using T0 rho0 1000; Cp0 4187; youngsModulus 1e9; poissonsRatio 0.35; epsilon0 1; f0 0.5; Pr 0.7; Tvap 273; Tbp 351; constantVolume false; } subModels { particleForces { sphereDrag; } injectionModels { model1 { type manualInjection; //coneNozzleInjection; massTotal 1.25e-4; parcelBasisType mass; SOI 0; positionsFile "reactingCloud1Positions"; U0 (0 0 0); //injectionMethod disc; //flowType constantVelocity; //UMag 0; //outerDiameter 2.0e-4; //innerDiameter 0; //duration 0.1; //position ( 0 0 -0.215 ); //direction ( 0 0 1 ); //parcelsPerSecond 900000; //flowRateProfile constant 1; //velocityProfile constant 0; //Cd constant 0.9; //thetaInner constant 0.0; //thetaOuter constant 10.0; sizeDistribution { type RosinRammler; RosinRammlerDistribution { minValue 1e-06; maxValue 0.00015; d 0.00015; n 3; } } } } injectionModel manualInjection; dispersionModel none; patchInteractionModel standardWallInteraction; heatTransferModel RanzMarshall; compositionModel singlePhaseMixture; phaseChangeModel liquidEvaporationBoil; surfaceFilmModel none; atomizationModel none; breakupModel ReitzDiwakar; // ReitzKHRT; stochasticCollisionModel none; radiation off; standardWallInteractionCoeffs { type rebound; } RanzMarshallCoeffs { BirdCorrection true; } singlePhaseMixtureCoeffs { phases ( liquid { C2H5OH 1; } ); } liquidEvaporationBoilCoeffs { enthalpyTransfer enthalpyDifference; activeLiquids ( C2H5OH ); } ReitzDiwakarCoeffs { solveOscillationEq yes; Cbag 6; Cb 0.785; Cstrip 0.5; Cs 10; } /* ReitzKHRTCoeffs { solveOscillationEq yes; B0 0.61; B1 40; Ctau 1; CRT 0.1; msLimit 0.2; WeberLimit 6; } TABCoeffs { y0 0; yDot0 0; Cmu 10; Comega 8; WeCrit 12; } */ } cloudFunctions {} // ************************************************************************* // Regards |
|
April 2, 2014, 09:29 |
|
#2 |
Member
Ananda Kannan
Join Date: Feb 2014
Location: Göteborg, Sweden
Posts: 55
Rep Power: 12 |
Hi!!
If I were you, I would first run the case using the debug build of openfoam. You would get a much better insight into the origin of the error. Regards ansubru |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
eddy dissipation model: combustion doesn't occur | roukaia | FLUENT | 2 | November 18, 2020 10:23 |
Mixture model + VOF? | blacksoil2012 | Fluent Multiphase | 2 | February 8, 2014 23:59 |
eddy dissipation model: combustion doesn't occur | roukaia | FLUENT | 0 | December 24, 2011 09:10 |
Problems bout CFD model of biomass gasification, Downdraft gasifier | wanglong | FLUENT | 2 | November 25, 2009 23:27 |
Reynolds Stress model in CFX vs Fluent | Tim | CFX | 1 | October 7, 2009 06:19 |