Computer Programs
NESC0877 K-FIX(3D).
K-FIX(3D), Transient 2 Phase Flow Hydrodynamic, X-Y-Z and Cylindrical Geometry, Eulerian Method
NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROBLEM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHOR, MATERIAL, CATEGORIES
[ top ]
[ top ]
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 |
|---|---|---|---|
| K-FIX(3D) | NESC0877/01 | Tested | 01-JUN-1982 |
Machines used:
| Package ID | Orig. computer | Test computer |
|---|---|---|
| NESC0877/01 | CDC 7600 | CDC 7600 |
[ top ]
3. DESCRIPTION OF PROBLEM OR FUNCTION
This package consists of two programs K-FIX(3D,FLX) which extend the transient, two-dimensional, two fluid program K-FIX (NESC Abstract 727) to perform three- dimensional calculations. The transient dynamics of three- dimensional, two-phase flow with interfacial exchange are calculated at all flow speeds. Each phase is described in terms of its own density, velocity, and temperature. The application is to flow in the annulus between two cylinders where the inner cylinder moves periodically perpendicular to its axis.
K-FIX(3D) is easily adaptable to a variety of two phase flow problems while K-FIX(3D,FLX) combines KFIX(3D), the three- dimensional version of the KFIX code, with the three-dimensional, elastic shell code FLX for application to a very specific class of problems. KFIX(3D,FLX) was developed specifically to calculate the coupled fluid structure dynamics of a light water reactor core support barrel under accident conditions. Motion may be induced by blowdown, prescribed displacement , or seismic action.
This package consists of two programs K-FIX(3D,FLX) which extend the transient, two-dimensional, two fluid program K-FIX (NESC Abstract 727) to perform three- dimensional calculations. The transient dynamics of three- dimensional, two-phase flow with interfacial exchange are calculated at all flow speeds. Each phase is described in terms of its own density, velocity, and temperature. The application is to flow in the annulus between two cylinders where the inner cylinder moves periodically perpendicular to its axis.
K-FIX(3D) is easily adaptable to a variety of two phase flow problems while K-FIX(3D,FLX) combines KFIX(3D), the three- dimensional version of the KFIX code, with the three-dimensional, elastic shell code FLX for application to a very specific class of problems. KFIX(3D,FLX) was developed specifically to calculate the coupled fluid structure dynamics of a light water reactor core support barrel under accident conditions. Motion may be induced by blowdown, prescribed displacement , or seismic action.
[ top ]
4. METHOD OF SOLUTION
In the K-FIX(3D), the six field equations for the two phases couple through mass, momentum, and energy exchange. The equations are solved using an Eulerian finite-difference technique that implicitly couples the rates of phase transitions, momentum, and energy exchange to determination of the pressure, density and velocity fields. The implicit solution is accomplished iteratively without linearizing the equations, thus eliminating the need for numerous derivative terms. With the three-dimensional K-FIX code calculations in Cartesian and cylindrical geometries can be performed. Obstacles built from the computing cells can be specified within the computing volume. In cylindrical geometry, calculations can be performed in the full 360 degrees or any angular segment. To enhance computing efficiency a new cell indexing scheme has been introduced, and computing time is reduced further by deletion of the viscous stress and heat conduction terms from the momentum and energy equations. FLX, which uses an explicit finite-difference solution algorithm to solve the shell equations, is explicitly coupled to the K-FIX(3D) fluid-dynamics program. A finite-difference numerical model is used for describing the dynamics of a core barrel.
In the K-FIX(3D), the six field equations for the two phases couple through mass, momentum, and energy exchange. The equations are solved using an Eulerian finite-difference technique that implicitly couples the rates of phase transitions, momentum, and energy exchange to determination of the pressure, density and velocity fields. The implicit solution is accomplished iteratively without linearizing the equations, thus eliminating the need for numerous derivative terms. With the three-dimensional K-FIX code calculations in Cartesian and cylindrical geometries can be performed. Obstacles built from the computing cells can be specified within the computing volume. In cylindrical geometry, calculations can be performed in the full 360 degrees or any angular segment. To enhance computing efficiency a new cell indexing scheme has been introduced, and computing time is reduced further by deletion of the viscous stress and heat conduction terms from the momentum and energy equations. FLX, which uses an explicit finite-difference solution algorithm to solve the shell equations, is explicitly coupled to the K-FIX(3D) fluid-dynamics program. A finite-difference numerical model is used for describing the dynamics of a core barrel.
[ top ]
[ top ]
[ top ]
[ top ]
[ top ]
NESC0877/01, included references:
- J.K. Dienes et al.:FLX: A Code for Coupled Fluid-Structure Analysis of Core Barrel
Dynamics. NUREG/CR-0959 LA-7927 R-4 (November 1979).
- W.C. Rivard and M.D. Torrey:
K-FIX: A Computer Program for Transient, Two-Dimensional,
Two-Fluid Flow. LA-NUREG-6623 NRC-4 (April 1977).
- W.C. Rivard and M.D. Torrey:
K-FIX: A Computer Program for Transient, Two Dimensional,
Two-Fluid Flow. LA-NUREG-6623, Suppl. II R4 (January 1979).
- Addendum to LA-NUREG-6623, Suppl. II
Corrections to the Listings of Three-Dimentional Code, Appendix
(3 pages)
- NESC Note 82-01 (October 9, 1981):
K-FIX(3D), NESC No. 877, Tape Description, Implementation,
Information, and Miscellaneous Subroutine Descriptions
(9 miscellaneous attachments to NESC Note, listed thereon)
[ top ]
11. MACHINE REQUIREMENTS
The K-FIX(3D) sample problem, excluding plotting and timing routines, took about 50000 (octal) words of a small core memory (SCM) and 65000 (octal) words of large core memory (LCM) storage. The K-FIX(3D, FLX) sample problem, excluding plotting and timing routines, took about 71300 (octal) words of SCM and 570000 (octal) words of LCM storage.
The K-FIX(3D) sample problem, excluding plotting and timing routines, took about 50000 (octal) words of a small core memory (SCM) and 65000 (octal) words of large core memory (LCM) storage. The K-FIX(3D, FLX) sample problem, excluding plotting and timing routines, took about 71300 (octal) words of SCM and 570000 (octal) words of LCM storage.
[ top ]
[ top ]
14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS
Memory
must be preset to zero when loading the programs and the parameter LCM=I must be added to the FTN card of K-FIX(3D,FLX) to allow proper indexing of very large arrays. Dependence on CDC system features and the local plotting and timing routines requires modification and incorporation of suitable replacement routines for implementation on different systems and in other computer environments.
Memory
must be preset to zero when loading the programs and the parameter LCM=I must be added to the FTN card of K-FIX(3D,FLX) to allow proper indexing of very large arrays. Dependence on CDC system features and the local plotting and timing routines requires modification and incorporation of suitable replacement routines for implementation on different systems and in other computer environments.
[ top ]
[ top ]
NESC0877/01
| File name | File description | Records |
|---|---|---|
| NESC0877_01.002 | K-FIX(3D) INFORMATION | 25 |
| NESC0877_01.003 | K-FIX(3D,FLX) SOURCE CARD IMAGES (FORTRAN-4) | 9826 |
| NESC0877_01.004 | K-FIX(3D,FLX) INPUT FOR SAMPLE CASE | 48 |
| NESC0877_01.005 | K-FIX(3D) SOURCE CARD IMAGES (FORTRAN-4) | 8342 |
| NESC0877_01.006 | K-FIX(3D) INPUT FOR SAMPLE CASE | 40 |
| NESC0877_01.007 | K-FIX(3D) SOURCE UPDATE (FORTRAN-4) | 3655 |
| NESC0877_01.008 | PERM3D SOURCE UPDATE (FORTRAN-4) | 115 |
| NESC0877_01.009 | ADDMAS SOURCE UPDATE (FORTRAN-4) | 81 |
| NESC0877_01.010 | ARRAYS SOURCE UPDATE (FORTRAN-4) | 129 |
| NESC0877_01.011 | K-FIX(3D) JCL | 29 |
| NESC0877_01.012 | DUMMY SUBROUTINES SOURCE (FORTRAN-4) | 45 |
| NESC0877_01.013 | K-FIX(3D,FLX) OUTPUT OF SAMPLE CASE | 1579 |
| NESC0877_01.014 | K-FIX(3D) OUTPUT OF SAMPLE CASE | 912 |
Keywords: LWR reactors, interfaces, pressure, reactor safety, shells, simulation, three-dimensional, two-phase flow, velocity.