Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiscale structures”

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
NASA NTRS2020

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗
DOE OSTI2024

Small scale tensile testing technique for measuring grain boundary strength of neutron-irradiated materials in focused ion beam systems

Irradiation induced defects on grain boundaries (GBs) can lead to GB strength degradation and cause intergranular cracking in structural materials, affecting the safe operation and service life of nuclear reactors. Establishing the quantitative relationship between the strength and character/microchemistry for individual GB in an effort of integrating experimental measurements and multiscale simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enhanced multiscale restriction-smoothed basis (MsRSB) preconditioning with applications to porous media flow and geomechanics

A novel method to enable application of the Multiscale Restricted Smoothed Basis (MsRSB) method to non M-matrices is presented. The original MsRSB method is enhanced with a filtering strategy enforcing M-matrix properties to enable the robust application of MsRSB as a preconditioner. Through applications to porous media flow and linear elastic geomechanics, the method is proven to be effective for scalar and vector problems with multipoint finite volume (FV) and finite element (FE) discretization schemes, respectively. Realistic complex (un)structured two- and three-dimensional test cases are considered to illustrate the method's performance.

58 GEOSCIENCES↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Understanding the Performance Gap between Polycrystalline and Single-Crystal Nickel-Rich Layered Oxide Cathodes

Singe-crystal (SC) nickel-rich layered oxide cathodes, composed of boundary-free particles with high tap density, offer significant advantages in volumetric energy density and mechanical strength compared with polycrystalline (PC) cathode materials. However, as the nickel content increases (≥80%), SC Ni-rich cathodes often suffer from faster performance degradation than PC cathodes of the same composition, and the underlying causes of this discrepancy remain poorly understood. Herein, we reveal the distinct Ni redox behaviors that govern the electrochemical performance of SC and PC Ni-rich cathodes using multiscale and operando characterization techniques. Our results indicate that the increasingly heterogeneous Ni oxidation process in SC cathodes leads to the additional irreversible oxygen redox activity that deteriorates both the mechanical and chemical structures. In contrast, PC cathodes, despite with more pronounced surface reconstruction, exhibit greater chemomechanical stability due to homogeneous redox reactions during charging. Consequently, we find that bulk degradation, more than surface reactions, ultimately leads to fast capacity decay of SC Ni-rich cathodes during cycling. In conclusion, this work offers a comprehensive view on the impact of Ni redox evolutions on the chemomechanical stability in Ni-rich layered oxide cathodes, providing new insights into the longstanding performance gap between SC and PC cathodes, and guiding the rational design of Ni-rich cathode architectures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Design nanoporous metal thin films via solid state interfacial dealloying

Thin-film solid-state interfacial dealloying (thin-film SSID) is an emerging technique to design nanoarchitecture thin films. The resulting controllable 3D bicontinuous nanostructure is promising for a range of applications including catalysis, sensing, and energy storage. Using a multiscale microscopy approach, we combine X-ray and electron nano-tomography to demonstrate that besides dense bicontinuous nanocomposites, thin-film SSID can create a very fine (5–15 nm) nanoporous structure. Not only is such a fine feature among one of the finest fabrications by metal-agent dealloying, but a multilayer thin-film design enables creating nanoporous films on a wider range of substrates for functional applications. Through multimodal synchrotron diffraction and spectroscopy analysis with which the materials’ chemical and structural evolution in this novel approach is characterized in details, we further deduce that the contribution of change in entropy should be considered to explain the phase evolution in metal-agent dealloying, in addition to the commonly used enthalpy term in prior studies. Here, the discussion is an important step leading towards better explaining the underlying design principles for controllable 3D nanoarchitecture, as well as exploring a wider range of elemental and substrate selections for new applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Kinetic characterization of 3D magnetic reconnection: A transformative step (Final Report)

We investigated magnetic reconnection in the Earth’s magnetopause, magnetotail, and the bow shock, by means of fully kinetic simulations and analysis of space data obtained by NASA’s Magnetospheric Multiscale mission. We elucidated signatures of particle energization in particle distribution functions in the vicinity of reconnection X-line, properties of reconnection in shock-driven turbulence, and wave excitations and nonlinear structures. We advanced the understanding of magnetic reconnection from the viewpoint of particle kinetics in the diffusion region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

Electron Diffusion Regions in Magnetotail Reconnection Under Varying Guide Fields

