|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Hafiz Ahmad
Join Date: May 2019
Location: Brunei
Posts: 19
Rep Power: 7 ![]() |
This is my first time testing a solver by using a significantly different geometry than the one provided in tutorials. Previously, I have fixed this issue (different solver) by changing the boundary conditions (the patch of U) of the test case. However, as of this solver, I am stuck. I'm not sure if it's because I have 2 inlets (an air inlet and a water inlet) or because of my boundary conditions. This is the out after I run twoPhaseEulerFoam:
Code:
/*---------------------------------------------------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 7 \\/ M anipulation | \*---------------------------------------------------------------------------*/ Build : 7-f1e51ce8d094 Exec : twoPhaseEulerFoam Date : Nov 15 2019 Time : 08:34:29 Host : "debian" PID : 7848 I/O : uncollated Case : /home/hafiz/OpenFOAM/hafiz-7/run/jetSpray nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 PIMPLE: No convergence criteria found PIMPLE: No corrector convergence criteria found Calculations will do 3 corrections PIMPLE: Operating solver in transient mode with 3 outer correctors Reading g Reading hRef Creating twoPhaseSystem Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } Calculating face flux field phi.air Selecting diameterModel for phase air: isothermal Selecting turbulence model type RAS Selecting RAS turbulence model mixtureKEpsilon RAS { RASModel mixtureKEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 1.92; Cp 0.25; sigmak 1; sigmaEps 1.3; } Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo eConst; equationOfState perfectFluid; specie specie; energy sensibleInternalEnergy; } Calculating face flux field phi.water Selecting diameterModel for phase water: constant Selecting turbulence model type RAS Selecting RAS turbulence model mixtureKEpsilon RAS { RASModel mixtureKEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 1.92; Cp 0.25; sigmak 1; sigmaEps 1.3; } Selecting default blending method: linear Selecting aspectRatioModel for (air in water): constant Selecting aspectRatioModel for (water in air): constant Selecting dragModel for (air and water): segregated Selecting swarmCorrection for (air and water): none Selecting dragModel for (air in water): SchillerNaumann Selecting swarmCorrection for (air in water): none Selecting dragModel for (water in air): SchillerNaumann Selecting swarmCorrection for (water in air): none Selecting virtualMassModel for (air in water): constantCoefficient Selecting virtualMassModel for (water in air): constantCoefficient Selecting heatTransferModel for (air in water): RanzMarshall Selecting heatTransferModel for (water in air): RanzMarshall Calculating field g.h Reading field p_rgh Creating field dpdt Creating field kinetic energy K No MRF models present No finite volume options present Courant Number mean: 0.0002 max: 0.0526316 Max Ur Courant Number = 0.088355 Calculating field DDtU1 and DDtU2 Starting time loop fieldAverage fieldAverage1: Starting averaging at time 0 Courant Number mean: 0.0002 max: 0.0526316 Max Ur Courant Number = 0.088355 Time = 1e-06 PIMPLE: Iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 1 Min(alpha.air) = 1 Max(alpha.air) = 1 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 1, Final residual = 1.71606e-09, No Iterations 1 smoothSolver: Solving for e.water, Initial residual = 1, Final residual = 2.2156e-16, No Iterations 1 min T.air 300 min T.water 300 GAMG: Solving for p_rgh, Initial residual = 0.996029, Final residual = 1.67269e-11, No Iterations 2 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 ? at tensorField.C:? #4 ? in "/home/hafiz/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/twoPhaseEulerFoam" #5 ? in "/home/hafiz/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/twoPhaseEulerFoam" #6 ? in "/home/hafiz/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/twoPhaseEulerFoam" #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #8 ? in "/home/hafiz/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/twoPhaseEulerFoam" Floating point exception Code:
checkMesh Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 25704 internal points: 0 faces: 50501 internal faces: 24799 cells: 12550 faces per cell: 6 boundary patches: 5 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 12550 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. ***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology waterInlet 1 4 ok (non-closed singly connected) airInlet 50 102 ok (non-closed singly connected) outlet 50 102 ok (non-closed singly connected) walls 25500 25704 ok (non-closed singly connected) defaultFaces 101 204 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (0 0 0) (0.25 0.05 0.0001) Mesh has 1 geometric (non-empty/wedge) directions (1 0 0) Mesh has 1 solution (non-empty) directions (1 0 0) All edges aligned with or perpendicular to non-empty directions. Boundary openness (7.92192e-20 1.45684e-18 5.48246e-15) OK. Max cell openness = 1.05879e-16 OK. Max aspect ratio = 1 OK. Minimum face area = 5e-08. Maximum face area = 1e-06. Face area magnitudes OK. Min volume = 5e-11. Max volume = 1e-10. Total volume = 1.25e-06. Cell volumes OK. Mesh non-orthogonality Max: 0 average: 0 Non-orthogonality check OK. Face pyramids OK. Max skewness = 2.74172e-13 OK. Coupled point location match (average 0) OK. Mesh OK. End U.air Code:
/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 7 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format binary; class volVectorField; object U.air; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { waterInlet { type fixedValue; value uniform (0 0 0); } airInlet { type fixedValue; value uniform (100 0 0); } outlet { type pressureInletOutletVelocity; phi phi.air; value $internalField; } walls { type fixedValue; value uniform (0 0 0); } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 7 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format binary; class volVectorField; object U.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { waterInlet { type fixedValue; value uniform (0 176.71 0); } airInlet { type fixedValue; value uniform (0 0 0); } outlet { type pressureInletOutletVelocity; phi phi.water; value $internalField;; } walls { type fixedValue; value uniform (0 0 0); } } // ************************************************************************* // PS: I am trying to simulate a water spray which comes from a small nozzle at the top left of the geometry and airflow comes from the inlet to further atomize the already present droplets. I know that lagrangian approach would be better than euler based solver but I read that some people have used it. PSS: I won't be replying for another 6-7 hours because it's already midnight where I am (just in case somebody is replies to this thread) |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries | NickG | OpenFOAM Installation | 3 | December 30, 2019 00:21 |
[blockMesh] blockMesh with double grading. | spwater | OpenFOAM Meshing & Mesh Conversion | 92 | January 12, 2019 09:00 |
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh | gschaider | OpenFOAM Community Contributions | 300 | October 29, 2014 18:00 |
OpenFOAM without MPI | kokizzu | OpenFOAM Installation | 4 | May 26, 2014 09:17 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 17:51 |