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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Blade loss transient dynamic analysis of turbomachinery

This paper reports on work completed to develop an analytical method for predicting the transient non-linear response of a complete aircraft engine system due to the loss of a fan blade, and to validate the analysis by comparing the results against actual blade loss test data. The solution, which is based on the component element method, accounts for rotor-to-casing rubs, high damping and rapid deceleration rates associated with the blade loss event. A comparison of test results and predicted response show good agreement except for an initial overshoot spike not observed in test. The method is effective for analysis of large systems.

Stallone, M. J.↗

NASA NTRS2020

Squash-Box Feasibility Driven Differential Dynamic Programming

Recently, Differential Dynamic Programming (DDP) and other similar algorithms have become the solvers of choice when performing non-linear Model Predictive Control (nMPC) with modern robotic devices. The reason is that they have a lower computational cost per iteration when compared with off-the-shelf Non-Linear Programming (NLP) solvers, which enables its online operation. However, they cannot handle constraints, and are known to have poor convergence capabilities. In this paper, we propose a method to solve the optimal control problem with control bounds through a squashing function (i.e., a sigmoid, which is bounded by construction). It has been shown that a naive use of squashing functions damage the convergence rate. To tackle this, we first propose to add a quadratic barrier that avoids the difficulty of the plateau produced by the sigmoid. Second, we add an outer loop that adapts both the sigmoid and the barrier; it makes the optimal control problem with the squashing function converge to the original control-bounded problem. To validate our method, we present simulation results for different types of platforms including a multi-rotor, a biped, a quadruped and a humanoid robot.

Navarro, Angel Santamaria↗

DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2006

An Agent-Based Model for Analyzing Control Policies and the Dynamic Service-Time Performance of a Capacity-Constrained Air Traffic Management Facility

Simple agent-based models may be useful for investigating air traffic control strategies as a precursory screening for more costly, higher fidelity simulation. Of concern is the ability of the models to capture the essence of the system and provide insight into system behavior in a timely manner and without breaking the bank. The method is put to the test with the development of a model to address situations where capacity is overburdened and potential for propagation of the resultant delay though later flights is possible via flight dependencies. The resultant model includes primitive representations of principal air traffic system attributes, namely system capacity, demand, airline schedules and strategy, and aircraft capability. It affords a venue to explore their interdependence in a time-dependent, dynamic system simulation. The scope of the research question and the carefully-chosen modeling fidelity did allow for the development of an agent-based model in short order. The model predicted non-linear behavior given certain initial conditions and system control strategies. Additionally, a combination of the model and dimensionless techniques borrowed from fluid systems was demonstrated that can predict the system s dynamic behavior across a wide range of parametric settings.

Conway, Sheila R.↗

DOE OSTI2022

Diurnal and Seasonal Dynamics of Solar-Induced Chlorophyll Fluorescence, Vegetation Indices, and Gross Primary Productivity in the Boreal Forest

Remote sensing of solar-induced chlorophyll fluorescence (SIF) provides a powerful proxy for gross primary productivity (GPP). It is particularly promising in boreal ecosystems where seasonal downregulation of photosynthesis occurs without significant changes in canopy structure or chlorophyll content. The use of SIF as a proxy for GPP is complicated by inherent non-linearities due to both physical (illumination effects) and ecophysiological (light use efficiencies) controls at fine spatial (tower/leaf) and temporal (half-hourly) scales. Here, to study the SIF-GPP relationship, we investigated the diurnal and seasonal dynamics of continuous tower-based measurements of SIF, GPP, and common vegetation indices at the Southern Old Black Spruce Site (SOBS) in Saskatchewan, CA over the course of two years. We find that SIF outperforms other vegetation indices as a proxy for GPP at all temporal scales but shows a non-linear relationship with GPP at a half-hourly resolution. At small temporal scales, SIF and GPP are predominantly driven by light and non-linearity between SIF and GPP is due to the light saturation of GPP. Averaged over daily and monthly scales, the relationship between SIF and GPP is linear due to a reduction in the observed PAR range. Seasonal changes in the light responses of SIF and GPP are driven by changes in light use efficiency which co-vary with changes in temperature, while illumination and canopy structure partially linearize the SIF-GPP relationship. Additionally, we find that the SIF-GPP relationship has a seasonal dependency. Our results help clarify the utility of SIF for estimating carbon assimilation in boreal forests.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1997

