Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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 415 records · Page 23

Holographic Kolmogorov-Sinai entropy and the quantum Lyapunov spectrum

In classical chaotic systems the entropy, averaged over initial phase space distributions, follows a universal behavior. While approaching thermal equilibrium it passes through a stage where it grows linearly, while the growth rate, the Kolmogorov-Sinai entropy (rate), is given by the sum over all positive Lyapunov exponents. A natural question is whether a similar relation is valid for quantum systems. We argue that the Maldacena-Shenker-Stanford bound on quantum Lyapunov exponents implies that the upper bound on the growth rate of the entropy, averaged over states in Hilbert space that evolve towards a thermal state with temperature T, should be given by πT times the thermal state’s von Neumann entropy. Strongly coupled, large N theories with black hole duals should saturate the bound. To test this we study a large number of isotropization processes of random, spatially homogeneous, far from equilibrium initial states in large N, $\mathcal{N}$ = 4 Super Yang Mills theory at strong coupling and compute the ensemble averaged growth rate of the dual black hole’s apparent horizon area. We find both an analogous behavior as in classical chaotic systems and numerical evidence that the conjectured bound on averaged entropy growth is saturated granted that the Lyapunov exponents are degenerate and given by λ i = ±2πT. This fits to the behavior of classical systems with plus/minus symmetric Lyapunov spectra, a symmetry which implies the validity of Liouville’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A method for characterization of multiple dynamic constitutive parameters of FRCs

We propose a method to measure multiple dynamic material constitutive parameters of unidirectional fiber reinforced composites (FRCs) in a single experiment. Dynamic short-beam shear (DSBS) experiments were performed on a modified Kolsky compression bar, with integration of high-speed imaging and digital image correlation (DIC). The unidirectional FRCs investigated were S-2 glass fiber reinforced matrix of TGDDM-Jeffamine® D230 with monoamine functionalized partially reacted substructures (mPRS) and commercially available SC-15. Analytical solutions of normal and shear strains of a composite beam were derived based on Timoshenko beam theory, assuming material to be transversely isotropic and have different moduli in tension and compression in each principle material orientation. Tensile and compressive moduli were inversely computed through monitoring normal strain slope when specimen was constantly loaded at a speed of ~7.3 m/s within a small deflection. Non-linear shear stress-strain behavior of the composite was described via Ramberg–Osgood equation. Finite element (FE) analysis was conducted in ABAQUS, simultaneously defining via user subroutine UMAT the transverse isotropy of material, bi-modulus constitutive model, and non-linear shear stress-strain relation. The method proposed in this work was validated by comparing strain distributions computed by FE model and DIC measurements. Comparing with traditional dynamic tensile, compressive, and shear experiments on FRCs, this method significantly simplifies the specimen preparation and design of complicated gripping fixtures for multiple experiments. Here, systematic errors resulting from variations of specimen geometry and dimension, loading direction, and instrumentation are reduced, thereby providing compatible data for numerical studies on impact behavior of composites.

36 MATERIALS SCIENCE↗

Fourier-informed knot placement schemes for B-spline approximation

Fitting B-splines to scientific data is especially challenging when the given data contain noise, jumps, or corners. Here, we describe how periodic data sets with these features can be efficiently approximated with B-splines by analyzing the Fourier spectrum of the data. Our method uses a collection of spectral filters to compute high-order derivatives, smoothed versions of noisy data, and the locations of jump discontinuities. Further, these quantities are then combined to choose knots that capture the qualitative features of the data, leading to accurate B-spline approximations with few knots. The method we introduce is direct and does not require any intermediate B-spline fitting before choosing the final knot distribution. Aside from fast Fourier transforms to transfer to and from Fourier space, the method runs in linear time with very little communication. We assess performance on several test cases in one and two dimensions, including data sets with jump discontinuities and noise. These tests show the method fits discontinuous data without spurious oscillations and remains effective in the presence of noise.

97 MATHEMATICS AND COMPUTING↗

Surface-Doped Zinc Gallate Colloidal Nanoparticles Exhibit pH-Dependent Radioluminescence with Enhancement in Acidic Media

