Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interfacial flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

AFM Slip Length Measurements for Water at Selected Phyllosilicate Surfaces

Most reported slip length measurements have been made at the surfaces of synthetic materials and modified synthetic materials. In contrast, few slip length measurements at the surface of unmodified natural mineral surfaces have been reported. In this regard, flow at the silica face surfaces of the phyllosilicate minerals, talc and mica, was considered. A slip boundary condition was expected at the nonpolar hydrophobic silica surface of talc leading to enhanced flow, and a no-slip boundary condition was expected at the hydrophilic silica surface of mica. Atomic force microscopy (AFM) slip length measurements were made at the talc and mica surfaces. The slip length results for the hydrophobic silica surface of talc were contrasted to the results for the hydrophilic silica surface of mica (no-slip flow). The results are discussed based on molecular dynamics simulations (MDS), as reported in the literature, and AFM images of surface nanobubbles. For nonpolar hydrophobic surfaces (such as talc), it is doubtful that the MDS interfacial water structure and the water exclusion zone (3.2 Å) account for the AFM slip flow with slip lengths as great as 95 nm. Rather, a better explanation for the AFM slip flow condition is based on reduced interfacial viscosity due to the presence of dissolved gas and the accommodation of pancake nanobubbles at the talc surface having a height dimension of magnitude similar to the slip length.

42 ENGINEERING↗
DOE OSTI2018

Multiphase flow studies for microscale hydrodynamics in the structured packed column

Post combustion carbon capture by solvent absorption in a structured packed column is a promising technology for mitigating greenhouse gas emissions. Computational fluid dynamic (CFD) modelling of such a column is a challenging multiscale problem due to the range in length and time scales. The microscale hydrodynamics play a key role in the overall column efficiency with the interfacial area significantly influencing the mass transfer between the gas and liquid phases. In this context, multiphase flow simulations using the volume of fluid (VOF) method in a representative elementary unit (REU) of the packed column can provide fundamental insights into the microscale hydrodynamics, such as, interfacial area and liquid holdup. The present paper systematically examines the impact of various factors (e.g., physical properties and contact angle) on the interfacial area. The results are compared with existing correlations and a scaling analysis is also performed. The solvent physical properties are characterized by the Kapitza number (Ka), a dimensionless number that depends only on fluid properties. So that the Ka number decreases with increasing viscosity. At a fixed liquid load, the interfacial area and liquid holdup are observed to increase with decreasing Ka number. The impact of contact angle (i.e., solid surface characteristics) is effectively investigated by modifying the wall boundary conditions. The interfacial area and liquid holdup are found to decrease with increasing contact angle. Subsequently, a phenomenological correlation for interfacial area is proposed that includes the impact of these parameters. Finally, this correlation may be used to predict the interfacial area for gas-liquid flow in a structured packing for rivulet to fully wetted flow regimes.

42 ENGINEERING↗
DOE OSTI2020

The Sensitivity of Hyporheic Exchange to Fractal Properties of Riverbeds

Hyporheic exchange in riverbeds is driven by current-bed topography interactions. Because riverbeds exhibit topographic roughness across scales, from individual grains to bedforms and bars, they can exhibit fractal patterns. This study analyzed the influence of fractal properties of riverbed topography on hyporheic exchange. A set of synthetic fractal riverbeds with different scaling statistics was used as inputs to sequentially coupled numerical simulations of turbulent channel flow and hyporheic flow. In the analysis, the maximum power spectrum (dune size) and the fractal dimension (topographic complexity) were considered as independent variables and we then investigated how interfacial fluxes and hyporheic travel times are functionally related to these variables. As the maximum power spectrum increases (i.e., dune height to flow depth ratio), the average interfacial flux increases logarithmically whereas it increases exponentially with an increase in fractal dimension. Hyporheic exchange is more sensitive to additional roughness (larger fractal dimensions) than to bedform size (larger maximum power). Our results imply that fractal properties of riverbeds are crucial to predicting hyporheic exchange. The predictive relationships we propose could be integrated with reduced complexity, large-scale models. Further, they can also be used to design artificial topographies that target hyporheic ecosystem services.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Single Bubble Drag Force Evaluation in Turbulent Flow based on DNS Results

