Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “exponential integrators”

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 307 records · Page 17

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗

Simulation Propulsion System and Trajectory Optimization

A number of new aircraft concepts have recently been proposed which tightly couple the propulsion system design and operation with the overall vehicle design and performance characteristics. These concepts include propulsion technology such as boundary layer ingestion, hybrid electric propulsion systems, distributed propulsion systems and variable cycle engines. Initial studies examining these concepts have typically used a traditional decoupled approach to aircraft design where the aerodynamics and propulsion designs are done a-priori and tabular data is used to provide inexpensive look ups to the trajectory ana-ysis. However the cost of generating the tabular data begins to grow exponentially when newer aircraft concepts require consideration of additional operational parameters such as multiple throttle settings, angle-of-attack effects on the propulsion system, or propulsion throttle setting effects on aerodynamics. This paper proposes a new modeling approach that eliminated the need to generate tabular data, instead allowing an expensive propulsion or aerodynamic analysis to be directly integrated into the trajectory analysis model and the entire design problem optimized in a fully coupled manner. The new method is demonstrated by implementing a canonical optimal control problem, the F-4 minimum time-to-climb trajectory optimization using three relatively new analysis tools: Open M-DAO, PyCycle and Pointer. Pycycle and Pointer both provide analytic derivatives and Open MDAO enables the two tools to be combined into a coupled model that can be run in an efficient parallel manner that helps to cost the increased cost of the more expensive propulsion analysis. Results generated with this model serve as a validation of the tightly coupled design method and guide future studies to examine aircraft concepts with more complex operational dependencies for the aerodynamic and propulsion models.

Optimization↗

New Pulsed Neutron Die-Away Experiments in Light Water

Thermal neutron scattering laws are important data for many nuclear science and engineering applications, especially criticality safety. Recently, pulsed-neutron die-away experiments have been proposed and used as an experiment to validate thermal neutron scattering laws. These experiments involve irradiating a target moderating material with pulses of neutrons from a neutron generator. The physics of thermal scattering greatly affect how the neutron population in the target exponentially decays via absorption and leakage. Herein, we present experimental results from a new pulsed-neutron die-away experiment of light water that was performed at Lawrence Livermore National Laboratory. The experiments were done with cylindrical targets of varying dimensions to modulate the experiment’s sensitivity to thermal scattering and to absorption. We compare the measured integral parameter of the experiment to simulations with MCNP6.2 ® and to past experiments in literature. The integral parameters that were more sensitive to TSLs showed a larger bias. This validation study is known to have errors because the model of the experiment lacks key features about the detectors and neutron generator. These features will be included in the final benchmark evaluation of the experiment which will be submitted to the International Criticality Safety Benchmark Evaluation Project.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Numerical solutions of Navier-Stokes equations for the structure of a trailing vortex

The structure and decay of a trailing vortex were analyzed during the numerical solutions of the full Navier-Stokes equations. Unsteady forms of the governing equations were recast in terms of circulation, vorticity, and stream function as dependent variables, and a second upwind finite difference scheme was used to integrate them with prescribed initial and boundary conditions. The boundary conditions at the outer edge and at the outflow section of the trailing vortex were considered. Different models of the flow were postulated, and solutions were obtained describing the development of the flow as integration proceeds in time. A parametric study was undertaken with a view to understanding the various phenomena that may possibly occur in the trailing vortex. Using the Hoffman and Joubert law of circulation at the inflow section, the results of this investigation were compared with experimental data for a Convair 990 wind model and a rectangular wing. With an exponentially decaying law of circulation at the inflow section and an adverse pressure gradient at the outer edge of the trailing vortex, solutions depict vortex bursting through the sudden expansion of the core and/or through the stagnation and consequent reversal of the flow on the axis. It was found that this bursting takes place at lower values of the swirl ratio as the Reynolds number increases.

Jain, A. C.↗

Hilbert-Space Fragmentation from Strict Confinement

