Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coalescence”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗
DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗
DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

In-situ imaging of nucleation and growth of superlattices from nanoscale colloidal nanoparticles

Nanoparticle-based superlattices have attracted extensive research efforts due to their versatile composition and packing-dependent optical, magnetic, electronic, catalytic, mechanical properties. Yet the experimental understanding of the thermodynamics and pathways of the crystallization processes for these systems, which govern their ultimate morphology, size, and packing structure, has for long remained underexplored. Here, this minireview highlights recent integration of liquid-phase transmission electron microscopy, which allows for real-space, real-time direct imaging of the crystallization processes at the single-particle level, with simulation and the conceptual framework established by pioneers such as Dr. A. A. Chernov. Different aspects of superlattice formation were revealed by direct imaging for the first time, such as the existence of prenucleation precursor in nonclassical crystallization, surface energy-dependent growth, and coalescence, allowing for the charting of phase coordinates and thermodynamic quantities at the nanoscale. We discuss the similarities and differences of crystal growth behaviors in atomic and nanoparticle systems, as well as engineering opportunities to fashion, shape, and achieve quality control in superlattice formation from nanoparticles, for next-generation optical and mechanical metamaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

In-situ TEM study of Kr ion irradiation tolerance of SiFeOC nanocomposite

Here, in this work, ion irradiation of polymer derived SiFeOC nanocomposite was carried out using 1.2 MeV Kr ions at room temperature and 600°C. The starting composite was composed of Fe 3 Si, SiC, SiOC, SiO 2 , and graphitic C. In-situ TEM investigations show uniform distribution of nano-crystalline Fe 3 Si and SiC phases in the amorphous SiOC matrix. During ion irradiation, the SiOC bulk microstructure and interfaces between Fe 3 Si or SiC crystallites and the SiOC matrix remain defect-free and demonstrate outstanding ion irradiation resistance. At room temperature, the crystalline domains are stable up to 2 dpa. At 600°C, Fe 3 Si crystallites are more stable than SiC; SiC crystallites are stable up to 4 dpa while the Fe 3 Si crystallites are stable up to 10 dpa. These crystallites also coalescence and amorphize simultaneously during ion irradiation. The exceptional tolerance to defect formation and irradiation of the SiFeOC nanocomposite provides important guidance to developing irradiation resistant fuels for advanced gas cooled reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Type IV failure in weldment of creep resistant ferritic alloys: II. Creep fracture and lifetime prediction

The creep rupture behavior (Type IV failure) for weldments of creep strength enhanced ferritic steel is numerically analyzed, using an integrated microstructure- and micromechanics-based finite element model. To account for the large microstructure gradients across weldments, a two-dimensional digital microstructure is constructed based on the actual observed microstructure of ferritic steel weldment by using the Voronoi-tessellation method. According to the fracture mechanism studies and literature experimental observations, the Type IV failure is identified as an intergranular creep fracture in the fine-grained or intercritical heated affected zone (FGHAZ or ICHAZ). In the present study, the following micromechanics model is employed to determine the micromechanical and microstructural origins for the failure process above, accounting for the underlying physical fracture mechanisms at different length scales, including nucleation of grain boundary cavities, their growth by competition of grain boundary diffusion and grain interior creep, viscous grain boundary sliding, and the emergence of microcracks by coalescence and their evolution to the ultimate failure. The methodology demonstrates the capabilities in modeling the Type IV failure and providing quantitative creep rupture lifetime prediction which shows an excellent agreement with long-term creep experimental data for creep strength enhanced ferritic steels and their weldments. In particular, the drop-off in time to rupture at high temperatures and low stress levels in the creep rupture curves is quantitatively predicted, and the transition of failure mechanisms from creep-controlled to diffusion-controlled creep fracture mechanism is illustrated.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Grain-boundary fracture mechanisms in Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolytes: When phase transformation acts as a temperature-dependent toughening mechanism

