Search NASA⌕ Search

SEARCH · Search NASA

Results for “Burgers”

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

Parallel-in-Time Solution of Scalar Nonlinear Conservation Laws

Here, we consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially nonoscillatory (WENO) reconstructions, and in time with high-order explicit Runge–Kutta methods. The solution of the global, discretized space-time problem is sought via a nonlinear iteration that uses a novel linearization strategy in cases of nondifferentiable equations. Under certain choices of discretization and algorithmic parameters, the nonlinear iteration coincides with Newton’s method, although, more generally, it is a preconditioned residual correction scheme. At each nonlinear iteration, the linearized problem takes the form of a certain discretization of a linear conservation law over the space-time domain in question. An approximate parallel-in-time solution of the linearized problem is computed with a single multigrid reduction-in-time (MGRIT) iteration; however, any other effective parallel-in-time method could be used in its place. The MGRIT iteration employs a novel coarse-grid operator that is a modified conservative semi-Lagrangian discretization and generalizes those we have developed previously for nonconservative scalar linear hyperbolic problems. Numerical tests are performed for the inviscid Burgers and Buckley–Leverett equations. For many test problems, the solver converges in just a handful of iterations with a convergence rate independent of mesh resolution, including problems with (interacting) shocks and rarefactions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Physics-Informed Active Learning With Simultaneous Weak-Form Latent Space Dynamics Identification

The parametric greedy latent space dynamics identification (gLaSDI) framework has demonstrated promising potential for accurate and efficient modeling of high-dimensional nonlinear physical systems. However, it remains challenging to handle noisy data. Here, to enhance robustness against noise, we incorporate the weak-form estimation of nonlinear dynamics (WENDy) into gLaSDI. In the proposed weak-form gLaSDI (WgLaSDI) framework, an autoencoder and WENDy are trained simultaneously to discover intrinsic nonlinear latent-space dynamics of high-dimensional data. Compared with the standard sparse identification of nonlinear dynamics (SINDy) employed in gLaSDI, WENDy enables variance reduction and robust latent space discovery, therefore leading to more accurate and efficient reduced-order modeling. Furthermore, the greedy physics-informed active learning in WgLaSDI enables adaptive sampling of optimal training data on the fly for enhanced modeling accuracy. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including viscous and inviscid Burgers' equations, time-dependent radial advection, and the Vlasov equation for plasma physics. With data that contains 5%–10% Gaussian white noise, WgLaSDI outperforms gLaSDI by orders of magnitude, achieving 1%–7% relative errors. Compared with the high-fidelity models, WgLaSDI achieves 121 to 1779x speed-up.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Simulation of the Impact of Point Defects and Edge Dislocations on X-Ray Diffraction in Hexagonal (Ni,Co) 1+2 x Ti 1– x O 3 Thin Films

In this work, a computer code for simulating high-resolution X-ray diffraction (XRD) data from disordered crystals with arbitrary spatial composition and local lattice parameters is developed. Simulated patterns are compared with the experimental data collected on a single phase, highly crystalline (Ni 0.42 Co 0.58 ) 2.22 Ti 0.39 O 3 solid-solution thin film exhibiting a large number of subsidiary minima. As a case study, the edge dislocations in hexagonal (Ni,Co) 3 O 3 thin films with Burgers vector parallel and perpendicular to the a-axis and c-axis, respectively, are modeled. No peak profiles are assumed and thus issues related to profile fitting are avoided. Both macroscopic features, such as film thickness, and atomic-scale structure, such as dislocations, are simultaneously modeled. Simulations are run on a desktop machine. Commonly applied database-based phase identification routines can result in wrong phase identification, unless intensities are properly modeled. Modeling tools of diffractometer manufacturers are compared with the present approach. An example of how simulated reciprocal lattice patterns can be used to choose relevant measurement geometries in defected thin films is given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

T$ \overline{T} $ deformation in SCFTs and integrable supersymmetric theories

We calculate the \( \mathcal{S} \) -multiplets for two-dimensional Euclidean \( \mathcal{N} \) = (0 , 2) and \( \mathcal{N} \) = (2 , 2) superconformal field theories under the T \( \overline{T} \) deformation at leading order of perturbation theory in the deformation coupling. Then, from these \( \mathcal{N} \) = (0 , 2) deformed multiplets, we calculate two- and three-point correlators. We show the \( \mathcal{N} \) = (0 , 2) chiral ring’s elements do not flow under the T \( \overline{T} \) deformation. Specializing to integrable supersymmetric seed theories, such as \( \mathcal{N} \) = (2 , 2) Landau-Ginzburg models, we use the thermodynamic Bethe ansatz to study the S-matrices and ground state energies. From both an S-matrix perspective and Melzer’s folding prescription, we show that the deformed ground state energy obeys the inviscid Burgers’ equation. Finally, we show that several indices independent of D -term perturbations including the Witten index, Cecotti-Fendley-Intriligator-Vafa index and elliptic genus do not flow under the T \( \overline{T} \) deformation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗
DOE OSTI2024