An Experimental Investigation on the Static Equilibria and Dynamics of Liquid Bridges

A liquid bridge is a volume of liquid held between two or more solid supports. In the case of small disk supports with a sharp edge, the contact line between the bridge and the support disk will be anchored along the edge of the disk. For these cases the solid presents a geometrical singularity and the contact angle is indeterminate within a given range. This dissertation presents research conducted on liquid bridges with anchored contact lines. The three major topics covered are: determining the role of support geometry on static equilibria, liquid bridge dynamical behavior, and forces exerted by a liquid bridge on a support structure. The work was primarily experimental and conducted in a "Plateau tank" that allowed for the simulation of equivalent low-gravity conditions. The main thrust of the experimental work involved the use of a high resolution optical measurement system for imaging the dynamic zone shape, measurement of the static and dynamic contact angles and non-invasive analysis of excited surface modes. The liquid bridge was manipulated by computer controlled linear actuators which allowed precise control over the physical characteristics of the bridge. Experiments have been carried out to locate a bifurcation point along the maximum volume axisymmetric stability margin. Below the critical slenderness the bifurcation from an axisymmetric to a stable nonaxisymmetric configuration is supercritical. However, above this critical slenderness, the bifurcation is subcritical. A series of experiments analyzed the effect on axisymmetric bridge stability by using support disks of different radii, The shape behavior as transition points were approached, as well as the limiting case of a vanishing support radius was investigated. Experiments were performed to determine the resonant frequencies of axisymmetric bridges subject to lateral vibrations. Anomolous results led to a series of experiments to characterize nonlinearities present in the dynamic bridge shape. Finally, an attempt was made to experimentally measure the force exerted by the bridge on the lower support disk. This was done through use of a force balance apparatus. Particular attention was paid to the behavior of the bridge as the minimum volume stability limit was approached.

Resnick, Andrew Howard↗

DOE OSTI2020

High-Fidelity Large-Signal Order Reduction Approach for Composite Load Model

With the increasing penetration of electronic loads and distributed energy resources, conventional load models cannot capture their dynamics. Therefore, a new comprehensive composite load model is developed by Western Electricity Coordinating Council (WECC). However, this model is a complex high-order non-linear system with multi-time-scale property, which poses challenges on power system studies with heavy computational burden. In order to reduce the model complexity, the authors firstly develop a large-signal order reduction (LSOR) method using singular perturbation theory. In this method, the fast dynamics are integrated into the slow ones to preserve transient characteristics of the former. Then, accuracy assessment conditions are proposed and embedded into the LSOR to improve and guarantee the accuracy of reduced-order model. Finally, the reduced-order WECC composite load model is derived by using the proposed algorithm. Overall, simulation results show that the reduced-order large-signal model significantly alleviates the computational burden while maintaining similar dynamic responses as the original composite load model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1985

The HEAO-2 Guest Investigator Program: Non-linear growth of instabilities in line-driven stellar winds

The linear instability of line-driven stellar winds to take proper account of the dynamical effect of scattered radiation were analyzed. It is found that: (1) the drag effect of the mean scattered radiation does greatly reduce the contribution of scattering lines to the instability at the very base of the wind, but the instability growth rate associated with such lines rapidly increases as the flow moves outward from the base, reaching more than 50% of the growth rate for pure absorption lines within a stellar radius of the surface, and eventually reaching 80% of that rate at large radii; (2) perturbations in the scattered radiation field may be important for the propagation of wind disturbances, but they have little effect on the wind instability; and (3) the contribution of strongly shadowed lines to the wind instability is often reduced compared to that of unshadowed lines, but their overall effect is not one of damping in the outer parts of the wind. It is concluded that, even when all scattering effects are taken into account, the bulk of the flow in a line-driven stellar wind is still highly unstable.

