Mark H. Anderson

A computational parameter study for the three-dimensional shock-bubble interaction

The morphology and time-dependent integral properties of the multifluid compressible flow resulting from the shock–bubble interaction in a gas environment are investigated using a series of three-dimensional multifluid-Eulerian simulations. The bubble consists of a spherical gas volume of radius 2.54 cm (128 grid points), which is accelerated by a planar shock wave. Fourteen scenarios are considered: four gas pairings, including Atwood numbers −0.8 A M ≤ 5.0. The data are queried at closely spaced time intervals to obtain the time-dependent volumetric compression, mean bubble fluid velocity, circulation and extent of mixing in the shocked-bubble flow. Scaling arguments based on various properties computed from one-dimensional gasdynamics are found to collapse the trends in these quantities successfully for fixed A . However, complex changes in the shock-wave refraction pattern introduce effects that do not scale across differing gas pairings, and for some scenarios with A > 0.2, three-dimensional (non-axisymmetric) effects become particularly significant in the total enstrophy at late times. A new model for the total velocity circulation is proposed, also based on properties derived from one-dimensional gasdynamics, which compares favourably with circulation data obtained from calculations, relative to existing models. The action of nonlinear-acoustic effects and primary and secondary vorticity production is depicted in sequenced visualizations of the density and vorticity fields, which indicate the significance of both secondary vorticity generation and turbulent effects, particularly for M > 2 and A > 0.2. Movies are available with the online version of the paper.

A diffusion layer model for steam condensation within the AP600 containment

Abstract Steam condensation plays a key role in removing heat from the atmosphere of the Westinghouse AP600 containment in case of a postulated accident. A model of steam condensation on containment surfaces under anticipated accident conditions is presented and validated against an extensive and sound database. Based on the diffusion layer theory and on the use of the heat/mass transfer analogy, one can deal with large temperature gradients across the gaseous boundary layer under high mass flux circumstances. The thermal resistance of the condensate film, as well as its wavy structure, have also been considered in this model. As compared to Anderson et al. (1998) (Experimental analysis of heat transfer within the AP600 containment under postulated accident conditions. Nucl. Eng. Des. (submitted)) experimental database, an average error lower than 15%, within the experimental confidence range, has demonstrated its remarkable accuracy. In particular, the model has shown a good response to the influence of primary variables in steam condensation (i.e. subcooling, noncondensable concentration and pressure), providing a mechanistic explanation for effects such as the presence of light noncondensable gas (i.e. helium as a simulant for hydrogen) in the gaseous mixture. In addition, the model has been contrasted against correlations used in safety analysis (i.e. Uchida, Tagami, Kataoka, etc.) and occasionally to Dehbi’s database. This cross-comparison has pointed out several shortcomings in the use of these correlations and has extended the model validation to other databases.

Experimental analysis of heat transfer within the AP600 containment under postulated accident conditions

Abstract The new AP600 reactor designed by Westinghouse uses a passive safety system relying on heat removal by condensation to keep the containment within the design limits of pressure and temperature. Even though some research has been done so far in this regard, there are some uncertainties concerning the behavior of the system under postulated accident conditions. In this paper, steam condensation onto the internal surfaces of the AP600 containment walls has been investigated in two scaled vessels with similar aspect ratios to the actual AP600. The heat transfer degradation in the presence of noncondensable gas has been analyzed for different noncondensable mixtures of air and helium (hydrogen simulant). Molar fractions of noncondensables/steam ranged from (0.4–4.0) and helium concentrations in the noncondensable mixture were 0–50% by volume. In addition, the effects of the bulk temperatures, the mass fraction of noncondensable/steam, the cold wall surface temperature, the pressure, noncondensable composition, and the inclination of the condensing surface were studied. It was found that the heat transfer coefficients ranged from 50 to 800 J s −1 K −1 m −2 with the highest for high wall temperatures at high pressure and low noncondensable molar fractions. The effect of a light gas (helium) in the noncondensable mixture were found to be negligible for concentrations less than approximately 35 molar percent but could result in stratification at higher concentrations. The complete study gives a large and relatively complete data base on condensation within a scaled AP600 containment structure, providing an invaluable set of data against which to validate models. In addition, specific areas requiring further investigation are summarized.