A High-Order Eulerian–Lagrangian Runge–Kutta Finite Volume (EL–RK–FV) Method for Scalar Nonlinear Conservation Laws

Abstract We present a class of high-order Eulerian–Lagrangian Runge–Kutta finite volume methods that can numerically solve Burgers’ equation with shock formations, which could be extended to general scalar conservation laws. Eulerian–Lagrangian (EL) and semi-Lagrangian (SL) methods have recently seen increased development and have become a staple for allowing large time-stepping sizes. Yet, maintaining relatively large time-stepping sizes post shock formation remains quite challenging. Our proposed scheme integrates the partial differential equation on a space-time region partitioned by linear approximations to the characteristics determined by the Rankine–Hugoniot jump condition. We trace the characteristics forward in time and present a merging procedure for the mesh cells to handle intersecting characteristics due to shocks. Following this partitioning, we write the equation in a time-differential form and evolve with Runge–Kutta methods in a method-of-lines fashion. High-resolution methods such as ENO and WENO-AO schemes are used for spatial reconstruction. Extension to higher dimensions is done via dimensional splitting. Numerical experiments demonstrate our scheme’s high-order accuracy and ability to sharply capture post-shock solutions with large time-stepping sizes.

Chen, Jiajie↗
DOE OSTI2025

Randomized Preconditioned Solvers for Strong Constraint 4D-Var Data Assimilation

The Strong Constraint 4D Variational (SC-4DVAR) data assimilation method is widely used in climate and weather applications. SC-4DVAR involves solving a minimization problem to compute the maximum a posteriori estimate, which we tackle using the Gauss-Newton method. The computation of the descent direction is expensive since it involves the solution of a large-scale and potentially ill-conditioned linear system, solved using the preconditioned conjugate gradient (PCG) method. Here, to address this cost, we efficiently construct scalable preconditioners using three different randomization techniques, which all rely on a certain low-rank structure involving the Gauss-Newton Hessian. The proposed techniques come with theoretical guarantees on the condition number, and at the same time, are amenable to parallelization. We also develop an adaptive approach to estimate the sketch size and choose between the reuse or recomputation of the preconditioner. We demonstrate the performance and effectiveness of our methodology on two representative model problems—the Burgers and barotropic vorticity equation—showing a drastic reduction in both the number of PCG iterations and the number of Gauss-Newton Hessian products after including the preconditioner construction cost.

Gauss-Newton↗
DOE OSTI2020

Revealing the Subsurface Basal $\langle$a$\rangle$ Dislocation Activity in Magnesium Through Lattice Rotation Analysis

A method was proposed in this study to reveal the subsurface basal dislocation activity in Mg-Y alloy and determine the corresponding Burgers vector. This is achieved by correlating the slip directions of dislocations to the lattice rotation represented by the {0001} pole figure. The identified basal slip system by this approach was verified by micro-Laue diffraction. This method can be applied as a complementary method to the conventional slip trace analysis to study the dislocation behavior of Mg alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

In Situ Micro-Pillar Compression to Examine Radiation-Induced Hardening Mechanisms of FeCrAl Alloys

Here, the effects of 5 MeV Fe 2+ ion irradiation at 300°C on the microstructure evolution and deformation behavior of a FeCrAl C26M alloy are presented. It has been found that dislocation loop density increases an order of magnitude from 1 dpa to 16 dpa irradiations, whereas, the dislocation loop size saturates with increasing damage. Micropillars, 600 nm in diameter and 1.3 µm in height, were fabricated and compressed inside grains with <001>, <011> and <111> crystallographic orientations, respectively. {112} <111> has been identified as the primary slip system in both unirradiated and irradiated alloy. The increase in yield stress after irradiation is observed with measurable variation along <001> and <011> vs. along <111>. By applying the Orowan dispersed barrier model, the increase of yield stress is found mainly due to the slip resistance of radiation generated defect loops. Detailed transmission electron microscopy (TEM) studies were performed to quantify the Burgers vector and the distribution of irradiation induced dislocations at elevated strains. It is revealed that localized shear instability is caused by avalanche slip events of ½<111> dislocations gliding out of tested pillars. Simultaneously, a large number of sessile/immobile <100> dislocations formed in the vicinity of slip band, leading to the hardening at elevated strains.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamic observation of dual-beam irradiated Fe and Fe-10Cr alloys at 435 °C

