Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear systems”

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 1,261 records · Page 70

Measurement bias in self-heating x-ray free electron laser experiments from diffraction studies of phase transformation in titanium

X-ray self-heating is a common by-product of X-ray Free Electron Laser (XFEL) techniques that can affect targets, optics, and other irradiated materials. Diagnosis of heating and induced changes in samples may be performed using the x-ray beam itself as a probe. However, the relationship between conditions created by and inferred from x-ray irradiation is unclear and may be highly dependent on the material system under consideration. Here, we report on a simple case study of a titanium foil irradiated, heated, and probed by a MHz XFEL pulse train at 18.1 keV delivered by the European XFEL using measured x-ray diffraction to determine temperature and finite element analysis to interpret the experimental data. We find a complex relationship between apparent temperatures and sample temperature distributions that must be accounted for to adequately interpret the data, including beam averaging effects, multivalued temperatures due to sample phase transitions, and jumps and gaps in the observable temperature near phase transformations. The results have implications for studies employing x-ray probing of systems with large temperature gradients, particularly where these gradients are produced by the beam itself. Finally, this study shows the potential complexity of studying nonlinear sample behavior, such as phase transformations, where biasing effects of temperature gradients can become paramount, precluding clear observation of true transformation conditions.

Crystallography↗

GFiRe—Gauge Field integrator for Reheating

We present a new numerical algorithm and code, GFiRe, for solving the nonlinear evolution of Abelian gauge fields coupled to complex scalar fields in homogeneous and isotropic spacetimes. We adopt a hybrid approach to solving the system: the spatial derivatives are discretized using standard Lattice Gauge Field Theory techniques, whereas the time evolution of the fields and scale factor is implemented with explicit, composite, symplectic integrators. An important property of our compound algorithm is that the discretized Gauss constraint is respected exactly, regardless of the order of the symplectic integrator. This remains true even when the background expansion is computed “self-consistently”; that is, when the expansion history is computed using spatial averaged components of the energy momentum tensor in the simulation volume. Hence, our code can also be used in cases where the fields dominate the energy density of the universe, for example, during reheating after inflation. We test the algorithm in scenarios of reheating where the inflaton is a complex scalar field with a potential ∝(2|Φ| 2 - v 2 ) 2 and is coupled to an Abelian gauge field. Tracing the evolution of the system through complex dynamics (including resonant excitation of fields, backreaction, formation of solitons, and changes in the equation of state) in a self-consistently expanding universe, we find the energy conservation violation (<10 -4 ) to be very stable and the Gauss constraint violation (<10 -6 ) to be dominated by differencing noise.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamics and transport in the boundary-driven dissipative Klein-Gordon chain

Motivated by experiments on chains of superconducting qubits, we consider the dynamics of a classical Klein-Gordon chain coupled to coherent driving and subject to dissipation solely at its boundaries. As the strength of the boundary driving is increased, this minimal classical model recovers the main features of the "dissipative phase transition" seen experimentally. Between the transmitting and non-transmitting regimes on either side of this transition (which support ballistic and diffusive energy transport respectively), we observe additional dynamical regimes of interest. These include a regime of superdiffusive energy transport at weaker driving strengths, together with a "dissipative breather crystal" regime at stronger driving strengths, which is characterized by emergent translation symmetry, ballistic energy transport, and coherent oscillations of a nonlinear normal mode. Here, we propose a non-local Lyapunov exponent as an experimentally measurable diagnostic of many-body chaos in this system, and more generally in open systems that are only coupled to an environment at their boundaries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Diffusive dynamics of critical fluctuations near the QCD critical point

A quantitatively reliable theoretical description of the dynamics of fluctuations in nonequilibrium is indispensable in the experimental search for the QCD critical point by means of ultrarelativistic heavy-ion collisions. Herein we consider the fluctuations of the net-baryon density which becomes the slow, critical mode near the critical point. Due to net-baryon number conservation the dynamics is described by the fluid dynamical diffusion equation, which we extend to contain a white noise stochastic current. Including nonlinear couplings from the 3d Ising model universality class in the free energy functional, we solve the fully interacting theory in a finite size system. We observe that purely Gaussian white noise generates non-Gaussian fluctuations, but finite size effects and exact net-baryon number conservation lead to significant deviations from the expected behavior in equilibrated systems. In particular the skewness shows a qualitative deviation from infinite volume expectations.With this benchmark established we study the real-time dynamics of the fluctuations. We recover the expected dynamical scaling behavior and observe retardation effects and the impact of critical slowing down near the pseudocritical temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deep Factorization Machine Learning for Disaggregation of Transmission Load Profiles with High Penetration of Behind-The-Meter Solar

