Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Non linear dynamics”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Chemoviscosity modeling for thermosetting resins - I

A new analytical model for chemoviscosity variation during cure of thermosetting resins was developed. This model is derived by modifying the widely used WLF (Williams-Landel-Ferry) Theory in polymer rheology. Major assumptions involved are that the rate of reaction is diffusion controlled and is linearly inversely proportional to the viscosity of the medium over the entire cure cycle. The resultant first order nonlinear differential equation is solved numerically, and the model predictions compare favorably with experimental data of EPON 828/Agent U obtained on a Rheometrics System 4 Rheometer. The model describes chemoviscosity up to a range of six orders of magnitude under isothermal curing conditions. The extremely non-linear chemoviscosity profile for a dynamic heating cure cycle is predicted as well. The model is also shown to predict changes of glass transition temperature for the thermosetting resin during cure. The physical significance of this prediction is unclear at the present time, however, and further research is required. From the chemoviscosity simulation point of view, the technique of establishing an analytical model as described here is easily applied to any thermosetting resin. The model thus obtained is used in real-time process controls for fabricating composite materials.

Hou, T. H.↗
NASA NTRS2010

Mixed Layer Sub-Mesoscale Parameterization - Part 1: Derivation and Assessment

Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,

Canuto, V. M.↗
NASA NTRS2015

Advanced Computational Modeling Approaches for Shock Response Prediction

Motivation: (1) The activation of pyroshock devices such as explosives, separation nuts, pin-pullers, etc. produces high frequency transient structural response, typically from few tens of Hz to several hundreds of kHz. (2) Lack of reliable analytical tools makes the prediction of appropriate design and qualification test levels a challenge. (3) In the past few decades, several attempts have been made to develop methodologies that predict the structural responses to shock environments. (4) Currently, there is no validated approach that is viable to predict shock environments overt the full frequency range (i.e., ~100 Hz to 10 kHz). Scope: (1) Model, analyze, and interpret space structural systems with complex interfaces and discontinuities, subjected to shock loads. (2) Assess the viability of a suite of numerical tools to simulate transient, non-linear solid mechanics and structural dynamics problems, such as shock wave propagation.

Derkevorkian, Armen↗
NASA NTRSDate not supplied

Validation of Ullage Collapse Due to Violent Lateral Slosh

Understanding the coupling between thermodynamics and propellant sloshing is important for cryogenic upper stages, especially during the lift-off phase and in missions including multiple engine restarts, as the fluid dynamic condition of the propellant is mainly influenced by various flight maneuvers. This paper describes the development and validation of a Volume of Fluid (VOF)-based CFD solver, Loci/STREAM-VOF, to the ullage collapse under violent later slosh conditions. Loci/STREAM-VOF is shown to be able to capture detailed slosh dynamics, including slosh frequency, non-linear slosh damping, and violent surface breakup due to forced acceleration. The comparison of surface shape and small drops formed during large amplitude slosh with experimental imagery shows excellent agreement throughout the sloshing period. Comparison of temperature histories in the ullage, at the gas-liquid interface, and inside liquid with the experiment conducted at JAXA also show very good agreement. The simulation points to a stratified temperature field inside the liquid near the interface, critical to the ullage pressure after the collapse. The predicted pressure collapse matches well with the experimental value due to sloshing. Parametric sensitivity studies suggest that high-fidelity Loci/STREM-VOF simulations can play an essential role in predicting the actual thermodynamics of ullage collapse when following the recommended simulation methodology.

H Q Yang↗
NASA NTRS2018

Modeling the Interactive Effects of Spruce Beetle Infestation and Climate on Subalpine Vegetation

