SEARCH · Engineering Papers

Results for “Dynamic Boundary”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Temperature and radiation effects on brittle versus ductile fracture behavior in miscible phase boundaries: insight from atomistic simulations

Temperature- and irradiation-assisted failure mechanisms in miscible phase boundaries are clarified via atomistic calculations. We first establish the temperature-dependent brittle-to-ductile transition in U–Zr miscible phase boundaries. Our results confirm that these boundaries are mostly brittle at low temperatures and ductile at elevated temperatures. We then investigate the changes induced by irradiation on the fracture mechanisms in such phase boundaries. The irradiation-induced defect accumulation follows a multi-stage process that starts with the accumulation of isolated small dislocation loops before transitioning to the saturation and growth of larger dislocation loops and end up with a reorganization into forest dislocations. The accumulation of loops is the primary feature to participate in the delineation between brittle and ductile interfacial fracture in irradiated phase boundaries. At low damage levels, radiation defect interactions with the crack tip are limited and U–Zr miscible boundaries fail through the emission of dislocations ahead of the crack tip followed by brittle cleavage in agreement with the classical Griffith’s criterion for crack stability. At higher damage levels, the failure mode transitions from brittle crack growth to ductile void growth. In this case, the microcrack tip is blunted by the high density of pre-existing, radiation-induced defects in the vicinity of the crack. This interaction leads to the development and growth of a cavity at the interface as opposed to interfacial cleavage.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Interference of electrical double layers: Confinement effects on structure, dynamics, and screening of ionic liquids

Ionic liquids are widely used as electrolytes in electronic devices in which they are subject to nanoconfinement within nanopores or nanofilms. Because the intrinsic width of an electrical double layer is on the order of several nanometers, nanoconfinement is expected to fundamentally alter the double layer properties. Furthermore, in confined systems, a large portion of the ions are interfacial, e.g., at the electrode interface, leading to significant deviations of electrostatic screening and ion dynamics as compared to bulk properties. In this article, we systematically investigate the interference between electrical double layers for nanoconfined ionic liquids and the resulting influence on the structure, dynamics, and screening behavior. We perform molecular dynamics simulations for the ionic liquids [BMIm + ][BF$^-_4$] and [BMIm + ][PF$^-_6$] confined between two flat electrodes at systematic separation distances between 1.5 nm and 4.5 nm for both conducting and insulating boundary conditions. We find that while ion dynamics is expectedly slower than in the bulk (by ~2 orders of magnitude), there is an unexpected non-linear trend with the confinement length that leads to a local maximum in dynamic rates at ~3.5–4.5 nm confinement. We show that this nonlinear trend is due to the ion correlation that arises from the interference between opposite double layers. We further evaluate confinement effects on the ion structure and capacitance and investigate the influence of electronic polarization of the ionic liquid on the resulting properties. This systematic evaluation of the connection between electrostatic screening and structure and dynamics of ionic liquids in confined systems is important for the fundamental understanding of electrochemical supercapacitors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Associating damage nucleation and distribution with grain boundary characteristics in Ta

Here, the effect of grain boundary (GB) characteristics (as represented by misorientation angle) on its dynamic response is investigated by sampling a large set of 185 GBs in Ta bi-crystals. These GBs encompass both ordered and disordered structures – typically referred to as metastable structures, in order to reflect the diversity in the local atomic structure observed in realistic GBs. The results reveal four distinct regimes in damage distribution and failure modes as a function of misorientation angle. The above trends also correlate with the variation in the corresponding void nucleation stress.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗

DOE OSTI2025

Molecular-dynamics study of diffusional creep in uranium mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. Yet, many aspects of its mechanical behavior and microstructural features are currently unknown. In this paper, molecular dynamics (MD) simulations are used to study UN's diffusional creep. Nanometer-sized polycrystals are used to simulate diffusional creep and to calculate an effective GB width. It is found that Nabarro-Herring creep is not dominant in the temperature range of 1700–2000 K and that the dominant diffusional creep mechanism is Coble creep with an activation energy of 2.28 ± 0.09 eV. A method is proposed to calculate the diffusional GB width and its temperature dependence in polycrystals. The effective GB width of UN is calculated as 2.69 ± 0.08 nm. This value fits very well with the prefactor of the phenomenological Coble creep formula. It is demonstrated that the most comprehensive thermal creep model for UN can be represented as the combination of our Coble creep model and the dislocation creep model proposed by Hayes et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Impedance Based Fixture Neutralization (FINE) Method to Replicate Field Environments in Laboratory Settings