The attractive mechanical properties and superior resistance to void-swelling make ferritic/martensitic alloys a promising structural material for advanced nuclear reactors. However, one anomaly that has intrigued researchers for more than 50 years is the proportion of two types of dislocation loops in Fe and Fe-Cr alloys with Burger vectors b=½<111> and b=<100>. Although the possible mechanisms responsible for the presence of <100> loops continue to be the subject of intense modeling studies, there remains incomplete experimental understanding of fundamental irradiation processes in Fe(Cr) alloys. Here, the dose dependence of the irradiation-induced microstructural evolution was examined from 0 to 20 displacement per atom (dpa) in high purity Fe and Fe-10Cr during simultaneous dual-beam (1 MeV Kr + 10 appm He/dpa) irradiation at 435 °C. We experimentally revealed that the mechanism for the formation of <100> loops may not follow the conventional simple dislocation reaction between two 1½<111> loops. Real-time dynamic formation and evolution of defects including black dot loops, loop coarsening, loop decoration, network dislocations, and cavities were demonstrated. Several results indicated that the addition of Cr and He could impede dislocation loop motion. The evolution of the defect size/density and relative fraction of ½<111> vs <100> loops were quantitatively summarized. With increasing dose, ½<111> loops became the dominant type of loop in both materials. Notably, <100> loops were predominantly observed near grain boundaries only for pure Fe, while arrays of nanoscale black dot defects composing the <100> loop strings were observed in plenty in Fe-10Cr.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding and quantifying electron beam effects during in situ TEM nanomechanical tensile testing on metal thin films

Transmission electron microscopy (TEM) imaging relies on high energy electrons for atomic scale resolution, however, the electrons themselves interact with and may alter the material being imaged. Using an in situ TEM MEMS-based nanomechanical testing technique, the effect of the electron beam (e-beam) on the deformation behavior of nanocrystalline Al and ultrafine-grained Au is investigated and quantified. We show that the e-beam enhances plastic deformation, leading to an increase in plastic strain rate and a decrease in true activation volume V* in Al (28 to 21b 3 , with b being the Burgers vector length). The e-beam has a much weaker effect on Au. The e-beam effect is not caused by knock-on damage, but rather an effective temperature increase due to additional atomic fluctuations provided by the e-beam. The effective temperature increase is larger for Al than Au. This e-beam effect does not change the deformation mechanisms, but instead accelerates the stress-driven, thermally activated plastic deformation. These experiments provide insight into the effects of the e-beam on plastic deformation in different metals and underscore the importance of understanding and quantifying these effects for proper interpretation of measured mechanical properties during in situ TEM experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomistic modeling of interface strengthening in Al-Si eutectic alloys

Al-Si cast alloys are usually composed of α-Al and Al-Si eutectic. Si flakes and Al matrix generally hold cube-on-cube orientation relationship with the primary interface (111) Al ∥(111) Si . Extensive experimental studies demonstrated that Si flakes cannot significantly improve mechanical properties of Al-Si cast alloys. We hypothesize that the weak strengthening effect associated with Si flakes might be attributed to thermomechanical properties of Al-Si interfaces besides their morphologies. To characterize Al-Si interfaces with a large lattice mismatch (> 30%), we proposed the quasi-coincident site lattice (Q-CSL) as reference lattice, and demonstrated that the Q-CSL Al-Si coherent interface has three characteristic coherent structures, one stable and low energy structure and two metastable and high energy structures. The translation vectors for the same type of coherent Q-CSL structures are consistent with three displacement shift complete (DSC) vectors. The two metastable structures can be obtained by shifting the low energy structure with three partial DSC vectors. Semi-coherent interface is composed of the low energy Q-CSL patches and three sets of interface misfit dislocations with Burgers vectors same as the DSC vectors. Atomistic simulations revealed that Al-Si interface exhibits low shear resistance. Ideal shear strength of the Q-CSL coherent interface is 110 MPa and semi-coherent interface is 20 MPa. The low shear resistance is attributed to the glide of interface misfit dislocations. Al-Si interface also exhibits low formation and migration energies of point defects. Owing to low shear strength and low formation and migration energies of point defects, interface sliding or shear readily happen under mechanical loading or during dislocation-interface interactions. Lattice dislocations can cross slip onto or climb along Al-Si interfaces. These reactions decrease the number of accumulated dislocation loops around Si flakes and promote nucleation and emission of lattice dislocations from Al-Si interfaces to matrix, consequently reduce the repulsive force on approaching dislocations and weaken Si flakes strengthening effect. In situ tension and compression tests in a scanning electron microscope reveal relatively weak strengthening effect due to Si flakes, consistent with the computed dislocation interaction with interfaces and shear behavior of interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kink mechanism in Cu/Nb nanolaminates explored by $\mathcal{in}$ $\mathcal{situ}$ pillar compression