As abnormal acidic pH symbolizes dysfunctions of cells, it is highly desirable to develop pH-sensitive luminescent materials for diagnosing disease and imaging-guided therapy using high-energy radiation. Herein, we explored near-infrared-emitting Cr-doped zinc gallate ZnGa 2 O 4 nanoparticles (NPs) in colloidal solutions with different pH levels under X-ray excitation. Ultrasmall NPs were synthesized via a facile hydrothermal method by controlling the addition of ammonium hydroxide precursor and reaction time, and structural characterization revealed Cr dopants on the surface of NPs. The synthesized NPs exhibited different photoluminescence and radioluminescence mechanisms, confirming the surface distribution of activators. It was observed that the colloidal NPs emit pH-dependent radioluminescence in a linear relationship, and the enhancement reached 4.6-fold when pH = 4 compared with the colloidal NPs in the neutral solution. This observation provides a strategy for developing new biomaterials by engineering activators on the nanoparticle surfaces for potential pH-sensitive imaging and imaging-guided therapy using high-energy radiation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cationic α-Diimine Nickel and Palladium Complexes Incorporating Phenanthrene Substituents: Highly Active Ethylene Polymerization Catalysts and Mechanistic Studies of syn/anti Isomerization

α-Diimine palladium complexes incorporating phenanthryl- and 6,7-dimethylphenanthrylimino groups have been synthesized and characterized. The (diimine)PdMeCl complexes prepared from 2,3-butanedione and acenaphthenequinone bearing the unsubstituted phenanthrylimino groups, 12a and 14a, respectively, exist as a mixture syn and anti isomers in a ca. 1:1 ratio. Separation and X-ray diffraction analysis of 14a-syn and 14a-anti isomers confirms the syn/anti assignments. The barrier to interconversion of 14a-syn and 14a-anti via ligand rotation, ΔG ‡ , was found to be 25.5 kcal/mol. The corresponding (diimine)PdMeCl complex prepared from acenaphthenequinone and incorporating the 6,7-dimethylphenanthrylimino group exists solely as the anti isomer, 14b, due to steric crowding which destabilizes the syn isomer. Analogous (diimine)NiBr2 complexes were prepared from 2,3-butanedione incorporating the phenanthrylimino group, 16a, and the 6,7-dimethylphenanthrylimino group, 16b. Nickel-catalyzed polymerizations of ethylene were carried out by activation of the dibromide complexes 16a and 16b using various aluminum alkyl activators. Complex 16a yields a bimodal distribution polymer, the low molecular weight fraction originating from the syn isomer and the high molecular weight fraction arising from the anti isomer. Polymerizations carried out by 16b yield only high molecular weight polymers with monomodal distributions due to the existence of a single isomer (anti) as the active catalyst. All polymers are linear or nearly so. All catalysts are highly active but catalysts derived from 16b are somewhat more active than 16a and exhibit turnover frequencies general over 10 6 and up to 5 × 10 6 per hour (40 °C, 27.2 atm ethylene, 15 min). Here, active palladium ethylene oligomerization catalysts were generated by conversion of the neutral methyl chloride complexes 14a and 14b to the cationic nitrile complexes 15a and 15b via halide abstraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

One-Parameter Analytical Derivation in Modified Budyko Framework for Unsteady-State Streamflow Elasticity in Humid Catchments

The changes in climate and catchment properties have altered the hydrological processes significantly at different spatiotemporal scales around the world. In particular for finer time scales, changes in water storage, which has been commonly neglected for long-term temporal scales, may play an important role on hydrological processes. Nevertheless, few studies addressed this question in modifying the Budyko framework, with most of them remaining a steady-state assumption on catchment characteristics. Here, in this work, we derive new analytical formulas of unsteady-state streamflow elasticity in a modified Budyko framework, incorporating both storage change and one specific parameter for catchment properties. We study 78 humid catchments in the USA with simulation data of daily hydrological processes from a probability-distribution-based hydrological model (the modified HyMOD). As indicated by results, the annual storage change ratio is linearly correlated with the annual aridity index, and this relationship can be used to estimate elasticity coefficients with our formulas. The estimated elasticity coefficients perform well in simulating the annual streamflow with the power-law model. For different catchments, variability of the unsteady-state elasticity is higher than that of the steady-state elasticity. Unsteady-state streamflow coefficients show significant linear correlation with catchment properties, such as the average slope, average elevation, and catchment area. This study provides a new analytical approach to investigate the interannual stability of catchments with varying climate and catchment properties.

54 ENVIRONMENTAL SCIENCES↗

Topological solitonic macromolecules