Heat Transfer and Fluid Flow Characteristics in Supercritical Pressure Water

A series of integral heat transfer measurements in a square annular flow passage was performed for bulk water temperatures of 175―400°C with upward mass velocities of 300 kg/m 2 s and 1000 kg/m 2 s and heat fluxes of 0, 200 kW/m 2 , and 440 kW/m 2 , all at a pressure of 25 MPa. Mean and turbulent velocities measured with a two-component laser Doppler velocimetry system along with simulations using the computational fluid dynamics (CFD) code FLUENT were used to explain the deterioration and enhancement of heat transfer in supercritical pressure water. At low mass velocities, the integral heat transfer measurements exhibited large localized wall temperature spikes that could not be accurately predicted with Nusselt correlations. Detailed mean and turbulent velocities along with FLUENT simulations show that buoyancy effects cause a significant reduction in turbulent quantities at a radial location similar to what is the law of the wall region for isothermal flow. At bulk temperatures near the pseudocritical temperature, high mass velocity integral heat transfer measurements exhibited an enhanced heat transfer with a magnitude dependent on the applied heat flux. Measured mean and turbulent velocities showed no noticeable changes under these conditions. FLUENT simulations show that the integrated effects of specific heat can be used to explain the observed effects. The experimentally measured heat transfer and local velocity data also serve as a database to compare existing CFD models, such as Reynolds-averaged Navier-Stokes (RANS) equations and possibly even large Eddy simulations (LES) and direct numerical simulations (DNS). Ultimately, these measurements will aid in the development of models that can accurately predict heat transfer to supercritical pressure water.

Shock-bubble interactions : Features of divergent shock-refraction geometry observed in experiments and simulations

The interaction of a planar shock wave with a spherical bubble in divergent shock-refraction geometry is studied here using shock tube experiments and numerical simulations. The particular case of a helium bubble in ambient air or nitrogen (A≈−0.8) is considered, for 1.4<M<3.0. Experimental planar laser diagnostics and three-dimensional multifluid Eulerian simulations clearly resolve features arising as a consequence of divergent shock refraction, including the formation of a long-lived primary vortex ring, as well as counter-rotating secondary and tertiary upstream vortex rings that appear at late times for M⩾2. Remarkable correspondence between experimental and numerical results is observed, which improves with increasing M, and three-dimensional effects are found to be relatively insignificant. Shocked-bubble velocities, length scales, and circulations extracted from simulations and experiments are used successfully to evaluate the usefulness of various analytical models, and characteristic dimensionle...

Experimental and numerical investigation of shock-induced distortion of a spherical gas inhomogeneity

Results are presented from a series of experiments and simulations, studying the interaction of a planar shock wave (2.0≤M≤5.0) with a discrete gas inhomogeneity. Experiments and computations confirm that the phenomenology of shock–bubble interactions is fundamentally altered by the changes in the Atwood number. In the case of a low Atwood number, the late time flow field is dominated by coherent vortical structures, whereas in the case of a high Atwood number, the shocked bubble is effectively reduced to a small core of compressed fluid, which trails behind a plume-like structure indicative of a well-developed mixing region.

Heat transfer in a supercritical fluid : Classification of heat transfer regimes

Because of the dramatic variation of physical properties with a modest change of temperature, no existing engineering correlation or models can accurately predict heat transfer of supercritical fluids. This paper seeks to classify the conditions where the existing models are applicable and to better understand these local heat transfer mechanisms. The first objective is the focus of this paper. FLUENT was employed to compute the wall temperatures for various heat flux and mass flux conditions and to be compared with experimental data. Because the model was developed for a wide range of flow conditions, it was necessary to make certain assumptions. The simulations showed a good agreement with high mass flux conditions, where buoyancy effects could be neglected. The FLUENT model, however, had difficulty predicting the localized low heat transfer rates seen in the combined condition of high heat flux and low mass flux. A new generalized parameter, dependent on the heat and mass flux, was developed to classify under which conditions this FLUENT standard model was applicable. This global Froude number can be used as the parameter to predict under which conditions the buoyancy effect will be dominant and lower heat transfer rates will occur.

