Search NASA⌕ Search

SEARCH · Search NASA

Results for “stress voiding”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Elastic Properties of Plasma-Exposed Tungsten Predicted by Molecular-Dynamics Simulations

We report results of systematic molecular-dynamics computations of the elastic properties of single-crystalline tungsten containing structural defects, voids and overpressurized He nanobubbles, related to plasma exposure of tungsten serving as a plasma- facing component (PFC) in nuclear fusion devices. Our computations reveal that the empty voids are centers of dilatation resulting in development of tensile stress in the tungsten matrix, whereas He-filled voids (nanobubbles) introduce compressive stress in the plasma-exposed tungsten. We find that the dependence of the elastic moduli of plasma-exposed tungsten, namely, the bulk, Young, and shear modulus, on its void fraction follows a universal exponential scaling relation. We also find that the elastic moduli of plasma-exposed tungsten soften substantially as a function of He content in the tungsten matrix, following an exponential scaling relation; this He-induced exponential softening is in addition to the softening caused in the matrix with increasing temperature. A systematic characterization of the dependence of the elastic moduli on the He bubble size reveals that He bubble growth affects significantly both the bulk modulus and the Poisson ratio of plasma-exposed tungsten, while its effect on the Young and shear moduli of the plasma-exposed material is weak. Our findings contribute directly to the development of a structure-properties database that is required for the predictive modeling of the dynamical response of PFCs in nuclear fusion devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Identification of mechanisms driving heterogeneous void growth in ductile aluminum

Void growth plays a central role in ductile fracture, yet the specific mechanisms that control this remain obscure. Classical models, such as those proposed by Rice and Tracey in 1969, are able to capture average rates of void growth, but cannot capture the heterogeneity of individual void growth. Building on recent work, the present study employs laboratory-based diffraction contrast tomography and in-situ x-ray computed tomography to investigate the effect of grain structure and other microstructural factors on void growth in an Al-2219 alloy. Crystal plasticity finite element (CP-FE) modeling is used alongside experimental data to evaluate the contributions of local mechanical states, grain orientation, grain size, and neighboring microstructural features. No strong linear relationships are found with any of the considered descriptors and void growth rate. Potential complex nonlinear relationships are explored with the use of a random forest regression model, which identifies initial void volume, void aspect ratio, local normal stress state, local shear stress state, and local equivalent plastic strain (EQPS) as features that most improve void growth rate predictions. The combination of these analyses suggests that these features should be prioritized to improve models of void growth.

Diffraction contrast tomography (DCT)↗
DOE OSTI2025

Morphological Heterogeneity Impact of Film Solid-State Cathode on Utilization and Fracture Dynamics

Structural heterogeneity in solid-state batteries can impact material utilization and fracture mechanisms. Crystallographically oriented lithium cobalt oxide film cathodes serve as a model electrode system for exploring how void distribution contributes to stress relief and build up during cycling. Real- and reciprocal-space operando and ex situ synchrotron based experiments are utilized to understand structural changes across multiple length scales contribute to stress generation and fracture. Nanotomography uncovers a depth-dependent porosity variation in the pristine electrode and highlights preferential fracture in regions of lower porosity during delithiation. Energy-dispersive X-ray diffraction and 3D X-ray absorption near-edge spectroscopy (XANES) reveal the underutilization of cathode material in these regions. 3D XANES also confirms preferential delithiation near the sub-grain boundaries. Chemo-mechanical modeling coupled with site-specific mechanical characterization demonstrate how stress accumulation in dense regions of the electrode leads to fracture and underutilization of active material. In conclusion, our findings reveal the importance of materials design to alleviate stress in small-volume changing cathodes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Advanced Materials and Manufacturing Technologies Nondestructive Examination Efforts at Idaho National Laboratory: Report of FY-24 Efforts

This report details FY-24 nondestructive examination (NDE) efforts at Idaho National Laboratory (INL) in support of the Advanced Materials and Manufacturing Technologies (AMMT) program. While the goal of this endeavor is to develop a multi-modal, multi-length scale workflow for nondestructive characterization of advanced manufactured (AM) nuclear reactor components, substantial development remains until this is a reality. In support of this effort X-ray computed tomography (XCT), X-ray diffraction (XRD), neutron computed tomography (nCT), neutron diffraction, lock-in thermography (LIT), multi-point lock-in thermography (MLIT), and positron annihilation spectroscopy (PAS) were all used on AM specimens to examine defects such as voids, porosity, and residual stress. In addition to summarizing the results of these NDE applications, recommendations for integrating these into a more comprehensive undertaking to promote NDE of engineering-scale components are also included.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Fracture characteristics and heat treatment of laser powder bed fusion additively manufactured GRCop-84 copper