Being ubiquitous, solitons have particle-like properties, exhibiting behaviour often associated with atoms. Bound solitons emulate dynamics of molecules, though solitonic analogues of polymeric materials have not been considered yet. Here we experimentally create and model soliton polymers, which we call “polyskyrmionomers”, built of atom-like individual solitons characterized by the topological invariant representing the skyrmion number. With the help of nonlinear optical imaging and numerical modelling based on minimizing the free energy, we reveal how topological point defects bind the solitonic quasi-atoms into polyskyrmionomers, featuring linear, branched, and other macromolecule-resembling architectures, as well as allowing for encoding data by spatial distributions of the skyrmion number. Application of oscillating electric fields activates diverse modes of locomotion and internal vibrations of these self-assembled soliton structures, which depend on symmetry of the solitonic macromolecules. Our findings suggest new designs of soliton meta matter, with a potential for the use in fundamental research and technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Integrated structured light architectures

Abstract The structural versatility of light underpins an outstanding collection of optical phenomena where both geometrical and topological states of light can dictate how matter will respond or display. Light possesses multiple degrees of freedom such as amplitude, and linear, spin angular, and orbital angular momenta, but the ability to adaptively engineer the spatio-temporal distribution of all these characteristics is primarily curtailed by technologies used to impose any desired structure to light. We demonstrate a laser architecture based on coherent beam combination offering integrated spatio-temporal field control and programmability, thereby presenting unique opportunities for generating light by design to exploit its topology.

42 ENGINEERING↗

A variational framework for residual-based adaptivity in neural PDE solvers and operator learning

Residual-based adaptive strategies are widely used in scientific machine learning yet remain largely heuristic. We introduce a variational framework that formalizes these methods through convex transformations of the residual, where different transformations correspond to distinct objective functionals. For instance, exponential weights target uniform error minimization, while linear weights recover quadratic error minimization. This perspective reveals adaptive weighting as a means of selecting sampling distributions that optimize a primal objective, directly linking discretization choices to error metrics. This principled approach yields three key benefits: it enables systematic design of adaptive schemes, reduces discretization error by lowering estimator variance, and enhances learning dynamics by improving gradient signal-to-noise ratio. Extending the framework to operator learning, we demonstrate substantial performance gains across diverse optimizers and architectures. Our results provide a theoretical perspective for residual-based adaptivity and establish a foundation for principled discretization and training.

97 MATHEMATICS AND COMPUTING↗

Data-Driven Exploration of Climate Attractor Manifolds For Long-Term Predictability

Focal Area: This white paper responds to Focal Area 3. We seek to gain insight into decadal-scale climate predictability by applying novel manifold-finding probabilistic AI techniques to the complex data produced by Earth System models (ESMs) such as E3SM. The associated portfolio of research activities leverages DOE’s asset mix of HPC platforms, climate expertise, climate simulation codes, and AI expertise. Science Challenge: Climate and climate models are dynamical systems exhibiting properties that are interpretable through chaos theory. The theory contains an important concept that is relevant to multi-decade-scale climate prediction: a chaotic attractor. While the space containing all the possible states of the Earth’s atmosphere and ocean, the possible weather, is large, the realized states tend to stay near the smaller-dimensioned attractor. This behavior is responsible for the “order behind the irregularity” [1] of climate phenomena. Climate change can be thought of as a change in the properties of the attractor, and predicting the climate over years to decades is equivalent to predicting how those properties will change. To date, the attractor has been a useful conceptual tool, but has not been amenable to direct characterization. A new development is the advent of efficient high-dimensional manifold-finding probabilistic AI techniques, which permit a data-driven characterization of the ESM attractor and its probability distribution over weather states. Such a characterization would result in a natural dimensional reduction — a “non-linear Principal Components Analysis (PCA) adapted to climate simulation data” — leading to important advances in scenario-based long-term climate prediction, long-term prediction of water cycle extremes, ESM verification, inter-model comparison, and process model development.

54 ENVIRONMENTAL SCIENCES↗

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermodynamics of the Inner Heliosheath

