|
[Sponsors] |
April 3, 2018, 09:02 |
Jet in cross-flow using sprayFoam
|
#1 |
Member
vishal
Join Date: Mar 2013
Posts: 73
Rep Power: 13 |
Hi...
I am trying to simulate jet in cross-flow using sprayFoam (water jet coming into air). I have used patchInjection and LISA atomization model. The phase change model, turbulence, and the heat transfer model are turned off. However, I am getting this error: #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib64/libc.so.6" #3 ? at makeBasicSprayParcelSubmodels.C:? #4 ? at ??:? #5 ? at ??:? #6 ? at ??:? #7 ? at ??:? #8 ? at ??:? #9 ? at ??:? #10 ? at ??:? #11 ? at ??:? #12 __libc_start_main in "/lib64/libc.so.6" #13 ? at ??:? /var/spool/slurmd/job01562/slurm_script: line 22: 120606 Floating point exception(core dumped) sprayFoam > log I understand it is something to do with division by zero . I have tried but unable to resolve the issue. Can anybody please tell me how to make this case work. I am attaching my case file here. I am using OF-2.4.0. |
|
April 3, 2018, 10:02 |
|
#2 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
I just had a look at your case without running it. You u file has no real outlet. You should have at least one boundary wit zeroGradient. Else you put material in your volume baut there is no chance for it to get out.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
April 3, 2018, 13:54 |
|
#3 |
Member
vishal
Join Date: Mar 2013
Posts: 73
Rep Power: 13 |
Hi Uwe, thanks for the suggestion. I tried it out but still not working. I also tried by increasing/reducing the number of parcels per second but of no help. Can you suggest me something, I am really struck here.
I have one more doubt, which injector in OF-2.2.2 represents a blob type injector? Thanks. |
|
April 4, 2018, 02:48 |
|
#4 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
Your case is a rather complicated one. If anybody should help you it would be a fine idea to explain something more:
- a sketch of the setup and what to be calculated - how are the properties of your case set? Where can one find which physical/chemical properties? - What is the base of your case? Is there a fundamental one which you modified? - How did you make safe, that the physical and chemical properties are correct and complete? What may be taken for safe because it is based on another well checked case? - Which modifications did you try? Which results did they have? - Do you write a protocol of your runs? I recommend this. Protocol a run, even if it crashes. Make sketches form the result or print them out.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
April 4, 2018, 11:15 |
|
#5 |
Member
vishal
Join Date: Mar 2013
Posts: 73
Rep Power: 13 |
Hello Uwe,
1) a sketch of the setup and what to be calculated click here link for the diagram of my case set-up. It's a jet in cross-flow simulation, where the cross-flow fluid is air and the jet fluid is water. I am trying to calculate the velocity and the droplet diameters (Sauter diameter and mean diameter) due to primary (LISA) and secondary atomization (ReitzDiwakar) using patchinjection model. My jet volume flow rate is 10e-6m3/s. 2) how are the properties of your case set? Where can one find which physical/chemical properties? As a first step, I am trying to run the simulation without combustion, turbulence or phase change. For the rest, U: both the inlets have fixed value, outlet zeroGradient (Changed after you suggested). p: zeroGradient at the inlet and 1e5 at the outlet T: both inlet have temperature of 300 k , oultet ZG N2 and O2: 0.766 and 0.234 at the cross-flow inlet and internal field. For jet inlet: zero Gradient was used. Ydefault: zeroGradient for everything. The rest variables. I have given some flag value as they are not being calculated. The physical and chemical properties are mainly copied from lagrangian/sprayFoam/aachenbomb. The sprayCloudProperties are copied from lagrangian/sprayFoam/aachenbomb tutorial and lagrangian/reactingParcelFoam/verticalchannel tutorial. 3)What is the base of your case? Is there a fundamental one which you modified? The base of my case is lagrangian/sprayFoam/aachenbomb tutorial. Here I have modified the injector type, primary atomization and geometry (blockMesh). Turned off combustion, heat transfer and turbulence. Also, I have changed the injection fluid from C7H16 to water. 4) How did you make safe, that the physical and chemical properties are correct and complete? What may be taken for safe because it is based on another well checked case? I consider the physical and chemical properties are correct and complete as they work fine for the aachenbomb tutorial. Moreover, I never found any one reporting any bug in this forum regarding the physical or chemical properties. 5) Which modifications did you try? Which results did they have? As I have mentioned in (3), I have changed the injection type, fluid, primary atomization and the numerical domain. The results obtained from the base case by running aachenbomb looks fine and correct. However, after making my changes, the simulation is blowing up. 6) Do you write a protocol of your runs? I recommend this. Protocol a run, even if it crashes. Make sketches form the result or print them out. Yes, I did write the protocol. They are attached in the .tar.gz file attached along with my case in the first post. If you have any more query, please let me know. Thanks. |
|
April 5, 2018, 01:45 |
|
#6 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
I'll have a closer look at it. Please have some patience.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
April 5, 2018, 03:22 |
|
#7 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
I used your model as you delivered it, only blockMesh before.
At least one scheme is missing: keyword div(((rho*nuEff)*dev2(T(grad(U))))) is undefined in dictionary ...at2_1/system/fvSchemes.divSchemes
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
April 5, 2018, 03:34 |
|
#8 |
Member
vishal
Join Date: Mar 2013
Posts: 73
Rep Power: 13 |
Thanks Uwe. I will incorporate this scheme and run it. Will get back to you with the result.
|
|
April 5, 2018, 12:20 |
|
#9 |
Member
vishal
Join Date: Mar 2013
Posts: 73
Rep Power: 13 |
Hi Uwe, I incorporated the div(((rho*nuEff)*dev2(T(grad(U))))) in fvScheme, the result is same. The simulation is blowing off after few time-steps.
Please click here for the log file of my simulation, along with the new case set-up. Thanks. |
|
April 6, 2018, 01:18 |
|
#10 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
I tried it again. Again, the case is not complete. I get in the first round:
Code:
--> FOAM FATAL ERROR: previous iteration field IOobject: volScalarField rho "C:/OpenFOAM/17.02/Z0Pilz-3.0.x/run/T2/0" not stored. Use field.storePrevIter() at start of iteration. From function GeometricField<Type, PatchField, GeoMesh>::prevIter() const in file /opt/OpenFOAM/OpenFOAM-3.0.x/src/OpenFOAM/lnInclude/GeometricField.C at line 839. FOAM aborting One last hint: Run the simulation without redirecting the log. Perhaps you don't have the last lines of the log stored, because the solver crashes irregularly. On the screen the log is complete. All in all, it seems to be a problem with the solver or one of it's libraries.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
Tags |
jet, lisa model, sprayfoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Strange Mean DPM Diameter along axial direction for droplet laden jet flow | HuifengGong | FLUENT | 2 | October 12, 2016 09:42 |
Laminar Wall Jet Flow over a Shallow Cavity | javadkz | FLUENT | 0 | October 3, 2016 03:41 |
Slot jet in cross flow. | Hamdy | Phoenics | 0 | March 5, 2006 05:01 |
flow simulation across a small fan | jane luo | Main CFD Forum | 15 | April 12, 2004 18:49 |
jet flow simulations | Nandu | CFX | 0 | September 1, 2002 12:30 |