Two-phase turbulent flows are widely encountered in the light water reactor engineering. The interaction between liquid phase and gas phase is bi-directional: the liquid turbulence influences the dispersion and distribution of the bubbles which in turn modulates the liquid turbulence with bubble-induced turbulence. Considering a single bubble, the effect of homogeneous turbulent flow is to modify the drag force compared to that in laminar flow. The bubble can modify the homogeneous turbulence through formation of vorticities behind the bubble. The correct description of the closure laws for the interfacial forces is of great importance in numerical simulation of bubbly flows. Various closure laws for the drag coefficient have been proposed in the literature in the last few decades. Among them, Tomiyama proposed a widely accepted correlation based on spherical bubbles in pure systems. While the development of closure model for two-phase laminar flow has reached a certain level of confidence, more studies and investigations are required for turbulent flows. In this paper, we want to numerically analyze the influence of homogeneous turbulent flow on the drag force experienced by a single bubble and compare its behavior with the laminar flow conditions based on Tomiyama's correlation. A set of parametric studies will be performed to justify the relationship between typical two-phase turbulent flow parameters and drag force. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Two-dimensional lift-up problem for a rigid porous bed

The present study analytically reinvestigates the two-dimensional lift-up problem for a rigid porous bed that was studied by Mei, Yeung, and Liu [“Lifting of a large object from a porous seabed,” J. Fluid Mech. 152, 203 (1985)]. Mei, Yeung, and Liu proposed a model that treats the bed as a rigid porous medium and performed relevant experiments. In their model, they assumed the gap flow comes from the periphery of the gap, and there is a shear layer in the porous medium; the flow in the gap is described by adhesion approximation [D. J. Acheson, Elementary Fluid Dynamics (Clarendon, Oxford, 1990), pp. 243-245.] and the pore flow by Darcy’s law, and the slip-flow condition proposed by Beavers and Joseph [“Boundary conditions at a naturally permeable wall,” J. Fluid Mech. 30, 197 (1967)] is applied to the bed interface. In this problem, however, the gap flow initially mainly comes from the porous bed, and the shear layer may not exist. Although later the shear effect becomes important, the empirical slip-flow condition might not physically respond to the shear effect, and the existence of the vertical velocity affects the situation so greatly that the slip-flow condition might not be appropriate. In contrast, the present study proposes a more general model for the problem, applying Stokes flow to the gap, the Brinkman equation to the porous medium, and Song and Huang’s [“Laminar poroelastic media flow,” J. Eng. Mech. 126, 358 (2000)] complete interfacial conditions to the bed interface. The exact solution to the problem is found and fits Mei’s experiments well. The breakout phenomenon is examined for different soil beds, mechanics that cannot be illustrated by Mei’s model are revealed, and the theoretical breakout times obtained using Mei’s model and our model are compared. The results show that the proposed model is more compatible with physics and provides results that are more precise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗
DOE OSTI2019

Recording interfacial currents on the subnanometer length and femtosecond time scale by terahertz emission

Electron dynamics at interfaces is a subject of great scientific interest and technological importance. Detailed understanding of such dynamics requires access to the angstrom length scale defining interfaces and the femtosecond time scale characterizing interfacial motion of electrons. In this context, the most precise and general way to remotely measure charge dynamics is through the transient current flow and the associated electromagnetic radiation. Here, we present quantitative measurements of interfacial currents on the subnanometer length and femtosecond time scale by recording the emitted terahertz radiation following ultrafast laser excitation. We apply this method to interlayer charge transfer in heterostructures of two transition metal dichalcogenide monolayers less than 0.7 nm apart. We find that charge relaxation and separation occur in less than 100 fs. This approach allows us to unambiguously determine the direction of current flow, to demonstrate a charge transfer efficiency of order unity, and to characterize saturation effects.

36 MATERIALS SCIENCE↗
NASA NTRS2008

Tungsten and Barium Transport in the Internal Plasma of Hollow Cathodes

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the flow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushedback to the emitter surface by the electric field and drag from the xenon ion flow. Thisbarium ion flux is sufficient to maintain a barium surface coverage at the downstream endgreater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length,so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollowcathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

Hollow cathodes↗
DOE OSTI2021

Shear-free mixing to achieve accurate temporospatial nanoscale kinetics through scanning-SAXS: ion-induced phase transition of dispersed cellulose nanocrystals

Time-resolved in situ characterization of well-defined mixing processes using small-angle X-ray scattering (SAXS) is usually challenging, especially if the process involves changes of material viscoelasticity. In specific, it can be difficult to create a continuous mixing experiment without shearing the material of interest; a desirable situation since shear flow both affects nanoscale structures and flow stability as well as resulting in unreliable time-resolved data. Here, we demonstrate a flow-focusing mixing device for in situ nanostructural characterization using scanning-SAXS. Given the interfacial tension and viscosity ratio between core and sheath fluids, the core material confined by sheath flows is completely detached from the walls and forms a zero-shear plug flow at the channel center, allowing for a trivial conversion of spatial coordinates to mixing times. With this technique, the time-resolved gel formation of dispersed cellulose nanocrystals (CNCs) was studied by mixing with a sodium chloride solution. It is observed how locally ordered regions, so called tactoids, are disrupted when the added monovalent ions affect the electrostatic interactions, which in turn leads to a loss of CNC alignment through enhanced rotary diffusion. The demonstrated flow-focusing scanning-SAXS technique can be used to unveil important kinetics during structural formation of nanocellulosic materials. However, the same technique is also applicable in many soft matter systems to provide new insights into the nanoscale dynamics during mixing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental methodology and theoretical framework in describing constrained plastic flow of FCC microscale tensile specimens

