Computer Programs
NESC0832 SOLA-DF.
last modified: 28-OCT-1991 | catalog | categories | new | search |

NESC0832 SOLA-DF.

SOLA-DF, Time-Dependent 2-D 2 Phase Flow, Eulerian Method with Various Boundary Conditions

top ]
1. NAME OR DESIGNATION OF PROGRAM:  SOLA-DF.
top ]
2. COMPUTERS

To submit a request, click below on the link of the version you wish to order. Rules for end-users are available here.

Program name Package id Status Status date
SOLA-DF NESC0832/01 Tested 01-OCT-1980
SOLA-DF NESC0832/02 Tested 28-OCT-1991

Machines used:

Package ID Orig. computer Test computer
NESC0832/01 CDC 7600 CDC 7600
NESC0832/02 DEC VAX 8810
top ]
3. DESCRIPTION OF PROBLEM OR FUNCTION

SOLA-DF is a numerical solution algorithm for gas-liquid mixture dynamics in two space dimensions and time. The two-phase system is described by a set of mixture equations plus a relation describing the relative flow of one phase  with respect to the other. The algorithm contains models to represent the interphase exchange rates of mass, momentum, and energy for water system mixtures.
top ]
4. METHOD OF SOLUTION

SOLA-DF uses a finite-difference technique, which includes the Implicit Continuous-fluid Eulerian (ICE) method as an optional feature, to solve the drift-flux equations of two phase flow in two-dimensional plane or cylindrical coordinates.
A wide variety of boundary conditions are possible, including free-slip, no-slip, specified inflow-outflow, continuative outflow,  periodic, and constant pressure boundaries. Implicit pressure gradients are used to eliminate the Courant condition as a stability requirements, making the code useful at all Mach numbers.
top ]
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

Simple equation of state, relative velocity model and phase transition model are assumed. When large phase transition rates are present, because of the decoupled calculation of phase change and the pressure iteration, care must be exercised to make sure the solution is independent of the time-step.
top ]
6. TYPICAL RUNNING TIME

The "grind" time, or time to advance the solution in one cell by one time-step, is typically 0.5 to 1.0 ms on a CDC7600. NESC executed the sample problem in less than 8 CP minutes on a CDC7600.
top ]
7. UNUSUAL FEATURES OF THE PROGRAM

SOLA-DF is designed to be easy to understand and modify so it can serve as a foundation for code development.
top ]
8. RELATED AND AUXILIARY PROGRAMS

The solution algorithm has evolved  from earlier SOLA series programs. The original SOLA was designed for problems involving a single, incompressible fluid in a fixed region. SOLA-SURF (NESC Abstract 651) is an extension which allows for the inclusion of free surfaces. SOLA-ICE (NESC Abstract 723) was designed to handle single-component, compressible fluids.
SOLA-DF  is a direct descendant of SOLA-ICE.
A variation of SOLA-DF is contained in SOLA-LOOP (NESC Abstrac 859), a network code, which handles systems of one-dimensional components  coupled through junctions. K-FIX (NESC Abstract 727) calculates the  transient dynamics of two-dimensional, two-phase flow with inter- facial exchange.
top ]
9. STATUS
Package ID Status date Status
NESC0832/01 01-OCT-1980 Tested at NEADB
NESC0832/02 28-OCT-1991 Screened
top ]
10. REFERENCES

- C.W. Hirt and N.C. Romero,
  Application of a Drift-Flux Model to Flashing in Straight Pipes,
  LA-6005-MS, July 1975.
NESC0832/01, included references:
- C.W. Hirt N.C. Romero, M.D. Torrey, and J.R. Travis:
  SOLA-DF: A Solution Algorithm for Nonequilibrium Two-Phase Flow,
  NUREG/CR-0690 (LA-7725-MS), (June 1979).
- R.M. Frank:
  J506 Advance Subroutine,
  J523 Draw Frame Subroutine,
  J563 Contour Plotting with Film Output Subroutine,
  LASL Program Library Write-up (January 1979).
- S.T. Bennion:
  J517 Draw Vector Subroutine,
  J526 Convert Coordinate Subroutine,
  J528 Define Graph Area Subroutine,
  LASL Program Library Write-up (January 1979).
- V. Gardiner:
  J508 Standardized Plot Subroutine,
  LASL Program Library Write-up (January 1979).
- SOLA-DF NESC No. 832.7600, Errata to LA-7725-MS Input Description,
  NESC Note 80-06 (July 9, 1979).
NESC0832/02, included references:
- C.W. Hirt N.C. Romero, M.D. Torrey, and J.R. Travis:
  SOLA-DF: A Solution Algorithm for Nonequilibrium Two-Phase Flow,
  NUREG/CR-0690 (LA-7725-MS), (June 1979).
- R.M. Frank:
  J506 Advance Subroutine,
  J523 Draw Frame Subroutine,
  J563 Contour Plotting with Film Output Subroutine,
  LASL Program Library Write-up (January 1979).
- S.T. Bennion:
  J517 Draw Vector Subroutine,
  J526 Convert Coordinate Subroutine,
  J528 Define Graph Area Subroutine,
  LASL Program Library Write-up (January 1979).
- V. Gardiner:
  J508 Standardized Plot Subroutine,
  LASL Program Library Write-up (January 1979).
- SOLA-DF NESC No. 832.7600, Errata to LA-7725-MS Input Description,
  NESC Note 80-06 (July 9, 1979).
top ]
11. MACHINE REQUIREMENTS

124000 (octal) words of memory are required for execution.
top ]
12. PROGRAMMING LANGUAGE(S) USED
Package ID Computer language
NESC0832/01 FORTRAN-IV
NESC0832/02 FORTRAN-IV
top ]
13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:  SCOPE and LTSS.
top ]
14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

SOLA-DF calls a number of LASL film plotting subroutines - ADV, CONTRJB, CONVRT, DGA, DRV, FRAME, LINCNT, and SPLOT. NESC supplied DUMMY routines for completeness. To duplicate the results in the report the ROUND option should be specified for the compiler. If ROUND is not specified, the results are printed at different time intervals.
top ]
15. NAME AND ESTABLISHMENT OF AUTHOR

            C.W. Hirth, N.C. Romero, M.D. Torrey,
            and J.R. Travis *
            Theoretical Division
            Los Alamos Scientific Laboratory
            P.O. Box 1663
            Los Alamos, New Mexico 87545
* Contact
top ]
16. MATERIAL AVAILABLE
NESC0832/01
File name File description Records
NESC0832_01.001 INFORMATION 5
NESC0832_01.002 SOURCE (F4,EBCDIC) 1401
NESC0832_01.003 INPUT DATA FOR SAMPLE PROBLEM 20
NESC0832_01.004 PRINTED OUTPUT OF SAMPLE PROBLEM 5908
NESC0832/02
File name File description Records
NESC0832_02.001 Information file 32
NESC0832_02.002 SOLA-DF source program (FORTRAN) 1422
NESC0832_02.003 SOLA-DF test case input data 20
NESC0832_02.004 SOLA-DF test case printed output 5908
top ]
17. CATEGORIES
  • H. Heat Transfer and Fluid Flow

Keywords: finite difference method, fluid flow, ice method, phase transformations, two-dimensional, two-phase flow.