Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential algebra”

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

Hidden zeros of the cosmological wavefunction

Motivated by the recent discovery of hidden zeros in particle and string amplitudes, we characterize zeros of individual graph contributions to the cosmological wavefunction of a scalar field theory. We demonstrate that these contributions factorize near these zeros for all tree graphs and provide evidence that this extends to loop graphs as well. We explicitly construct polytopal realizations of the relevant graph associahedra and show that the cosmological zeros have natural geometric and physical interpretations. As a byproduct, we establish an equivalence between the wavefunction coefficients of chain graphs and flat-space Tr(ϕ 3 ) amplitudes, enabling us to leverage the cosmological zeros to uncover the recently discovered hidden zeros of colored amplitudes.

Cosmological models↗
NASA NTRS1990

Calculation of turbulent flow in complex geometries with a second-moment closure model

A full Reynolds stress turbulence model is applied to calculate the flow in various channels with complex cross-section. The model is shown to be more robust and to produce more consistent results than algebraic stress models. The relationship between various versions of the latter is provided. It is shown that differential or algebraic stress models derived from a direct application of a popular linear approximation for the pressure-strain terms in the Reynolds stress equations will always underpredict the secondary motion, and hence the accompanying distortion of other flow properties. A simple approximation which leads to better predictions is presented and related to previous practice.

Demuren, A. O.↗
DOE OSTI2023

Ground vehicle lane-keeping assistance system via differential flatness output feedback control and algebraic derivative estimation

Vehicle run-off-road is one of the most frequent and fatal traffic accidents in the United States. Various lane-keeping assistance (LKA) systems have been developed in the last decade to help drivers stay on the road. Most of them are built upon linear driver–vehicle–road (DVR) models and treat road curvature as a disturbance. Albeit effective, their control performance would degrade if road curvature varies rapidly. Here, this paper proposes a novel nonlinear DVR model by integrating a driver steering model into vehicle–road kinematics, which explicitly considers road curvature. Particularly, this nonlinear DVR system has been proven to be differentially flat, and a flatness-based LKA system is designed. Additionally, Model-Free Control is introduced to compensate for system modeling errors. Hardware-in-the-loop simulations and diver-in-the-loop experiments validate the proposed control framework and demonstrate the performance enhancement with respect to a representative linear robust LKA system.

42 ENGINEERING↗
DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗
DOE OSTI2020

Chiral algebra, localization, modularity, surface defects, and all that

We study the 2D vertex operator algebra (VOA) construction in 4D N = 2 superconformal field theories on S 3 × S 1 , focusing on both old puzzles and new observations. The VOA lives on a two-torus T 2 ⊂ S 3 × S 1 , it is 1 2 Z -graded, and this torus is equipped with the natural choice of spin structure (1,0) for the Z + 1 2 -graded operators, corresponding to the NS sector vacuum character. By analyzing the possible refinements of the Schur index that preserves the VOA, we find that it admits discrete deformations, which allows access to the remaining spin structures (1,1), (0,1), and (0,0), of which the latter two involve the inclusion of a particular surface defect. For Lagrangian theories, we perform the detailed analysis: we describe the natural supersymmetric background, perform localization, and derive the gauged symplectic boson action on a torus in any spin structure. In the absence of flavor fugacities, the 2D and 4D path integrals precisely match, including the Casimir factors. We further analyze the 2D theory: we identify its integration cycle and the two-point functions and interpret flavor holonomies as screening charges in the VOA. Next, we make some observations about modularity; the T-transformation acts on our four partition functions and lifts to a large diffeomorphism on S 3 × S 1 . More interestingly, we generalize the four partition functions on the torus to an infinite family labeled by both the spin structure and the integration cycle inside the complexified maximal torus of the gauge group. Members of this family transform into one another under the full modular group, and we confirm the recent observation that the S-transform of the Schur index in Lagrangian theories exhibits logarithmic behavior. Finally, we comment on how locally our background reproduces the Ω-background.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Labeled trees and the efficient computation of derivations

The effective parallel symbolic computation of operators under composition is discussed. Examples include differential operators under composition and vector fields under the Lie bracket. Data structures consisting of formal linear combinations of rooted labeled trees are discussed. A multiplication on rooted labeled trees is defined, thereby making the set of these data structures into an associative algebra. An algebra homomorphism is defined from the original algebra of operators into this algebra of trees. An algebra homomorphism from the algebra of trees into the algebra of differential operators is then described. The cancellation which occurs when noncommuting operators are expressed in terms of commuting ones occurs naturally when the operators are represented using this data structure. This leads to an algorithm which, for operators which are derivations, speeds up the computation exponentially in the degree of the operator. It is shown that the algebra of trees leads naturally to a parallel version of the algorithm.

Grossman, Robert↗
DOE OSTI2022

Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers

In recent years, the SUite of Nonlinear and DIfferential/ALgebraic equation Solvers (SUNDIALS) has been redesigned to better enable the use of application-specific and third-party algebraic solvers and data structures. Throughout this work, we have adhered to specific guiding principles that minimized the impact to current users while providing maximum flexibility for later evolution of solvers and data structures. The redesign was done through the addition of new linear and nonlinear solvers classes, enhancements to the vector class, and the creation of modern Fortran interfaces. The vast majority of this work has been performed “behind-the-scenes,” with minimal changes to the user interface and no reduction in solver capabilities or performance. These changes allow SUNDIALS users to more easily utilize external solver libraries and create highly customized solvers, enabling greater flexibility on extreme-scale, heterogeneous computational architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Toward performance-portable PETSc for GPU-based exascale systems

