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 559 records · Page 31

A cosmic ray driven instability

The interaction between energetic charged particles and thermal plasma which forms the basis of diffusive shock acceleration leads also to interesting dynamical phenomena. For a compressional mode propagating in a system with homogeneous energetic particle pressure it is well known that friction with the energetic particles leads to damping. The linear theory of this effect has been analyzed in detail by Ptuskin. Not so obvious is that a non-uniform energetic particle pressure can addition amplify compressional disturbances. If the pressure gradient is sufficiently steep this growth can dominate the frictional damping and lead to an instability. It is important to not that this effect results from the collective nature of the interaction between the energetic particles and the gas and is not connected with the Parker instability, nor with the resonant amplification of Alfven waves.

Dorfi, E. A.↗

Three axis rotational maneuver and vibration stabilization of elastic spacecraft

A control law for three-axis rotational maneuvers of a spacecraft beam-tip body configuration based on non-linear inversion and modal velocity feedback is presented. A decoupling attitude control law is presented such that in the closed-loop system the attitude angles of the spacecraft are independently controlled, using the control moments acting on the space vehicle. This controller asymptotically decouples the flexible dynamics from the rigid one and also allows the decomposition of the elastic dynamics into two subsystems representing the transverse deflections of the beam in two orthogonal planes. These low-order subsystems are used for the derivation of a modal velocity feedback stabilizer using the force and moment actuators at the end body. Simulation results are presented to show the capability of the controller.

Singh, Sahjendra N.↗

Organized motions underlying turbulent shear flows

The objective of this project is to determine the nature and significance of the organized motions underlying turbulent shear flow. There is considerable experimental evidence for the existence of such motions. In particular, one consistently observes longitudinal streaks with a spacing of about 100 in wall units in the near-wall region of wall-bounded shear flows. Recently, an analysis based on the direct resonance mechanism has predicted the appearance of streaks with precisely such a spacing. Also, the minimum channel simulations of Jimenez and Moin have given a strong dynamical significance to that spanwise length scale. They have shown that turbulent-like flows can not be maintained when the spanwise wavelength of the motion is constrained to be less than about that critical number. A critical review of the direct resonance ideas and the non-linear theory of Benney and Gustavsson is presented first. It is shown how this leads to the later mean flow-first harmonic theory of Benney. Finally, we note that a different type of analysis has led to the prediction streaks with a similar spacing. This latter approach consists of looking for optimum fields and directly provides deep insights into why a particular structure or a particular scale should be preferred.

Waleffe, F.↗

An adaptive identification and control scheme for large space structures

A unified identification and control scheme capable of achieving space at form performance objectives under nominal or failure conditions is described. Preliminary results are also presented, showing that the methodology offers much promise for effective robust control of large space structures. The control method is a multivariable, adaptive, output predictive controller called Model Predictive Control (MPC). MPC uses a state space model and input reference trajectories of set or tracking points to adaptively generate optimum commands. For a fixed model, MPC processes commands with great efficiency, and is also highly robust. A key feature of MPC is its ability to control either nonminimum phase or open loop unstable systems. As an output controller, MPC does not explicitly require full state feedback, as do most multivariable (e.g., Linear Quadratic) methods. Its features are very useful in LSS operations, as they allow non-collocated actuators and sensors. The identification scheme is based on canonical variate analysis (CVA) of input and output data. The CVA technique is particularly suited for the measurement and identification of structural dynamic processes - that is, unsteady transient or dynamically interacting processes such as between aerodynamics and structural deformation - from short, noisy data. CVA is structured so that the identification can be done in real or near real time, using computationally stable algorithms. Modeling LSS dynamics in 1-g laboratories has always been a major impediment not only to understanding their behavior in orbit, but also to controlling it. In cases where the theoretical model is not confirmed, current methods provide few clues concerning additional dynamical relationships that are not included in the theoretical models. CVA needs no a priori model data, or structure; all statistically significant dynamical states are determined using natural, entropy-based methods. Heretofore, a major limitation in applying adaptive controllers to LSS problems has been the requirement for accurate models. We feel that MPC coupled with the CVA identification scheme will resolve this limitation. The results are based on research performed under an SBIR grant, administered by the Air Force Office of Scientific Research.

Carroll, J. V.↗