CORROSION OF 316L STAINLESS STEEL ALLOY AND HASTELLOY-N SUPERALLOY IN MOLTEN EUTECTIC LiF-NaF-KF SALT AND INTERACTION WITH GRAPHITE

Abstract Molten FLiNaK salt [46.5%LiF-11.5%NaF-42%KF (mol%)] has been proposed for use as a secondary reactor coolant and medium for transfer of high-temperature process heat from nuclear reactors to chemical plants. Two alloys—Hastelloy-N superalloy (Hastelloy-N) and Type 316L stainless steel alloy (316L steel)—were exposed to molten FLiNaK salt in a 316L steel crucible under argon cover gas at 850·C for 1000 h. Graphite was also introduced into the test with the goal of studying the corrosion behavior of relevant reactor material combinations. The results show that corrosion of 316L steel occurred primarily through surface depletion of Cr. Contrarily, Hastelloy-N experienced weight gain due to the electrochemical plating of corrosion products, Fe and Cr, derived from the 316L steel crucible. The graphite sample enhanced the corrosion of the 316L steel sample and crucible, which induced the formation of (Cr,Fe)7C3 and (Mo,Cr,Fe)2C carbides on the surface of graphite. These carbide formations were attributed to the nonelectric transfer between 316L steel and graphite. Besides reducing the availability of chromium to plate, the presence of graphite did not change the basic corrosion of the 316L steel and plating process of Hastelloy-N.

In Situ Measurements of Spectral Emissivity of Materials for Very High Temperature Reactors

Abstract An experimental facility for in situ measurements of high-temperature spectral emissivity of materials in environments of interest to the gas-cooled very high temperature reactor (VHTR) has been developed. The facility is capable of measuring emissivities of seven materials in a single experiment, thereby enhancing the accuracy in measurements due to even minor systemic variations in temperatures and environments. The system consists of a cylindrical silicon carbide (SiC) block with seven sample cavities and a deep blackbody cavity, a detailed optical system, and a Fourier transform infrared spectrometer. The reliability of the facility has been confirmed by comparing measured spectral emissivities of SiC, boron nitride, and alumina (Al2O3) at 600°C against those reported in literature. The spectral emissivities of two candidate alloys for VHTR, INCONEL® alloy 617 (INCONEL is a registered trademark of the Special Metals Corporation group of companies) and SA508 steel, in air environment at 700°C were measured.

Impact of Corrosion Test Container Material in Molten Fluorides

The effects of crucible material choice on alloy corrosion rates in immersion tests in molten LiF–NaF–KF (46.5–11.5-42 mol. %) salt held at 850 °C for 500 hrs are described. Four crucible materials were studied. Molten salt exposures of Incoloy-800H in graphite, Ni, Incoloy-800H, and pyrolytic boron nitride (PyBN) crucibles all led to weight-loss in the Incoloy-800H coupons. Alloy weight loss was ~30 times higher in the graphite and Ni crucibles in comparison to the Incoloy-800H and PyBN crucibles. It is hypothesized galvanic coupling between the alloy coupons and crucible materials contributed to the higher corrosion rates. Alloy salt immersion in graphite and Ni crucibles had similar weight-loss hypothesized to occur due to the rate limiting out diffusion of Cr in the alloys to the surface where it reacts with and dissolves into the molten salt, followed by the reduction of Cr from solution at the molten salt and graphite/Ni interfaces. As a result, both the graphite and the Ni crucibles provided sinks for the Cr, in the formation of a Ni–Cr alloy in the case of the Ni crucible, and Cr carbide in the case of the graphite crucible.