CFDCode FLUTAN
Basics:
FLUTAN is a highly vectorized computer code for steady and transient 3d fluiddynamic and thermalhydraulic analyses in Cartesian or cylinder coordinates. It was developed in order to simulate single phase flows with small compressibility. The conservation equations for mass, momentum, energy, concentration of a contaminant, and turbulence quantities are solved. Material property libraries are available for the fluids water, air, sodium, lead, and leadbismuth.
Numerics:
The equations are discretized on a structured grid by a finite volume method. A staggered grid is used for the velocities. The discretization of the diffusive terms is performed by a central difference scheme. A first order upwind or one of two second order upwind methods (QUICK and LECUSSO) can be chosen for the convective terms. Several Poisson solvers are available for pressure calculation by means of the ICE technique, e.g. the CRESOR solver. A first order implicit Eulermethod is used for time discretization.
Newly developed and implemented numerical features are achieving a larger efficiency and an improved accuracy: One is an explicitly treated static local grid refinement method for the Cartesian grid. The other one is a new method to use body fitted coordinates in a code basing on a staggered grid: the conservation equations are transformed from the Cartesian to a general curvilinear system by keeping the physical Cartesian velocity components as dependent variables and by defining three Cartesian velocity components on every cell surface.


Calculated instantaneous temperature field on the surfaces of a sump cooling experiment. View from below on the fluid domain. The blue areas represent the cooler surface temperatures, the straggly red and yellow surface represents the temperature distribution on the heater plate. (Groetzbach et al. 2001) 
Physical models:
Several turbulence models are available in FLUTAN basing on transport equations for some turbulence quantities. The most important one is the Turbulence Model for Buoyant Flows (TMBF) which consists of a first order keps model in a lowReynolds number formulation and a second order fiveequations turbulent heat flux model. The TMBF contains a number of models which were developed for variable Prandtl number applications on the basis of our Direct Numerical Simulation results from the TURBITcode. In several benchmarks it turned out that the TMBF in its current development status is a powerful tool at least for forced and mixed convection even for liquid metal flows. The model development activities concentrate on achieving a wider application range of the more sophisticated turbulence models.
Applying a porous body concept with meshcell dependent porosities and directiondependent surface permeabilities allows to model effectively internal structures with scales below the grid widths and to mimic internal and external walls at cell surfaces.
Special thermal boundary conditions are available like a heat exchanger model, a 1d wall model, and a model for heat radiation from solid surfaces. A 3d heat conduction model for solid structures was developed for simulating their internal nonuniform transport of heat. The structure temperatures are discretized on a separate grid on which the heat conduction equation is solved.
Pre and Postprocessing:
These tools were undergoing a permanent extension and modification. By developing adequate interfaces, interactive geometry specification and meshing was possible; a TCL/TK based user interface allowed to handle the code system and to control interactively running jobs; several output interfaces allow to analyse the results by some graphics tools, e.g. it is possible to animate even 3d presentations of data by the AVS software.


Application fields: FLUTAN was carefully tested and verified. It successfully participated in several benchmarks, some with blind predictions for some experiments. It was already intensively used in the nuclear field to
 investigate the decay heat removal from pooltype LMFBRs,
 to analyse the cooling of PWR containment shells by natural convection and radiation,
 to interpret the singlephase longterm cooling of a core melt in an optional sump cooling concept of the EPR (see upper figure above),
 and to predict the flow and cooling behaviour in ADS targets and in the related MEGAPIE target module (see lower figure above).
References to the code
 Bottoni, M., Willerding, G., Advanced solution algorithms for transient multidimensional thermohydraulic flow problems in complex geometries with the program COMMIX2/KfK. Nucl. Engng & Design 100 (1987) pp.35165.
 Borgwaldt, H., CRESOR, A Robust Vectorized Poisson Solver Implemented in the COMMIX2(V) ThermalHydraulic Code. Int. Conf. On Supercomputing in Nuclear Applications (SNA 90). Mito City, Japan, March 1216 (1990) pp. 346351.
 Guenther, C., Fortgeschrittene UpwindDifferenzenVerfahren zur numerischen Loesung der KonvektionsDiffusionsgleichung. Habilitation, Universitaet Karlsruhe (1992) KfK 4697
 Willerding, G., Baumann, W., FLUTAN 2.0 Input Specifications. Forschungszentrum Karlsruhe (1996) FZKA 5712.
 Ammann, T., Entwicklung eines impliziten Verfahrens zur lokalen Verfeinerung eines strukturierten Gitters; Dissertation, Universitaet Karlsruhe (1997) FZKA 5864
 Carteciano L.N., Weinberg D., Mueller U., Development and Analysis of a Turbulence Model for Buoyant Flows (96kB). Proc. 4th World Conf. on Experimental Heat Transfer, Fluid Mechanics and Thermodynamics, Bruxelles, Belgium, June 26 (1997) 3 Pisa Edition ETS, pp. 133946
 Jin, X., Rechenverfahren zur Diskretisierung von Stroemungen in komplexer Geometrie mittels koerperangepasster Gitter. Dissertation, Universitaet Karlsruhe (2001) FZKA 6596
 Carteciano, L., Groetzbach, G., Validation of turbulence models for a free hot sodium jet with different buoyancy flow regimes using the computer code FLUTAN; Forschungszentrum Karlsruhe, FZKA 6600 (2003)
References to some applications:
 Weinberg, D.; Rust, K.; Hoffmann, H., Overview report of RAMONANEPTUN program on passive decay heat removal. Forschungszentrum Karlsruhe, FZKA5667 (1996)
 Baumann, W., Carteciano, L., Weinberg, D., Thermal propagation effects in a vertical turbulent flow behind a jet block  A benchmark exercise. J. of Hydraulics Research 35 (1997) pp. 843864
 Cheng, X., Mueller, U., Turbulent Natural Convection Coupled with Thermal Radiation in Large Vertical Channels with Asymmetric Heating. Int. J. Heat Mass Transfer 41 (1998) pp. 16811692
 Groetzbach, G., Carteciano, L.N., Dorr, B., Jin, X., Analysis of an LWR Sump Cooling Concept (684kB). 29th IAHR Congress, Special Seminar Advances in Industrial Hydraulics and Applications to Energy Production, Beijing, China, September 1721 (2001)
 Groetzbach, G., Carteciano, L. N., Dorr, B., Analysis of the integral heat transfer characteristics of the MEGAPIE target (1033kB). Jahrestagung Kerntechnik 2002, Kerntechnische Gesellschaft e.V., INFORUM Bonn (2002) pp. 559 – 563
For further references on the FLUTAN code and its applications consider the databank of FZK publications and search for the publications of the authors listed above.
Contact person:
G. Groetzbach, email: groetzbach at kit.edu