Kinetic structures of electron diffusion regions (EDRs) under finite guide fields in magnetotail reconnection are reported. The EDRs with guide fields 0.14–0.5 (in unit of the reconnecting component) are detected by the Magnetospheric Multiscale spacecraft. The key new features include the following: (1) cold inflowing electrons accelerated along the guide field and demagnetized at the magnetic field minimum while remaining a coherent population with a low perpendicular temperature, (2) wave fluctuations generating strong perpendicular electron flows followed by alternating parallel flows inside the reconnecting current sheet under an intermediate guide field, and (3) gyrophase bunched electrons with high parallel speeds leaving the X‐line region. The normalized reconnection rates for the three EDRs range from 0.05 to 0.3. The measurements reveal that finite guide fields introduce new mechanisms to break the electron frozen‐in condition.

Chen, L.‐J.↗
NASA NTRS2021

Direct Multipoint Observations Capturing the Reformation of a Supercritical Fast Magnetosonic Shock

Using multipoint Magnetospheric Multiscale (MMS) observations in an unusual string-of-pearls configuration, we examine in detail observations of the reformation of a fast magnetosonic shock observed on the upstream edge of a foreshock transient structure upstream of Earth’s bow shock. The four MMS spacecraft were separated by several hundred km, comparable to suprathermal ion gyro-radius scales or several ion inertial lengths. At least half of the shock reformation cycle was observed, with a new shock ramp rising up out of the “foot” region of the original shock ramp. Using the multipoint observations, we convert the observed time-series data into distance along the shock normal in the shock’s rest frame. That conversion allows for a unique study of the relative spatial scales of the shock’s various features, including the shock’s growth rate, and how they evolve during the reformation cycle. Analysis indicates that: the growth rate increases during reformation, electron-scale physics play an important role in the shock reformation, and energy conversion processes also undergo the same cyclical periodicity as reformation. Strong, thin electron-kinetic-scale current sheets and large-amplitude electrostatic and electromagnetic waves are reported. Evidence is also presented of nonlinear wave decay from electromagnetic whistler-mode “lion roars” to electrostatic solitary waves in the downstream plasma regime. Results highlight the critical cross-scale coupling between electron-kinetic-and ion-kinetic-scale processes and details of the nature of nonstationarity, shock-front reformation at collisionless, fast magnetosonic shocks.

D L Turner↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Microstrain screening towards defect-less layered transition metal oxide cathodes

Microstrain and the associated surface-to-bulk propagation of structural defects are known to be major roadblocks to developing high-energy and long-life batteries. However, the origin and effects of microstrain during the synthesis of battery materials remain largely unknown. Here, in this study, we perform microstrain screening during real-time and realistic synthesis of sodium layered oxide cathodes. Evidence gathered from multiscale in situ synchrotron X-ray diffraction and microscopy characterization collectively reveals that the spatial distribution of transition metals within individual precursor particles strongly governs the nanoscale phase transformation, local charge heterogeneity and accumulation of microstrain during synthesis. This unexpected dominance of transition metals results in a counterintuitive outward propagation of defect nucleation and growth. These insights direct a more rational synthesis route to reduce the microstrain and crystallographic defects within the bulk lattice, leading to significantly improved structural stability. The present work on microstrain screening represents a critical step towards synthesis-by-design of defect-less battery materials.

25 ENERGY STORAGE↗
DOE OSTI2021

Multiscale operando X-ray investigations provide insights into electro-chemo-mechanical behavior of lithium intercalation cathodes

The electrochemical performance and cycle life of lithium-ion batteries (LIBs) depend on the electrochemical, chemical, and mechanical behavior of electrodes and electrolytes. Despite extensive studies conducted previously, challenges exist to decouple these behaviors, capture the evolution of electro-chemo-mechanical behavior in realistic conditions, and correlate atomic-scale stress evolution to micro-scale bulk mechanical degradation. Here, we report multiscale operando techniques to investigate polydisperse battery electrodes by integrating volume-averaged quantitative synchrotron X-ray scattering with high-resolution transmission X-ray microscopy (TXM). The former provides us information spanning a wide spatial range, from Angstrom-level atomic structures to micrometer-level particle scales, while the latter provides time-resolved 2D images of the particles during cycling. The complementarity of the two operando techniques is demonstrated by an over-lithiation test of LiCoO 2 electrodes, where particles crack and eventually pulverize. Additionally, the techniques are applied to study LiCoO 2 cycling stability from 3.0 V to 4.5 V. Operando X-ray scattering result shows nanometer-scale features keep forming in LiCoO 2 electrodes during cycling, resulting in an increased projected area observed by the TXM experiment. The formation of such features is inhibited by a polymer coating on the electrode, leading to vastly improved cycling stability. The polymer coating alleviates LiCoO 2 surface deterioration, reduces side product generation, and inhibits LiCoO 2 particles volume expansion during the cycling test. These operando multimodal X-ray techniques presented herein thus offer a novel, multiscale diagnostic modality for studying existing and emerging battery materials, aiding the development of next-generation LIBs.