Transient Vibration Prediction for Rotors on Ball Bearings Using Load-dependent Non-linear Bearing Stiffness

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus an accurate rotordynamic transient analysis requires bearing forces to be determined at each step of the transient solution. Analyses have been carried out to show the effect of accurate bearing transient forces (accounting for non-linear speed and load dependent bearing stiffness) as compared to conventional use of average rolling-element bearing stiffness. Bearing forces were calculated by COBRA-AHS (Computer Optimized Ball and Roller Bearing Analysis - Advanced High Speed) and supplied to the rotordynamics code ARDS (Analysis of Rotor Dynamic Systems) for accurate simulation of rotor transient behavior. COBRA-AHS is a fast-running 5 degree-of-freedom computer code able to calculate high speed rolling-element bearing load-displacement data for radial and angular contact ball bearings and also for cylindrical and tapered roller beatings. Results show that use of nonlinear bearing characteristics is essential for accurate prediction of rotordynamic behavior.

Fleming, David P.↗

Analytical integration of the tractions induced by non-singular dislocations on an arbitrary shaped triangular quadratic element

Here, an analytical model is proposed to evaluate the nodal force induced by a segment of dislocation upon an arbitrary shaped triangular element. This calculation is required in hybrid methods that associate dislocation dynamics to boundary or finite element to solve simultaneously the evolution of large ensembles of dislocations with complex boundary conditions. Nodal forces are defined as the triple integration of the unbalanced traction field induced by a straight dislocation upon the surface of the element. Following our previous approach (Queyreau et al 2014 Modelling Simul. Mater. Sci. Eng. 22 035004) on a simpler geometry and in the case of linear isotropic elasticity, triple integrals are solved by sequences of integration by parts that exhibit recurrence relations. The traction field is defined and finite everywhere even at the core of dislocations, thanks to the use of the non-singular stress expression formulated by Cai et al (2006 J. Mech. Phys. Solids 54 561–587). The nodal force expressions can be used when considering both a single convolution or double convolution of the Green's function with the core distribution. A solution is also proposed for the case of a semi-infinite segment through the study of the asymptotic behavior of the analytical expressions. The proposed approach is exact and very computationally efficient. The choice of an arbitrary shaped triangular element and quadratic shape functions allow the consideration of complex geometries and comply with automatic meshing procedures. These analytical expressions could also be employed to estimate dislocation interactions with interfaces.

42 ENGINEERING↗

Potential Non-Linearities in the High Latitude Circulation and Ozone Response to Stratospheric Aerosol Injection

The impacts of Stratospheric Aerosol Injection (SAI) on the atmosphere and surface climate depend on when and where the sulfate aerosol precursors are injected, as well as on how much surface cooling is to be achieved. We use a set of CESM2(WACCM6) SAI simulations achieving three different levels of global mean surface cooling and demonstrate that unlike some direct surface climate impacts driven by the reflection of solar radiation by sulfate aerosols, the SAI-induced changes in the high latitude circulation and ozone are more complex and could be non-linear. This manifests in our simulations by disproportionally larger Antarctic springtime ozone loss, significantly larger intra-ensemble spread of the Arctic stratospheric jet and ozone responses, and non-linear impacts on the extratropical modes of surface climate variability under the strongest-cooling SAI scenario compared to the weakest one. These potential non-linearities may add to uncertainties in projections of regional surface impacts under SAI.

54 ENVIRONMENTAL SCIENCES↗

Focused wave interactions with floating structures: a blind comparative study

The paper presents results from the Collaborative Computational Project in Wave Structure Interaction (CCP-WSI) Blind Test Series 2. Without prior access to the physical data, participants, with numerical methods ranging from low-fidelity linear models to fully non-linear Navier-Stokes (NS) solvers, simulate the interaction between focused wave events and two separate, taut-moored, floating structures: a hemispherical-bottomed cylinder and a cylinder with a moonpool. The ‘blind’ numerical predictions for heave, surge, pitch and mooring load, are compared against physical measurements. Dynamic time warping is used to quantify the predictive capability of participating methods. In general, NS solvers and hybrid methods give more accurate predictions; however, heave amplitude is predicted reasonably well by all methods; and a WEC-Sim implementation, with CFD-informed viscous terms, demonstrates comparable predictive capability to even the stronger NS solvers. Large variations in the solutions are observed (even among similar methods), highlighting a need for standardisation in the numerical modelling of WSI problems.

42 ENGINEERING↗

Tidal Prediction