Rybicki, G. B.↗

DOE OSTI2023

Progress with applications of three-ion ICRF scenarios for fusion research: A review

The viability of magnetic confinement fusion as an energy source depends on achieving the high ion temperatures required for D-T fusion. Among the available techniques, plasma heating with waves in the ion cyclotron range of frequencies (ICRF) is a prominent method for bulk ion heating in fusion plasmas. Furthermore, a detailed understanding of the non-linear physics of alpha heating and the complex impact of MeV-range fast ions on plasma dynamics becomes progressively more important. This paper provides a comprehensive overview of recent developments with the three-ion ICRF scenarios on Alcator C-Mod, ASDEX Upgrade and JET tokamaks. The results demonstrate the flexibility of these novel scenarios for heating bulk ions in D-T ≈ 50%-50% plasmas and efficient generation of MeV-range fast ions in multi-ion species plasmas. In conclusion, several key results relevant for ITER and future fusion reactors are highlighted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1979

Stochastic modeling of the time-averaged equations for climate dynamics

Two analyses on climate dynamics are presented, both based on a simplified set of barotropic equations for representing large scale non-linear atmospheric circulation characteristics. The statistical properties of the set of equations were investigated for small values of a parameter k sub o, corresponding to the physical case of large zonal westerlies and relatively weak eddy motions. The effect of seasonal type forcing on the solution of these equations was also studied and some preliminary numerical results are presented.

Spreiter, J. R.↗

NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

DOE OSTI2023

Social network structure and the spread of complex contagions from a population genetics perspective

Ideas, behaviors, and opinions spread through social networks. If the probability of spreading to a new individual is a non-linear function of the fraction of the individuals’ affected neighbors, such a spreading process becomes a “complex contagion”. This non-linearity does not typically appear with physically spreading infections, but instead can emerge when the concept that is spreading is subject to game theoretical considerations (e.g. for choices of strategy or behavior) or psychological effects such as social reinforcement and other forms of peer influence (e.g. for ideas, preferences, or opinions). Here we study how the stochastic dynamics of such complex contagions are affected by the underlying network structure. Motivated by simulations of complex contagions on real social networks, we present a framework for analyzing the statistics of contagions with arbitrary non-linear adoption probabilities based on the mathematical tools of population genetics. The central idea is to use an effective lower-dimensional diffusion process to approximate the statistics of the contagion. This leads to a tradeoff between the effects of ”selection” (microscopic tendencies for an idea to spread or die out), random drift, and network structure. Our framework illustrates intuitively several key properties of complex contagions: stronger community structure and network sparsity can significantly enhance the spread, while broad degree distributions dampen the effect of selection compared to random drift. Finally, we show that some structural features can exhibit critical values that demarcate regimes where global contagions become possible for networks of arbitrary size. Our results draw parallels between the competition of genes in a population and memes in a world of minds and ideas. Our tools provide insight into the spread of information, behaviors, and ideas via social influence, and highlight the role of macroscopic network structure in determining their fate.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1988

Real-time optical image processing techniques

Nonlinear real-time optical processing on spatial pulse frequency modulation has been pursued through the analysis, design, and fabrication of pulse frequency modulated halftone screens and the modification of micro-channel spatial light modulators (MSLMs). Micro-channel spatial light modulators are modified via the Fabry-Perot method to achieve the high gamma operation required for non-linear operation. Real-time nonlinear processing was performed using the halftone screen and MSLM. The experiments showed the effectiveness of the thresholding and also showed the needs of higher SBP for image processing. The Hughes LCLV has been characterized and found to yield high gamma (about 1.7) when operated in low frequency and low bias mode. Cascading of two LCLVs should also provide enough gamma for nonlinear processing. In this case, the SBP of the LCLV is sufficient but the uniformity of the LCLV needs improvement. These include image correlation, computer generation of holograms, pseudo-color image encoding for image enhancement, and associative-retrieval in neural processing. The discovery of the only known optical method for dynamic range compression of an input image in real-time by using GaAs photorefractive crystals is reported. Finally, a new architecture for non-linear multiple sensory, neural processing has been suggested.