Garnet-type, solid electrolytes, such as Li 7 La 3 Zr 2 O 12 (LLZO), are a promising alternative to liquid electrolytes for lithium-metal batteries. However, such solid-electrolyte materials frequently exhibit undesirable lithium (Li) metal plating and fracture along grain boundaries. In this work, we employ atomistic simulations to investigate the mechanisms and key fracture properties associated with intergranular fracture along one such boundary. Our results show that, in the case of a Σ5 (310) grain boundary, this boundary exhibits brittle fracture behavior, i.e. the absence of dislocation activity ahead of the propagating crack tip, accompanied with a decrease in work of separation, peak stress, and maximum stress intensity factor as the temperature increases from 300 K to 1500 K. As the crack propagates, we predict two temperature-dependent Li clustering regimes. For temperatures at or below 900 K, Li tends to cluster in the bulk region away from the crack plane driven by a void-coalescence mechanism concomitant a simultaneous cubic-to-tetragonal phase transition. The tetragonalization of LLZO in this temperature regime acts as an emerging toughening mechanism. At higher temperatures, this phase transition mechanism is suppressed leading to a more uniform distribution of Li throughout the grain-boundary system and lower fracture properties as compared to lower temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The application of synchrotron micro-computed tomography to characterize the three-dimensional microstructure in irradiated nuclear fuel

This research demonstrated the first time that synchrotron X-ray micro-computed tomography (μ-CT) was utilized to analyze the three-dimensional (3-D) morphology of the microstructure of irradiated nuclear fuel. Both absorption contrast and propagation-based phase contrast enhanced μ-CT techniques were evaluated. Porosity and phase region volumes were determined in neutron irradiated uranium 10 wt.% zirconium (U-10Zr) using propagation-based phase contrast enhanced μ-CT. Volumes were determined by image processing software and algorithms, with the porosity volume subsequently compared to the ASTM E562 method. The 3-D porosity of a small region (~ 8 x 10 5 µm 3 ) located in the intermediate fuel redistribution zone that was removed via a focused ion beam was determined to be 7.2% via image processing. A detailed 3-D analysis of the porosity with respect to volume and morphology was provided. Five growth stages of pore evolution were observed including: nucleation, growth, coalescence, interconnected porosity, and extended/interconnected porosity. Four distinct phase regions were identified, including the pores, via both absorption contrast and propagation-based phase contrast enhanced μ-CT. A quantitative comparison of fuel swelling between this localized microscopic U-10Zr fuel region, previously irradiated in the Fast Flux Test Facility, and historic U-10Zr fuels irradiated in the Experimental Breeder Reactor-II (EBR-II) was conducted. It was determined that for similar burnups, the localized swelling in this microscopic fuel region was 7.7%, which was approximately one-half of the value of the macroscopic swelling reported in EBR-II irradiated fuel. With the advent of new specimen preparation tools for irradiated materials, as well as advances in characterization techniques, researchers can now obtain previously inaccessible microstructural insight into irradiated fuels. This possibility naturally lends itself to improved fuel performance models for the nuclear community.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Evolution of dislocation loops and voids in post-irradiation annealed ThO 2 : A combined in-situ TEM and cluster dynamics investigation