The ever-growing integration of distributed energy resources (DERs), especially behind-the-meter (BTM) solar generations, poses imperative operational challenges to system operators such as regional transmission organizations (RTOs). It is important for RTOs to effectively and accurately extract actual load profiles at the transmission level for a single node with significant BTM solar injection. This paper first illustrates the necessity of disaggregating the daily actual load profile of a single node. Furthermore, by segmenting nodes with selected timeseries features, nodes with significant BTM solar generation are identified. Lastly, a bi-level framework is proposed, comprising reference node disaggregation and DeepFM nodal disaggregation, aimed at disaggregating the nodal load profiles from which system operators require more information. By adopting a hybrid Deep Factorization Machine (DeepFM) model, the model achieve accurate results by extracting both linear and nonlinear relations between nodes in the same region and the zonal load and nodal load profile. To overcome the lack of ground truth, this paper segments the load profile into daytime, nighttime, and zero-crossing points and utilizes the latter two for evaluation purposes. The proposed disaggregation procedure is validated using real world, minute-level, normalized, and anonymized nodal data in the PJM service territory.

42 ENGINEERING↗

Correlation-Aided Robust Decentralized Dynamic State Estimation of Power Systems with Unknown Control Inputs

This paper proposes a correlation-aided robust adaptive unscented Kalman filter for power system decentralized dynamic state estimation with unknown inputs, termed as robust AUKF-UI. The temporal and spatial correlations among the unknown inputs are used to derive a vector auto-regressive (VAR) model in an adaptive manner. This VAR model is further integrated together with state transition and measurement models for joint state and unknown inputs estimation. This allows taking into account the implicit cross-correlations between the states and the unknown inputs. As a result, the rank requirement for unknown input vector estimation is relaxed and the local generator frequency measurement is not required. The temporal correlations of time series innovation vectors, predicted state and input vectors are also leveraged by the robust AUKFUI to detect, identify and process bad data. Without these correlations, it is very challenging to address bad data with unknown inputs. Simulation results carried out on the IEEE 39-bus system demonstrate that the proposed robust AUKF-UI achieves much better results than other methods in the presence of low measurement redundancy, strong nonlinearity, and bad data.

Zhao, Junbo↗

Pretest Modeling A Spent Nuclear Fuel Seismic Shake Test

The U.S. Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a dry canister storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental U.S. to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. This paper describes the pretest predictions made with LS-DYNA models of the mockup storage systems. The test will use two dry storage system configurations, a mockup vertical concrete cask (VCC) and a mockup horizontal storage module (HSM). The test will use a production-quality canister and basket. Within the canister will be four instrumented fuel assemblies with fuel rods containing surrogate mass and 28 instrumented dummy assemblies that are intended to match the mass and outer dimension of a fuel assembly. The finite element models include models of the VCC and HSM on the shake table to calculate the system level dynamic responses and separate single fuel assembly models to calculate stress and strain in fuel assembly components. Both types of models include nonlinear behavior like rod-to-rod contact and the ability for VCC’s to rock and slide. This paper presents the expected response of the VCC, HSM, and fuel assemblies to the shake table testing that is planned to start in April of 2024. The earthquake conditions represent seismic hazards in the 2,000-to-20,000-year return period range. The test data is expected to confirm the expectation that fuel rod cladding will remain intact, fuel assembly structural components like guide tubes will remain intact, and no significant VCC sliding or tipping will occur in the range of conditions to be tested with the shake table.

Klymyshyn, Nicholas A.↗

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Indexing habitat opportunity for juvenile anadromous fishes in tidal-fluvial wetland systems

In intertidal wetlands, annual habitat use by migrating juvenile salmon is limited by both physical and biological constraints. Physical limits are set by inundation patterns and water quality parameters, while biological bounds are determined by species specific migration timing and seasonal residency behaviors. We developed a new metric called the fish habitat opportunity index (FHOI) to quantify the total time per calendar year that salmon could access intertidal sites, based on periods of inundation and juvenile salmon seasonal migration. The index was applied using data collected at a recently reconnected tidal freshwater wetland area during 2006-2009. We found the wetland to be accessible to different salmon species for different lengths of time and periods of the year. On average, chum salmon (Oncorhynchus keta) had remarkably consistent opportunities to use marsh habitat over the 4 years sampled, while opportunities were more variable for Chinook (O. tshawytscha), and coho salmon (O. kisutch). Inundation varied nonlinearly with tidal height, and fish had access to the productive marsh edge for approximately 40-50% of the time during periods of access. The FHOI is widely applicable to wetland systems throughout the tidal fluvial continuum, and provides a metric for assessing fish habitat opportunity that is more realistic than simple wetted area calculations. Finally, the FHOI affords a method of a priori comparison among potential restoration sites to help resource managers predict the likely relative benefit to juvenile salmon.