Laser Powder Bed Fusion (LPBF) of Glenn Research Copper 84 (GRCop-84), a Cr 2 Nb (8 at. % Cr, 4 at. % Nb) precipitation hardened alloy, produces a fully dense high conductivity alloy with a yield strength of 500 MPa and ultimate tensile strength (UTS) of 740 MPa, superior to other competing copper alloys, and 20% elongation at fracture for material stressed perpendicular to the build direction. The high thermal stability of the Cr 2 Nb precipitate in the copper matrix reduces coarsening and maintains a 300 MPa yield, 520 MPa UTS and 26% elongation after a 900 °C, 5-h heat treatment, while a 3 h 450 °C heat treatment increases yield to 810 MPa, UTS to 970 MPa with 9% elongation at fracture, for samples stressed perpendicular to the build direction. Tensile strength anisotropy based on print direction was attributed to internal stress and columnar grain formation. Void nucleation during tensile fracture was initiated by brittle fracture of precipitate particles within the copper matrix. Fracture cusps contain matching precipitate fragment geometry on opposing sides located near the cusp center in at least 80% of fracture cusps. An optimal precipitate size of 100 nm is predicted for maximum tensile strength from precipitates on fracture surfaces, while tensile testing with varying heat treatments shows maximum strength with 100 nm and smaller precipitates. Cr2Nb precipitates are shown to transition between polycrystalline and monocrystalline structures at high temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
DOE OSTI2024

Shock-induced twinning/detwinning and spall failure in Cu–Ta nanolaminates at atomic scales

Here, this study provides new insights into the role of interfaces on the deformation and failure mechanisms in shock-loaded Cu–Ta–Cu trilayer system. The thickness of the Ta layer, piston velocities, and shock pulse durations were varied to explore the impact of impedance mismatch and loading conditions on spallation behavior and twin formation. It was found that the interfaces play a crucial role in the dynamic response of these multilayered systems since secondary reflection waves generated at the interfaces significantly affected the peak stress and pressure profiles, influencing void nucleation and failure modes. In the trilayer systems, failure predominantly occurred at interfaces and within the Ta layer, with void nucleation sites and twinning behavior being markedly different compared to single-crystal Cu and Ta. Increasing the Ta layer thickness modified the wave interactions, leading to different failure locations. Higher piston velocities were associated with increased spall strength by enhancing wave interactions and void formation, particularly at the interfaces and within the Ta layer, under specific configurations. Additionally, shorter shock pulse durations facilitated earlier initiation of the release fan, reducing twin formation and altering the failure dynamics by accelerating twin annihilation and pressure release.

36 MATERIALS SCIENCE↗
DOE OSTI2026

TEM Approaches for Microstructure-Informed Prediction of Mechanical Properties in Structural Alloys

Predicting the mechanical performance of structural alloys from their evolving microstructure remains a major challenge in materials science, particularly for nuclear structural materials, where irradiation-induced defects span multiple types and length scales and interact through complex mechanisms. The dispersed barrier hardening (DBH) [1] and Friedel–Kroupa–Hirsch (FKH) [2,3] models have been widely used to evaluate the hardening contributions of individual obstacles and to estimate tensile strength from quantified microstructures; however, when multiple size-dependent obstacles coexist and evolve, predicting temperature-dependent tensile strength becomes significantly more complex, and a fully consistent hardening model is still lacking. Transmission electron microscopy (TEM) plays a central role in refining hardening models and enabling predictive assessments of tensile strength evolution by providing quantitative characterization of dislocations, irradiation-induced defects (e.g., dislocation loops and cavities), precipitates, and grain structure (Fig. 1.). These experimentally measured defect densities are incorporated into physically based hardening models with size- and shape- dependent obstacle strengths [4], using root-sum-square superposition for obstacles of comparable strength and linear superposition for dissimilar ones [5]. In addition, recent advances in TEM [6-8], including high-resolution imaging, 4D-STEM strain mapping, EDS/EELS elemental analysis, and flash-polishing-based TEM specimen preparation and extraction-replica methods (Fig. 2), further improve the accuracy of microstructural quantification. By comparison with prior studies as well as our own results, we show that when TEM-derived microstructural information is carefully integrated with physically grounded hardening models, yield strength (or irradiation-induced hardening) measured at room temperature can be predicted with good quantitative agreement across multiple alloy classes. In-situ TEM combined with high-temperature mechanical testing represents an important next step for refining hardening models by directly probing dislocation–obstacle interactions across varying irradiation doses and temperatures [9]. Because the barrier strength factor (α) depends on both temperature and obstacle size, it should not be treated as a constant fitting parameter; rather, it must be explicitly evaluated to achieve physically meaningful predictions of mechanical behaviour at operating temperatures. This presentation therefore discusses why all strengthening contributions (e.g., Peierls stress, solid-solution strengthening, voids, bubbles, dislocation loops, dislocation lines, and grain boundaries) must be considered collectively, why appropriate superposition methods are essential when obstacles possess different barrier strength factors, how hardness measurements can be meaningfully related to tensile properties, and how TEM-derived microstructural information can be systematically incorporated into hardening models. More broadly, it outlines a pathway toward microstructure-informed prediction of mechanical properties and supports the goal of establishing science-based tools for evaluating structural materials in extreme environments [10].