The effect of isochronal annealing on the evolution of dislocation loop and void population in proton irradiated ThO 2 has been investigated. Post-irradiation annealing in other actinide oxides like UO 2 shows significant loop coarsening. ThO 2 samples were irradiated with 2 MeV protons up to a dose of 0.1 dpa at 600°C. Post-irradiation isochronal annealing was performed at 600, 800, 1000 and 1100°C for 1 h at each temperature using in-situ TEM. Only faulted 1/3<111> type dislocation loops were observed, and their sizes and distribution were characterized. Further, the population of self-interstitial atom (SIA) dislocation loops did not show any significant growth and coarsening. Additionally, nanometric voids were observed at annealing temperatures of 1000 and 1100°C. Using cluster dynamics (CD), we have studied the nucleation and growth of point defects and defect clusters, i.e., SIA prismatic dislocation loops and nanometric and sub-nanometric voids in proton irradiated ThO 2 . The CD model was further utilized to predict the growth and coarsening of loops and voids during isochronal annealing at the experimental and higher temperatures. The model did not predict significant SIA loop growth which closely corresponds to the TEM observations. CD predicted SIA loop coarsening is insignificant even at high annealing temperature of 1500°C because the model only considers the growth of defect clusters by absorption of like point defects, i.e., SIA loops absorb interstitials and voids absorb vacancies, and cannot account for their migration and coalescence due to elastic interaction. The CD model also predicts the evolution of nanometric voids having mean size within the error bounds of TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Evolution of extended defects in UO 2 during high temperature annealing

The evolution of thermophysical characteristics of uranium oxide (UO 2 ) under irradiation is a critical aspect of nuclear fuel performance of light water reactors. This study examines the effect of high-temperature annealing at 1000 °C - 1600 °C on defect evolution in krypton irradiated UO 2 at room temperature. Transmission electron microscopy analysis reveals that the size of irradiation-induced bubbles and dislocation loops increases but their density decreases with annealing temperatures. Perfect dislocation loops with = <110> dominate, regardless of the annealing temperature. A rate theory model for extended defect evolution is applied to conclude that current understanding of dynamical loop and bubble evolution at high-temperature annealing is governed by defect coalescence.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Impact of grain boundary and surface diffusion on predicted fission gas bubble behavior and release in UO 2 fuel

In this work, we quantify the impact of grain boundary (GB) and surface diffusion on fission gas bubble evolution and fission gas release in UO 2 nuclear fuel using simulations with a hybrid phase field/cluster dynamics model. Here, we begin with a comprehensive literature review of uranium vacancy and xenon atom diffusivity in UO 2 through the bulk, along GBs, and along surfaces. In our model we represent fast GB and surface diffusion using a heterogeneous diffusivity that is a function of the order parameters that represent bubbles and grains. We find that the GB diffusivity directly impacts the rate of gas release via GB transport, and that the GB diffusivity is likely below 104 times the lower value from Olander and van Uffelen. We also find that the surface diffusivity impacts bubble coalescence and mobility, and that the bubble surface diffusivity is likely below 10 -4 times the value from Zhou and Olander.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Microstructurally informed synchrotron x-ray analysis revealing helium defect transitions in ultrafine grained tungsten

The formation of insoluble gaseous defects in materials due to nuclear transmutation or ion implantation involves the diffusion of impurity atoms to form atomic defect clusters that coalesce into bubbles or cavities and ultimately degrade the material properties. Transmission electron microscopy (TEM) is limited in its ability to resolve sub-nanometer gas clusters whereas X-ray diffraction (XRD) provides information pertaining to local atomic changes. Here, in this study, helium (He) implanted ultrafine grained tungsten is explored through a multimodal defect characterization campaign combining TEM-informed Small Angle X-ray Scattering (SAXS) analysis, XRD lattice parameter measurements, and nanoscale He cluster quantification from a region of reciprocal space accessible via Wide Angle X-ray Scattering (WAXS). Moderate elevated temperature implantations are shown to produce high concentrations of sub-nanoscale He clusters and small, homogeneously distributed cavities, which collectively are linked to lattice expansion and further substantiated through complementary atomistic simulations. Increased implantation temperatures encourage the diffusion of these defects to the grain boundaries (GBs), leading to lattice relaxation and the growth of large GB cavities manifesting as bimodal size distributions in the SAXS analysis. Overall, our results demonstrate the utility of multimodal synchrotron X-ray analysis in bridging the gap between microscale He cavity quantification and atomic-scale defect analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nanostructured ligament and fiber Al–doped Li 7 La 3 Zr 2 O 12 scaffolds to mediate cathode-electrolyte interface chemistry