Tides are the most predictable of oceanographic phenomena, due both to the simplicity and predictability of the astronomical forcing and to the near linearity of the ocean's dynamical response. In the classical and simplest scenario, tidal prediction is based on harmonic analysis of past measurements at a fixed location. Limits to predictability arise because isolated astronomical spectral lines are broadened into "cusps" of incoherent energy, for example through interactions with non-tidal flows. Tidal prediction at locations without past measurements has historically been a major challenge, but, owing to near-global observations of modern satellite altimeters, the empirical harmonic approach now yields reasonably accurate predictions throughout most of the open ocean. Advances in numerical modeling and data assimilation allow these predictions to be refined (especially in shallow seas where observations remain insufficient to directly constrain tidal wave structure) and extended to include tidal currents. We review recent progress in the development of global and regional-scale tidal prediction capabilities, summarize accuracy of available charts, and briefly consider outstanding issues. Satellite altimetry has also helped unravel the global tidal energy budget and has clarified the role of internal tides as a sink of tidal energy. We summarize these results, and then turn to the challenging problem of predicting internal tides. To the extent that low modes remain coherent with the surface tide, elevations can be directly mapped. We review current efforts in this direction, which are already producing charts with predictive capability. However, internal tides (especially higher modes) interact much more strongly with lower-frequency ocean flows, so a significant fraction of this tidal signal is intermittent and incoherent. We close with a brief review of ongoing efforts to model global tides in combination with wind-forced ocean motions.

Harmonic Constants↗

Transition Zones at the Changing Coastal Terrestrial‐Aquatic Interface

Coastal soils are a significant but highly uncertain component of global biogeochemical cycles. These systems experience spatial and temporal variability in biogeochemical processes, driven by marsh‐to‐upland gradients and hydrological fluctuations. These fluctuations make it difficult to understand and predict biogeochemical processes in these highly dynamic systems. We studied coastal soil biogeochemistry and its variability (a) at regional scales and (b) across transects from upland forest to marsh, in two contrasting regions—Lake Erie, a freshwater lacustrine system, and Chesapeake Bay, a saltwater estuarine system. Salinity‐related analytes were a key source of variability in soil biogeochemistry, not just in the saltwater system, but surprisingly, also in the freshwater system. We had hypothesized linear trends in biogeochemical parameters along the TAI—however, contrary to expectations, transition soils were not consistently intermediate between upland and marsh endmembers; the non‐monotonic trends of C, P, Fe along our transects suggest that these do not behave as expected and may be difficult to model and predict—thus these are key analytes to study in our regions. Rapidly changing soil factors across coastal gradients (e.g., Ca, K, CEC, and TS) may act as precursors to ecosystem shifts. Our comprehensive soil characterization represents a snapshot of a single timepoint of surface soils and provides essential data for mechanistic modeling of ecosystem dynamics across coastal transects.