wetland restoration, estuary, habitat opportunity,↗

A General Decentralized Dynamic State Estimation With Synchronous Generator Magnetic Saturation

We report the saturation of the nonlinear magnetic circuit of synchronous generators is often neglected when performing Kalman filter-based dynamic state estimation (DSE), yielding significant estimation bias. This letter addresses this problem and proposes a generalized DSE framework to handle magnetic saturation. Moreover, this letter derives a state initialization procedure that improves the Kalman filter tracking speed. The framework is flexible in dealing with different saturation functions and generator models. Numerical results on the Texas 2000-bus system verify the effectiveness of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Analog Systems for Edge Optimization

Over the past decade, analog computing has the subject of substantial research interest providing a path toward improved computational efficiency in the post-Dennard era. Analog matrix vector multiplication (MVM) accelerators provide a popular approach given the ubiquity of MVM operations in numerous applications. However, historically analog computing systems can struggle with applications requiring high precision due to the inherent susceptibility of these systems to analog non-idealities. Therefore, prior work on analog systems has focused either on applications known to be tolerant of limited precision (e.g., neural network inference), or using expensive techniques to emulate high-precision using many analog MVM operations. In this work, we propose an alternative approach. Motivated by recent advances in inexact nonlinear solvers and optimizers, we explore the potential of co-designing optimization algorithms which can take full advantage of the fundamentally inexact analog MVM operations. To enable these co-designed algorithms we also develop a general mathematical theory of the precision and energy efficiency of analog operations, and a new system architecture for tightly-coupled analog and digital computation. Finally, we examine the applicability of analog computing to a wider class of symmetric positive definite systems and find potential in using analog operations as a sparse approximate inverse preconditioner. With these core innovations, this project provides a path toward effectively implementing optimization algorithms on power-constrained autonomous and semi-autonomous systems.

97 MATHEMATICS AND COMPUTING↗

A Fast Temporal Decomposition Procedure for Long-Horizon Nonlinear Dynamic Programming

We propose a fast temporal decomposition procedure for solving long-horizon nonlinear dynamic programs. The core of the procedure is sequential quadratic programming (SQP) that utilizes a differentiable exact augmented Lagrangian as the merit function. Within each SQP iteration, we approximately solve the Newton system using an overlapping temporal decomposition strategy. We show that the approximate search direction is still a descent direction of the augmented Lagrangian provided the overlap size and penalty parameters are suitably chosen, which allows us to establish the global convergence. Moreover, we show that a unit step size is accepted locally for the approximate search direction and further establish a uniform, local linear convergence over stages. This local convergence rate matches the rate of the recent Schwarz scheme (Na et al. 2022). However, the Schwarz scheme has to solve nonlinear subproblems to optimality in each iteration, whereas we only perform a single Newton step instead. Numerical experiments validate our theories and demonstrate the superiority of our method.

97 MATHEMATICS AND COMPUTING↗

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗

Efficacy of vibro-impact energy harvesting absorbers on controlling dynamical systems under vortex-induced vibrations and base excitation

Here this study investigates improving the efficacy of an energy harvesting absorber's ability to control a structure under vortex-induced vibrations, base excitation, and a combination of the two by including mechanical amplitude stoppers. The nonlinear reduced-order model is developed through modifying trilinear spring models to represent the impact forces, a modified van der Pol oscillator to represent the forcing due to the vortex-induced vibrations and using the Euler-Lagrange principle to express the equations of motion. It is seen that a soft stopper stiffness and a 5mm gap performs the most effectively of increasing the power generated from the absorber while still greatly reducing the primary structure's amplitude. By changing the stopper's location towards the middle of the energy harvesting absorber, the large effects of the impact forces are reduced and improves the efficacy of medium and hard stopper stiffnesses to generate near the amount of power the soft stopper does, while greatly improving the control of the primary structure. When the system is under combined loadings, the large oscillations of the synchronization region cause the effective configuration to be that of a 27.5 mm gap with soft stiffnesses. The results shows that medium stiffness stoppers with small gaps generate large aperiodic regions due to the high impact force. When the oscillations are close to the stoppers, the beating phenomenon is observed and is not overpowered by the vibro-impact force.

16 TIDAL AND WAVE POWER↗

