CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

expandableInterFoam

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

Like Tree1Likes
  • 1 Post By linch

Reply
 
LinkBack Thread Tools Display Modes
Old   April 27, 2012, 07:28
Default expandableInterFoam
  #1
Senior Member
 
Illya Shevchuk
Join Date: Aug 2009
Location: Darmstadt, Germany
Posts: 174
Rep Power: 8
linch is on a distinguished road
Hi guys,

I'm working on a modification of the interFoam solver that considers isobaric expansibility of both phases.

So my modified VOF and pressure correction equations look like:
\frac{\partial{\alpha_1}}{\partial{t}}+\nabla\cdot \left( \alpha_1 \vec{u}_1 \right) =Su_{\alpha_1}
\nabla \cdot \left( \left( \frac{1}{A_D} \right)_f \nabla_f p \right)=\nabla \cdot \vec{u}^* + Su_p
where Su_p=\nabla \cdot \vec{u} is velocity divergence computed from the continuity equation.

The problem with this formulation is, that it is not conservative and requires consistent and accurate calculation of the source terms Su_{\alpha_1} and Su_p. For now it works fine within in each phase, but I still have considerable mass imbalance in the interface region.

The second formulation I try to implement is a conservative one. VOF equation works well:
\frac{\partial{(\alpha_1 \rho_1)}}{\partial{t}}+\nabla\cdot \left( \alpha_1 \rho_1 \vec{u}_1 \right) =0

But I get pressure-velocity-oscillations if I solve following pressure correction equation:
\nabla \cdot \left( \left( \frac{1}{A_D} \right)_f \rho_f \nabla_f p \right)=\nabla \cdot \left( \rho \vec{u} \right)^* + Su_p
with Su_p=\frac{\partial \rho}{\partial t}

It looks like the classical chessboard problem (cartesian mesh), which I thought shouldn't be there because of the Rhie&Chow-like discretization of pressure. And that's actually my question to you: why does it happen and how could I prevent it?

Can anyone help me with a smart advice?

Best regards,
Ilya
Attached Images
File Type: jpg p_rgh.jpg (30.3 KB, 34 views)
Bernhard likes this.
linch is offline   Reply With Quote

Old   May 2, 2012, 09:30
Default
  #2
Senior Member
 
Illya Shevchuk
Join Date: Aug 2009
Location: Darmstadt, Germany
Posts: 174
Rep Power: 8
linch is on a distinguished road
Update: reducing the convergence criteria (i.e. from 1e-08 to 1e-12 for p_rgh) reduces the pressure oscillations dramatically, but still I'm interested to know where do the oscillations come from?
linch is offline   Reply With Quote

Reply

Tags
compressible, interfoam, pressure correction, vof

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



All times are GMT -4. The time now is 00:13.