We report Nano metallic laminates (NMLs) exhibit different failure modes depending on the loading conditions due to their mechanical anisotropies. Kinking is a typical failure mode in many NMLs compressed along a layer-parallel direction. However, a detailed description of the microstructure evolution during kink band (KB) formation and an in-depth understanding of the formation mechanisms are lacking. In this work, the KB process is investigated in Cu/Nb NMLs by in situ micro pillar compression in the scanning electron microscope (SEM) along a layer-parallel direction. Post-mortem S/TEM and transmission Kikuchi diffraction (TKD) analyses show that kink banding leads to significant microstructure changes characterized by an accumulation of geometrically necessary dislocations (GNDs) and of tilt geometrically necessary boundaries (GNBs) near KB boundaries (KBBs). The distinct microstructure evolution implies that KB formation is facilitated by the inhomogeneous microstructures resulting in constrained deformation modes. Specifically, dislocations active on slip planes nearly parallel to the interfaces make a major contribution to kink evolution after the onset of kinking. Once layer-parallel slip systems are activated, preexisting lattice dislocations and dislocations nucleating from interfaces will accumulate as GNDs near KBBs via the stochastic storage of lattice dislocations that have certain Burgers vectors. GNDs can further transform into GNBs via cross-slip and climb driven processes near the KBB. Furthermore, GNBs near KBBs can grow by incorporating more GNDs or by coalescence to accommodate the KB evolution. We further hypothesize that microstructural perturbations and their ensuing stresses can initiate KB formation in Cu/Nb NMLs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic modeling of surface and grain boundary dislocation nucleation in FCC metals

Dislocation nucleation plays a critical role in the plastic deformation of crystalline materials. However, it is challenging to predict the active mode and associated rate of dislocation nucleation under typical experimental loading conditions through molecular dynamics simulation due to timescale limitations. In this report we use the free-end nudged elastic band method to determine the activation energies and activation volumes of dislocation nucleation in four typical face-centered cubic metals of Au, Al, Cu and Ni. We focus on the representative processes of surface and grain boundary dislocation nucleation. The atomistically determined activation volumes of these dislocation nucleation processes are larger than 10b3 (with b being the Burgers vector length) under typical experimental loading conditions. These results are compared with experimentally measured activation volumes in ultrafine-grained and nanocrystalline metals, thereby providing mechanistic insight into their rate-controlling deformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The complexity of { 11 2 ¯ 1 } twin interface structure and energetics in HCP materials

Here, the quest for a deeper understanding of the mechanical properties of HCP alloys has been an ongoing effort. Due to the lack of a sufficient number of slip systems, twin nucleation and migration originate as a substantial plastic deformation mechanism in HCP metals. The paper uncovers the energetics and atomistic mechanism of the {$11\overline{2}1$} < $11\overline{2}6$ > which is one of the most prominent twin modes in HCP materials that has not been fully understood. In view of corrugated {$11\overline{2}1$} interfaces, the positioning of lattice/motif atoms needs to be clarified. First, we establish the initial twin boundary (TB) structure unambiguously via the lattice offsets resulting in volume invariancy and correct relative positioning of the twin variants. Secondly, the TB migration mechanism is determined from a crystallographic (zero net shuffle, $Σ\vec{s}$ = 0) and energetic standpoints (minimum energy path). The Generalized Planar Fault Energy (GPFE) curves are calculated for Ti, Hf, Mg, and Zr. Such GPFE curves for this twin mode are calculated for the first time, and point to unusually high unstable energies (300-850 mJ/m 2 ), and low migration barriers (<60 mJ/m 2 ) for twinning accompanied with large shuffle magnitudes, |$\vec{s}$|, comparable to the Burger’s vector (|$\vec{s}$| ≅ $b_e$ > 0.44Å)

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal plasticity model of BCC metals from large-scale MD simulations

