Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear dynamical 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,711 records · Page 95

PLEXUS: A Pattern-Oriented Runtime System Architecture for Resilient Extreme-Scale High-Performance Computing Systems

For high-performance computing (HPC) system designers and users, meeting the myriad challenges of next-generation exascale supercomputing systems requires rethinking their approach to application and system software design. Among these challenges, providing resiliency and stability to the scientific applications in the presence of high fault rates requires new approaches to software architecture and design. As HPC systems become increasingly complex, they require intricate solutions for detection and mitigation for various modes of faults and errors that occur in these large-scale systems, as well as solutions for failure recovery. These resiliency solutions often interact with and affect other system properties, including application scalability, power and energy efficiency. Therefore, resilience solutions for HPC systems must be thoughtfully engineered and deployed.In previous work, we developed the concept of resilience design patterns, which consist of templated solutions based on well-established techniques for detection, mitigation and recovery. In this paper, we use these patterns as the foundation to propose new approaches to designing runtime systems for HPC systems. The instantiation of these patterns within a runtime system enables flexible and adaptable end-to-end resiliency solutions for HPC environments. The paper describes the architecture of the runtime system, named Plexus, and the strategies for dynamically composing and adapting pattern instances under runtime control. This runtime-based approach enables actively balancing the cost-benefit trade-off between performance overhead and protection coverage of the resilience solutions. Based on a prototype implementation of PLEXUS, we demonstrate the resiliency and performance gains achieved by the pattern-based runtime system for a parallel linear solver application.

Hukerikar, Saurabh↗

A Hybrid Fuel Cell and Battery Storage Power Management for Grid-Interactive EV Charging Station

With the increasing adoption of renewable energy sources in grid-interactive Electric Vehicle (EV) charging stations, the role of energy storage systems has become critical. While large energy storage systems have mitigated the intermittency of renewable energy, integrating multi-source energy management with prioritized charging can further enhance the reliability of charging stations (CS). This paper presents a decentralized energy management (DEM) approach combining battery energy storage (BES) and fuel cell (FC) systems using a rule-based line resistance correction droop (LRCD) control technique. The proposed droop control dynamically adjusts the gain to balance the state-of-charge (SoC) of the BES, enhancing power support longevity and improving battery life under varying capacity conditions by reducing current stress. Additionally, the paper addresses the challenges of using fuel cells in linear regions to optimize efficiency and manage various charging scenarios. The CS integrates unity power factor grid interaction, and power support for auxiliary loads, maintaining harmonic distortion within 5% during grid islanding. The approach evaluates DC bus voltage regulation under various scenarios of PV array power fluctuations and dynamic load variations, in both grid-connected and standalone operations. In conclusion, the proposed control strategy is validated on a laboratory prototype through various dynamic load variation and grid islanding scenarios.

