CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM

MULES advection scheme for compressible flows?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 2, 2012, 10:35
Default MULES advection scheme for compressible flows?
  #1
New Member
 
Pharg Mandadapu
Join Date: Jul 2011
Posts: 16
Rep Power: 5
pharg_yrartibra is on a distinguished road
The MULES advection scheme is described in the report "Drop impact onto a liquid layer of finite thickness: dynamics of the cavity evolution", as well as in the report "A Coupled Pressure Based Solution Algorithm Based on the Volume-Of-Fluid Approach for Two or More Immiscible Fluids" where it is improved to work better for three or more phases.

In the first report, the flow is considered to be incompressible. I don't know about the second report because I can't access it right now (I'm getting a 404 message) but I guess that it also treats the flow as incompressible. But what about if you have a compressible flow? Can you modify the MULES scheme to cope with compressible flows as well? Is it described somewhere how you can do that and how the modified version of MULES for compressible flows would work? Thanks in advance.

Last edited by pharg_yrartibra; August 22, 2012 at 20:41.
pharg_yrartibra is offline   Reply With Quote

Old   August 9, 2012, 17:35
Default
  #2
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
Well, until I could read in both papers, no description of MULES is given far than naming it. MULES is based in Flux Corrected Transport (Zalezak's work) so that you could search about FCT in FVM for compressible flows and then do some changes in FOAM's code.

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   August 9, 2012, 19:27
Default
  #3
New Member
 
Pharg Mandadapu
Join Date: Jul 2011
Posts: 16
Rep Power: 5
pharg_yrartibra is on a distinguished road
Thanks, I found the reports Flux-corrected transport, I: SHASTA, a fluid transport algorithm that works by Jay P. Boris and David L. Book and Fully multidimensional flux-corrected transport algorithms for fluids by Steven T. Zalesak.

Is there any difference between the method used in Zalesak's paper and the method used in interFoam (which I guess is MULES)? Also, do you know if these methods conserves mass?
pharg_yrartibra is offline   Reply With Quote

Old   August 9, 2012, 22:31
Default
  #4
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
The principal conceptual difference between Zalesak's FCT and MULES we've found is the calculations of lambdas, their are the factors which control the amount of anti-diffusion used in the method. In Zalesak's work these are constant in each time-step while in MULES the calculation is iterative, but I don't know which is the loop finishing criterion. As is coded I think the method is less diffusive than FCT, but at the same time it could be more unstable. The other difference is the implementation, even though FCT is theoretically extensible to multi-dimension, that's not show in the paper. MULES is fully multidimensional, working by faces, the code is very tricky to understand. Some parts are not clear for me, I asked in this forum without answer, for example:

MULES::explicitSolve dimensions inconsistence

maybe it is proper to ask directly in FOAM Mantis.

There is a paper from Rudman where he shows the use of FCT in VOF problems, I think it could be treated as a "bridge" between Zalesak's work and the modern treatment found in FOAM MULES

Volume-Tracking Methods for Interfacial Flow Calculations

he compares it with SLIC and other reconstruction methods typically used in VOF, there aren't much differences. In fact in the last OpenFOAM workshop there was a lot of activity around VOF working with CICSAM (which had been dropped off from vanilla FOAM) and other schemes.

Respect to the conservativeness, MULES conserves volume, which in incompressible solvers implies mass conservation, but not in compressible ones. The other objective is to preserve boundedness. Maybe you can ask to Kent Wardle, he is working in interFoam related solver and dealing with the conservation:

multiphaseInterFoam

Hope this helps, and please post your findings to share the knowledge.

Regards
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Reply

Tags
compressible flow, mules

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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
Total variation diminishing scheme for advection or advection dispersion equation cooljd Main CFD Forum 7 July 27, 2012 09:32
Advection Scheme on an C-Grid Mr. Matti Main CFD Forum 0 April 14, 2010 13:19
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 12:16
TVD scheme applied to axsymetric nozzle flows Zhou Hua Main CFD Forum 0 May 7, 1999 01:10


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