Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Controlling quantum many-body dynamics in driven Rydberg atom arrays

Dynamic stabilization of an array Large-scale systems comprising one-dimensional chains and two-dimensional arrays of excited atoms held in a programmable optical lattice are a powerful platform with which to simulate emergent phenomena. Bluvstein et al. built an array of up to 200 Rydberg atoms and subjected the system to periodic excitation. Under such driven excitation, they found that the array of atoms stabilized, freezing periodically into what looked like time crystals. Understanding and controlling the dynamic interactions in quantum many-body systems lies at the heart of contemporary condensed matter physics and the exotic phenomena that can occur. Science , this issue p. 1355

Science & Technology - Other Topics↗
DOE OSTI2020

Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport

The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗
DOE OSTI2021

A field–particle correlation analysis of a perpendicular magnetized collisionless shock

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tailoring plasticity mechanisms in compositionally graded hierarchical steels fabricated using additive manufacturing

Abstract While there exists in nature abundant examples of materials with site-specific gradients in microstructures and properties, engineers and designers have traditionally used monolithic materials with discrete properties. Now, however, additive manufacturing (AM) offers the possibility of creating structures that mimic some aspects of nature. One example that has attracted attention in the recent years is the hierarchical structure in bamboo. The hierarchical architecture in bamboo is characterized by spatial gradients in properties and microstructures and is well suited to accommodate and survive complex stress states, severe mechanical forces, and large deformations. While AM has been used routinely to fabricate functionally graded materials, this study distinguishes itself by leveraging AM and physical metallurgy concepts to trigger cascading deformation in a single sample. Specifically, we have been successful in using AM to fabricate steel with unique spatial hierarchies in structure and property to emulate the structure and deformation mechanisms in natural materials. This study shows an improvement in the strength and ductility of the nature-inspired “hierarchical steel” compared with conventional cast stainless steels. In situ characterization proves that this improvement is due to the sequential activation of multiple deformation mechanisms namely twinning, transformation-induced plasticity, and dislocation-based plasticity. While significantly higher strengths can be achieved by refining the chemical and processing technique, this study sets the stage to achieve the paradigm of using AM to fabricate structures which emulate the flexibility in mechanical properties of natural materials and are able to adapt to in-service conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A framework for discrete optimization of stellarator coils

Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Operator dynamics in Floquet many-body systems

We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading. Published by the American Physical Society 2025

Yoshimura, Takato (ORCID:0000000309159846)↗
DOE OSTI2023

Frequency-bin photonic quantum information

Discrete frequency modes, or bins, present a blend of opportunities and challenges for photonic quantum information processing. Frequency-bin-encoded photons are readily generated by integrated quantum light sources, naturally high-dimensional, stable in optical fiber, and massively parallelizable in a single spatial mode. Yet quantum operations on frequency-bin states require coherent and controllable multifrequency interference, making them significantly more challenging to manipulate than more traditional spatial degrees of freedom. In this mini-review, we describe recent developments that have transformed these challenges and propelled frequency bins forward. Focusing on sources, manipulation schemes, and detection approaches, we introduce the basics of frequency-bin encoding, summarize the state of the art, and speculate on the field’s next phases. Given the combined progress in integrated photonics, high-fidelity quantum gates, and proof-of-principle demonstrations, frequency-bin quantum information is poised to emerge from the lab and leave its mark on practical quantum information processing—particularly in networking where frequency bins offer unique tools for multiplexing, interconnects, and high-dimensional communications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Tracer-field cross-correlations with k -nearest neighbour distributions

ABSTRACT In astronomy and cosmology significant effort is devoted to characterizing and understanding spatial cross-correlations between points – e.g galaxy positions, high energy neutrino arrival directions, X-ray and AGN sources, and continuous fields – e.g. weak lensing meiand Cosmic Microwave Background maps. Recently, we introduced the k-nearest neighbour (kNN) formalism to better characterize the clustering of discrete (point) data sets. Here, we extend it to the point – field cross-correlations analysis. It combines kNN measurements of the point data set with measurements of the field smoothed at many scales. The resulting statistics are sensitive to all orders in the joint clustering of the points and the field. We demonstrate that this approach, unlike the 2-pt cross-correlation, can measure the statistical dependence of two data sets even when there are no linear (Gaussian) correlations between them. We further demonstrate that this framework is far more effective than the two point function in detecting cross-correlations when the continuous field is contaminated by high levels of noise. For a particularly high level of noise, the cross-correlation between haloes and the underlying matter field in a cosmological simulation, between 10 h−1 Mpc and 30 h−1 Mpc, is detected at >5σ significance using the technique presented here, when the two-point cross-correlation significance is ∼1σ. Finally, we show that kNN cross-correlations of haloes and the matter field can be well modelled on quasi-linear scales using the Hybrid Effective Field Theory (HEFT) framework, with the same set of bias parameters as are used for 2-pt cross-correlations. The substantial improvement in the statistical power of detecting cross-correlations using this method makes it a promising tool for various cosmological applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Rapid prototyping of arbitrary 2D and 3D wireframe DNA origami