25 ENERGY STORAGE↗
DOE OSTI2024

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗
DOE OSTI2020

Solitary Magnetic Structures at Quasi–Parallel Collisionless Shocks: Formation

Solitary magnetic structures known as SLAMS (short large-amplitude magnetic structures) have been considered as essential elements of collisionless shocks with quasi-parallel geometries. Yet the physics underlying their formation remains an open question. In this paper, we use measurements from the magnetospheric multiscale mission combined with fully kinetic simulations to study the formation of SLAMS. We find that gyro-resonance between solar wind ions and right-hand circularly polarized electromagnetic waves results in magnetic field amplification. Gyro-trapping by the growing magnetic field builds up the plasma density that further enhances the current and field. The solitary nature of SLAMS stems from a beat-like magnetic field envelope where the maximum sets the initial location for nonlinear growth. Our results present a conceptual advance on SLAMS, and may shed new light on the open question of magnetic field amplification at astrophysical shocks.

58 GEOSCIENCES↗
NASA NTRS2018

Simulation of Stochastic Mud-Crack Damage Formation in an Environmental Barrier Coating

The integrated Finite Element Analysis–Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button–sized disk–shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300oC to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.

residual tensile stress↗
NASA NTRS2023

Reconnection in a Pinch

A recently published analysis of current sheets has updated the classic Harris 1D static solution by considering multiple classes of charged particle trajectories in a generalized and dynamic current sheet. It uses a 1D PIC simulation to describe dynamic pinching and bifurcation of the current sheet. These 1D results strongly suggest that properties of the inflowing plasma, including the plasma beta, have an important effect on the equilibrium thickness of the pinched current sheet. Since 1D studies cannot describe magnetic reconnection, the time appears right to carry such 1D studies over to 2D or 3D simulations to explore current sheet thickness effects on reconnection. The Magnetospheric Multiscale Mission (MMS), with its well-resolved multipoint measurements of collisionless plasma and fields, has found that collisionless reconnection is accompanied by non-adiabatic motions of electrons that only occur in magnetic structures with a narrow scale comparable to electron inertial lengths (d e ). The recent 1D studies suggest that a plasma pinch to such scales may only occur for inflowing magnetized plasmas with relatively low plasma beta. We conclude that a parametric exploration of simulated and observed reconnection inflow conditions, particularly plasma beta, should shed light on the enablement of reconnection in collisionless plasmas.

reconnection↗
DOE OSTI2023

High-fidelity multiscale model development for accelerated fuel qualification for high discharge burnup [Poster]

Reactor vendors seek economic benefits associated with increasing nuclear fuel service lifetime in the existing light-water reactors fleet. FFRD phenomena represent a major safety concern which still needs to be addressed. Formation of high burnup structure (HBS) in conjunction with LOCA can lead to relocating fuel to escape the fuel pin and get dispersed into the primary coolant system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

On Simulating the Formation of Structured, Crystalline Systems via Non-classical Pathways

Observations in crystal growth and assembly from recent in situ methods suggest alternative, non-classical crystallization pathways play an important role in the determination of the micro- and meso- structures in crystalline systems. These processes display parallels that cross-cut multiple disciplines investigating crystallization across four orders of magnitude in size scales and widely differing environments, hinting that alternative crystal growth pathways may be a fundamental scheme in natural crystal formation. Using a system of short-range attractive microbeads, we demonstrate that the addition of a small concentration of sub-species incommensurate with the lattice spacing of the dominant species results in a stark change in crystal size and morphology. These changes are attributed to the presence of fleeting, amorphous-like configurations of beads that ultimately change the melting and growth dynamics in preferred directions. From these real-time observations, we hypothesize the amorphous mineral precursors present in biological mineralized tissues undergo similar non-classical crystallization processes resulting in the complex structures found in biomineralization.

36 MATERIALS SCIENCE↗