Test article vibration test responses will vary from one test laboratory to another due to differences in shaker and fixture dynamic characteristics. This is also seen between the field and laboratory due to the differing dynamic characteristics between field assembly loading and boundary conditions compared to the laboratory configuration. This work introduces a technique called Impedance Fixture Neutralization which customizes input forces to cause consistent responses for the Device Under Test across different vibration testing conditions. The customized force neutralizes the dynamic variations between configurations for the Device Under Test. The device responses can be replicated in several situations: when using a different fixture with the same attachment points, when the test article mounting location on the fixture changes, when the force location changes, or any combination of these situations. Impedance Fixture Neutralization uses the uncoupled dynamic characteristics of the test article and excitation fixtures to customize the input thereby causing the same test article responses between two mounting configurations (either field to laboratory or two different laboratory configurations). The application of the technique is shown using an analytical model of a two-beam system and analytically using experimental FRFs from plate and frame component hardware characterization tests. In both cases a device under test is attached to two test fixtures wherein the dynamic differences in the configurations are neutralized.

Laboratory↗

DOE OSTI2023

High-temperature chromium diffusion in austenitic stainless steel: Ab initio molecular dynamics simulations

Chromium self-diffusion through stainless steel (SS) matrix and along grain boundaries is an important mechanism controlling SS structural materials corrosion. For this work, Cr diffusion in austenitic SS was simulated using canonical ab initio molecular dynamics with realistic models of type-316 SS bulk, with and without Cr vacancies, and a low-energy Σ3 twin boundary typically observed at active corrosion sites. Cr self-diffusion coefficients at 750 and 850 °C calculated using Einstein’s diffusion equation are 4.2 × 10 –6 and 8.1 × 10 –6 Å 2 ps –1 in pristine bulk, 3.8 × 10 –3 and 5.5 × 10 –3 Å 2 ps –1 in bulk including Cr vacancies, and 9.5 × 10 –2 and 1.0 × 10 –1 Å 2 ps –1 at a Σ3[1 1 1]60° twin boundary.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Permeation of niobium through grain boundaries in copper

Low mutual solubility is expected to avert interdiffusion in layered composites of phase-separating metals. However, we show that Nb diffuses through polycrystalline Cu—despite the minimal bulk solubility of these elements—due to short circuit transport of Nb along certain grain boundaries in Cu. Atomistic modeling demonstrates that Nb-permeable Cu grain boundaries exhibit negative enthalpy of mixing of Nb, resulting in an enthalpically-stabilized solution of highly mobile Nb atoms that easily diffuse through the boundary. By contrast, Mo, which also has minimal solubility with Cu, has positive enthalpy of mixing at Nb-permeable Cu grain boundaries and does not diffuse through them. Furthermore, our findings suggest material selection and grain boundary engineering as pathways for designing improved diffusion barriers and thermally stable laminate composites.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Permeation of niobium through grain boundaries in copper

36 MATERIALS SCIENCE↗

DOE OSTI2025

Fundamental microscopic properties as predictors of large-scale quantities of interest: Validation through grain boundary energy trends

Correlations between fundamental microscopic properties computable from first principles, which we term canonical properties, and complex large-scale quantities of interest (QoIs) provide an avenue to predictive materials discovery. Here, we propose that such correlations can be efficiently discovered through simulations utilizing approximate interatomic potentials (IPs), which serve as an ensemble of “synthetic materials”. As a proof of principle we build a regression model relating canonical properties to the symmetric tilt grain boundary (GB) energy curves in face-centered cubic crystals, characterized by the scaling factor in the universal lattice matching model of Runnels et al. (2016), which we take to be our QoI. Our analysis recovers known correlations of GB energy to other properties and discovers new ones. We also demonstrate, using available density functional theory (DFT) GB energy data, that the regression model constructed from IP data is consistent with DFT results, confirming the assumption that the IPs and DFT belong to same statistical pool and thereby validating the approach. Regression models constructed in this fashion can be used to predict large-scale QoIs based on first-principles data and provide a general method for training IPs for QoIs beyond the scope of first-principles calculations.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Deformation accommodating periodic computational domain for a uniform velocity gradient

Many multiscale methods for granular materials use periodic computational domains to consider particle scale interactions and then calculate the stress to drive the continuum scale calculations. For problems involving large material deformations, the computation domain often needs to be reinitialized because of the distortion, causing the loss of the history information of the system. This work introduces an algorithm to accommodate a large deformation of the material while maintaining the computational domain cuboid to avoid the domain reinitialization during the computation. The algorithm uses a rotating frame of reference, in which the velocity gradient can be represented by an upper triangular matrix. The deformation caused by the upper triangular matrix is treated by the image system implied in the periodicity to maintain the computational domain cuboid. The effect from the rotation of the reference frame is considered using the inertial forces. Finally, simulations of simple and pure shear motions are carried out to illustrate the algorithm.