In the subalpine zone of the Rocky Mountains, climate change is predicted to result in an increase in the frequency and severity of spruce beetle outbreaks. Climate change itself may affect vegetation, potentially leading to changes in species composition. The direct and indirect effects of climate and disturbances on forest composition, biomass, and dynamics open the possibility for non-linear ecosystem responses. Modeling studies allow for the study of the interaction of these effects and their impact on the forest system. University of Virginia Forest Model Enhanced (UVAFME), an individual-based gap model that simulates forest dynamics and characteristics, is updated with a spruce beetle subroutine that calculates the probability for beetle infestation and potential mortalityof each tree on a plot. The updated model is then run with multiple scenarios that combine beetle infestation with current or altered climate at sites across the southern Rocky Mountains. Results show that spruce beetle infestations acted to facilitate competition with invading lower-elevation species, resulting in an increase in the biomass of historically lower elevation species and a further decline in Engelmann spruce biomass than occurred with solely bark beetle disturbance or solely climate change. We also found an initial enhancing effect between spruce beetle infestation and climate change; however, by the end of 100 yr of climate change and potential beetle infestation, climate had a dampening effect on spruce beetle infestation, through loss of host trees. These results are an important step in understanding the possible futures for vegetation of the Rocky Mountains as well as for spruce forests across the western United States and Canada.

Foster, Adrianna C.↗
NASA NTRS2004

Overview of the GRC Stirling Convertor System Dynamic Model

A Stirling Convertor System Dynamic Model has been developed at the Glenn Research Center for controls, dynamics, and systems development of free-piston convertor power systems. It models the Stirling cycle thermodynamics, heat flow, gas, mechanical, and mounting dynamics, the linear alternator, and the controller. The model's scope extends from the thermal energy input to thermal, mechanical dynamics, and electrical energy out, allowing one to study complex system interactions among subsystems. The model is a non-linear time-domain model containing sub-cycle dynamics, allowing it to simulate transient and dynamic phenomena that other models cannot. The model details and capability are discussed.

Lewandowski, Edward J.↗
DOE OSTI2025

Linear Discriminant Analysis-Based Machine Learning and All-Atom Molecular Dynamics Simulations for Probing Electro-Osmotic Transport in Cationic-Polyelectrolyte-Brush-Grafted Nanochannels

Deciphering the correct mechanisms governing certain phenomena in polyelectrolyte (PE) brush grafted systems, revealed through atomistic simulations, is an extremely challenging problem. In a recent study, our all-atom molecular dynamics (MD) simulations revealed a non-linearly large electroosmotic (EOS) flow (in the presence of an applied electric field) in nanochannels grafted with PMETAC [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] brushes. Given the lack of any formal procedure that would have directed us to identify the correct factors responsible for such an occurrence, we needed to spend several months and devote significant analyses to unravel the involved mechanisms. In this paper, we propose a Linear Discriminant Analysis (LDA) based Machine Learning (ML) approach to address this gap. At first, we obtain data on certain basic features from the all-atom MD data. These basic features represent the number of atoms of certain species around one atom of another (or same) species. Here, we obtain such data on basic features for a reference case (case of an EOS flow in PMETAC-brush-grafted nanochannels with a smaller electric field) and a perturbed case (case of an EOS flow in PMETAC-brush-grafted nanochannels with a larger electric field) in bins in which the nanochannel half height has been divided into. These datasets are high-dimensional dataset, to which the LDA is applied. This leads to the projection of the data (between the reference and the perturbed states) in a highly separated form on a 1D line. From such LDA calculations, we are able to identify the relative importance of the different basic features in ensuring this separation of the data (between the reference and the perturbed states) on the 1D line. This relative importance of the different basic features is quantified as “importance scores” for the different features, which in turn tell us what to study and where to study. Such knowledge enables us to rapidly identify the key factors responsible for the non-linearly large EOS transport in PMETAC-brush-grafted nanochannels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Solving differential‐algebraic equations in power system dynamic analysis with quantum computing

