Assessment of CTF boiling transition and critical heat flux modeling capabilities using the OECD/NRC BFBT and PSBT benchmark databases

The need to refine models for best-estimate calculations, based on good-quality experimental data, has been expressed in many recent meetings in the field of nuclear applications. The modeling needs arising in this respect should not be limited to the currently available macroscopic methods but should be extended to next-generation analysis techniques that focus on more microscopic processes. One of the most valuable databases identified for the thermalhydraulics modeling was developed by the Nuclear Power Engineering Corporation (NUPEC), Japan. From 1987 to 1995, NUPEC performed steady-state and transient critical power and departure from nucleate boiling (DNB) test series based on the equivalent full-size mock-ups. Considering the reliability not only of the measured data, but also other relevant parameters such as the system pressure, inlet sub-cooling and rod surface temperature, these test series supplied the first substantial database for the development of truly mechanistic and consistent models for boiling transition and critical heat flux. Over the last few years the Pennsylvania State University (PSU) under the sponsorship of the U.S. Nuclear Regulatory Commission (NRC) has prepared, organized, conducted and summarized the OECD/NRC Full-size Fine-mesh Bundle Tests (BFBT) Benchmark. The international benchmark activities have been conducted in cooperation with the Nuclear Energy Agency/Organization for Economic Co-operation and Development (NEA/OECD) and Japan Nuclear Energy Safety (JNES) organization, Japan. Consequently, the JNES has made available the Boiling Water Reactor (BWR) NUPEC database for the purposes of the benchmark. Based on the success of the OECD/NRC BFBT benchmark the JNES has decided to release also the data based on the NUPEC Pressurized Water Reactor (PWR) subchannel and bundle tests for another follow-up international benchmark entitled OECD/NRC PWR Subchannel and Bundle Tests (PSBT) benchmark. This paper presents an application of the joint Penn State University/Technical University of Madrid (UPM) version of the well-known subchannel code COBRA-TF, namely CTF, to the critical power and departure from nucleate boiling (DNB) exercises of the OECD/NRC BFBT and PSBT benchmarks

Nonlocal electron heat flux revisited

A known nonlocal model of electron heat flux, applying for (scale length/thermal ion-electron mean-free path) of order Z)1/2(e*/T)312, ionization number Z, large, and e*~ 6.5 T (the energy of electrons carrying most of the flux), is reconsidered. The large e*/T ratio simplifies the complete formalism. A simple flux formula, exact for both smooth and steep profiles, is given. Thermoelectric effects and other models are discussed.

Self-consistent, nonlocal electron heat flux at arbitrary ion charge number

A single, nonlocal expression for the electron heat flux, which closely reproduces known results at high and low ion charge number 2, and “exact” results for the local limit at all 2, is derived by solving the kinetic equation in a narrow, tail-energy range. The solution involves asymptotic expansions of Bessel functions of large argument, and (Z-dependent)order above or below it, corresponding to the possible parabolic or hyperbolic character of the kinetic equation; velocity space diffusion in self-scattering is treated similarly to isotropic thermalization of tail energies in large Z analyses. The scale length H characterizing nonlocal effects varies with Z, suggesting an equal dependence of any ad hoc flux limiter. The model is valid for all H above the mean-free path for thermal electrons.

Self-similar expansion of laser plasmas with nonlocal heat flux

A previous hydrodynamic model of the expansion of a laser-produced plasma, using classical (Spitzer) heat flux, is reconsidered with a nonlocal heat flux model. The nonlocal law is shown to be valid beyond the range of validity of the classical law, breaking down ultimately, however, in agreement with recent predictions.

Nonlocal electron heat-flux

Electron thermal conduction in a not quite collisional unmagnetlzed plasma is analysed. The failure of classical results for temperature scale-length up to 100 times larger than thermal mean-free-path for electron scattering, and large ion-charge number Z , is discussed. Recent results from a nonlocal model of conduction at large Z are reviewed. Closed form expressions for Braginskii's coefficients a ,/3 , y for Z =0(1) are derived. An extension of the nonlocal model for Z =0(1) is discussed.