Lin, Yan-Ru [ORNL] (ORCID:0000000339991473)↗
DOE OSTI2025

Attribution of heterogeneous stress distributions in low-grain polycrystals under conditions leading to damage

In high-purity polycrystalline metallic materials, voids tend to favor grain boundaries as nucleation sites due to the elevated stress states produced by granular interactions and the weakened grain boundary from the relative atomic disorder. To quantify the key factors of this elevated stress state, simple compression of a small multi-grain cylinder of body-centered cubic tantalum was simulated using a single crystal plasticity model that incorporates non-Schmid effects. Four increasingly complex synthetic microstructures were created to tractably incorporate grain boundary interactions, and a statistically significant number of combinations were performed by varying the initial crystallographic orientations of the microstructure. Most of these simulations produce the maximum von Mises stress on a grain boundary and less frequently at the multi-grain junctions. To build a statistical model for the maximum von Mises stress at the grain boundary, physically based features that could contribute to the elevated stress state were selected. Then, a learning algorithm based on information theory was used to identify which of these features contributed the most information to the data set. The identified features include a grain’s propensity to accommodate both elastic and plastic deformations and their directional components. The misalignment of the direction of each grain’s mechanical response was found to be strongly correlated to the magnitude of the stress near the grain boundary. For all of the synthetic microstructures, the statistical models produce a residual distribution that is nearly Gaussian with a variance of, at most, 10% of the prior distribution. The successful performance of the statistical model implies the correct identification of the physical features that cause severe stress localization in polycrystalline materials. The statistical models constructed here can be used to formulate a physically motivated void nucleation model which is sensitive to a microstructure’s propensity to produce elevated stress states. As a result, these statistical models also enable the design of material microstructures, in which the crystallographic orientation is chosen to resist void nucleation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Role of grain boundary crystallography on void growth in FCC metals

The growth of an isolated spherical void at and near a grain boundary or a triple junction is studied by means of full-field crystal plasticity simulations with explicit representation of the void. To examine the fundamental aspects of void/grain boundary interactions in these cases, a three-dimensional large-strain elasto-viscoplastic fast-Fourier transform (LS-EVP-FFT) model with axisymmetric tensile periodic boundary conditions was developed. We employ this model to reveal the role of crystallography (grain orientations, grain boundary misorientation and inclination angle with respect to the largest component of the applied stress) and type of loading on void growth. It was found that grain orientation was the primary crystallographic feature governing void growth in strain-rate controlled simulations. For strains small enough to not significantly alter the local geometry, the overall growth of the void in bicrystal boundaries was not strongly affected by grain boundary inclination angle. The relative amount of void growth experienced within each grain in a bicrystal or tri-crystal was strongly dependent on the grain boundary inclination angle. We found a subset of orientations that exhibit faster void growth in bicrystal boundaries than inside the constituent single crystals. Last, we show that when placed primarily in one grain but near the grain boundary, void growth was similar to void growth when placed at the grain boundary. Furthermore, the growing void remained in the original grain and did not grow into the other grain, regardless of crystallographic character.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Effects of sample bias on adhesion of magnetron sputtered Cr coatings on SiC

Swelling of SiC at 300 C due to in-service neutron irradiation causes tensile residual stresses in coatings which are expected to adversely affect the performance of coated SiC composite fuel cladding for light water reactors. Matching the coating swelling with the substrate, a solution common for thermal expansion, is not practical in the case of neutron irradiation. Biasing samples during magnetron sputtering deposition induces compressive residual stress which may counteract this. In this study, chromium coatings were deposited on SiC by DC magnetron sputtering with no external heating at bias voltages of –50V, –75V, and –100V. The effects of the bias voltage on morphology, residual stress, microstrain, texture, and adhesion are shown. The low deposition temperature resulted in the coating microstructure evolution following an energetic particle bombardment dominated trend. Additionally, at the two lower bias voltages knock-on implantation dominated increasing the residual stress and microstrain while at the highest bias voltage, thermal spike migration allowed for defect relaxation. When the knock-on induced compressive residual stress exceeded 0.8 GPa microcrack formation in the SiC substrate decreased coating adhesion. While no microcracks formed at the lowest bias voltage, insufficient atomic mobility during coating growth lead to voids forming in the coating. A balance is needed to form void-free coatings that have high compressive residual stress.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Void coalescence and ductile failure in IN718 investigated via high-energy synchrotron X-ray tomography and diffraction

