 |
 |
 |
|
|
|
[Sponsors] | ||||
Fluctuating Liquid Penetration for Wedge sprayFoam |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 21, 2013, 05:17
|
Fluctuating Liquid Penetration for Wedge sprayFoam
|
#1 |
|
New Member
Ah Cai
Join Date: Dec 2010
Location: Malaysia
Posts: 16
Rep Power: 15  |
Hi all,
I am running a 1 degree angle wedge case using sprayFoam (OF-2.2.x) for diesel spray simulation. The injector is placed along the axis-symmetrical line. When I extract the Liquid Penetration Length (LPL), the LPL is found to be fluctuating. Refer to the attachments. However, vapour penetration length (VPL) does not face this problem. When LPL is extracted using dieselFoam (OF-2.0.x), the fluctuation did not happen. So, I am suspecting it is a coding problem in sprayFoam 2.2.x. I am still new to OpenFOAM thus solving this code problem would be a tough challenge for me. Any help with this is greatly appreciated! Attached is the sprayCloudProperties and blockMeshDict if needed. Thanks. Ong |
|
|
|
|
|
December 21, 2013, 11:07
|
|
#2 |
|
New Member
RJ HO
Join Date: Dec 2012
Posts: 21
Rep Power: 13 |
HI,
I have a question that is not relevant to yours. Do you happen to have any problem with total injected mass? I have been having problems with incorrect injected mass on wedge cases. Regards RJ |
|
|
|
|
|
|
January 6, 2014, 02:45
|
|
#3 |
|
New Member
Join Date: Aug 2013
Posts: 4
Rep Power: 12 |
Hi..
I am new to OF and trying to run a 5 degree angle wedge case using sprayFoam But I get following Error: BlockMeshDict is attacjed for refernce. Thanks In Advance: Error is:: Code:
Solving cloud sprayCloud
Cloud: sprayCloud
Current number of parcels = 6
Current mass in system = 7.94663e-10 /** Expected Mass is 8.33e-8 kg for 5 deg.**//
Linear momentum = (-2.3317e+09 0 2.94384e+09)
|Linear momentum| = 3.7554e+09
Linear kinetic energy = 3.06861e+28
Rotational kinetic energy = 0
model1:
number of parcels added = 6
mass introduced = 7.95718e-10
Parcels absorbed into film = 0
New film detached parcels = 0
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Temperature min/max = 320.367, 371.071
Mass transfer phase change = 1.055e-12
D10, D32, Dmax (mu) = 2.28102e-12, 8.01824e-12, 1089.75
Liquid penetration 95% mass (m) = 0.000159027
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 void Foam::multiply<Foam::Tensor<double> >(Foam::Field<Foam::Tensor<double> >&, Foam::UList<double> const&, Foam::UList<Foam::Tensor<double> > const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::compressible::RASModels::kEpsilon::divDevRhoReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#6
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/sprayFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/sprayFoam"
Floating point exception (core dumped)
aniket@aniket-Inspiron-N5110:~/OpenFOAM/aniket-2.2.1/run/tutorials/lagrangian/sprayFoam/aachenBomb$
Last edited by wyldckat; August 27, 2017 at 13:11. Reason: Added [CODE][/CODE] markers |
|
|
|
|
|
|
January 29, 2014, 09:08
|
|
#4 | |
|
New Member
Ah Cai
Join Date: Dec 2010
Location: Malaysia
Posts: 16
Rep Power: 15 |
Quote:
Yes. If im not mistaken i was also having a problem with incorrect mass injected. I have stopped using OF22x and revert to an older version, i.e. OF20x which seems to have less problem for my wedge case. Sorry but I think I might not be much help in OF22x~ Cheers |
||
|
|
|
||
|
January 29, 2014, 09:10
|
|
#5 | |
|
New Member
Ah Cai
Join Date: Dec 2010
Location: Malaysia
Posts: 16
Rep Power: 15 |
Quote:
Did you run checkMesh?? Did it show an error? if no, then most likely is not your mesh problem, might be your settings problem~ Last edited by wyldckat; August 27, 2017 at 13:13. Reason: Added [CODE][/CODE] markers |
||
|
|
|
||
|
January 29, 2014, 17:14
|
same error different case
|
#6 |
|
New Member
david mckelvey
Join Date: Nov 2013
Posts: 14
Rep Power: 12 |
hi
i got the same error message, I'm running a transient, compressible gas pipe flow simulation, using sonicFoam i get the following error after 0.01482 sec with a time step of 1e-05, according to checkMesh my mesh is ok, I'm guessing its something to do with the settings, but I'm not sure what exactly Code:
Time = 0.01483 Courant Number mean: -1.84532e+55 max: 5.6429e+58 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 void Foam::multiply<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so" #4 Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so" #5 Foam::compressible::RASModels::RNGkEpsilon::divDevRhoReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so" #6 in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/sonicFoam" #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #8 in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/sonicFoam" Floating point exception (core dumped) any advice would be welcome David Last edited by wyldckat; August 27, 2017 at 13:13. Reason: Added [CODE][/CODE] markers |
|
|
|
|
|
|
January 29, 2014, 18:00
|
|
#7 |
|
Senior Member
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 740
Rep Power: 28 |
When you say your injector is on the axis, do you mean directly on the axis? Which direction is it pointing? And how are you extracting the penetration length?
|
|
|
|
|
|
|
February 7, 2014, 10:19
|
|
#8 | |
|
New Member
Join Date: Aug 2013
Posts: 4
Rep Power: 12 |
Quote:
First, pardon me for too late reply,, I ran checkMesh ..and it showed: (note: here BC used is cyclic) Code:
Create time
Create polyMesh for time = 0
Time = 0
Mesh stats
points: 174023
faces: 510081
internal faces: 494419
cells: 168100
faces per cell: 5.97561
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0
Overall number of cells of each type:
hexahedra: 164000
prisms: 4100
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0
Checking topology...
Boundary definition OK.
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
front 4100 4242 ok (non-closed singly connected)
back 4100 4242 ok (non-closed singly connected)
tankWall 4100 4242 ok (non-closed singly connected)
inlet 1681 1723 ok (non-closed singly connected)
outlet 1681 1723 ok (non-closed singly connected)
axis 0 0 ok (empty)
Checking geometry...
Overall domain bounding box (0 -0.00349 0) (0.099939 0.00349 1)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (5.27124e-15 -6.19357e-17 9.43227e-18) OK.
***High aspect ratio cells found, Max aspect ratio: 5.98308e+194, number of cells 168100
<<Writing 168100 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 5.06068e-09. Maximum face area = 2.43902e-05. Face area magnitudes OK.
Min volume = 1.66667e-300. Max volume = 2e-300. Total volume = 3.34833e-295. Cell volumes OK.
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 acos in "/lib/x86_64-linux-gnu/libm.so.6"
#4 Foam::polyMesh::checkFaceOrthogonality(Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, bool, bool, Foam::HashSet<int, Foam::Hash<int> >*) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/checkMesh"
#6
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/checkMesh"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/checkMesh"
Floating point exception (core dumped)
New blockMesh for cyclic BC is attached here Last edited by wyldckat; August 27, 2017 at 13:14. Reason: Added [CODE][/CODE] markers |
||
|
|
|
||
|
March 19, 2014, 09:59
|
|
#9 | |
|
New Member
Join Date: Nov 2013
Posts: 1
Rep Power: 0 |
Quote:
you need to put neighbour patch as front and back instead of axis. Last edited by wyldckat; August 27, 2017 at 13:14. Reason: Added [CODE][/CODE] markers |
||
|
|
|
||
|
January 5, 2015, 15:34
|
|
#10 | |
|
New Member
Ehsan
Join Date: Oct 2014
Location: Iran
Posts: 2
Rep Power: 0 |
Quote:
my project is exactly similiar to your work and I am running a 1 degree angle wedge case using sprayFoam (OF-2.3.x) for diesel spray simulation. The injector is placed along the axis-symmetrical line. can you help me how to extract the Liquid Penetration Length (LPL) and vapour penetration length (VPL)? thanks. |
||
|
|
|
||
|
July 26, 2015, 07:18
|
plot liquid penetration length
|
#11 | |
|
Member
alvaro
Join Date: Apr 2015
Posts: 33
Rep Power: 11 |
Quote:
I want display the parameteres LPL and VPL in paraview, and then extract the values to plot (without have to write a post-processing function if it's possilbe). Could you explain me how did you get this values? Thanks. PD: I'm using sprayFoam in OF v2.3.0 Last edited by alvariten; July 26, 2015 at 08:52. |
||
|
|
|
||
|
June 8, 2016, 08:30
|
Spray length or spray penetration
|
#12 |
|
New Member
Join Date: Jan 2016
Posts: 2
Rep Power: 0 |
Hello,
Can anybody please help me in figuring out on how to extract spray length or spray penetration from sprayFoam results. I have simulated the case using foamJob utility to extract log files. There is a context in log file "Liquid penetration 95% (m)" what does it mean? I have also used foamVTK to generate VTK output. Should I load VTK using paraview and extract the parcel data?. I am very confused rigth now. Can someone please help me to figure out the solution to my problem. Thanks in advance. Yasho |
|
|
|
|
|
|
July 8, 2017, 07:57
|
|
#13 | |
|
Member
Sadegh Ebadi
Join Date: Apr 2015
Posts: 75
Rep Power: 11 |
Quote:
I have the same problem, may anyone explain how can I extract Liquid Penetration Length (LPL) and vapour penetration length (VPL)? |
||
|
|
|
||
|
August 27, 2017, 13:23
|
|
#14 | ||
|
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128  |
Greetings to all, although I'm coming rather late to look into this.
@omid20110: I finally taken a look at the PM you sent me. In it you mentioned: Quote:
Quote:
The best I can do right now is wave my finger at this wiki page: https://openfoamwiki.net/index.php/Contrib_dieselFoam - and ask: Is this the type of case you are running? Either way, the best would be to let us know which tutorial case can be used as a reference to study it and so that I or anyone else can figure out how to use it. Best regards, Bruno
__________________
|
|||
|
|
|
|||
|
October 13, 2017, 22:38
|
|
#15 | |
|
Member
Sadegh Ebadi
Join Date: Apr 2015
Posts: 75
Rep Power: 11 |
Quote:
Thanks for your reply. what do you mean by this: "...which tutorial case can be used as a reference..."? |
||
|
|
|
||
|
March 22, 2018, 02:53
|
|
#16 | |
|
New Member
Suchit
Join Date: Jul 2015
Location: India
Posts: 5
Rep Power: 10 |
Hello everybody,
I have gone through this thread and many other threads on this forum. But still, I think the issue is open and there is no specific information about how to plot vapor penetration length. According to one paper which I read recently, its definition is Quote:
I could extract LPL and SMD from the log file. But what is the variable thats associated with VPL? or should I modify the solver so that it ill be written in the log file?. Till now, I have tried utilities like sample and probe. But I dont think they will give the answer. Any help in this will be surely useful. |
||
|
|
|
||
|
| Tags |
| diesel, liquid penetration length, spray simulation, sprayfoam |
| Thread Tools | Search this Thread |
| Display Modes | |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Extracting Liquid Penetration Length | ahcai007 | OpenFOAM Running, Solving & CFD | 15 | April 18, 2022 09:01 |
| Liquid penetration length and lift-off length | lunar1984 | STAR-CD | 3 | October 5, 2017 05:00 |
| sprayFoam crashes | lukasfischer | OpenFOAM Running, Solving & CFD | 3 | July 14, 2013 11:08 |
| ploting liquid penetration | mooyii | OpenFOAM | 2 | September 9, 2011 18:24 |
| Need some help for total liquid fraction | linus | FLUENT | 0 | December 19, 2006 03:29 |
3Likes
Linear Mode