We derive annual sky maps of the proton temperature in the inner heliosheath (IHS), and track their temporal evolution over the years 2009–2016 of Interstellar Boundary Explorer observations. Other associated thermodynamic parameters also determined are the density, kappa (the parameter that characterizes kappa distributions), temperature rate, polytropic index, and entropy. We exploit the theory of kappa distributions and their connection with polytropes, to (i) express a new polytropic quantity Π that remains invariant along streamlines where temperature and density may vary, (ii) parameterize the proton flux in terms of the Π invariant and kappa, and (iii) derive the temperature and density, respectively, from the slope and intercept of the linear relationship between kappa and logarithm of Π. We find the following thermodynamic characteristics: (1) temperature sky maps and histograms shifted to their lowest values in 2012 and their highest in 2015; (2) temperature negatively correlated with density, reflecting the subisothermal polytropic behavior; (3) temperature positively correlated with kappa, revealing characteristics of the mechanism responsible for generating kappa distributions; (4) processes in IHS are subisothermal tending toward isobaric, consistent with previously published results; (5) linear relationship between kappa and polytropic indices, revealing characteristics of the particle potential energy; and (6) entropy positively correlated with polytropic index, aligned with the underlying theory that entropy increases toward the isothermal state where the kappa distribution reduces to the Maxwell–Boltzmann description.

79 ASTRONOMY AND ASTROPHYSICS↗

Discovery of False Data Injection Attacks on Power Grid Frequency Controllers with Reinforcement Learning [Poster]

While inverter-based DER (distributed energy resources) are instrumental to integrating renewable energy into the power grid, they reduce the grid's mechanical inertia, thereby increasing the risk of frequency instabilities. To compensate for frequency instability risks, the grid must also undergo a transformation to include digital technologies that allow for two-way communication between the utility and customers. The current and future state of the power grid allows for building a cleaner energy landscape. However, the grid may also become vulnerable to novel cyber threats. To preemptively protect the power grid against elaborate cyber-attacks, we propose to discover potential threats via reinforcement learning. In this work, the focus is on studying false data injection attacks that target the control logic of frequency controllers. We show that a reinforcement learning agent can successfully discover how to best inject false data into linear droop controllers.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Taylor-Expansion-Based Robust Power Flow in Unbalanced Distribution Systems: A Hybrid Data-Aided Method

Traditional power flow methods often adopt certain assumptions designed for passive balanced distribution systems, thus lacking practicality for unbalanced operation. moreover, their computation accuracy and efficiency are heavily subject to unknown errors and bad data in measurements or prediction data of distributed energy resources (ders). to address these issues, this paper proposes a hybrid data-aided robust power flow algorithm in unbalanced distribution systems, which combines taylor series expansion knowledge with a data-driven regression technique. the proposed method initiates a linearization power flow model to derive an explicitly analytical solution by modified taylor expansion. to mitigate the approximation loss that surges due to the der integration and bad data, we further develop a data-aided robust support vector regression approach to estimate the errors efficiently. comparative analysis in the 13-bus and 123-bus ieee unbalanced feeders shows that the proposed hybrid algorithm achieves superior computational efficiency, with guaranteed accuracy and robustness against outliers.

data-driven↗

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model↗

Optimizing design and dispatch of a renewable energy system with combined heat and power

We embellish a mixed-integer program that prescribes a set of renewable energy, conventional generation, and storage technologies to procure, along with a corresponding dispatch strategy. Specifically, we add combined heat and power to this set. The model minimizes fixed and operational costs less incentives for the use of various technologies, subject to a series of component interoperability and system-wide constraints. The resulting mixed-integer linear program contains hundreds of thousands of variables and constraints. We demonstrate how to efficiently formulate and solve the corresponding instances such that we produce near-optimal solutions in minutes. A previous rendition of the model required hours of solution time for the same instances.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Data-Driven Linear Parameter-Varying Modeling and Control of Flexible Loads for Grid Services: Preprint

Flexible loads have great potential to improve the power grid's flexibility and stability. To effectively control large ensembles of heterogeneous loads, reliable models thereof are required. The paper presents a data-driven modeling and control approach to manage flexible loads for providing grid services. We leverage a linear parameter-varying autoregressive moving average (LPV-ARMA) model to describe the aggregate load response, where the parameters in the model are used to capture external environment impacts (e.g., weather). A gain-scheduling feedback controller is then developed to adapt to environmental variations. This data-driven approach can be easily applied to different types of loads in various environmental conditions. In addition to the ensemble controller, distributed load controllers are designed to deliver grid services, while maintaining the quality of service of inherent load tasks. We demonstrate the work on the IEEE 37-node distribution system for real-time power regulation services through control of thermostatically controlled loads.

control↗