42 ENGINEERING↗

DOE OSTI2021

Competing magnetic interactions and emergent phase diagrams in double perovskite Y 2 Ni x Co 1–x MnO 6

We present a comparative study of double perovskites Y 2 Ni x Co 1–x MnO 6 for x = 1, 0.5, and 0. The polycrystalline samples of Y 2 Ni x Co 1–x MnO 6 with space group P2 1 /n were synthesized via sol-gel technique. X-ray photoelectron spectroscopy (XPS) confirms the presence of majority Ni 2+ /Mn 4+ in Y 2 NiMnO 6 (YNMO) and a mixed valence state Ni 2+/3+ /Co 2+/3+ /Mn 4+/3+ in co-doped Y 2 Ni 0.5 Co 0.5 MnO 6 (YNCMO). The temperature dependent magnetization results suggest homogeneous substitution of Ni ions with Co ions with minimal contribution of 3+ cationic magnetic interactions in YNCMO. The superexchange ferromagnetic (FM) interaction has a major contribution to the magnetism of YNMO. The strength of antiferromagnetic (AFM) coupling resulting from the antiphase boundaries in Y 2 Ni x Co 1–x MnO 6 increases with the increase in the Co concentration. The AC susceptibility (AC-χ) study further validates the presence of higher density of antiphase boundaries, which is reflected by the dynamics of domain walls in YNCMO. The isothermal magnetic entropy change (ΔS M ) as functions of temperature and magnetic field is exploited to assess the stabilization of different magnetic phases. Finally, the findings of ΔS M (T,µ 0 H) lead to the proposed new magnetic phase diagrams for Y 2 Ni x Co 1–x MnO 6 for x = 1 and 0.5, in comparison with the previously established phase diagram of Y 2 CoMnO 6 .

36 MATERIALS SCIENCE↗

DOE OSTI2021

Hot deformation behavior of an industrially cast large grained low density austenitic steel

Hot compression testing was performed on specimens of an Fe–30Mn–9Al–1Si-0.9C-0.5Mo wt.% low density steel which were extracted from an industrial casting with large grain size (320 ± 50 μm). Grain boundaries act as nucleation sites for dynamic recrystallization in the present material and the large grain size results in fewer nucleation sites, which delays dynamic recrystallization to larger strains and/or higher temperatures and causes higher rates of work hardening relative to similar materials but with smaller grain sizes. Here, a relatively high activation energy for hot deformation of 470 ± 90 kJ mol –1 was determined from a Zener-Hollomon analysis of the flow stresses at different temperatures and strain rates, attributed in part to the large grain size.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Mechanisms of Three-Dimensional Solid-Phase Epitaxial Crystallization of Strontium Titanate

Strontium titanate (SrTiO 3 , STO) is a complex metal oxide with a cubic perovskite crystal structure. Due to its easily described and understood crystal structure in the cubic phase, STO is an ideal model system for exploring the mechanistic details of solid-phase epitaxy (SPE) in complex oxides. SPE is a crystallization approach that aims to guide crystal growth at low homologous temperatures to achieve targeted microstructures. Beyond planar thin films, SPE can also exploit the addition of a chemically inert, noncrystallizing, amorphous obstacle in the path of crystallization to generate complex three-dimensional structures. The introduction of this mask fundamentally alters the SPE process, inducing a transition from two- to three-dimensional geometries and from vertical to lateral crystal growth under the influence of the crystal/mask/amorphous boundary. Using a combination of molecular dynamics simulations and experiments, we identify several unique phenomena in the nanoscale growth behaviors in both conventional (unmasked) and masked SPE. Examining conventional SPE of STO, we find that crystallization at the interface is strongly correlated to, and potentially driven by, density fluctuations in the region of the amorphous STO near the crystalline/amorphous interface with a strong facet dependence. In the masked case, we find that the crystalline growth front becomes nonplanar near contact with the mask. We also observe a minimum vertical growth requirement prior to lateral crystallization. Both phenomena depend on the relative bulk and interfacial free energies of the three-phase (crystal/mask/amorphous) system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Global peak water limit of future groundwater withdrawals