Brewster angle-cavity ringdown spectroscopy for low temperature plasma measurements in multiphases

Here we report on the development of a Brewster angle-cavity ringdown spectroscopy (BA-CRDS) system for low temperature plasma diagnostics. The system can measure gas species in solutions, with a detection limit (minimum detectable absorbance) of 9.1 × 10 -5 , which is equivalent to a detection limit of 0.04 parts per billion for measuring OH radicals in water at 308 nm. With higher reflectivity ringdown mirrors and improved design of a Brewster angle cell, the detection limit can potentially be up to 10 -6 or lower. In this exploratory study, the absorption cross sections of HgBr 2 and H 2 O 2 in the aqueous phase at 256 nm are measured to be (1.8 ± 0.1) × 10 -18 cm 2 and (5.2 ± 0.5) × 10 -20 cm 2 , respectively. Furthermore, temporal profiles of absorbance from distilled water, HgBr 2 , and H 2 O 2 solutions when interacting with a helium atmospheric plasma jet are individually characterized at different plasma powers, gas flow rates, and/or solute concentrations. The observed linear temporal profiles of absorbance from the plasma-interacted water suggest formation of H 2 O 2 from plasma-generated OH radicals, while the nonlinear temporal profiles from the plasma-treated HgBr 2 solutions reveal possible removal of HgBr 2 by OH radicals. Our results demonstrate that the new BA-CRDS system is a powerful tool for quantification of reactive plasma species in multiphases or other complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling active nematics via the nematic locking principle

Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.

Mitchell, Kevin A. [Univ. of California, Merced, C↗

Computation and consequences of high order amplitude- and parameter-dependent tune shifts in storage rings for high precision measurements

Nonlinear effects of the various electric field and magnetic field components of storage rings to confine the particles and bend their trajectory can cause substantial amplitude-dependent tune shifts within the beam. Furthermore, tune shifts are often sensitive to variations of system parameters, e.g. total particle momentum offsets $δp$. Such amplitude- and parameter-dependent tune shifts influence the dynamics and stability of a beam in particle storage rings. Thus, it is critical for high precision measurements to analyze and understand these influences. On this basis, we present in this paper normal form methods for the calculation of high order amplitude and system parameter dependencies of the horizontal and vertical tunes in storage rings using the differential algebra (DA) framework within COSY INFINITY. A storage ring is simulated using COSY INFINITY to generate a DA Poincaré return map describing the transverse phase space behavior after each revolution in the storage ring. The map is expanded around the parameter-dependent closed orbit of the system before transforming the resulting map into normal form coordinates to extract the high order tune dependencies on the phase space amplitude and variation in the system parameters. As a specific example, a storage ring similar to the Storage Ring of the Muon g-2 Experiment at Fermilab (E989) is investigated.

43 PARTICLE ACCELERATORS↗

Deep Learning enabled spectral energy conversion for in situ exposure measurements

A detector-specific deep learning (DL) approach is presented for spectra-to-exposure conversion using large-format sodium iodide (NaI(Tl)) detectors deployed for in situ environmental radiation measurements in emergency response scenarios. Accurate determination of exposure from NaI spectra is challenging due to poor energy resolution, partial energy absorption, and the strong sensitivity of traditionally deployed analytical conversion methods to calibrated source geometry and pre-deployment assumptions. Here, to address these limitations, a multi-layer perceptron model was trained on a hybrid in situ /Monte Carlo dataset constructed to span a broad range of photon energies, spatial extents, and realistic deployment variability, representative of general in situ emergency response conditions. The DL model was evaluated against commonly fielded analytical approaches under matched simulation conditions, including a single-factor method, a G-function method, and a modeled pressurized ion chamber (PIC) baseline. This study was intentionally computational in scope to enable controlled, like-for-like comparisons between conversion techniques while minimizing confounding real-world variability. Comparison to the modeled PIC provides contextual benchmarking and is not intended as a field inter-comparison with deployed instruments. Across the evaluated 20 keV to 3 MeV energy range, the DL approach consistently exhibited higher accuracy and reduced variance relative to the analytical methods against a deterministically calculated exposure. This may indicate improved robustness to spectral complexity without reliance on source-, geometric-, or spectral region-specific optimization. While results do not represent real-world validation, the presented work demonstrates that deep learning may effectively learn the nonlinear detector response-to-exposure relationship for asymmetric NaI(Tl) detectors and offers a promising pathway for improving in situ exposure estimation using spectroscopic systems already integrated into initial real-time emergency response operations.

61 RADIATION PROTECTION AND DOSIMETRY↗