Assessment of CTF boiling transition and critical heat flux modeling capabilities using the OECD/NRC BFBT and PSBT benchmark databases

Over the last few years, the Pennsylvania State University (PSU) under the sponsorship of the US Nuclear Regulatory Commission (NRC) has prepared, organized, conducted, and summarized two international benchmarks based on the NUPEC data—the OECD/NRC Full-Size Fine-Mesh Bundle Test (BFBT) Benchmark and the OECD/NRC PWR Sub-Channel and Bundle Test (PSBT) Benchmark. The benchmarks’ activities have been conducted in cooperation with the Nuclear Energy Agency/Organization for Economic Co-operation and Development (NEA/OECD) and the Japan Nuclear Energy Safety (JNES) Organization. This paper presents an application of the joint Penn State University/Technical University of Madrid (UPM) version of the well-known sub-channel code COBRA-TF (Coolant Boiling in Rod Array-Two Fluid), namely, CTF, to the steady state critical power and departure from nucleate boiling (DNB) exercises of the OECD/NRC BFBT and PSBT benchmarks. The goal is two-fold: firstly, to assess these models and to examine their strengths and weaknesses; and secondly, to identify the areas for improvement.

Ultra high temperature latent heat energy storage and thermophotovoltaic energy conversion

A conceptual energy storage system design that utilizes ultra high temperature phase change materials is presented. In this system, the energy is stored in the form of latent heat and converted to electricity upon demand by TPV (thermophotovoltaic) cells. Silicon is considered in this study as PCM (phase change material) due to its extremely high latent heat (1800 J/g or 500 Wh/kg), melting point (1410 C), thermal conductivity (~25 W/mK), low cost (less than $2/kg or $4/kWh) and abundance on earth. The proposed system enables an enormous thermal energy storage density of ~1 MWh/m3, which is 10e20 times higher than that of lead-acid batteries, 2e6 times than that of Li-ion batteries and 5e10 times than that of the current state of the art LHTES systems utilized in CSP (concentrated solar power) applications. The discharge efficiency of the system is ultimately determined by the TPV converter, which theoretically can exceed 50%. However, realistic discharge efficiencies utilizing single junction TPV cells are in the range of 20e45%, depending on the semiconductor bandgap and quality, and the photon recycling efficiency. This concept has the potential to achieve output electric energy densities in the range of 200-450 kWhe/m3, which is comparable to the best performing state of the art Lithium-ion batteries.

Multiple solutions and numerical analysis to the dynamic and stationary models coupling a delayed energy balance model involving latent heat and discontinuous albedo with a deep ocean

We study a climatologically important interaction of two of the main components of the geophysical system by adding an energy balance model for the averaged atmospheric temperature as dynamic boundary condition to a diagnostic ocean model having an additional spatial dimension. In this work, we give deeper insight than previous papers in the literature, mainly with respect to the 1990 pioneering model by Watts and Morantine. We are taking into consideration the latent heat for the two phase ocean as well as a possible delayed term. Non-uniqueness for the initial boundary value problem, uniqueness under a non-degeneracy condition and the existence of multiple stationary solutions are proved here. These multiplicity results suggest that an S-shaped bifurcation diagram should be expected to occur in this class of models generalizing previous energy balance models. The numerical method applied to the model is based on a finite volume scheme with nonlinear weighted essentially non-oscillatory reconstruction and Runge–Kutta total variation diminishing for time integration.

Evapotranspiration of grapevine trained to a gable trellis system under netting and black plastic mulching