Small scale mechanical testing (SSMT) allows for the quantification of the mechanical responses for microscale and nanoscale material volumes. It has significant potential in improving our fundamental understanding of plastic deformation and in providing an engineering solution for testing complex and/or critical materials. Importantly, previous studies have revealed that plastic intermittency, manifested as strain bursts, are characteristic of the plastic flow behavior at the small length scale. In this work, interfacial constraint has been utilized to suppress strain bursts in order to stabilize plasticity. A new parameter, blocked volume ratio (BVR), is introduced to quantify the degree of constraint for a microscale specimen due to an interface. Here, we showed that the plastic flow behavior in FCC Ni-based metals is highly dependent on the BVR of the specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Stability, intermittency and universal Thorpe length distribution in a laboratory turbulent stratified shear flow

Stratified shear flows occur in many geophysical contexts, from oceanic overflows and river estuaries to wind-driven thermocline layers. In this study, we explore a turbulent wall-bounded shear flow of lighter miscible fluid into a quiescent fluid of higher density with a range of Richardson numbers $0.05\lesssim Ri\lesssim 1$ . In order to find a stability parameter that allows close comparison with linear theory and with idealized experiments and numerics, we investigate different definitions of $Ri$ . We find that a gradient Richardson number defined on fluid interface sections where there is no overturning at or adjacent to the maximum density gradient position provides an excellent stability parameter, which captures the Miles–Howard linear stability criterion. For small $Ri$ the flow exhibits robust Kelvin–Helmholtz instability, whereas for larger $Ri$ interfacial overturning is more intermittent with less frequent Kelvin–Helmholtz events and emerging Holmboe wave instability consistent with a thicker velocity layer compared with the density layer. We compute the perturbed fraction of interface as a quantitative measure of the flow intermittency, which is approximately 1 for the smallest $Ri$ but decreases rapidly as $Ri$ increases, consistent with linear theory. For the perturbed regions, we use the Thorpe scale to characterize the overturning properties of these flows. The probability distribution of the non-zero Thorpe length yields a universal exponential form, suggesting that much of the overturning results from increasingly intermittent Kelvin–Helmholtz instability events. Finally, the distribution of turbulent kinetic energy, conditioned on the intermittency fraction, has a similar form, suggesting an explanation for the universal scaling collapse of the Thorpe length distribution.

58 GEOSCIENCES↗
DOE OSTI2021

Are Contact Angle Measurements Useful for Oxide-Coated Liquid Metals?

This work establishes that static contact angles for gallium-based liquid metals have no utility despite the continued and common use of such angles in the literature. In the presence of oxygen, these metals rapidly form a thin (~1–3 nm) surface oxide “skin” that adheres to many surfaces and mechanically impedes its flow. This property is problematic for contact angle measurements, which presume the ability of liquids to flow freely to adopt shapes that minimize the interfacial energy. We show here that advancing angles for a metal are always high (>140°)—even on substrates to which it adheres—because the solid native oxide must rupture in tension to advance the contact line. The advancing angle for the metal depends subtly on the substrate surface chemistry but does not vary strongly with hydrophobicity of the substrate. During receding measurements, the metal droplet initially sags as the liquid withdraws from the “sac” formed by the skin and thus the contact area with the substrate initially increases despite its volumetric recession. Furthermore, the oxide pins at the perimeter of the deflated “sac” on all the surfaces are tested, except for certain rough surfaces. With additional withdrawal of the liquid metal, the pinned angle gets smaller until eventually the oxide “sac” collapses. Thus, static contact angles can be manipulated mechanically from 0° to >140° due to hysteresis and are therefore uninformative. We also provide recommendations and best practices for wetting experiments, which may find use in applications that use these alloys such as soft electronics, composites, and microfluidics.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Influence of convection on microstructure