Ductile failure through the growth and coalescence of voids is of particular relevance for many engineering materials. Yet, the lack of experimental measurements of the mechanical state of the material at an appropriate length scale has limited further understanding of the regime of ductile failure. In this study, local grain-scale experimental measurements are used to determine their relevance for describing and predicting ductile failure. Additive manufactured materials, due to the presence of inherent porosity, as well as the ability to tailor internal porosity, provides a promising avenue to study ductile failure. Selective laser melting is used to manufacture a specially designed specimen with two large, internal voids, in addition to the natural porosity, which is characteristic to the production process. The initial porosity and its evolution upon tensile loading are characterized via micro tomography. Several locations of void coalescence are captured in the sample, revealing activity of multiple modes of failure. Finite element simulations, with a simplified J2 plasticity model, instantiated with the initial void structure, is deployed, in which geometric localizations of stress concentration corresponded with experimentally observed sites of coalescence but was inadequate in capturing the appropriate failure mechanism. Furthermore, from the experimental results of the far-field high-energy diffraction microscopy, the heterogeneous micromechanical state is identified and tracked around voids due to the local grain interactions. These experiments determine narrow bands of low stress triaxiality, at the onset of failure, which highlight the path of coalescence through intervoid shearing. Diffraction spot spreading analysis aided characterization of intragranular plasticity and strain heterogeneity, which can be coupled with high triaxiality to form the conditions inducive for coalescence through intervoid necking. In this study, the experimentally determined grain-scale description of the micromechanical state provides a physical basis that can accurately capture both the sites and the mechanism of void coalescence at the onset of ductile failure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Orientation-selected micro-pillar compression of additively manufactured 316L stainless steels: Comparison of as-manufactured, annealed, and proton-irradiated variants

In this work, irradiation response and deformation mechanisms of additively manufactured (AM) 316L stainless steel were studied by atomic scale characterization and micro-pillar compression. The AM 316L stainless steels were fabricated by direct energy deposition, a laser-based additive manufacturing process. Irradiation with 2 MeV protons at 360°C was performed to create ~1.8 displacements-per-atom (dpa) damage in AM 316L. Deformation behaviors of the as-manufactured, annealed, and proton-irradiated variants were studied, focusing on the effects of manufacturing-induced pores, residual stress, and irradiation-introduced defects (dislocation loops and voids). Micro-pillars were prepared from grains of pre-selected orientation, avoiding contributions of grain boundaries and allowing determination of resolved shear stress on {111} glide planes. Transmission electron microscopy was used to characterize the pre- and post-deformation microstructure. It was found that in the as-manufactured alloy variant, moving dislocations were the major deformation carrier, with noticeable blocking by fabrication-induced pores, In the annealed variant, hardness was reduced, and deformation was also accomplished by dislocation gliding. In the proton-irradiated variant, significant twinning was observed. Comparing measured resolved shear stress and predicted critical stress for dislocation dissociation, we conclude that irradiation hardening became high enough to activate twinning. Therefore, the deformation mechanism changes from dislocation gliding to twinning. The study is important for both processing optimization and performance evaluation of AM alloys for reactor applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗
DOE OSTI2024

Catch bond kinetics are instrumental to cohesion of fire ant rafts under load

Dynamic networks composed of constituents that break and reform bonds reversibly are ubiquitous in nature owing to their modular architectures that enable functions like energy dissipation, self-healing, and even activity. While bond breaking depends only on the current configuration of attachment in these networks, reattachment depends also on the proximity of constituents. Therefore, dynamic networks composed of macroscale constituents (not benefited by the secondary interactions cohering analogous networks composed of molecular-scale constituents) must rely on primary bonds for cohesion and self-repair. Toward understanding how such macroscale networks might adaptively achieve this, we explore the uniaxial tensile response of 2D rafts composed of interlinked fire ants ( S. invicta ). Through experiments and discrete numerical modeling, we find that ant rafts adaptively stabilize their bonded ant-to-ant interactions in response to tensile strains, indicating catch bond dynamics. Consequently, low-strain rates that should theoretically induce creep mechanics of these rafts instead induce elastic-like response. Our results suggest that this force-stabilization delays dissolution of the rafts and improves toughness. Nevertheless, above 35 % strain low cohesion and stress localization cause nucleation and growth of voids whose coalescence patterns result from force-stabilization. These voids mitigate structural repair until initial raft densities are restored and ants can reconnect across defects. However mechanical recovery of ant rafts during cyclic loading suggests that—even upon reinstatement of initial densities—ants exhibit slower repair kinetics if they were recently loaded at faster strain rates. These results exemplify fire ants’ status as active agents capable of memory-driven, stimuli-response for potential inspiration of adaptive structural materials.

Science & Technology - Other Topics↗
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↗