Abstract Power system dynamics are generally modeled by high dimensional non‐linear differential‐algebraic equations (DAEs) given a large number of components forming the network. These DAEs' complexity can grow exponentially due to the increasing penetration of distributed energy resources, whereas their computation time becomes sensitive due to the increasing interconnection of the power grid with other energy systems. This paper demonstrates the use of quantum computing algorithms to solve DAEs for power system dynamic analysis. We leverage a symbolic programming framework to equivalently convert the power system's DAEs into ordinary differential equations (ODEs) using index reduction methods and then encode their data into qubits using amplitude encoding. The system non‐linearity is captured by Hamiltonian simulation with truncated Taylor expansion so that state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can solve the power system's DAEs accurately with a computational complexity polynomial in the logarithm of the system dimension. We also illustrate the use of recent advanced tools in scientific machine learning for implementing complex computing concepts, that is, Taylor expansion, DAEs/ODEs transformation, and quantum computing solver with abstract representation for power engineering applications.

computational complexity↗
DOE OSTI2022

The DREAM Endstation at the Linac Coherent Light Source

Free-electron lasers (FEL), with their ultrashort pulses, ultrahigh intensities, and high repetition rates at short wavelength, have provided new approaches to Atomic and Molecular Optical Science. One such approach is following the birth of a photo electron to observe ion dynamics on an ultrafast timescale. Such an approach presents the opportunity to decipher the photon-initiated structural dynamics of an isolated atomic and molecular species. It is a fundamental step towards understanding single- and non-linear multi-photon processes and coherent electron dynamics in atoms and molecules, ultimately leading to coherent control following FEL research breakthroughs in pulse shaping and polarization control. A key aspect for exploring photoinduced quantum phenomena is visualizing the collective motion of electrons and nuclei in a single reaction process, as dynamics in atoms/ions proceed at femtosecond (10 –15 s) timescales while electronic dynamics take place in the attosecond timescale (10 –18 s). Here, we report on the design of a Dynamic Reaction Microscope (DREAM) endstation located at the second interaction point of the Time-Resolved Molecular and Optical (TMO) instrument at the Linac Coherent Light Source (LCLS) capable of following the photon–matter interactions by detecting ions and electrons in coincidence. The DREAM endstation takes advantage of the pulse properties and high repetition rate of LCLS-II to perform gas-phase soft X-ray experiments in a wide spectrum of scientific domains. With its design ability to detect multi-ions and electrons in coincidence while operating in step with the high repetition rate of LCLS-II, the DREAM endstation takes advantage of the inherent momentum conservation of reaction product ions with participating electrons to reconstruct the original X-ray photon–matter interactions. In this report, we outline in detail the design of the DREAM endstation and its functionality, with scientific opportunities enabled by this state-of-the-art instrument.

47 OTHER INSTRUMENTATION↗
NASA NTRS1996

Guidance of Nonlinear Nonminimum-Phase Dynamic Systems

The research work has advanced the inversion-based guidance theory for: systems with non-hyperbolic internal dynamics; systems with parameter jumps; and systems where a redesign of the output trajectory is desired. A technique to achieve output tracking for nonminimum phase linear systems with non-hyperbolic and near non-hyperbolic internal dynamics was developed. This approach integrated stable inversion techniques, that achieve exact-tracking, with approximation techniques, that modify the internal dynamics to achieve desirable performance. Such modification of the internal dynamics was used (a) to remove non-hyperbolicity which is an obstruction to applying stable inversion techniques and (b) to reduce large preactuation times needed to apply stable inversion for near non-hyperbolic cases. The method was applied to an example helicopter hover control problem with near non-hyperbolic internal dynamics for illustrating the trade-off between exact tracking and reduction of preactuation time. Future work will extend these results to guidance of nonlinear non-hyperbolic systems. The exact output tracking problem for systems with parameter jumps was considered. Necessary and sufficient conditions were derived for the elimination of switching-introduced output transient. While previous works had studied this problem by developing a regulator that maintains exact tracking through parameter jumps (switches), such techniques are, however, only applicable to minimum-phase systems. In contrast, our approach is also applicable to nonminimum-phase systems and leads to bounded but possibly non-causal solutions. In addition, for the case when the reference trajectories are generated by an exosystem, we developed an exact-tracking controller which could be written in a feedback form. As in standard regulator theory, we also obtained a linear map from the states of the exosystem to the desired system state, which was defined via a matrix differential equation.