Wireframe DNA origami assemblies can now be programmed automatically from the top-down using simple wireframe target geometries, or meshes, in 2D and 3D, using either rigid, six-helix bundle (6HB) or more compliant, two-helix bundle (DX) edges. While these assemblies have numerous applications in nanoscale materials fabrication due to their nanoscale spatial addressability and high degree of customization, no easy-to-use graphical user interface software yet exists to deploy these algorithmic approaches within a single, standalone interface. Further, top-down sequence design of 3D DX-based objects previously enabled by DAEDALUS was limited to discrete edge lengths and uniform vertex angles, limiting the scope of objects that can be designed. Here, we introduce the open-source software package ATHENA with a graphical user interface that automatically renders single-stranded DNA scaffold routing and staple strand sequences for any target wireframe DNA origami using DX or 6HB edges, including irregular, asymmetric DX-based polyhedra with variable edge lengths and vertices demonstrated experimentally, which significantly expands the set of possible 3D DNA-based assemblies that can be designed. ATHENA also enables external editing of sequences using caDNAno, demonstrated using asymmetric nanoscale positioning of gold nanoparticles, as well as providing atomic-level models for molecular dynamics, coarse-grained dynamics with oxDNA, and other computational chemistry simulation approaches.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Can strong substorm-associated MeV electron injections be an important cause of large radiation belt enhancements?

It has become well-established that strong outer radiation belt enhancements are due to wave-driven electron energization by whistler-mode chorus waves. However, in this study, we examine strong MeV electron injections on 10 July 2019 and find substantial evidence that such injections may be a crucial contributor to outer radiation belt enhancement events. For such an examination, it is essential to precisely separate temporal flux changes from spatial variations observed as Van Allen Probes move along their orbits. Employing a new “hourly snapshot” analysis approach, we discover unprecedented details of electron flux evolutions that suggest that for this event, the outer belt enhancement was not continuous but instead intermittent, mostly composed of 4 large discrete injection-driven flux increases. The injections appear as sharp flux increases when observed near apogee. Otherwise, by comparing hourly snapshots for different times, we infer injections and infer temporally stable fluxes between injections, despite strong and continuous chorus emission. The fast and intermittent electron flux growth successively extending earthwards implies cumulative outer belt enhancement via a series of repetitive inward transport associated with injection-induced electric fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Evaluation of riverbed magnetic susceptibility for mapping biogeochemical hot spots in groundwater‐impacted rivers

Abstract Redox hot spots occurring as metal‐rich anoxic groundwater discharges through oxic wetland and river sediments commonly result in the formation of iron (Fe) oxide precipitates. These redox‐sensitive precipitates influence the release of nutrients and metals to surface water and can act as ‘contaminant sponges’ by absorbing toxic compounds. We explore the feasibility of a non‐invasive, high‐resolution magnetic susceptibility (MS) technique to efficiently map the spatial variations of magnetic Fe oxide precipitates in the shallow bed of three rivers impacted by anoxic groundwater discharge. Laboratory analyses on Mashpee River (MA, USA) sediments demonstrate the sensitivity of MS to sediment Fe concentrations. Field surveys in the Mashpee and Quashnet rivers (MA, USA) reveal several discrete high MS zones, which are associated with likely anoxic groundwater discharge as evaluated by riverbed temperature, vertical head gradient, and groundwater chemistry measurements. In the East River (CO, USA), widespread cobbles/rocks exhibit high background MS from geological ferrimagnetic minerals, thereby obscuring the relatively small enhancement of MS from groundwater induced Fe oxide precipitates. Our study suggests that, in settings with low geological sources of magnetic minerals such as lowland rivers and wetlands, MS may serve as a complementary tool to temperature methods for efficiently mapping Fe oxide accumulation zones due to anoxic groundwater discharges that may function as biogeochemical hot spots and water quality control points in gaining systems.

Wang, Chen↗
DOE OSTI2022

Numerical Analysis of Combustion Dynamics in a Full-Scale Rotating Detonation Rocket Engine using Large Eddy Simulations

