Viscoplastic self-consistent formulation as generalized material model for solid mechanics applications
Not Available
SEARCH · Search NASA
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.
Not Available
Abstract Recent findings suggest that ions are strongly correlated in atmospheric pressure plasmas if the ionization fraction is sufficiently high ( ≳ 10 − 5 ). A consequence is that ionization causes disorder-induced heating (DIH), which triggers a significant rise in ion temperature on a picosecond timescale. This is followed by a rise in the neutral gas temperature on a longer timescale of up to nanoseconds due to ion–neutral temperature relaxation. The sequence of DIH and ion–neutral temperature relaxation suggests a new mechanism for ultrafast neutral gas heating. Previous work considered only the case of an instantaneous ionization pulse, whereas the ionization pulse extends over nanoseconds in many experiments. Here, molecular dynamics simulations are used to analyze the evolution of ion and neutral gas temperatures for a gradual ionization over several nanoseconds. The results are compared with published experimental results from a nanosecond pulsed discharge, showing good agreement with a measurement of fast neutral gas heating.
With the recent advances in additive manufacturing (AM), long-range periodic lattice assemblies are being developed for vibration and shock mitigation components in aerospace and military applications with unique geometric and topological structures. There has been extensive work in understanding the static properties associated with varying topology of these lattice architectures, but there is almost no understanding of microstructural affects in such structures under high-strain rate dynamic loading conditions. Here we report the shock behavior of lattices with varying intrinsic grain structures achieved by post process annealing. High resolution 316L stainless steel lattices were 3D printed by a laser-powder bed fusion machine and characterized by computed tomography. Subsequent annealing resulted in stress-relieved and recrystallized lattices. Overall the lattices had strong cubic texture aligning with the x-, y- and z-directions of the build with a preference outside the build direction (z). The recrystallized sample had more equiaxed polygonal grains and a layer of BCC ferrite at the surface of the structure approximately 1 grain thick. Upon dynamic compression the as-deposited lattice showed steady compaction behavior while the heat-treated lattices exhibit negative velocity behavior indicative of failure. We attribute this to the stiffer BCC ferrite in the annealed lattices becoming damaged and fragmenting during compression.
Cathode-directed streamer evolution in near atmospheric air is modeled in 3D pin-to-plane geometries using a 3D kinetic Particle-In-Cell (PIC) code that simulates particle-particle collisions via the Direct Simulation Monte Carlo (DSMC) method. Due to the computational challenges associated with a complete 360° volumetric domain, a practical alternative was achieved using a wedge domain and a range of azimuthal angles was explored (5°, 15°, 30°, and 45°) to study possible effects on the streamer growth and propagation due to the finite wedge angle. A DC voltage of 6 kV is administered to a hemispherical anode of radius 100 μm, with a planar cathode held at ground potential, generating an over-volted state with an electric field of 4 MV/m across a 1500 μm gap. The domain is seeded with an initial ion and electron density of 10 18 m -3 at 1 eV temperature confined to a spherical region of radius 100 μm centered at the tip of the anode. The air chemistry model [1] includes standard Townsend breakdown mechanisms (electron-neutral elastic, excitation, ionization, attachment, and detachment collision chemistry and secondary electron emission) as well as streamer mechanisms (photoionization and ion-neutral collisions) via tracking excited state neutrals which can then either quench via collisions or spontaneously emit a photon based on specific Einstein-A coefficients [2, 3]. In this work, positive streamer dynamics are formally quantified for each wedge angle in terms of electron velocity and density as temporal functions of coordinates r, Φ, and z. Applying a random plasma seed for each simulation, particles of interest are tracked with near femtosecond temporal resolution out to 1.4 ns and spatially binned. This process is repeated six times and results are averaged. Prior 2D studies have shown that the reduced electric field, E/n, can significantly impact streamer evolution [4]. We extend the analysis to 3D wedge geometries, to limit computational costs, and examine the wedge angle’s effect on streamer branching, propagation, and velocity. Results indicate that the smallest wedge angle that produced an acceptably converged solution is 30°. The potential effects that a mesh, when under-resolved with respect to the Debye length, can impart on streamer dynamics and numerical heating were not investigated, and we explicitly state here that the smallest cell size was approximately 10 times the minimum λ D in the streamer channel at late times. This constraint on cell size was the result of computational limitations on total mesh count.
The Real Time Nuclear Activation Detector (RTNAD) array at NIF measures the distribution of 14 MeV neutrons emitted by deuterium-tritium (DT) fueled inertial confinement fusion implosions. The uniformity of the neutron distribution is an important indication of implosion symmetry and DT shell integrity. The array consists of 48 LaBr 3 (Ce) crystal gamma-ray spectrometers mounted outside the NIF target chamber, which continuously monitor the slow decay of the 909 keV gamma-ray line from activated 89 Zr located in Zr cups surrounding each crystal. The measured decay rate dramatically increases during a DT implosion in proportion to the number of 14 MeV neutrons striking each Zr cup. The neutrons produce activated 89 Zr through an (n, 2n) reaction on 90 Zr, which is insensitive to low energy neutrons. The neutron flux along the detector line-of-sight at shot time is determined by extrapolating the fitted 909 keV decay curve back to shot time. Automatic analysis algorithms were developed to handle the non-stop data stream. The large number of detectors and the high statistical accuracy of the array enable the spherical harmonic modes of the neutron angular distribution to be measured up to L ≤ 4 to provide a better understanding of implosion dynamics. In addition, these data combined with measurements of the down-scattered neutrons can be used to derive fuel areal density distributions. This paper will describe the RTNAD hardware and analysis procedures.
The synthesis of binary nanocrystal superlattices (BNSLs) enables the targeted integration of orthogonal physical properties, such as photoluminescence and magnetism, into a single superstructure, unlocking a vast design space for multifunctional materials. However, the formation mechanism of BNSLs remains poorly understood, restricting the prediction of the structure and properties of superlattices. Here we use a combination of in situ scattering and molecular simulation to elucidate the self-assembly of two common BNSLs (AlB 2 and NaZn 13 ) through emulsion templating. Our self-assembly experiments reveal that no intermediate structures precede the formation of the final binary phases, indicating that their formation proceeds through classical nucleation. Using simulations, we find that, despite the formation of AlB 2 and NaZn 13 typically being attributed to entropy, their self-assembly is most consistent with the nanocrystals possessing short-range interparticle attraction, which we find can accelerate nucleation kinetics in BNSLs. Here, we also find homogeneous, classical nucleation in simulations, corroborating our experiments. These results establish a robust correspondence between experiment and theory, paving the way towards prediction of BNSLs.
We present a computationally efficient time-domain model capable of simulating parametric resonances in a floating body in waves. The model assumes all wave forces to be linear, but the inertia and restoring forces acting on the body are expanded to second order in body motions. The simulation speed on a standard computer is approximately 40 times faster than real time. The model is applied to a soft-moored floating axisymmetric body which absorbs energy through heave, but is otherwise free to move in six degrees of freedom. Under certain conditions, we show that the body responds parametrically with large amplitudes not only in surge and pitch, but also in sway, roll, and yaw, provided it is given some small initial displacement in one of these out-of-plane modes. The predictions are confirmed by simulations using state-of-the-art nonlinear Froude-Krylov and computational fluid dynamics models.
Spatial-Temporal Graph Learning (ST-GL) is a prominent research area due to its unique capability to effectively learn real-world graphs. Applications of ST-GL pose stringent and various demands on not only real-time inference with low energy cost and high ac- curacy but also fast training. Unfortunately, as Moore’s Law approaches its limits and ST-GL model complexity drastically grows, the gap between digital hardware’s computational power and ST- GL application demands is widening. In response, this paper introduces Nature-GL, a nature-powered graph learning paradigm that exploits the principle of entropy increase to advance graph learning. In particular, Nature-GL transforms both the training and inference of real-valued ST-GL into electron-speed natural anneal- ing processes of a parameterized dynamical system that represents the target graphs. Experimental results across four real-world ap- plications with six datasets demonstrate that Nature-GL achieves orders-of-magnitude speedups in both training and inference, delivering higher accuracy compared to Graph Neural Networks.
Benchmarks are a cornerstone of modern machine learning, enabling reproducibility, comparison, and scientific progress. However, AI benchmarks are increasingly complex, requiring dynamic, AI-focused workflows. Rapid evolution in model architectures, scale, datasets, and deployment contexts makes evaluation a moving target. Large language models often memorize static benchmarks, causing a gap between benchmark results and real-world performance. Beyond traditional static benchmarks, continuous adaptive benchmarking frameworks are needed to align scientific assessment with deployment risks. This calls for skills and education in AI Benchmark Carpentry. From our experience with MLCommons, educational initiatives, and programs like the DOE's Trillion Parameter Consortium, key barriers include high resource demands, limited access to specialized hardware, lack of benchmark design expertise, and uncertainty in relating results to application domains. Current benchmarks often emphasize peak performance on top-tier hardware, offering limited guidance for diverse, real-world scenarios. Benchmarking must become dynamic, incorporating evolving models, updated data, and heterogeneous platforms while maintaining transparency, reproducibility, and interpretability. Democratization requires both technical innovation and systematic education across levels, building sustained expertise in benchmark design and use. Benchmarks should support application-relevant comparisons, enabling informed, context-sensitive decisions. Dynamic, inclusive benchmarking will ensure evaluation keeps pace with AI evolution and supports responsible, reproducible, and accessible AI deployment. Community efforts can provide a foundation for AI Benchmark Carpentry.
We present an ab initio free energy model derived from a fully relativistic density functional theory (DFT) electronic structure with dynamic magnetism for δ -plutonium (face-centered cubic, fcc). The DFT model is extended with orbital-orbital interaction in a parameter free orbital polarization (OP) mechanism consistent with previous modeling of plutonium. Gibbs free energy is built from components associated with the temperature dependence of the electronic structure and the corresponding electronic entropy, lattice vibrations within an anharmonic lattice dynamics model, and dynamical fluctuations of the magnetization density, i.e. magnetic fluctuations. The fluctuation model consists of transverse and longitudinal modes driven by temperature induced excitations of the DFT + OP electronic structure. The ab initio model thus incorporates fluctuating states beyond the electronic ground state. Thanks to the dynamic magnetism, the theory predicts excellent thermodynamic properties and a Gibbs free energy in accord with CALPHAD and semi-empirical modeling developed from the thermodynamic observables. The magnetic fluctuations further explain anomalous behaviors of the thermal expansion in plutonium. Specifically, a thermal expansion for the δ -plutonium system turning from positive to negative at temperatures above room temperature, a tendency for gallium to reduce and remove the negative thermal expansion depending on composition, and a positive thermal expansion for the high temperature ϵ phase.
Classical hydrodynamics is a remarkably versatile description of the coarse-grained behaviour of many-particle systems once local equilibrium has been established. The form of the hydrodynamical equations is determined primarily by the conserved quantities present in a system. Some quantum spin chains are known to possess, even in the simplest cases, a greatly expanded set of conservation laws, and recent work suggests that these laws strongly modify collective spin dynamics, even at high temperature. In this work, by probing the dynamical exponent of the one-dimensional Heisenberg antiferromagnet KCuF 3 with neutron scattering, we find evidence that the spin dynamics are well described by the dynamical exponent z = 3/2, which is consistent with the recent theoretical conjecture that the dynamics of this quantum system are described by the Kardar–Parisi–Zhang universality class. This observation shows that low-energy inelastic neutron scattering at moderate temperatures can reveal the details of emergent quantum fluid properties like those arising in non-Fermi liquids in higher dimensions.
We present the first study of the quark mass dependence of the heavy quark momentum and spatial diffusion coefficients using lattice QCD with light dynamical quarks corresponding to a pion mass of 320 MeV. We find that, for the temperature range 195 MeV < 𝑇 < 293 MeV, the spatial diffusion coefficients of the charm and bottom quarks are smaller than those obtained in phenomenological models that describe the 𝑝 𝑇 spectra and elliptic flow of open heavy flavor hadrons.
We present high-temperature thermodynamic properties for graphite from first-principles anharmonic theory. The ab initio electronic structure is obtained from density-functional theory coupled to a lattice dynamics method that is used to model anharmonic lattice vibrations. This combined approach produces free energies and specific heats for graphite that compare well with available experiments and results from models that empirically represent experimental data, such as CALPHAD. We show that anharmonic theory for the phonons is essential for accurate thermodynamic quantities above about 1000 K.
Molecular electronic devices require precise control over the flow of current in single molecules. However, the electron transport properties of single molecules critically depend on dynamic molecular conformations in nanoscale junctions. Here, in this work, we report a unique strategy for controlling molecular conductance using shape-persistent molecules. Chemically diverse, charged ladder molecules, synthesized via a one-pot multicomponent ladderization strategy, show a molecular conductance (d[log( G/G 0 )]/d x ≈ -0.1 nm -1 ) that is nearly independent of junction displacement, in stark contrast to the nanogap-dependent conductance (d[log( G/G 0 )]/d x ≈ -7 nm -1 ) observed for non-ladder analogues. Ladder molecules show an unusually narrow distribution of molecular conductance during dynamic junction displacement, which is attributed to the shape-persistent backbone and restricted rotation of terminal anchor groups. These principles are further extended to a butterfly-like molecule, thereby demonstrating the strategy's generality for achieving gap-independent conductance. Overall, our work provides important avenues for controlling molecular conductance using shape-persistent molecules. Achieving robust and controllable conductance in single-molecule junctions is challenging due to the dynamic nature of molecular conformations that fluctuate over operational timescales. A strategy using shape-persistent molecules has now been developed that demonstrates nearly junction-displacement-independent conductance, providing a stable solution for single-molecule electronic properties.
Additively manufactured (AM) metals offer the potential for customizable, cost-effective components, but qualification and certification are crucial. Key to this process is understanding pore dynamics under stress, typically analyzed using micro-computed tomography. This study introduces laboratory-scale “betatron” x-rays from laser wakefield acceleration as a high-throughput alternative for x-ray tomography of advanced materials, such as AM AlSi10Mg alloys. Coupled with 3D finite element modeling, this method provides detailed insights into stress-porosity interactions. The approach delivers high-resolution scans, revealing that pore shape and local triaxiality significantly influence fracture dynamics, supporting advanced material characterization. This work also demonstrates the potential and versatility of laser-betatron x-ray μCT for generating large datasets to accelerate our understanding of the stochastic, process-specific nature of pore formation in AM alloys.
Computing uses energy. At the bare minimum, erasing information in a computer increases the entropy. Landauer has calculated %7E k B T log(2) Joules is dissipated per bit of energy erased. While the success of Moores law has allowed increasing computing power and efficiency for many years, these improvements are coming to an end. This project asks if there is a way to continue those gains by circumventing Landauer through reversible computing. We explore a new reversible computing paradigm, asynchronous ballistic reversible computing or ABRC. The ballistic nature of data in ABRC matches well with superconductivity which provides a low-loss environment and a quantized bit encoding the fluxon. We discuss both these and our development of a superconducting fabrication process at Sandia. We describe a fully reversible 1-bit memory cell based on fluxon dynamics. Building on this model, we propose several other gates which may also offer reversible operation.
The rapidly growing industry of crop biostimulants leverages the application of plant growth promoting rhizobacteria (PGPR) to promote plant growth and health. However, introducing nonnative rhizobacteria may impact other aspects of ecosystem functioning and have legacy effects; these potential consequences are largely unexplored. Nontarget consequences of PGPR may include changes in resident microbiomes, nutrient cycling, pollinator services, functioning of other herbivores, disease suppression, and organic matter persistence. Importantly, we lack knowledge of whether these ecosystem effects may manifest in adjacent ecosystems. The introduced PGPR can leave a functional legacy whether they persist in the community or not. Legacy effects include shifts in resident microbiomes and their temporal dynamics, horizontal transfer of genes from the PGPR to resident taxa, and changes in resident functional groups and interaction networks. Ecosystem functions may be affected by legacies PGPR leave following niche construction, such as when PGPR alter soil pH that in turn alters biogeochemical cycling rates. Here, we highlight new research directions to elucidate how introduced PGPR impact resident microbiomes and ecosystem functions and their capacity for legacy effects.
Incorporating 2,3‐dihydrofuran (DHF) into dicyclopentadiene (DCPD) enables frontally polymerizable, deconstructible thermosets via acid hydrolysis. At 5–20 mol% DHF and low initiator concentration (~100 ppm), DHF inhibits front propagation, lowers the glass transition temperature ( T g ) by up to 20°C, and quenches the front, especially at 20% DHF. This paper demonstrates that increasing the initiator concentration to 250–1000 ppm overcomes inhibition, enhances front speeds by ~50%, and suppresses quenching at the expense of pot‐life. Above 5 mol% DHF, an emergence of a second exotherm is observed, attributed to increased rate dependence on the inhibition of the initiator by DHF. Void formation from volatile DHF is suppressed by applying 110 kPa nitrogen, though at the cost of reduced front speed. Dynamic mechanical analysis shows that higher initiator concentration produced more highly cross‐linked polymers, with a ~20% increase in storage modulus at 25°C and a 20°C increase in T g .