Devasia, Santosh↗
DOE OSTI2023

$\text{GPLaSDI}$: Gaussian Process-based interpretable Latent Space Dynamics Identification through deep autoencoder

Numerically solving partial differential equations (PDEs) can be challenging and computationally expensive. This has led to the development of reduced-order models (ROMs) that are accurate but faster than full order models (FOMs). Recently, machine learning advances have enabled the creation of non-linear projection methods, such as Latent Space Dynamics Identification (LaSDI). LaSDI maps full-order PDE solutions to a latent space using autoencoders and learns the system of ODEs governing the latent space dynamics. By interpolating and solving the ODE system in the reduced latent space, fast and accurate ROM predictions can be made by feeding the predicted latent space dynamics into the decoder. In this paper, we introduce GPLaSDI, a novel LaSDI-based framework that relies on Gaussian process (GP) for latent space ODE interpolations. Using GPs offers two significant advantages. First, it enables the quantification of uncertainty over the ROM predictions. Second, leveraging this prediction uncertainty allows for efficient adaptive training through a greedy selection of additional training data points. This approach does not require prior knowledge of the underlying PDEs. Consequently, GPLaSDI is inherently non-intrusive and can be applied to problems without a known PDE or its residual. Here we demonstrate the effectiveness of our approach on the Burgers equation, Vlasov equation for plasma physics, and a rising thermal bubble problem. Our proposed method achieves between 200 and 100,000 times speed-up, with up to 7% relative error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

An energy-dependent electro-thermal response model of CUORE cryogenic calorimeter

The Cryogenic Underground Observatory for Rare Events (CUORE) is the most sensitive experiment searching for neutrinoless double-beta decay (0νββ) in 130 Te. CUORE uses a cryogenic array of 988 TeO 2 calorimeters operated at ~10 mK with a total mass of 741 kg. To further increase the sensitivity, the detector response must be well understood. Here, we present a non-linear thermal model for the CUORE experiment on a detector-by-detector basis. We have examined both equilibrium and dynamic electro-thermal models of detectors by numerically fitting non-linear differential equations to the detector data of a subset of CUORE channels which are well characterized and representative of all channels. Further, we demonstrate that the hot-electron effect and electric-field dependence of resistance in NTD-Ge thermistors alone are inadequate to describe our detectors' energy-dependent pulse shapes. We introduce an empirical second-order correction factor in the exponential temperature dependence of the thermistor, which produces excellent agreement with energy-dependent pulse shape data up to 6 MeV. We also present a noise analysis using the fitted thermal parameters and show that the intrinsic thermal noise is negligible compared to the observed noise for our detectors.

47 OTHER INSTRUMENTATION↗
NASA NTRS2012

Development of a High Fidelity Dynamic Module of the Advanced Resistive Exercise Device (ARED) Using Adams