In this work, we study one-dimensional spin-1/2 models in which strict confinement of Ising domain walls leads to the fragmentation of Hilbert space into exponentially many disconnected subspaces. Whereas most previous works emphasize dipole moment conservation as an essential ingredient for such fragmentation, we instead require two commuting U(1) conserved quantities associated with the total domain-wall number and the total magnetization. The latter arises naturally from the confinement of domain walls. Remarkably, while some connected components of the Hilbert space thermalize, others are integrable by Bethe ansatz. We further demonstrate how this Hilbert-space fragmentation pattern arises perturbatively in the confining limit of $\mathbb{Z}_2$ gauge theory coupled to fermionic matter, leading to a hierarchy of timescales for motion of the fermions. This model can be realized experimentally in two complementary settings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Generation of Near-Equipartition Magnetic Fields in Turbulent Collisionless Plasmas

The mechanisms that generate “seed” magnetic fields in our Universe and that amplify them throughout cosmic time remain poorly understood. By means of fully kinetic particle-in-cell simulations of turbulent, initially unmagnetized plasmas, here we study the genesis of magnetic fields via the Weibel instability and follow their dynamo growth up to near-equipartition levels. In the kinematic stage of the dynamo, we find that the rms magnetic field strength grows exponentially with rate γ B ≃0.4⁢u rms /L, where L/2⁢π is the driving scale and u rms is the rms turbulent velocity. In the saturated stage, the magnetic field energy reaches about half of the turbulent kinetic energy. Here, magnetic field growth is balanced by dissipation via reconnection, as revealed by the appearance of plasmoid chains. At saturation, the integral-scale wave number of the magnetic spectrum approaches k int ≃12⁢π/L. Our results show that turbulence—induced by, e.g., the gravitational buildup of galaxies and galaxy clusters—can magnetize collisionless plasmas with large-scale near-equipartition fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Streaming Readout and Data-Stream Processing With ERSAP

With the exponential growth in the volume and complexity of data generated at high-energy physics and nuclear physics research facilities, there is an imperative demand for innovative strategies to process this data in real or near-real-time. Given the surge in the requirement for high-performance computing, it becomes pivotal to reassess the adaptability of current data processing architectures in integrating new technologies and managing streaming data. This paper introduces the ERSAP framework, a modern solution that synergizes flow-based programming with the reactive actor model, paving the way for distributed, reactive, and high performance in data stream processing applications. Additionally, we unveil a novel algorithm focused on time-based clustering and event identification in data streams. The efficacy of this approach is further exemplified through the data-stream processing outcomes obtained from the recent beam tests of the EIC prototype calorimeter at DESY.

Vardan, Gyurjyan↗

Hybrid learning techniques for scientific data reduction with performance guarantees

The research initiatives supported by the U.S. Department of Energy (DOE) Grant DE-SC0022265 are fundamentally aimed at pioneering advanced machine learning (ML) techniques for scientific data compression within high-performance computing (HPC) environments. This comprehensive body of work addresses the critical challenge posed by the exponential growth of data generated by scientific simulations in domains such as fusion energy, climate modeling, and computational fluid dynamics (CFD). A core objective is to develop compression algorithms that achieve substantial data reduction—often by orders of magnitude—while rigorously ensuring the fidelity of both the primary data (PD) and scientifically crucial derived quantities of interest (QoI). The methodologies deployed under this grant integrate sophisticated deep learning architectures, prominently featuring autoencoders, advanced generative models like conditional diffusion, and hybrid learning techniques. Key innovations include the development of Guaranteed Autoencoders (GAE) and the Guaranteed Conditional Diffusion with Tensor Correction (GCDTC) framework, which provide explicit, instance-level error bounds on reconstructed data. Furthermore, specialized strategies such as nonlinear constraint satisfaction are employed to preserve the integrity of QoI, a vital requirement for the trustworthiness of downstream scientific analyses. This research also focuses on the design and implementation of scalable, GPU-accelerated software pipelines that seamlessly integrate into existing HPC workflows, ensuring both computational efficiency and practical applicability. The CAESAR framework, for example, unifies foundation and generative models to create an adaptive and efficient compression solution for spatio-temporal scientific data. Collectively, these efforts represent a significant advancement in mitigating the scientific data deluge, enabling more effective data management, accelerated scientific discovery, and optimized utilization of HPC resources.