Accurate crystal plasticity models that faithfully capture the behavior of single crystals under a wide range of loading conditions, such as loading direction, strain rate, and temperature, are still lacking. Here we introduce a novel approach in which a crystal plasticity (CP) model is informed directly from and calibrated to large-scale quantum-accurate MD simulations in which single crystal BCC Ta serves as a testbed material. By analyzing our large set of MD simulations several key insights are obtained leading us to modify constitutive assumptions in order to address deficiencies of existing CP models. Importantly, we observe that the standard notion of fixed slip systems – pairs of slip directions and slip planes – is inadequate for describing high-rate plasticity in BCC tantalum at room temperature. Instead, pencil glide defined as dislocation motion in the maximum resolved shear stress planes (MRSSP) of each Burgers vector is fully consistent with our MD simulation data while providing significant simplifications of the constitutive relations. Herein, our resulting new CP model closely matches the behavior of single crystals observed in high-rate MD simulations while being fully consistent with lower rate experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-pressure phase transitions in a laser directed energy deposited Fe-33Cu Alloy

Here, additively manufactured Fe-Cu alloys contain both equilibrium face-centered cubic (FCC) and metastable body-centered cubic (BCC) crystal structure Cu precipitates depending on their size. However, the stability of these nanoscale precipitates under extreme conditions such as high pressures has not been reported. This study investigates the phase transformations and microstructural stability of laser directed energy deposition (DED-LB) made Fe67Cu33 alloy (nominal composition in at.%) under high static pressure deformation using in-situ synchrotron X-ray diffraction under pressure, postmortem high-resolution scanning transmission electron microscopy (HR-STEM), and molecular dynamics (MD) simulations. In-situ XRD results reveal a reversible phase transformation from BCC to the hexagonal close-packed (HCP) structure in the Fe grains at an onset pressure of 16.4 GPa, significantly higher than reported in the literature for pure Fe. Although no high-pressure phase transition was observed in the FCC Cu grains through XRD, HR-STEM analysis uncovers a phase transition to the HCP structure in nanoscale metastable BCC Cu precipitates within the BCC Fe matrix. After decompression, the Fe matrix reverted back to the BCC structure with periodic lath martensite, while regions of the nanoscale BCC Cu precipitates retained the metastable HCP structure. MD simulations support the BCC → HCP transition in the embedded nanoscale coherent Cu precipitates, consistent with the classical Burgers mechanism. Thus, by leveraging in-situ XRD observation, postmortem high-resolution S/TEM microscopy, and MD simulations, this study offers profound insights into the distinctive phase transformations induced by high pressure in DED-LB Fe67Cu33 alloy, which is distinguished by its hierarchical microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Distorted dislocation cores and asymmetric glide resistances in titanium

The determination of Critical Resolved Shear Stress (CRSS) in titanium for basal, prismatic, and pyramidal slip-planes without empirical constants is presented by combining Density Functional Theory (DFT) and anisotropic elasticity. A new mechanism leading to tension-compression (T-C) asymmetry in the CRSS levels has been revealed for the first time. The conditions for this asymmetry are established, involving a complex interplay between the dislocation core-structure and core-advance behavior. The three conditions for T-C asymmetry that need to be simultaneously satisfied can be summarized as: (1) a medium stacking fault width, d, (3 < d/b F < 10, where b F is the magnitude of the Burgers vector of the full dislocation), (2) an asymmetric core-structure of the extended dislocation (ξ 1 ≠ ξ 2 , where ξ 1 and ξ 2 are the core-widths of the first and second partials, respectively), and (3) intermittent motion of the partials (Δd/b F ≠ 0, where Δd is the magnitude of fluctuation in stacking fault width during intermittent motion). Pyramidal-slip in titanium satisfies all three conditions, resulting in significant T-C asymmetry. The CRSS theory considers a Wigner-Seitz (W-S) cell based domain area assigned to each lattice site for the calculations of core-energies accurately capturing the slip-plane lattice. The W-S based approach is essential due to the lower symmetry of the HCP crystal circumventing potential errors associated with the one-dimensional atomic-row approximation. Dislocation core structures are obtained for all the slip-systems in titanium showing significant distortion of the disregistry profile governing the core-advance behavior. The CRSS values predicted from the theory show agreement with the experimental CRSS levels reported in the literature.

Dislocation core↗