NASA's Digital Astronaut Project (DAP) implements well-vetted computational models to predict and assess spaceflight health and performance risks, and enhance countermeasure development. DAP provides expertise and computation tools to its research customers for model development, integration, or analysis. DAP is currently supporting the NASA Exercise Physiology and Countermeasures (ExPC) project by integrating their biomechanical models of specific exercise movements with dynamic models of the devices on which the exercises were performed. This presentation focuses on the development of a high fidelity dynamic module of the Advanced Resistive Exercise Device (ARED) on board the ISS. The ARED module, illustrated in the figure below, was developed using the Adams (MSC Santa Ana, California) simulation package. The Adams package provides the capabilities to perform multi rigid body, flexible body, and mixed dynamic analyses of complex mechanisms. These capabilities were applied to accurately simulate: Inertial and mass properties of the device such as the vibration isolation system (VIS) effects and other ARED components, Non-linear joint friction effects, The gas law dynamics of the vacuum cylinders and VIS components using custom written differential state equations, The ARED flywheel dynamics, including torque limiting clutch. Design data from the JSC ARED Engineering team was utilized in developing the model. This included solid modeling geometry files, component/system specifications, engineering reports and available data sets. The Adams ARED module is importable into LifeMOD (Life Modeler, Inc., San Clemente, CA) for biomechanical analyses of different resistive exercises such as squat and dead-lift. Using motion capture data from ground test subjects, the ExPC developed biomechanical exercise models in LifeMOD. The Adams ARED device module was then integrated with the exercise subject model into one integrated dynamic model. This presentation will describe the development of the Adams ARED module including its capabilities, limitations, and assumptions. Preliminary results, validation activities, and a practical application of the module to inform the relative effect of the flywheels on exercise will be discussed.

Humphreys, B. T.↗
NASA NTRS2001

The Use of Non-Standard Devices in Finite Element Analysis

A general mathematical description of the response behavior of thin-skin pneumatic envelopes and many other membrane and cable structures produces under-constrained systems that pose severe difficulties to analysis. These systems are mobile, and the general mathematical description exposes the mobility. Yet the response behavior of special under-constrained structures under special loadings can be accurately predicted using a constrained mathematical description. The static response behavior of systems that are infinitesimally mobile, such as a non-slack membrane subtended from a rigid or elastic boundary frame, can be easily analyzed using such general mathematical description as afforded by the non-linear, finite element method using an implicit solution scheme if the incremental uploading is guided through a suitable path. Similarly, if such structures are assembled with structural lack of fit that provides suitable self-stress, then dynamic response behavior can be predicted by the non-linear, finite element method and an implicit solution scheme. An explicit solution scheme is available for evolution problems. Such scheme can be used via the method of dynamic relaxation to obtain the solution to a static problem. In some sense, pneumatic envelopes and many other compliant structures can be said to have destiny under a specified loading system. What that means to the analyst is that what happens on the evolution path of the solution is irrelevant as long as equilibrium is achieved at destiny under full load and that the equilibrium is stable in the vicinity of that load. The purpose of this paper is to alert practitioners to the fact that non-standard procedures in finite element analysis are useful and can be legitimate although they burden their users with the requirement to use special caution. Some interesting findings that are useful to the US Scientific Balloon Program and that could not be obtained without non-standard techniques are presented.

Schur, Willi W.↗
DOE OSTI2024

Optimal Control of SOEC-Based Hydrogen Production Systems for Demand Response Using Deep Reinforcement Learning in Smart Grids

Solid oxide electrolysis cell (SOEC) hydrogen production technology can range in size from small, appliance-size equipment to large-scale, central production facilities that can be tied directly to renewable or non-greenhouse-gas-emitting forms of electricity production, making it an ideal resource for demand response (DR). The SOEC hydrogen production system is a complex integrated system that encompasses fluid dynamics, electrical dynamics, and electrochemical and thermal dynamics, all of which involve non-linearity and non-convexity. Proper control of the SOEC hydrogen production system is crucial to enable its participation in the DR program. Here, to overcome the difficulty of designing an explicit control law for such nonlinear systems with nonconvex optimization features in DR applications, deep reinforcement learning (DRL) is explored to achieve the optimal control of the SOEC system for DR participation. Specifically, a twin delayed deterministic policy gradient (TD3) control framework is applied to achieve optimal response performance during DR events by considering power tracking error and hydrogen production efficiency with a suitable reward function. Two case studies with grid connections for tracking different DR commands were investigated. The first case study involved operating conditions reaching the boundaries, while the second involved operating conditions within the boundaries. The results showed that the proposed DRL-based control for SOEC can track the DR signal in a timely manner while maintaining high energy efficiency.