97 MATHEMATICS AND COMPUTING↗

Final report- UFL - RAPIDS2: A SciDAC Institute for Computer Science, Data, and Artificial Intelligence

The research initiatives supported by the U.S. Department of Energy (DOE) Grant DE-SC0022265 are fundamentally aimed at pioneering advanced machine learning (ML) techniques for scientific data compression within high-performance computing (HPC) environments. This comprehensive body of work addresses the critical challenge posed by the exponential growth of data generated by scientific simulations in domains such as fusion energy, climate modeling, and computational fluid dynamics (CFD). A core objective is to develop compression algorithms that achieve substantial data reduction—often by orders of magnitude—while rigorously ensuring the fidelity of both the primary data (PD) and scientifically crucial derived quantities of interest (QoI). The methodologies deployed under this grant integrate sophisticated deep learning architectures, prominently featuring autoencoders, advanced generative models like conditional diffusion, and hybrid learning techniques. Key innovations include the development of Guaranteed Autoencoders (GAE) and the Guaranteed Conditional Diffusion with Tensor Correction (GCDTC) framework, which provide explicit, instance-level error bounds on reconstructed data. Furthermore, specialized strategies such as nonlinear constraint satisfaction are employed to preserve the integrity of QoI, a vital requirement for the trustworthiness of downstream scientific analyses. This research also focuses on the design and implementation of scalable, GPU-accelerated software pipelines that seamlessly integrate into existing HPC workflows, ensuring both computational efficiency and practical applicability. The CAESAR framework, for example, unifies foundation and generative models to create an adaptive and efficient compression solution for spatio-temporal scientific data. Collectively, these efforts represent a significant advancement in mitigating the scientific data deluge, enabling more effective data management, accelerated scientific discovery, and optimized utilization of HPC resources.

97 MATHEMATICS AND COMPUTING↗

Suzaku observation of IGR J16318-4848

We report on the first Suzaku observation of IGR J16318-4848, the most extreme example of a new group of highly absorbed X-ray binaries that have recently been discovered by the International Gamma-Ray Astrophysics Laboratory (INTEGRAL). The Suzaku observation was carried out between 2006 August 14 and 17, with a net exposure time of 97 ks. The average X-ray spectrum of the source can be well described (chi-square (sub red)= 0.99) with a continuum model typical for neutron stars i.e., a strongly absorbed power law continuum with a photon index of 0.676(42) and an exponential cutoff at 20.5(6) keY. The absorbing column is N(sub H) = 1.95(3) X 10(exp 24)/square cm. Consistent with earlier work, strong fluorescent emission lines of Fe K-alpha, Fe K-beta, and Ni K-alpha are observed. Despite the large N(sub H), no Compton shoulder is seen in the lines, arguing for a non-spherical and inhomogeneous absorber. Seen at an average 5-60 keV absorbed flux of 3.4 x 10(exp -10) erg/square cm/second, the source exhibits significant variability on timescales of hours.

Barragan, Laura↗

Bioreactor Development for CO2-Based In Situ Resource Utilization Manufacturing