In eutectic growth, as the solid phases grow they reject atoms to the liquid. This results in a variation of melt composition along the solid/liquid interface. In the past, mass transfer in eutectic solidification, in the absence of convection, was considered to be governed only by the diffusion induced by compositional gradients. However, mass transfer can also be generated by a temperature gradient. This is called thermotransport, thermomigration, thermal diffusion or the Soret effect. A theoretical model of the influence of the Soret effect on the growth of eutectic alloys is presented. A differential equation describing the compositional field near the interface during unidirectional solidification of a binary eutectic alloy was formulated by including the contributions of both compositional and thermal gradients in the liquid. A steady-state solution of the differential equation was obtained by applying appropriate boundary conditions and accounting for heat flow in the melt. Following that, the average interfacial composition was converted to a variation of undercooling at the interface, and consequently to microstructural parameters. The results obtained show that thermotransport can, under certain circumstances, be a parameter of paramount importance.

Wilcox, William R.↗
DOE OSTI2023

INVESTIGATION OF ROCK-FLUID INTERACTIONS USING GEOMATERIAL MICROFLUIDICS

Illite and Illite-Smectite clays are reportedly the Caney Shale's most prevalent clay minerals. The swelling and fines migration of clays may occur when water-based fluids contact the clay minerals. This may result in the reduction of pore space and permeability, and wettability alteration. Therefore, we quantified the swelling potentials, fines migration and wettability alteration of Illite and Illite-Smectite clays upon exposure to model brines and a produced water, and a model oil and a crude oil. We measured density, viscosity, pH, TDS, and interfacial tension as they are important in multiphase flow and displacement characteristics. Geomaterial microfluidics is recently being used to study physicochemical interactions of solid-fluid systems. It facilitates visualizing the surface behavior upon exposure to various fluids. This study includes preparation and characterization of Illite-Smectite clay-coated glass capillary tubes. Illite and Illite-Smectite clay-coated geomaterial microfluidic surfaces were prepared to investigate the effects of first contact fluid, brine salinity, and aging on wettability. In this work, advancing and receding contact angles of model brine-air, model oil-air, produced water-air, crude oil-air, model brine-model oil systems, and produced water-crude oil systems have been measured in the untreated, and clay-coated microfluidic channels of 1000 µm width and 50 µm depth. The aging effect was studied by retaining the non-aqueous and aqueous solutions in the corresponding microfluidic channels for approximately 24 h in the case of channels first contacted by aqueous and non-aqueous solutions, respectively. The advancing and receding contact angles of Caney shale-crude oil-produced water systems were also measured. The experimental findings indicate the swelling potential, fines migration, and wettability alteration of Illite and Illite-Smectite clays. The results signify the influence of first contact fluid, brine salinity, and aging on wettability of the untreated and clay-coated microfluidic surfaces. The experiments also demonstrate the wetting nature of the Caney shale samples from reservoir and non-reservoir zones.

02 PETROLEUM↗
DOE OSTI2016

Results from the assessment of COBRA-IE multidimensional flow capability with selected canonical problems

A systematic assessment of the capabilities of thermal-hydraulic codes such as RELAP5-3D and COBRA-IE to handle multi-dimensional flow is ongoing. The assessment is being performed in hierarchical fashion by isolating simple processes first and building complexity along the way. The first stage in the process has been to develop simple and fundamental canonical problems which challenge key aspects of the multidimensional flow process while removing other complex processes such as interfacial mass and heat transfer. Results of this work have been presented for RELAP5-3D at the NURETH-16 conference [1], [2], [3]. The purpose of this paper is to illustrate similar results obtained using COBRA-IE. The results are presented for two canonical problems. The first pertains to the simulation of Faraday waves (CP06) in an annulus. The second problem compares the solution of 2D single-phase, inviscid, frictionless, incompressible, irrotational flow (CP07) to the potential flow solution. Results from COBRA-IE are compared with similar results obtained with RELAP5-3D which were previously published in the 2015 NURETH papers series on this topic. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1999

Condensation of Forced Convection Two-Phase Flow in a Miniature Tube

A physical/mathematical model of annular film condensation at the inlet of a miniature tube has been developed. In the model, the liquid flow is coupled with the vapor flow along the liquid-vapor interface through the interfacial temperature, heat flux, shear stress, and pressure jump conditions due to surface tension effects. The model predicts the shape of the liquid-vapor interface along the condenser and leads to the conclusion that there is complete condensation at a certain distance from the condenser inlet. The numerical results show that complete condensation of the incoming vapor is possible at comparatively low heat loads and that this is a special case of a more general condensation regime with two-phase bubbly flow downstream of the initial annular film condensation region. Observations from the flow visualization experiment confirm the existence and qualitative features of annular film condensation leading to the complete condensation phenomenon in a small diameter (3.25 mm) circular tube condenser.

Begg, E.↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