CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Cavitating in Foam14

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 26, 2007, 08:06
Default Does anybody know about cavita
  #1
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 7
zjucfd is on a distinguished road
Does anybody know about cavitating solver in Foam1.4? I first run the tutorials case of Nozzle2D, it works well. Then I attempted to try the flow around a hyrofoil, but got bad results. The pressure vary a lot, beyond my imagination, even below zero at some spots. The only reason I can figure out is that the solver introduced the barotropic model, which allows the liquid to be compressible. Could anybody give me a hint that what can I do. Thanks a lot!
zjucfd is offline   Reply With Quote

Old   June 26, 2007, 08:53
Default Hey, what is the pressure r
  #2
Senior Member
 
Markus Hartinger
Join Date: Mar 2009
Posts: 102
Rep Power: 7
hartinger is on a distinguished road
Hey,

what is the pressure range you're expecting from experiments or analytical solutions? What are you getting?

Does the flow field look reasonable? experiments?

Mesh good enough? Timestep appropriate? Is the solution smooth?

cavitatingFoam can produce and handle negative pressure. Sure, it is not physical?

regards
markus
hartinger is offline   Reply With Quote

Old   June 26, 2007, 11:02
Default Hi, I let the outlet pressu
  #3
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 7
zjucfd is on a distinguished road
Hi,

I let the outlet pressure fixed at 1.e5pa, and the inlet velocity 0.6m/s. In order to run a no-cavitation case as a test, the saturation vapour pressure is set to a very low value, ensuring that no cavitation occurs.

But the obtained flow field is very odd. The maximum value of velocity reaches about 50m/s, and the pressure reaches 5pa. so, it looks unreasonable at all.

What's going wrong with my calculation? It's due to the boundary conditions? of which I'm not very sure. Can you recommend me some materials about this cavitation model in Foam1.4?

thanks,
best regards!
Scofield
zjucfd is offline   Reply With Quote

Old   June 26, 2007, 12:00
Default to check whether your boundary
  #4
Senior Member
 
Markus Hartinger
Join Date: Mar 2009
Posts: 102
Rep Power: 7
hartinger is on a distinguished road
to check whether your boundary conditions are correct, apply them to the nozzle2D case. with such a low speed you shoudn't have any problems
hartinger is offline   Reply With Quote

Old   June 26, 2007, 12:20
Default Dear Hartinger, I have sent
  #5
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 7
zjucfd is on a distinguished road
Dear Hartinger,

I have sent you an e_mail with my case files enclosed. It's very kind of you to give me some advice on it. I think my case is different from nozzle2D, which is a pressure driven flow.
zjucfd is offline   Reply With Quote

Old   June 27, 2007, 12:09
Default Does anybody can take a look a
  #6
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 7
zjucfd is on a distinguished road
Does anybody can take a look at my case files? I suppose that the inlet pressure is fixed at 1e5pa, and velocity about 1m/s. Thanks a lot!

zjucfd is offline   Reply With Quote

Old   June 27, 2007, 21:54
Default sorry, the mesh file is too bi
  #7
New Member
 
Connie Scofield
Join Date: Mar 2009
Posts: 12
Rep Power: 7
zjucfd is on a distinguished road
sorry, the mesh file is too big to upload.
0.tar.gz
constant.tar.gz
system.tar.gz
zjucfd is offline   Reply With Quote

Old   July 6, 2007, 04:20
Default I tried to perform cavitatingF
  #8
Member
 
Ruben I. Mukhamadeev
Join Date: Mar 2009
Location: Obninsk, Kaluga reg., Russian Federation
Posts: 69
Rep Power: 7
benru is on a distinguished road
I tried to perform cavitatingFoam from tutorial (nozzle2D) - and calcs was stopped with errors :

Exec : cavitatingFoam . nozzle2D
Date : Jul 06 2007
Time : 12:18:28
Host : r117-1
PID : 6746
Root : /home/rmukhamadeev/OpenFOAM/rmukhamadeev-1.4/run/tutorials/cavitatingFoam
Case : nozzle2D
Nprocs : 1
Create time

Create mesh for time = 0

Reading thermodynamicProperties

Reading transportProperties

Reading field p

Creating compressibilityModel

Selecting compressibility model linear
Reading field U

Reading/calculating face flux field phiv

Reading/calculating face flux field phi

phiv Courant Number mean: 0 max: 0 acoustic max: 0.298593

Starting time loop

phiv Courant Number mean: 0 max: 0 acoustic max: 0.298593
deltaT = 1.19999e-10
Time = 1.19999e-10

smoothSolver: Solving for rho, Initial residual = 0.50711, Final residual = 0.254106, No Iterations 1000
max-min rho: 832.754 828.236
max-min gamma: 0 0
#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xb7f6d420]
#3 Foam::polyMesh::calcDirections() const
#4 Foam::polyMesh::directions() const
#5 Foam::fvMatrix<foam::vector<double> >::solve(Foam::Istream&)
#6 Foam::lduMatrix::solverPerformance Foam::solve<foam::vector<double> >(Foam::tmp<foam::fvmatrix<foam::vector<double> > > const&)
#7 main
#8 __libc_start_main
#9 __gxx_personality_v0 at ../sysdeps/i386/elf/start.S:122

Could anybody point a cause ? Need help !

Ruben
benru is offline   Reply With Quote

Old   July 6, 2007, 05:59
Default Thank you, will try Ruben
  #9
Member
 
Ruben I. Mukhamadeev
Join Date: Mar 2009
Location: Obninsk, Kaluga reg., Russian Federation
Posts: 69
Rep Power: 7
benru is on a distinguished road
Thank you, will try

Ruben
benru is offline   Reply With Quote

Old   August 9, 2007, 02:30
Default Hello all, I am trying to u
  #10
New Member
 
Ido Silverman
Join Date: Mar 2009
Posts: 13
Rep Power: 7
idosil is on a distinguished road
Hello all,

I am trying to understand the cavitatingFoam solver and have a few questions:

1. in UEqn - the viscosity term appears twice, once implicit with the laplacian operator and a second time explicit as div of (muf grad U). Should not each term be multiplied by 0.5?

2. in pEqn - I understand that both p.boundaryField().updateCoeffs() and U.correctBoundaryConditions() fix recalculate boundary condition for these fields. what is the difference?

3. in pEqn - rUA and HbyA are recalculated at each inner iteration but the UEqn equation matrixs were build only once outside of this loop. does the matrixs implicitly change when the boundary or fields are recalculated?

Thanks, Ido
idosil is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to set cavitation number for cavitating foil Karim FLUENT 2 June 28, 2009 21:05
cavitating flow John Main CFD Forum 2 May 4, 2008 17:26
To modify the interFoam solver for cavitating zjucfd OpenFOAM Running, Solving & CFD 5 August 26, 2007 08:03
Cavitating Flows CFdtoy Main CFD Forum 0 July 20, 2004 16:51
2D cavitating hydrofoil Therus CFX 1 February 26, 2004 11:14


All times are GMT -4. The time now is 10:02.