Sustainable long-duration manned missions on both the Moon and Mars will require in situ resource utilization (ISRU). Carbon dioxide (CO2) has great potential as a harvestable resource, making up 95% of the atmosphere on Mars and being produced as respiratory waste in spacecraft and future planetary habitats. Through ISRU, biomanufacturing has the capability to produce a near limitless array of products from local space resources. Here, a CO2-based ISRU recombinant protein producing bioreactor and associated biomanufacturing organisms were designed to produce a highly stable carbonic anhydrase (CA). Candidate organisms were selected by growth characterization on acetate and formic acid, carbon substrates that are synthesized via electrochemical conversion of CO2. To improve growth on the CO2 producing substrate formic acid and for direct integration of ISRU CO2, a synthetic Calvin-Benson-Bassam cycle was designed for use in Escherichia coli. Multiplex genetic modification in E. coli was facilitated by a tailored CRISPR/Cas9 and λ red recombineering two-vector system. For expression of CA, a blue light regulated T7 promoter was employed for dynamic and small molecule free induction. Efficient bioproduction through a fed-batch exponential feeding strategy was determined via mass balance analysis from ISRU substrates to biomass and CA yield. Flux balance analysis was used to model ISRU substrate metabolism and metabolic pathway engineering in candidate organisms under cultivation strategy conditions for both metabolism reconstruction and pathway design optimization. Finally, a small-scale, disposable bag bioreactor for use in the NASA Bioculture System infrastructure was designed to enable CO2-based CA biomanufacturing in reduced-gravity environments.

biomanufacturing↗

Two-dimensional materials for bio-realistic neuronal computing networks

Two-dimensional (2D) van der Waals materials have found broad utility in a diverse range of applications including electronics, optoelectronics, renewable energy, and quantum information technologies. Meanwhile, exponentially growing digital data coupled with the ubiquity of artificial intelligence algorithms have generated significant interest in edge neuromorphic computing as an alternative to centralized cloud computing. The drive to incorporate neuroscience principles into computing hardware is motivated by the low power consumption, parallel processing, and reconfigurability of the human brain. The diverse library of 2D materials with atomic-level thicknesses, exceptional electrostatic tunability, and integration versatility is particularly well-suited for realizing bio-realistic synaptic and neuronal functionality. Here, we summarize past and present work in this field and outline the frontier challenges that have not yet been overcome. Here we also delineate potential solutions and suggest that the neuroscience principles of criticality and synchrony have the potential to inspire breakthrough applications of 2D materials in neuronal computing networks.

36 MATERIALS SCIENCE↗

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗

The Stokes line width and uncertainty relations

For a function given by contour integral the two types (conventions) of asymptotic representations are considered: the usual representation by asymptotic series in inverse powers of large parameters and the special division of contour integral in contributions of high and low saddle points. It is shown that the width of the recessive term formation zone (Stokes strip) in the second convention is determined by uncertainty relation and is much less than the zone width in the first convention. The reasons of such a difference is clarified. The results of the work are useful for understanding of formation region of the exponentially small process arising on the background of the strong one.

Nikishov, A. I.↗

4U 1909+07: A Hidden Pearl

We present a detailed spectral and timing analysis of the High Mass X-ray Binary (HMXB) 4U 1909+07 with INTEGRAL and RXTE. 4U 1909+07 is a persistent accreting X-ray pulsar with a period of approximately 605 s. The period changes erratically consistent with a random walk expected for a wind accreting system. INTEGRAL detects the source with an average of 2.4 cps (corresponding to 15 mCrab), but sometimes exhibits flaring activity up to 50 cps (i.e. 300 mCrab). The strongly energy dependent pulse profile shows a double peaked structure at low energies and only a single narrow peak at energies above 20 keV. The phase averaged spectrum is well described by a powerlaw modified at higher energies by an exponential cutoff and photoelectric absorption at low energies. In addition at 6.4 keV a strong iron fluorescence line and at lower energies a black body component are present. We performed phase resolved spectroscopy to study the pulse phase dependence of the spectral parameters: while most spectral parameters are constant within uncertainties, the blackbody normalization and the cutoff folding energy vary strongly with phase.

Kreykenbohm, Ingo↗

Bioreactor Development for CO2-Based In Situ Resource Utilization Manufacturing