08 HYDROGEN↗
DOE OSTI2023

Detailed Simulations of a Three-Stage Supercritical Carbon Dioxide Axial Compressor with a Focus on the Shrouded Stator Cavities Flow

This paper describes the findings of detailed simulations performed to investigate the impact of seal teeth cavity leakage flow on the aerodynamic and thermal performance of a three-stage supercritical CO 2 axial compressor. The study compares a shrouded stator configuration (with cavities) to a cantilevered stator configuration (without cavities) to highlight their differences. High-fidelity computational fluid dynamics simulations were performed using non-linear harmonic (NLH) and mixing plane assumptions, considering various possible rotor/stator interface configurations for mixing plane calculations. The key performance parameters for each case were compared, and the best-performing configuration selected for further analysis. The individual stage performance parameters are also examined and compared between the cantilevered and shrouded configurations. It was observed that in the shrouded case, the leakage flow enters the cavity downstream of the stator trailing edge and gets entrained into the primary flow upstream of the stator, leading to boundary layer changes at the hub and degradation of stator and downstream rotor performance. Vortical flow structures were also observed in the stator wells, which tended to change the flow angles around the region, thereby affecting mixing and velocity distribution, resulting in a slight deterioration of compressor performance. Additionally, the study examines windage heating due to shear work from rotating walls, including the seal teeth surface. The amount of shear work done on the leakage flow and the corresponding rise in fluid temperature were quantified, tabulated, and further compared with a simple analytical model, showing good agreement between them and, hence, validating the numerical approach used.

42 ENGINEERING↗
NASA NTRS1983

Space Shuttle propulsion parameter estimation using optimal estimation techniques, volume 1

The mathematical developments and their computer program implementation for the Space Shuttle propulsion parameter estimation project are summarized. The estimation approach chosen is the extended Kalman filtering with a modified Bryson-Frazier smoother. Its use here is motivated by the objective of obtaining better estimates than those available from filtering and to eliminate the lag associated with filtering. The estimation technique uses as the dynamical process the six degree equations-of-motion resulting in twelve state vector elements. In addition to these are mass and solid propellant burn depth as the ""system'' state elements. The ""parameter'' state elements can include aerodynamic coefficient, inertia, center-of-gravity, atmospheric wind, etc. deviations from referenced values. Propulsion parameter state elements have been included not as options just discussed but as the main parameter states to be estimated. The mathematical developments were completed for all these parameters. Since the systems dynamics and measurement processes are non-linear functions of the states, the mathematical developments are taken up almost entirely by the linearization of these equations as required by the estimation algorithms.

Source record↗
NASA NTRS1986

Effects of orography on planetary scale flow

The earth's orography is composed of a wide variety of scales, each contributing to the spectrum of atmospheric motions. A well studied subject (originating with Charney and Eliassen) is the direct forcing of planetary scale waves by the planetary scale orography: primarily the Tibetan plateau and the Rockies. However, because of the non-linear terms in the equations of dynamic meteorology, even the smallest scales of mountain induced flow can contribute to the planetary scale if the amplitude of the small scale disturbance is sufficintly large. Two possible mechanisms for this are illustrated. First, preferentially located lee cyclones can force planetary waves by their meridional transport of heat and momentum (Hansen and Chen). Recent theories are helping to explain the phenomena of lee cyclogenesis (e.g., Smith, 1984, J.A.S.). Second, mesoscale mountain wave and severe downslope wind phenomena produce such a large local drag, that planetary scale waves can be produced. The mechanism of upscale transfer is easy to understand in this case as the standing planetary scale wave has a wavelength which depends on the mean structure of the atmosphere, and not on the width of the mountain (just as in small scale lee wave theory). An example of a theoretical description of a severe wind flow with very large drag is shown.

Smith, R. B.↗