Patel, Kaizad F. [Pacific Northwest National Labor↗

Guideline requirements for serviceable spacecraft grasping/berthing/docking interfaces based on simulations and flight experience (survey paper)

As space vehicles and structures become larger and more complex, the development of systems to assist humans in assembling, operating, maintaining, and performing space rescue or retrieval of these vehicles and structures becomes increasingly important. With the diversity of international spacecraft, both manned and unmanned, planned to be in orbit in the future, a set of guidelines for berthing and docking subsystems is mandatory if servicing, resupply, and retrieval is to become practical on an international level. Successful interaction between these space systems and ground and/or space-based humans requires standardized and effective operational interface designs, particularly with respect to space grasping/berthing/docking interface mechanisms. This paper defines the spacecraft mechanical interfaces necessary to create a standard dynamic envelope for joining two free-flying spacecraft in a 'hard' berth or dock with each other in space. A review was made of past space flights and dynamics simulations dating back to 1962 to obtain necessary parameters and their values for successful manually controlled and autonomous spacecraft docking/berthing. The various spacecraft docking/berthing mechanisms and concepts are illustrated along with their dynamic capture and impact tolerances including maximum contact velocity along the approach axis and in the y-z plane; capture linear misalignment tolerances; and maximum capture roll, pitch, and yaw angles. From this data sets of recommended guidelines parameters were developed for autonomous and manual impact docking tolerances, non-impact grasping/berthing tolerances (end effectors), berthing contact conditions, and alignment tolerances after rigidizing. Also, detailed requirements were developed for mechanical design interface features, as well as latching, unlatching, and separation tolerances. This data was drafted in the form of a proposed ANSI Standard guideline, reviewed, and added to by members of the committee representing several spacecraft manufacturers, NASA, and the USAF, and a consensus was reached.

Thompson, Allen B.↗

Simulation and Analysis of Technology and Operational Procedures to Reduce the Combined Effects of Emissions and Contrails

The development and evaluation of concepts and technology to support future air traffic management systems require a hierarchy of models ranging from real-time simulations to extensive field evaluations. Air traffic simulation models such as Airspace Concept Evaluation System, Center Tracon Automation System, Future Air traffic management Concept Evaluation Tool and others are used to design air traffic systems balancing the conflicting objectives of maximizing safety, meeting future demands for airports and airspace and increase efficiency of traffic flows in the presence of uncertain weather. The impact of aviation emissions and contrails on climate imposes another constraint on the design of aircraft and aviation operations. The understanding of the complex interaction between physical climate system, the carbon and other greenhouse gas emissions and aviation activity can be improved by the development of integrated assessment models that include emission and climate models together with air traffic simulations. The modeling of aircraft emissions and their interaction with each other to change the concentration levels of different gasses in the atmosphere and the resulting impact of the radiative forcing on the equilibrium of the Earth's atmosphere is complex and requires the use of coupled atmosphere-ocean general circulation models together with three-dimensional models of carbon cycle and chemistry of other non-CO2 greenhouse gases. These models are computationally intensive and unsuitable for studies involving the generation of multiple scenarios. Simple emission and climate models, based on the input-output relations of linear systems, capture the fundamental emission to climate impact behavior by careful selection of key variables and their dynamics. The impact of various greenhouse gases depends on the total concentration, effect per unit change in atmospheric concentration and the spatial distribution of the gas. All these quantities are influenced by the lifetime of the gas. The impact of a greenhouse gas depends on the interval of assessment, which may vary from a few decades to a few centuries. Climate metrics are aimed at providing a common scale to compare different greenhouse gases. If the metrics are to be used as a tool in developing and evaluating aviation operations, they should be transparent and easy to apply. Global Warming Potential and Aggregate Global Temperature Potential are some of the commonly used metrics. This paper integrates a national-level air traffic simulation and optimization capability with simple climate models and carbon cycle models, and climate metrics to assess the impact of aviation on climate. The capability brings together metrics, which are useful in aviation operations together with metrics used in climate studies. The capability can be used to make trade-offs between extra fuel cost and reduction in climate impact. There is considerable uncertainty in our understanding of the radiative forcing associated with emissions and contrails. The parameters in the simulation can be used to evaluate the effect of various uncertainties in emission models and contrails. It can also be used to evaluate the impact of different decision horizons. Alternatively, the optimization results from the simulation can be used as inputs to other tools that monetize global climate impacts like the FAA's Aviation Environmental Portfolio Management Tool for Impacts.

impact of aviation on the environment↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

Development of a non-linear simulation for generic hypersonic vehicles - ASUHS1

A nonlinear simulation is developed to model the longitudinal motion of a vehicle in hypersonic flight. The equations of motion pertinent to this study are presented. Analytic expressions for the aerodynamic forces acting on a hypersonic vehicle which were obtained from Newtonian Impact Theory are further developed. The control surface forces are further examined to incorporate vehicle elastic motion. The purpose is to establish feasible equations of motion which combine rigid body, elastic, and aeropropulsive dynamics for use in nonlinear simulations. The software package SIMULINK is used to implement the simulation. Also discussed are issues needing additional attention and potential problems associated with the implementation (with proposed solutions).

Salas, Juan↗

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NaGaS 2 : An Elusive Layered Compound with Dynamic Water Absorption and Wide‐Ranging Ion‐Exchange Properties

Abstract Most ternary sulfides belonging to the MGaS 2 structure‐type have been known for many years and are well‐characterized. Surprisingly, there have been no reports of the NaGaS 2 composition, which contains Na, a monovalent cation slightly larger in size than Li, found in LiGaS 2 , a compound known for its non‐linear optical properties. Now it is demonstrated for the first time that the unique reversible water absorption in NaGaS 2 has resulted in its absence from previous reports owing to difficulties encountered when characterizing this compound by SC XRD. The layered structure of this compound coupled with uniquely easy migration of water molecules between the layers allows for ion exchange with 3d and 5f metal cations. Some cations, for example, Ni 2+ , facilitate exfoliation of the layers, providing a facile synthetic route to a new class of 2D chalcogenide materials and furthermore demonstrating that NaGaS 2 can readily uptake uranyl species from aqueous solutions.

Klepov, Vladislav V.↗

Combining 3-Momentum and Kinetic Energy on Galilei/Newton Spacetime

Without the mass-energy equivalence available on Minkowski spacetime M, it is not possible on 4-dimensional non-relativistic Galilei/Newton spacetime G to combine 3-momentum and total mass-energy in a single tensor object. However, given a fiducial frame, it is possible to combine 3-momentum and kinetic energy into a linear form (particle) or (1,1) tensor (continuum) in a manner that exhibits increased unity of classical mechanics on flat relativistic and non-relativistic spacetimes M and G. As on M, for a material continuum on G, the first law of thermodynamics can be considered a consequence of a unified dynamical law for energy-momentum rather than an independent postulate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