Sustainable long-duration manned missions on both the Moon and Mars will require in situ resource utilization (ISRU). Carbon dioxide (CO2) has great potential as a harvestable resource, making up 95% of the atmosphere on Mars and being produced as respiratory waste in spacecraft and future planetary habitats. Through ISRU, biomanufacturing has the capability to produce a near limitless array of products from local space resources, which include pharmaceuticals, bioplastics, chemical feedstocks, and industrial enzymes. Here, a CO2-based ISRU recombinant protein bioreactor and associated biomanufacturing organisms were designed to produce a highly stable carbonic anhydrase (CA). Initial work characterized candidate organisms for growth on acetate and formic acid, carbon substrates that can be synthesized via electrochemical conversion of CO2. To improve growth on the CO2 producing substrate formic acid and for direct integration of ISRU CO2, a synthetic Calvin-Benson-Bassam cycle was designed for use in Cyberlindnera jadinii and Escherichia coli. Genetic modifications in E. coli will be facilitated by a tailored CRISPR/Cas9 and λ red recombineering two-vector system. For expression of CA, a blue light regulated T7 promoter was employed for dynamic and small molecule free induction. Efficient bioproduction through a fed-batch exponential feeding strategy was determined via mass balance calculations from ISRU substrates to biomass and CA yield. Flux balance analysis was used to model ISRU substrate metabolism and metabolic pathway engineering in candidate organisms under cultivation strategy conditions for both metabolism reconstruction and pathway design optimization. Finally, a small-scale, disposable bag bioreactor concept for use in the NASA Bioculture System infrastructure was designed to enable CO2-based CA production in reduced-gravity environments.

Biomanufacturing, Pathway Engineering, Flux Balanc↗

Nonperturbative negative geometries: amplitudes at strong coupling and the amplituhedron

The amplituhedron determines scattering amplitudes in planar N = 4 super Yang-Mills by a single “positive geometry” in the space of kinematic and loop variables. We study a closely related definition of the amplituhedron for the simplest case of four-particle scattering, given as a sum over complementary “negative geometries”, which provides a natural geometric understanding of the exponentiation of infrared (IR) divergences, as well as a new geometric definition of an IR finite observable F(g, z) — dually interpreted as the expectation value of the null polygonal Wilson loop with a single Lagrangian insertion — which is directly determined by these negative geometries. This provides a long-sought direct link between canonical forms for positive (negative) geometries, and a completely IR finite post-loop-integration observable depending on a single kinematical variable z, from which the cusp anomalous dimension Γ cusp (g) can also be straightforwardly obtained. We study an especially simple class of negative geometries at all loop orders, associated with a “tree” structure in the negativity conditions, for which the contributions to F(g, z) and Γ cusp can easily be determined by an interesting non-linear differential equation immediately following from the combinatorics of negative geometries. This lets us compute these “tree” contributions to F(g, z) and Γ cusp for all values of the ‘t Hooft coupling. The result for Γ cusp remarkably shares all main qualitative characteristics of the known exact results obtained using integrability.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Advanced interactive display formats for terminal area traffic control

The basic design considerations for perspective Air Traffic Control displays are described. A software framework has been developed for manual viewing parameter setting (MVPS) in preparation for continued, ongoing developments on automated viewing parameter setting (AVPS) schemes. The MVPS system is based on indirect manipulation of the viewing parameters. Requests for changes in viewing parameter setting are entered manually by the operator by moving viewing parameter manipulation pointers on the screen. The motion of these pointers, which are an integral part of the 3-D scene, is limited to the boundaries of screen. This arrangement has been chosen, in order to preserve the correspondence between the new and the old viewing parameter setting, a feature which contributes to preventing spatial disorientation of the operator. For all viewing operations, e.g. rotation, translation and ranging, the actual change is executed automatically by the system, through gradual transitions with an exponentially damped, sinusoidal velocity profile, in this work referred to as 'slewing' motions. The slewing functions, which eliminate discontinuities in the viewing parameter changes, are designed primarily for enhancing the operator's impression that he, or she, is dealing with an actually existing physical system, rather than an abstract computer generated scene. Current, ongoing efforts deal with the development of automated viewing parameter setting schemes. These schemes employ an optimization strategy, aimed at identifying the best possible vantage point, from which the Air Traffic Control scene can be viewed, for a given traffic situation.

Grunwald, Arthur J.↗