Over the past 50 years, humans have extracted the Earth’s groundwater stocks at a steep rate, largely to fuel global agro-economic development. Given society’s growing reliance on groundwater, we explore ‘peak water limits’ to investigate whether, when and where humanity might reach peak groundwater extraction. Using an integrated global model of the coupled human–Earth system, we simulate groundwater withdrawals across 235 water basins under 900 future scenarios of global change over the twenty-first century. Here we find that global non-renewable groundwater withdrawals exhibit a distinct peak-and-decline signature, comparable to historical observations of other depletable resources (for example, minerals), in nearly all (98%) scenarios, peaking on average at 625 km 3 yr -1 around mid-century, followed by a decline through 2100. The peak and decline occur in about one-third (82) of basins, including 21 that may have already peaked, exposing about half (44%) of the global population to groundwater stress. Most of these basins are in countries with the highest current extraction rates, including the United States, Mexico, Pakistan, India, China, Saudi Arabia and Iran. Furthermore, these groundwater-dependent basins will probably face increasing costs of groundwater and food production, suggesting important implications for global agricultural trade and a diminished role for groundwater in meeting global water demands during the twenty-first century.

58 GEOSCIENCES↗

DOE OSTI2023

Implications of phonon anisotropy on thermal conductivity of fluorite oxides

Abstract Fluorite oxides are attractive ionic compounds for a range of applications with critical thermal management requirements. In view of recent reports alluding to anisotropic thermal conductivity in this face-centered cubic crystalline systems, we perform a detailed analysis of the impact of direction-dependent phonon group velocities and lifetimes on the thermal transport of fluorite oxides. We demonstrate that the bulk thermal conductivity of this class of materials remains isotropic despite notable anisotropy in phonon lifetime and group velocity. However, breaking the symmetry of the phonon lifetime under external stimuli including boundary scattering present in nonequilibrium molecular dynamics simulations of finite size simulation cell gives rise to apparent thermal conductivity anisotropy. We observe that for accurate determination of thermal conductivity, it is important to consider phonon properties not only along high symmetry directions commonly measured in inelastic neutron or x-ray scattering experiments but also of those along lower symmetry. Our results suggests that certain low symmetry directions have a larger contribution to thermal conductivity compared to high symmetry ones.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Phase transitions in high-purity zirconium under dynamic compression

In this work, we present results from ramp compression experiments on high-purity Zr that show the α → ω, ω → β, as well as reverse β → ω phase transitions. Simulations with a multi-phase equation of state and phenomenological kinetic model match the experimental wave profiles well. While the dynamic α ! ! transition occurs ~ 9 GPa above the equilibrium phase boundary, the ω → β transition occurs within 0.9 GPa of equilibrium. We estimate that the dynamic compression path intersects the equilibrium ω – β line at P = 29:2 GPa, and T = 490 K. The thermodynamic path in the interior of the sample lies ~ 100 K above the isentrope at the point of the ω → β transition. Approximately half of this dissipative temperature rise is due to plastic work, and half is due to the non-equilibrium α → ω transition. The inferred rate of the α → ω transition is several orders of magnitude higher than that measured in dynamic diamond anvil cell (DDAC) experiments in an overlapping pressure range. We discuss a model for the influence of shear stress on the nucleation rate. We find that the shear stress s ji has the same effect on the nucleation rate as a pressure increase δP = cϵ ij s ji /(ΔV/V); where c is a geometric constant ~ 1 and, ϵ ij are the transformation shear strains. The small fractional volume change ΔV/V ≈ 0:1 at the α → ω transition amplifies the effect of shear stress, and we estimate that for this case δP is in the range of several GPa. Correcting our transition rate to a hydrostatic rate brings it approximately into line with the DDAC results, suggesting that shear stress plays a significant role in the transformation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Multi-Scale Thermal Analysis for Design of SiC-Based Medium Voltage Motor Drive: Preprint

This paper describes a multi-scale thermal analysis approach for the design of an air-cooled 1.7-kV SiC MOSFET-based medium-voltage variable-speed motor drive. The National Renewable Energy Laboratory supported the project objectives by providing thermal management design support to the project partner The Ohio State University (OSU). The scope of the models and limited resources required efficient and flexible thermal models to be developed. Two modeling techniques are described that significantly reduced model run time and enabled more complex models to be run faster while retaining needed accuracy. The first technique uses the effectiveness-NTU method to extract convection boundary conditions from a computational fluid dynamics (CFD) model that can be applied to a fast-running finite element analysis model. The second is a porous media technique that enables system-level CFD simulations that incorporate effects from heat exchangers (e.g., pin fin heat sinks) that run in a fraction of the time required for fully resolved CFD simulations. The multi-scale approach to the thermal analysis enabled fast and accurate simulation for the converter design ranging from the die level up to the full system with 36 submodules. The modeling results were validated against experimental data from system tests performed by OSU.

DIRECT ENERGY CONVERSION↗