In this study, scaffold structures of electrospun aluminum–substituted lithium lanthanum zirconate Li 7 La 3 Zr 2 O 12 (Al-LLZO) were synthesized and used as an additive in a LiNi 0.6 Mn 0.2 Co 0.2 O 2 composite cathode. The scaffolds were crystalized in the cubic phase after calcination at 700 °C. The Al-LLZO scaffold morphology was dependent on the precursor formulation (aqueous and dimethylformamide. The aqueous precursors resulted in scaffolds of densely coalesced ligaments, whereas the dimethylformamide precursors resulted in high–aspect ratio nanofiber scaffolds. The long-term cycling stability and rate performance of the cells were found to depend on the Al-LLZO scaffold morphology. The uniformly dispersed Al-LLZO fibers resulted in a more stable cathode electrolyte interface formation through the reduced decomposition of the LiPF 6 salt during cycling, resulting in a better high-rate and long-term cycling performance.

25 ENERGY STORAGE↗
DOE OSTI2026

Mechanistic origin of solvent-dependent thermal stability in sodiated Sn anodes for sodium-ion batteries

Understanding the thermal stability of high-energy density alloy anodes is critical for the safe deployment of sodium-ion batteries (SIBs). Here, accelerating rate calorimetry (ARC), post-mortem characterizations, and density functional theory (DFT) calculations are combined to understand the thermal reactivity of fully sodiated Sn, Sn-hard carbon (HC) blends, and HC anodes in carbonate- and ether-based electrolytes. ARC measurements show that propylene carbonate (PC) causes earlier self-heating rate (SHR) onset and higher reactivity than tetraethylene glycol dimethyl ether (TEGDME), indicating inferior thermal stability. Sodiated Sn exhibits better thermal stability than sodiated HC, while Sn-HC blends show intermediate behavior that improves with increasing Sn content. Post-ARC analyses reveal desodiation of Na15Sn4 to metallic Sn with particle coalescence, whereas Sn-HC blends and HC retain finer morphologies. PC promotes Sn oxidation to SnO, while TEGDME suppresses oxide formation; NaPF6-containing electrolytes additionally form NaF. DFT calculations show that PC adsorption lowers Na extraction energy and enhances interfacial electronic interactions, facilitating Na release and reductive decomposition. These results establish a direct correlation between solvent-dependent reaction pathways and thermal stability in SIB alloy anodes.

Accelerating rate calorimetry↗
DOE OSTI2021

The influence of preexisting host rock damage on fault network localization

The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.

58 GEOSCIENCES↗
DOE OSTI2020

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗
DOE OSTI2023

Achieving superior mechanical properties: Tailoring multicomponent microstructure in AISI 9254 spring steel through a two-stage Q&P process and nanoscale carbide integration

In the pursuit of lightweight, durable steel, we have successfully developed a multicomponent structure in AISI 9254 spring steel using a two-stage quenching and partitioning (Q&P) process. The primary objective of this process was to engineer an optimized microstructure consisting of nanobainite, martensite, and nano-carbides. Utilizing the insights gained from the results of advanced techniques such as X-ray diffraction (XRD), scanning electron microscopy (SEM), transmission electron microscopy (TEM), and atom probe tomography (APT) performed on the as-received AISI 9254 spring steel, we refined the quenching and partitioning (Q&P) path, leading to the successful establishment of a bainitic transformation for superior mechanical properties. Our tensile tests revealed a high yield strength (≈ 1600 ± 25 MPa) and ultimate tensile strength (≈ 1850 ± 50 MPa), along with considerable elongation (≈ 11.15 ± 0.25%). We also identified that pre-formed martensite lath defects and high silicon content play crucial roles during the Q&P process, preventing carbide coalescence and increasing strain-hardening capacity. Finally, this study demonstrates the potential of a Q&P process to generate high-strength, ductile steel for automotive and aerospace applications.

36 MATERIALS SCIENCE↗