Large eddy simulations (LES) using detailed chemistry and leveraging adaptive mesh refinement (AMR) are performed to gain insights into the combustion dynamics within a full-scale methane-oxygen non-premixed rotating detonation rocket engine (RDRE) employing impinging discrete injection schemes. In particular, a comparative analysis of two operating conditions corresponding to the same global equivalence ratio but different mass flow rates is carried out to investigate the resultant impact on detonation wave characteristics and RDRE global performance. Multiple co-rotating detonation waves with spatially-distributed wave structure and preferential alignment with the inner wall of the annulus (due to asymmetry in fuel distribution) are encountered under both conditions. Both cases exhibit pre-detonation deflagrative burning in the fill region, while one of the cases shows higher susceptibility to backflow into the feed plenums due to lower plenum pressures. Furthermore, heat release analysis shows that the thrust obtained from the RDRE is closely linked to the distribution of total heat release between detonative and deflagrative combustion. On the other hand, combustion efficiency is associated with the fraction of heat release occurring in fuel-rich versus fuel-lean regions within the RDRE.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Gradient-based surface nuclear magnetic resonance for groundwater investigation

In medical magnetic resonance imaging, spatial localization (imaging) is based upon the application of controlled magnetic field gradients on top of the main magnetic field to spatially modulate the frequency and/or phase of the nuclear magnetic resonance (NMR) signal across the volume of investigation. In this work, we have applied similar physical principles to produce controlled magnetic field gradients during surface NMR-based groundwater investigations. In this approach, a gradient pulse of variable amplitude or duration is applied immediately after the excitation pulse to cause predictable phase encoding of the NMR signal as a function of depth. This approach is also applicable to emerging surface NMR detection methods that use a prepolarization field with fast nonadiabatic turn-off to generate detectable NMR signals from the shallow subsurface. In this case, the gradient pulse is applied after terminating the prepolarization field and provides a heretofore unavailable means of localizing the NMR response as a function of depth. The application of gradients can also be combined with tip-angle-based modulation to yield higher imaging resolution than can be achieved through either gradient- or tip-angle-based imaging alone. We implemented this new gradient-based capability into a surface NMR gradient generation accessory that is compatible with the GMR-Flex instrument and developed surface NMR-specific forward modeling and linear inverse models. We validated the accuracy of this novel gradient-based sNMR technology using computer simulations, experiments using a small pool filled with a discrete layer of bulk water, and field experiments at well-characterized groundwater test sites along Ebey Island, WA, and Larned, KS. The gradient-based sNMR imaging observations were compared with high-resolution direct push NMR results observed at these sites. The results of computer simulations and field experiments indicate improvements in both the detection (signal-to-noise ratio) and spatial resolution of shallow subsurface water content using gradient-based surface NMR, compared with traditional surface NMR imaging methods.

Geochemistry & Geophysics↗
DOE OSTI2025

A Particle-in-cell Method for Plasmas with A Generalized Momentum Formulation, Part III: A family of Gauge Conserving Methods

In this paper, we introduce a new family of spatially co-located field solvers for particle-in-cell applications which evolve the potential formulation of Maxwell’s equations under the Lorenz gauge. Our recent work [2] introduced the concept of time-consistency, which connects charge conservation to the preservation of the gauge at the semi-discrete level. It will be shown that there exists a large family of time discretizations which satisfy this property. Additionally, it will be further shown that for large classes of time marching methods, the satisfaction of the gauge condition automatically implies the satisfaction of Gauss’s law for electricity, with the potential formulation ensuring that that Gauss’s law for magnetism is satisfied by definition. We focus on popular time marching methods including centered differences, backward differences, and diagonally-implicit Runge-Kutta methods, which are coupled to a spectral discretization in space. We demonstrate the theory by testing the methods on a relativistic Weibel instability and a drifting cloud of electrons.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Exploring the MeV sky with a combined coded mask and Compton telescope: the Galactic Explorer with a Coded aperture mask Compton telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds upon the heritage of past and current missions improving the sensitivity and, very importantly, the angular resolution. This consists in combining a Compton telescope and a coded-mask telescope. We delineate the motivation for such a concept and we define the scientific goals for such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which, in turn, enables identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Three-state nematicity in the triangular lattice antiferromagnet Fe 1/3 NbS 2

Nematic order is the breaking of rotational symmetry in the presence of translational invariance. While originally defined in the context of liquid crystals, the concept of nematic order has arisen in crystalline matter with discrete rotational symmetry, most prominently in the tetragonal Fe-based superconductors where the parent state is four-fold symmetric. In this case the nematic director takes on only two directions, and the order parameter in such 'Ising-nematic' systems is a simple scalar. In this study, using a spatially resolved optical polarimetry technique, we show that a qualitatively distinct nematic state arises in the triangular lattice antiferromagnet Fe 1/3 NbS 2 . The crucial difference is that the nematic order on the triangular lattice is a Z 3 or three-state Potts-nematic order parameter. As a consequence, the anisotropy axes of response functions such as the resistivity tensor can be continuously reoriented by external perturbations. This discovery lays the groundwork for devices that exploit analogies with nematic liquid crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