The evapotranspiration (ET c) of a table grape vineyard (Vitis vinifera, cv. Red Globe) trained to a gable trellis under netting and black plastic mulching was determined under semiarid conditions in the central Ebro River Valley during 2007 and 2008. The netting was made of high-density polyethylene (pores of 12 mm2) and was placed just above the ground canopy about 2.2 m above soil surface. Black plastic mulching was used to minimize soil evaporation. The surface renewal method was used to obtain values of sensible heat flux (H) from high-frequency temperature readings. Later, latent heat flux (LE) values were obtained by solving the energy balance equation. For the May–October period, seasonal ET c was about 843 mm in 2007 and 787 mm in 2008. The experimental weekly crop coefficients (K cexp) fluctuated between 0.64 and 1.2. These values represent crop coefficients adjusted to take into account the reduction in ET c caused by the netting and the black plastic mulching. Average K cexp values during mid- and end-season stages were 0.79 and 0.98, respectively. End-season K cexp was higher due to combination of factors related to the precipitation and low ET o conditions that are typical in this region during fall. Estimated crop coefficients using the Allen et al. (1998) approach adjusting for the effects of the netting and black plastic mulching (K cFAO) showed a good agreement with the experimental K cexp values.

Evapotranspiration and crop coefficients of rice (Oryza sativa L.) under sprinkler irrigation in a semiarid climate determined by the surface renewal method

The evapotranspiration (ETc) of sprinkler-irrigated rice was determined for the semiarid conditions of NE Spain during 2001, 2002 and 2003. The surface renewal method, after calibration against the eddy covariance method, was used to obtain values of sensible heat flux (H) from high-frequency temperature readings. Latent heat flux values were obtained by solving the energy balance equation. Finally, lysimeter measurements were used to validate the evapotranspiration values obtained with the surface renewal method. Seasonal rice evapotranspiration was about 750–800 mm. Average daily ETc for mid-season (from 90 to 130 days after sowing) was 5.1, 4.5 and 6.1 mm day−1 for 2001, 2002 and 2003, respectively. The experimental weekly crop coefficients fluctuated in the range of 0.83–1.20 for 2001, 0.81–1.03 for 2002 and 0.84–1.15 for 2003. The total growing season was about 150–160 days. In average, the crop coefficients for the initial (Kcini), mid-season (Kcmid) and late-season stages (Kcend) were 0.92, 1.06 and 1.03, respectively, the length of these stages being about 55, 45 and 25 days, respectively.

Ion charge number and flux saturation effects in the corona of a laser-irradiated pellet