The Portable Extensible Toolkit for Scientific computation (PETSc) library delivers scalable solvers for nonlinear time-dependent differential and algebraic equations and for numerical optimization. The PETSc design for performance portability addresses fundamental GPU accelerator challenges and stresses flexibility and extensibility by separating the programming model used by the application from that used by the library, and it enables application developers to use their preferred programming model, such as Kokkos, RAJA, SYCL, HIP, CUDA, or OpenCL, on upcoming exascale systems. Furthermore, a blueprint for using GPUs from PETSc-based codes is provided, and case studies emphasize the flexibility and high performance achieved on current GPU-based systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2022

Physics-Informed Sparse Gaussian Process for Probabilistic Stability Analysis of Large-Scale Power System with Dynamic PVs and Loads

This work proposes a physics-informed sparse Gaussian process (SGP) for probabilistic stability assessment of large-scale power systems in the presence of uncertain dynamic PVs and loads. The differential and algebraic equations considering uncertainties from dynamic PVs and loads are reformulated to a nonlinear mapping relationship that allows the application of SGP. Thanks to the nonparametric characteristic of Gaussian process, the proposed framework does not require distributions of uncertain inputs and this distinguishes it from existing approaches. As the original Gaussian process is not scalable to large-scale systems with high dimensional uncertain inputs, this paper develops the SGP with a stochastic variational inference technique. It leads to approximately two orders of complex reduction. A data pre-processing step is also introduced to tackle the coexistence of stable and unstable cases by sample clustering and constructing separate SGPs. The probabilistic transient stability index is analyzed to assess system stability under different uncertain dynamics loads and PVs. Comparisons are performed with the sampling-based, the polynomial chaos expansion-based, and traditional Gaussian process-based methods on the modified IEEE 118-bus and Texas 2000-bus systems under various scenarios, including different levels of uncertainties and the existence of nonlinear correlations among dynamic PVs. The impacts of data quality and quantity issues are also investigated. It is shown that the proposed SGP achieves significantly improved computational efficiency while maintaining high accuracy with a limited number of data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

PETSc/TAO developments for GPU-based early exascale systems

The Portable Extensible Toolkit for Scientific Computation (PETSc) library provides scalable solvers for nonlinear time-dependent differential and algebraic equations and for numerical optimization via the Toolkit for Advanced Optimization (TAO). PETSc is used in dozens of scientific fields and is an important building block for many simulation codes. During the U.S. Department of Energy’s Exascale Computing Project, the PETSc team has made substantial efforts to enable efficient utilization of the massive fine-grain parallelism present within exascale compute nodes and to enable performance portability across exascale architectures. We recap some of the challenges that designers of numerical libraries face in such an endeavor, and then discuss the many developments we have made, which include the addition of new GPU backends, features supporting efficient on-device matrix assembly, better support for asynchronicity and GPU kernel concurrency, and new communication infrastructure. In conclusion, we evaluate the performance of these developments on some pre-exascale systems as well as the early exascale systems Frontier and Aurora, using compute kernel, communication layer, solver, and mini-application benchmark studies, and then close with a few observations drawn from our experiences on the tension between portable performance and other goals of numerical libraries.

Exascale Computing Project (ECP)↗
DOE OSTI2023

High-Performance Electron Sources: Numerical Methods and Beam Dynamics at the Precision Frontier

Electron sources have a wide range of applications and there are many stakeholders that express continuing need for improvements and performance enhancements. Whether we consider ultra-cold, high-brightness, high-charge or high-average current source needs, there are some common themes from the point of view of the beam dynamics involved. These are the following ones: ability to model accurately the emission processes, including the presence of often complicated cathode and other boundary surfaces with a wide range of spatial scales; ability to deal accurately and efficiently with a large number of particles interacting pair-wise, including the stochastic part of these interactions with a wide range of spatial scales; and ability to propagate the particle distributions in time, including collisions with a wide range of temporal scales. These tasks require high precision and accuracy since the goal is usually generation, transport and preservation of very high-quality beams. This grant addressed one of the rem

43 PARTICLE ACCELERATORS↗
NASA NTRS1965

Universal variables.

Universal variable parameters and formulas used for two-body problem, differential correction and various orbits and perturbations

ORBIT CALCULATION↗
NASA NTRS1972

Equations of motion in the linear approximation.

Attempt to develop a gauge-invariant theory of the motion of singularities in an n-dimensional Riemannian space. A gauge-invariant theory is developed for a nongeodetic particle represented by a time-like world line in the linearized field theory. The solution of the basic algebraic and differential equations is to be constructed from the source line, the future null cones emanating from it, a scale factor, and nothing else, and the leading term of the solution is the Schwarzschild field. A class of possible solutions is examined, and it is shown that this class contains just one acceptable member - namely, a particle in constant acceleration.

Robinson, I.↗
NASA NTRS1977

Some comments on series solutions

The use of power series and truncated power series in the MACSYMA system for algebraic manipulation is illustrated. Algebraic and differential equations are solved using Taylor series or asymptotic series. Deficiencies of the current scheme are noted, and remedies suggested.

Fateman, R. J.↗