Khalid, Mohd [Oak Ridge National Laboratory (ORNL)↗

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗

ELAS: A general-purpose computer program for the equilibrium problems of linear structures. Volume 2: Documentation of the program

A general purpose digital computer program for the in-core solution of linear equilibrium problems of structural mechanics is documented. The program requires minimum input for the description of the problem. The solution is obtained by means of the displacement method and the finite element technique. Almost any geometry and structure may be handled because of the availability of linear, triangular, quadrilateral, tetrahedral, hexahedral, conical, triangular torus, and quadrilateral torus elements. The assumption of piecewise linear deflection distribution insures monotonic convergence of the deflections from the stiffer side with decreasing mesh size. The stresses are provided by the best-fit strain tensors in the least squares at the mesh points where the deflections are given. The selection of local coordinate systems whenever necessary is automatic. The core memory is used by means of dynamic memory allocation, an optional mesh-point relabelling scheme and imposition of the boundary conditions during the assembly time.

Utku, S.↗

SIR-B ocean-wave enhancement with fast Fourier transform techniques

Shuttle Imaging Radar (SIR-B) imagery is Fourier filtered to remove the estimated system-transfer function, reduce speckle noise, and produce ocean scenes with a gray scale that is proportional to wave height. The SIR-B system response to speckled scenes of uniform surfaces yields an estimate of the stationary wavenumber response of the imaging radar, modeled by the 15 even terms of an eighth-order two-dimensional polynomial. Speckle can also be used to estimate the dynamic wavenumber response of the system due to surface motion during the aperture synthesis period, modeled with a single adaptive parameter describing an exponential correlation along track. A Fourier filter can then be devised to correct for the wavenumber response of the remote sensor and scene correlation, with subsequent subtraction of an estimate of the speckle noise component. A linearized velocity bunching model, combined with a surface tilt and hydrodynamic model, is incorporated in the Fourier filter to derive estimates of wave height from the radar intensities corresponding to individual picture elements.

Tilley, David G.↗

High-resolution Wave Propagation Method for Stratified Flows

The implementation of the multidimensional f-waves Riemann solver for the time-dependent, three-dimensional, nonhydrostatic, meso- and microscale atmospheric flows is described in detail. The Riemann solver employs flux-based wave decomposition (f-waves) for the calculation of Godunov fluxes in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The scheme incorporates the source term due to gravity without introducing discretization errors which is an important property in the context of atmospheric flows. The resulting flow solver is conservative, accurate, stable, and well-balanced. The implementation of the solver is evaluated using benchmark test cases for atmospheric dynamics.

Riemann problem↗

Stability Assessment and Tuning of an Adaptively Augmented Classical Controller for Launch Vehicle Flight Control

Recently, a robust and practical adaptive control scheme for launch vehicles [ [1] has been introduced. It augments a classical controller with a real-time loop-gain adaptation, and it is therefore called Adaptive Augmentation Control (AAC). The loop-gain will be increased from the nominal design when the tracking error between the (filtered) output and the (filtered) command trajectory is large; whereas it will be decreased when excitation of flex or sloshing modes are detected. There is a need to determine the range and rate of the loop-gain adaptation in order to retain (exponential) stability, which is critical in vehicle operation, and to develop some theoretically based heuristic tuning methods for the adaptive law gain parameters. The classical launch vehicle flight controller design technics are based on gain-scheduling, whereby the launch vehicle dynamics model is linearized at selected operating points along the nominal tracking command trajectory, and Linear Time-Invariant (LTI) controller design techniques are employed to ensure asymptotic stability of the tracking error dynamics, typically by meeting some prescribed Gain Margin (GM) and Phase Margin (PM) specifications. The controller gains at the design points are then scheduled, tuned and sometimes interpolated to achieve good performance and stability robustness under external disturbances (e.g. winds) and structural perturbations (e.g. vehicle modeling errors). While the GM does give a bound for loop-gain variation without losing stability, it is for constant dispersions of the loop-gain because the GM is based on frequency-domain analysis, which is applicable only for LTI systems. The real-time adaptive loop-gain variation of the AAC effectively renders the closed-loop system a time-varying system, for which it is well-known that the LTI system stability criterion is neither necessary nor sufficient when applying to a Linear Time-Varying (LTV) system in a frozen-time fashion. Therefore, a generalized stability metric for time-varying loop=gain perturbations is needed for the AAC.

VanZwieten, Tannen↗

Finite element appendage equations for hybrid coordinate dynamic analysis.

The increasingly common practice of idealizing a spacecraft as a collection of interconnected rigid bodies to some of which are attached linearly elastic flexible appendages leads to equations of motion expressed in terms of a combination of discrete coordinates describing the arbitrary rotational motions of the rigid bodies and distributed or modal coordinates describing the small, time-varying deformations of the appendages: such a formulation is said to employ a hybrid system of coordinates. In the present paper the existing literature is extended to provide hybrid coordinate equations of motion for a finite element model of a flexible appendage attached to a rigid base undergoing unrestricted motions and some of the advantages of the finite element approach are noted. Transformations to the modal coordinates appropriate for the general case are provided.

Likins, P. W.↗

Numerical methods for multibody systems

This article gives a brief summary of some results obtained by Nasser on modeling and simulation of inequality problems in multibody dynamics. In particular, the augmented Lagrangian method discussed here is applied to a constrained motion problem with impulsive inequality constraints. A fundamental characteristic of the multibody dynamics problem is the lack of global convexity of its Lagrangian. The problem is transformed into a convex analysis problem by localization (piecewise linearization), where the augmented Lagrangian has been successfully used. A model test problem is considered and a set of numerical experiments is presented.

Glowinski, Roland↗

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗

Ultrafast temporal phase-resolved nonlinear optical spectroscopy in the molecular frame

In an ultrafast nonlinear optical interaction, the electric field of the emitted nonlinear signal provides direct access to the induced nonlinear transient polarization or transient currents and thus carries signatures of ultrafast dynamics in a medium. Measurement of the electric field of such signals offers sensitive observables to track ultrafast electron dynamics in various systems. In this work, we resolve the real-time phase of the electric field of a femtosecond third-order nonlinear optical signal in the molecular frame. The electric field emitted from impulsively pre-aligned gas-phase molecules at room temperature, in a degenerate four-wave mixing scheme, is measured using a spectral interferometry technique. The nonlinear signal is measured around a rotational revival to extract its molecular-frame angle dependence from pump-probe time-delay scans. By comparing these measurements for two linear molecules, carbon dioxide and nitrogen, we show that the measured second-order phase parameter (temporal chirp) of the signal is sensitive to the valence electronic symmetry of the molecules, whereas the amplitude of the signal does not show such sensitivity. We compare measurements to theoretical calculations of the chirp observable in the molecular frame. This work is an important step towards using electric field measurements in nonlinear optical spectroscopy to study ultrafast dynamics of electronically excited molecules in the molecular frame.

47 OTHER INSTRUMENTATION↗

Signatures of an energetic charged body streaming in a plasma

A charged body moving in a plasma can excite a variety of linear and nonlinear waves in the form of trailing wakes, fore-wake shocks, and precursor solitons. These structures can further interact with the background plasma to create secondary effects that can serve as signatures of the passage of the charged body. Using particle-in-cell simulations, we carry out a basic investigation of the dynamics of a plasma system that is being traversed by an energetic charged body. Using two different shapes of this charged source, namely, an idealized infinite length planar source and a two-dimensional thin rectangular source, we examine the differences in the nature of the excited wave structures and their consequent impact on the background plasma. Our simulations reveal interesting features such as the dependence of the precursor speeds on the total charge on the driving source, local particle trapping, and energization of the trapped particles in various regions along the traversal path leading to the formation of energetic charged beamlets. So our basic findings could find practical applications such as in analyzing the trajectories of charged objects like space debris orbiting in the ionosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, one can assume that the total temperature is constant. Therefore, the energy equations can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed in this paper. It is shown that the system is strictly hyperbolic and well posed for the initial-value problem. Boundary conditions are described such that the linearized system is well posed. The hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗

An application of high authority/low authority control and positivity

Control Dynamics Company (CDy), in conjunction with NASA Marshall Space Flight Center (MSFC), has supported the U.S. Air Force Wright Aeronautical Laboratory (AFWAL) in conducting an investigation of the implementation of several DOD controls techniques. These techniques are to provide vibration suppression and precise attitude control for flexible space structures. AFWAL issued a contract to Control Dynamics to perform this work under the Active Control Technique Evaluation for Spacecraft (ACES) Program. The High Authority Control/Low Authority Control (HAC/LAC) and Positivity controls techniques, which were cultivated under the DARPA Active Control of Space Structures (ACOSS) Program, were applied to a structural model of the NASA/MSFC Ground Test Facility ACES configuration. The control systems design were accomplished and linear post-analyses of the closed-loop systems are provided. The control system designs take into account effects of sampling and delay in the control computer. Nonlinear simulation runs were used to verify the control system designs and implementations in the facility control computers. Finally, test results are given to verify operations of the control systems in the test facility.

Seltzer, S. M.↗

Dynamic loading of spur gears with linear or parabolic tooth profile modification

A computer simulation was conducted to investigate the effects of both linear and parabolic tooth profile modification on the dynamic response of low-contact-ratio spur gears. The effect of the total amount of modification and the length of the modification zone were studied at various loads and speeds to find the optimal profile modification for minimal dynamic loading. Design charts consisting of normalized maximum dynamic load curves were generated for gear systems operated at various loads and with different tooth profile modification. An optimum profile modification can be determined from these design charts to minimize the dynamic loads of spur gear systems.

Lin, Hsiang Hsi↗

Recursive Newton-Euler formulation of manipulator dynamics

A recursive Newton-Euler procedure is presented for the formulation and solution of manipulator dynamical equations. The procedure includes rotational and translational joints and a topological tree. This model was verified analytically using a planar two-link manipulator. Also, the model was tested numerically against the Walker-Orin model using the Shuttle Remote Manipulator System data. The hinge accelerations obtained from both models were identical. The computational requirements of the model vary linearly with the number of joints. The computational efficiency of this method exceeds that of Walker-Orin methods. This procedure may be viewed as a considerable generalization of Armstrong's method. A six-by-six formulation is adopted which enhances both the computational efficiency and simplicity of the model.

Nasser, M. G.↗

Dynamic loading of spur gears with linear or parabolic tooth profile modifications

A computer simulation was conducted to investigate the effects of both linear and parabolic tooth profile modification on the dynamic response of low-contact-ratio spur gears. The effect of the total amount of modification and the length of the modification zone were studied at various loads and speeds to find the optimal profile modification for minimal dynamic loading. Design charts consisting of normalized maximum dynamic load curves were generated for gear systems operated at various loads and with different tooth profile modification. An optimum profile modification can be determined from these design charts to minimize the dynamic loads of spur gear systems.

Lin, H. H.↗

Hybrid state-space self-tuning control of uncertain linear systems

The paper presents a hybrid state-space self-tuner using a new dual-rate sampling scheme for digital adaptive control of continuous-time uncertain linear systems. A state-space-based recursive least-squares algorithm, together with a variable forgetting factor, is used for direct estimations of both the equivalent discrete-time uncertain linear system parameters and the associated discrete-time state of a continuous-time uncertain linear system from the sampled input and output data. An analogue optimal regional pole-placement design method is used for designing an optimal observer-based analogue controller. A suboptimal observer-based digital controller is then designed from the designed analogue controller using digital redesign technique. To enhance the robustness of parameter identification and state estimation algorithms, a dynamic bound for a class of uncertain bilinear parameters and a fast-rate digital controller are developed at each fast-sampling period. Also, to accommodate computation loads and computation delay for developing the advanced hybrid self-tuner, the designed analogue controller and observer gains are both updated at each slow-sampling period. This control technique has been successfully applied to benchmark control problems.

Shieh, L. S.↗