The quasisteady structure of the corona of a laser-irradiated pellet is completely determined for arbitrary Z, (ion charge number} and re/ra (ratio of critical and ablation radii), and for heat-flux saturation factor/above approximately 0.04. The ion-to-electron temperature ratio at rc grows sensibly with Z,; all other quantities depend weakly and nonmonotonically on Z,. For rc /ra close to unity, and all Z, of interest (Z, < 47}, the flow is subsonic at rc. For a given laser power W, flux saturation may decrease (low/) or increase (high/) the ablation pressure Pa relative to the value obtained when saturation is not considered; in some cases a decrease in/with W fixed increases Pa. For intermediate^ ~0.1), Pa cc (W/r* )2/3 p\n\pc = critical density), independently of rc/ra; for/~0.6, Pa «s larger by a factor of about [rc/raf13. For rjra > 1.2 roughly, the mass ablation rate is C{Z,) [{m/kZ.f^Kr^Pl) l,\ independent of pc and/, and barely dependent on Z,(m, is ion mass; k, Boltzmann's constant; K, conductivity coefficient; and C, a tabulated function).

Nonlocal electron heat relaxation in a plasma shock at arbitrary ionization number

A recently obtained nonlocal expression for the electron heat flux valid for arbitrary ionization numbers Z is used to study the structure of a plane shock wave in a fully ionized plasma. Nonlocal effects are only important in the foot of the electronic preheating region, where the electron temperature gradient is the steepest. The results are quantified as a function of a characteristic Knudsen number of that region. This work also generalizes to arbitrary values of Z previous results on plasma shock wave structure.

Metodología de simulación numérica de la transferencia de calor por ebullición en película en torno a superficies esféricas refrigeradas por líquido saturado en convección libre: aplicación a reactores nucleares de lecho de bolas (PBR)

En la presente Tesis se realizó un análisis numérico, usando el código comercial Ansys-Fluent, de la refrigeración de una bola de combustible de un reactor de lecho de bolas (PBR, por sus siglas en inglés), ante un escenario de emergencia en el cual el núcleo sea desensamblado y las bolas se dejen caer en una piscina de agua, donde el mecanismo de transferencia de calor inicialmente sería la ebullición en película, implicando la convección y la radiación al fluido. Previamente se realizaron pruebas de validación, comparando los resultados numéricos con datos experimentales disponibles en la literatura para tres geometrías diferentes, lo cual permitió seleccionar los esquemas y modelos numéricos con mejor precisión y menor costo computacional. Una vez identificada la metodología numérica, todas las pruebas de validación fueron ampliamente satisfactorias, encontrándose muy buena concordancia en el flujo de calor promedio con los datos experimentales. Durante estas pruebas de validación se lograron caracterizar numéricamente algunos parámetros importantes en la ebullición en película con los cuales existen ciertos niveles de incertidumbre, como son el factor de acoplamiento entre convección y radiación, y el factor de corrección del calor latente de vaporización. El análisis térmico de la refrigeración de la bola del reactor por ebullición en película mostró que la misma se enfría, a pesar del calor de decaimiento, con una temperatura superficial de la bola que desciende de forma oscilatoria, debido al comportamiento inestable de la película de vapor. Sin embargo, la temperatura de esta superficie tiene una buena uniformidad, notándose que las áreas mejor y peor refrigeradas están localizadas en la parte superior de la bola. Se observó la formación de múltiples domos de vapor en diferentes posiciones circunferenciales, lo cual causa que el área más caliente de la superficie se localice donde se forman los domos más grandes. La separación entre los domos de vapor fue consistente con la teoría hidrodinámica, con la adición de que la separación entre domos se reduce a medida que evolucionan y crecen, debido a la curvatura de la superficie. ABSTRACT A numerical cooling analysis of a PBR fuel pebble, after an emergency scenario in which the nucleus disassembly is made and the pebbles are dropped into a water pool, transmitting heat by film boiling, involving convection and radiation to the fluid, is carried out in this Thesis. First, were performed validation tests comparing the numerical results with experimental works available for three different geometries, which allowed the selection of numerical models and schemes with better precision and lower computational cost. Once identified the numerical methodology, all validation tests were widely satisfactory, finding very good agreement with experimental works in average heat flux. During these validation tests were achieved numerically characterize some important parameters in film boiling with which there are certain levels of uncertainty, such as the coupling factor between convection and radiation, and the correction factor of the latent heat of vaporization. The thermal analysis of pebble cooling by film boiling shows that despite its decay heat, cooling occurs, with pebble surface temperature descending from an oscillatory manner, due to the instability of the vapor film. However, the temperature of this surface has a good uniformity, noting that the best and worst refrigerated area is located at the top of the pebble. The formation of multiple vapor domes at different circumferential positions is observed, which cause that the hottest area of the surface was located where biggest vapor domes were formed. The separation between vapor domes was consistent with the hydrodynamic theory, with the addition that the separation is reduced as the vapor dome evolves and grows, due to the surface curvature.

Bubble growth by injection of gas into viscous liquids in cylindrical and conical tubes.

The laminar low Mach number flow of a gas in a tube is analyzed for very small and very large values of the inlet-to-wall temperature ratio. When this ratio tends to zero, pressure forces confine the cold gas to a thin core around the axis of the tube. This core is neatly bounded by an ablation front that consumes it at a finite distance from the tube inlet. When the temperature ratio tends to infinity, the temperature of the gas increases smoothly from the wall to the axis of the tube and the shear stress and heat flux are positive at the wall despite the fact that the viscosity and thermal conductivity of the gas scaled with their inlet values tend to zero at the wall.