Liu, Hua-Kuang↗

DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A sui generis whipping-instability-based self-sequencing multi-monodisperse 2D spray from an anisotropic microfluidic liquid jet device

Well-defined aerosols pave the way for versatile basic and applied research. Here, we demonstrate a unique whipping instability that generates from a high-aspect-ratio microfluidic device resulting in a unique steady-state gas-focused whipping jet (WJ) without any need for electrification. This WJ device emanates a multi-monodisperse whipping spray jet with a two-dimensional (2D) profile. We demonstrate this phenomenon based on various fluidic parameters theoretically and experimentally. The 2D WJ’s unique behavior is derived using analytical fluid dynamics to explain jet diameter, whipping regime, and spreading angle. The phenomenon is further characterized experimentally by measuring the angle with respect to the flow rate, the distances between droplets, the droplet shapes, and the reproducibility of these parameters. We also explain the precise fabrication of such inexpensive devices. Lastly, we highlight these devices’ potential use as sample environments in versatile applications ranging from cryoelectron microscopy over mass spectrometry to drug formulation and structural studies at X-ray free-electron lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Benchmark calculations of radiative forcing by greenhouse gases

Changes in concentrations of greenhouse gases lead to changes in radiative fluxes throughout the atmosphere. The value of this change, the instantaneous radiative forcing, varies across climate models, due partly to differences in the distribution of clouds, humidity, and temperature across models, and partly due to errors introduced by approximate treatments of radiative transfer. This paper describes an experiment within the Radiative Forcing Model Intercomparision Project that uses benchmark calculations made with line–by–line models to identify parameterization error in the representation of absorption and emission by greenhouse gases. Clear–sky instantaneous forcing by greenhouse gases is computed using a set of 100 profiles, selected from a re–analysis of present–day conditions, that represent the global annual mean forcing from pre–industrial times to the present–day with sampling errors of less than 0.01 W m –2 . Six contributing line–by–line models agree in their estimate of this forcing to within 0.025 W m –2 while even recently–developed parameterizations have typical errors four or more times larger, suggesting both that the samples reveal true differences among line–by–line models and that parameterization error will be readily identifiable. Agreement among line–by–line models is better in the longwave than in the shortwave where differing treatments of the water vapor vapor continuum affect estimates of forcing by carbon dioxide and methane. The impacts of clouds on instantaneous radiative forcing are estimated from climate model simulations. Here, the adjustment due to stratospheric temperature change by assuming fixed dynamical heating. Adjustments are large only for ozone and for carbon dioxide, for which stratospheric cooling introduces modest non–linearity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Dependence of the impurity transport on the dominant turbulent regime in ELM-y H-mode discharges on the DIII-D tokamak

Laser blow-off injections of aluminum and tungsten have been performed on the DIII-D tokamak to investigate the variation of impurity transport in a set of dedicated ion and electron heating scans with a fixed value of the external torque. The particle transport is quantified via the Bayesian inference method, which, constrained by a combination of a charge exchange recombination spectroscopy, soft X-ray measurements, and VUV spectroscopy provides a detailed uncertainty quantification of transport coefficients. Contrasting discharge phases with a dominant electron and ion heating reveal 3-fold drop in the impurity confinement time and order of magnitude increase in midradius impurity diffusion, when additional electron heating is applied. Further, the calculated stationary aluminum density profiles reverse from peaked in electron heated to hollow in the ion heated case, following a similar trend as electron and carbon density. Comparable values of a core diffusion have been observed for W and Al ions, while differences in the propagation dynamics of these impurities are attributed to pedestal and edge transport. Modeling of the core transport with non-linear gyrokinetics code CGYRO, significantly underpredicts the magnitude of the variation in Al transport. Diffusion increases 3-times steeper with additional electron heat flux and 10-times lower diffusion is observed in ion heated case than predicted by the modeling. Here, the CGYRO model quantitatively matches the increase of the Al diffusion when approaching the linear threshold for the transition from the ion temperature gradient (ITG) to trapped electron mode (TEM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