Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear modeling”

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 361 records · Page 20

Extended MHD modeling of nonlinear instabilities in fusion and space plasmas

A number of different sub-projects where pursued within this DOE early career project. The primary focus was on using fully nonlinear, curvilinear, extended MHD simulations of instabilities with applications to fusion and space plasmas. In particular, we performed comprehensive studies of the dynamics of the double tearing mode in different regimes and confi gurations, using Cartesian and cyclindrical geometry and investigating both linear and non-linear dynamics. In addition to traditional extended MHD involving Hall term and electron pressure gradient, we also employed a new multi-fluid moment model, which shows great promise to incorporate kinetic effects, in particular off-diagonal elements of the pressure tensor, in a fluid model, which is naturally computationally much cheaper than fully kinetic particle or Vlasov simulations. We used our Vlasov code for detailed studies of how weak collisions effect plasma echos. In addition, we have played an important supporting role working with the PPPL theory group around Will Fox and Amitava Bhattacharjee on providing simulation support for HED plasma experiments performed at high-powered laser facilities like OMEGA-EP in Rochester, NY. This project has support a great number of computational advances in our fluid and kinetic plasma models, and has been crucial to winning multiple INCITE computer time awards that supported our computational modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HALO: A full-orbit model of nonlinear interaction of fast particles with eigenmodes

HALO (HAgis LOcust) solves the initial value Vlasov-Maxwell problem perturbatively for application to certain nonlinear wave-particle problems in tokamak plasmas. It uses the same basic approach as the HAGIS code for wave evolution but is built on the LOCUST-GPU full-orbit code for the solution of the Hamiltonian fast particle motion in cylindrical coordinates. The wave amplitude and particle evolution include all finite Larmour radius effects. We describe and benchmark the currently implemented Alfvén eigenmode workflow, demonstrating correct particle motion, linear and nonlinear power transfer. Here, the formulation and numerical scheme are sufficiently general as to allow easy future implementation of different kinds of eigenmodes, such as modes close to the ion-cyclotron frequency. The code can model multiple eigenmodes and multiple fast ion species simultaneously, and supports the general form of the equilibrium distribution in constants of motion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nonlinear instability of a forced baroclinic Rossby wave

The nonlinear instability of forced baroclinic Rossby waves to finite-amplitude perturbations was investigated using a simple two-layer quasi-geostrophic weakly nonlinear model and a low-order spectral model, developed for this purpose. The results were applied to a study of the interaction of planetary-scale stationary eddies with synoptic-scale transient eddies. A comparison of the energetics of the weakly nonlinear solution with that seen in observational analysis and in GCM experiments indicate three basic properties of the interaction between the two types of eddies: (1) eddy available potential energy is transferred from the planetary-scale stationary eddies to the synoptic-scale transient eddies, (2) eddy kinetic energy is transferred in the opposite direction, and (3) the eddy available potential energy transfer dominates the transfer of eddy kinetic energy. It is concluded that the simple weakly nonlinear model used in these studies does indeed capture much of the fundamental behavior of wave-wave interactions seen in the atmosphere.

Feldstein, Steven B.↗

Nonlinear Shell Modeling of Thin Membranes with Emphasis on Structural Wrinkling

Thin solar sail membranes of very large span are being envisioned for near-term space missions. One major design issue that is inherent to these very flexible structures is the formation of wrinkling patterns. Structural wrinkles may deteriorate a solar sail's performance and, in certain cases, structural integrity. In this paper, a geometrically nonlinear, updated Lagrangian shell formulation is employed using the ABAQUS finite element code to simulate the formation of wrinkled deformations in thin-film membranes. The restrictive assumptions of true membranes, i.e. Tension Field theory (TF), are not invoked. Two effective modeling strategies are introduced to facilitate convergent solutions of wrinkled equilibrium states. Several numerical studies are carried out, and the results are compared with recent experimental data. Good agreement is observed between the numerical simulations and experimental data.

Tessler, Alexander↗

Relationships between nonlinear normal modes and response to random inputs

The ability to model nonlinear structures subject to random excitation is of key importance in designing hypersonic aircraft and other advanced aerospace vehicles. When a structure is linear, superposition can be used to construct its response to a known spectrum in terms of its linear modes. Superposition does not hold for a nonlinear system, but several works have shown that a system's dynamics can still be understood qualitatively in terms of its nonlinear normal modes (NNMs). Here, this work investigates the connection between a structure's undamped nonlinear normal modes and the spectrum of its response to high amplitude random forcing. Two examples are investigated: a spring-mass system and a clamped-clamped beam modeled within a geometrically nonlinear finite element package. In both cases, an intimate connection is observed between the smeared peaks in the response spectrum and the frequency-energy dependence of the nonlinear normal modes. In order to understand the role of coupling between the underlying linear modes, reduced order models with and without modal coupling terms are used to separate the effect of each NNM's backbone from the nonlinear couplings that give rise to internal resonances. In the cases shown here, uncoupled, single-degree-of-freedom nonlinear models are found to predict major features in the response with reasonable accuracy; a highly inexpensive approximation such as this could be useful in design and optimization studies. More importantly, the results show that a reduced order model can be expected to give accurate results only if it is also capable of accurately predicting the frequency-energy dependence of the nonlinear modes that are excited.

42 ENGINEERING↗

On the Importance of the Dynamics of Discretizations

It has been realized recently that the discrete maps resulting from numerical discretizations of differential equations can possess asymptotic dynamical behavior quite different from that of the original systems. This is the case not only for systems of Ordinary Differential Equations (ODEs) but in a more complicated manner for Partial Differential Equations (PDEs) used to model complex physics. The impact of the modified dynamics may be mild and even not observed for some numerical methods. For other classes of discretizations the impact may be pronounced, but not always obvious depending on the nonlinear model equations, the time steps, the grid spacings and the initial conditions. Non-convergence or convergence to periodic solutions might be easily recognizable but convergence to incorrect but plausible solutions may not be so obvious - even for discretized parameters within the linearized stability constraint. Based on our past four years of research, we will illustrate some of the pathology of the dynamics of discretizations, its possible impact and the usage of these schemes for model nonlinear ODEs, convection-diffusion equations and grid adaptations.

Sweby, Peter K.↗

Nonlinear modal resonances in low-gravity slosh-spacecraft systems

Nonlinear models of low gravity slosh, when coupled to spacecraft vibrations, predict intense nonlinear eigenfrequency shifts at zero gravity. These nonlinear frequency shifts are due to internal quadratic and cubic resonances between fluid slosh modes and spacecraft vibration modes. Their existence has been verified experimentally, and they cannot be correctly modeled by approximate, uncoupled nonlinear models, such as pendulum mechanical analogs. These predictions mean that linear slosh assumptions for spacecraft vibration models can be invalid, and may lead to degraded control system stability and performance. However, a complete nonlinear modal analysis will predict the correct dynamic behavior. This paper presents the analytical basis for these results, and discusses the effect of internal resonances on the nonlinear coupled response at zero gravity.

Peterson, Lee D.↗

Theory of multiple quantum dot formation in strained-layer heteroepitaxy

We develop a theory for the experimentally observed formation of multiple quantum dots (QDs) in strained-layer heteroepitaxy based on surface morphological stability analysis of a coherently strained epitaxial thin film on a crystalline substrate. Using a fully nonlinear model of surface morphological evolution that accounts for a wetting potential contribution to the epitaxial film's free energy as well as surface diffusional anisotropy, we demonstrate the formation of multiple QD patterns in self-consistent dynamical simulations of the evolution of the epitaxial film surface perturbed from its planar state. The simulation predictions are supported by weakly nonlinear analysis of the epitaxial film surface morphological stability. We find that, in addition to the Stranski-Krastanow instability, long-wavelength perturbations from the planar film surface morphology can trigger a nonlinear instability, resulting in the splitting of a single QD into multiple QDs of smaller sizes, and predict the critical wavelength of the film surface perturbation for the onset of the nonlinear tip-splitting instability. The theory provides a fundamental interpretation for the observations of “QD pairs” or “double QDs” and other multiple QDs reported in experimental studies of epitaxial growth of semiconductor strained layers and sets the stage for precise engineering of tunable-size nanoscale surface features in strained-layer heteroepitaxy by exploiting film surface nonlinear, pattern forming phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development and validation of a general purpose linearization program for rigid aircraft models

A FORTRAN program that provides the user with a powerful and flexible tool for the linearization of aircraft models is discussed. The program LINEAR numerically determines a linear systems model using nonlinear equations of motion and a user-supplied, nonlinear aerodynamic model. The system model determined by LINEAR consists of matrices for both the state and observation equations. The program has been designed to allow easy selection and definition of the state, control, and observation variables to be used in a particular model. Also, included in the report is a comparison of linear and nonlinear models for a high performance aircraft.

Duke, E. L.↗

Development and validation of a general purpose linearization program for rigid aircraft models

This paper discusses a FORTRAN program that provides the user with a powerful and flexible tool for the linearization of aircraft models. The program LINEAR numerically determines a linear systems model using nonlinear equations of motion and a user-supplied, nonlinear aerodynamic model. The system model determined by LINEAR consists of matrices for both the state and observation equations. The program has been designed to allow easy selection and definition of the state, control, and observation variables to be used in a particular model. Also, included in the report is a comparison of linear and nonlinear models for a high-performance aircraft.

Duke, E. L.↗

Updating Lunar Subsurface Geotechnical Properties Based Upon Nonlinear Winkler Modeling of Apollo and Lunokhod Measurements

Returning human and robotic explorers to the lunar surface has identified a need for a renewed understanding of lunar geotechnical properties related to landing on and exploring the lunar surface. We review fundamental lunar geotechnical properties such as bulk density, cohesion, friction angle and shear strength, and apply new data derived from Apollo and Lunokhod geotechnical measurements to update the understanding of lunar regolith behavior. The analysis derives a new, nonlinear relationship of lunar geotechnical properties as a function of depth.

Lunar Geotechnical engineering↗

Modeling the nonlinear hysteretic response in DAE experiments of Berea sandstone: A case-study

Dynamic acousto-elasticity (DAE) allows probing the instantaneous state of a material while the latter slowly and periodically is changed by an external, dynamic source. In DAE investigations of geo-materials, hysteresis of the material's modulus defect displays intriguing features which have not yet been interpreted in terms of any specific mechanism occurring at atomic or mesoscale. Here, experimental results on dry Berea sandstone, which is the rock type best investigated by means of a DAE technique, are analyzed in terms of three rheological models providing simplified representations of mechanisms involving dislocations interacting with point defects which are distributed along the dislocations' core or glide planes, and microcracks with finite stiffness in compression. Constitutive relations linking macroscopic strain and stress are derived. From the latter, the modulus defect associated to each mechanism is recovered. These models are employed to construct a composite one which is capable of reproducing several of the main features observed in the experimental data. The limitations of the present approach and, possibly, of the current implementation of DAE are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Precession Index and a Nonlinear Energy Balance Climate Model

A simple nonlinear energy balance climate model yields a precession index-like term in the temperature. Despite its importance in the geologic record, the precession index e sin (Omega)S, where e is the Earth's orbital eccentricity and (Omega)S is the Sun's perigee in the geocentric frame, is not present in the insolation at the top of the atmosphere. Hence there is no one-for-one mapping of 23,000 and 19,000 year periodicities from the insolation to the paleoclimate record; a nonlinear climate model is needed to produce these long periods. A nonlinear energy balance climate model with radiative terms of form T n, where T is surface temperature and n less than 1, does produce e sin (omega)S terms in temperature; the e sin (omega)S terms are called Seversmith psychroterms. Without feedback mechanisms, the model achieves extreme values of 0.64 K at the maximum orbital eccentricity of 0.06, cooling one hemisphere while simultaneously warming the other; the hemisphere over which perihelion occurs is the cooler. In other words, the nonlinear energy balance model produces long-term cooling in the northern hemisphere when the Sun's perihelion is near northern summer solstice and long-term warming in the northern hemisphere when the aphelion is near northern summer solstice. (This behavior is similar to the inertialess gray body which radiates like T 4, but the amplitude is much lower for the energy balance model because of its thermal inertia.) This seemingly paradoxical behavior works against the standard Milankovitch model, which requires cool northern summers (Sun far from Earth in northern summer) to build up northern ice sheets, so that if the standard model is correct it must be more efficient than previously thought. Alternatively, the new mechanism could possibly be dominant and indicate southern hemisphere control of the northern ice sheets, wherein the southern oceans undergo a long-term cooling when the Sun is far from the Earth during northern summer. The cold water eventually flows north, cooling the northern hemisphere. This might explain why the northern oceans lag the southern ones when it comes to orbital forcing.

Rubincam, David↗

A nonlinear steady model for moist hydrostatic mountain waves

The dynamics of hydrostatic gravity waves generated by the passage of a steady, stably stratified, moist flow over a two-dimensional topography is considered. Coriolis effects are neglected. The cloud region is determined by the dynamics, and within that region the Brunt-Vaisala frequency takes on a value smaller than the outside value. In both the dry and cloudy regions the Brunt-Vaisala frequency is constant with height. The moist layer is considered to be either next to the mountain or at midlevels and to be deep enough so that an entire cloud forms in that layer. The nonlinearity in the flow and lower boundary affects the dynamics of these waves and wave drag. The latter is found to depend upon: (1) the location of the moist layer with respect to the ground, (2) the amount of moisture, (3) the degree of nonlinearity and (4) the departure from symmetry in the bottom topography.

Barcilon, A.↗

Closed-Loop Evaluation of an Integrated Failure Identification and Fault Tolerant Control System for a Transport Aircraft

Formal robustness analysis of aircraft control upset prevention and recovery systems could play an important role in their validation and ultimate certification. Such systems developed for failure detection, identification, and reconfiguration, as well as upset recovery, need to be evaluated over broad regions of the flight envelope or under extreme flight conditions, and should include various sources of uncertainty. To apply formal robustness analysis, formulation of linear fractional transformation (LFT) models of complex parameter-dependent systems is required, which represent system uncertainty due to parameter uncertainty and actuator faults. This paper describes a detailed LFT model formulation procedure from the nonlinear model of a transport aircraft by using a preliminary LFT modeling software tool developed at the NASA Langley Research Center, which utilizes a matrix-based computational approach. The closed-loop system is evaluated over the entire flight envelope based on the generated LFT model which can cover nonlinear dynamics. The robustness analysis results of the closed-loop fault tolerant control system of a transport aircraft are presented. A reliable flight envelope (safe flight regime) is also calculated from the robust performance analysis results, over which the closed-loop system can achieve the desired performance of command tracking and failure detection.

Shin, Jong-Yeob↗

Recurrence in truncated Boussinesq models for nonlinear waves in shallow water

The rapid spatial recurrence of weakly nonlinear and weakly dispersive progressive shallow-water waves is examined using a numerical integration technique on the discretized and truncated form of the Boussinesq equations. This study primarily examines recurrence in wave fields with Ursell number O(1) and characterizes the sensitivity of recurrence to initial spectral shape and number of allowed frequency modes. It is shown that the rapid spatial recurrence is not an inherent property of the considered Boussinesq systems for evolution distances of 10-50 wavelengths. The main result of the study is that highly truncated Boussinesq models of resonant shallow-water ocean surface gravity waves predict rapid multiple recurrence cycles, but that this is an artifact dependent on the number of allowed modes. For initial conditions consisting of essentially all energy concentrated in a single mode, damping of the recurrence cycles increases as the number of low-power background modes increases. When more than 32 modes are allowed, the recurrence behavior is relatively insensitive to the number of allowed modes.

Elgar, Steve↗

Lidar-Enhanced Wind Turbine Control: Past, Present, and Future

The main challenges in harvesting energy from the wind arise from the unknown incoming turbulent wind field. Balancing the competing interests of reduction in structural loads and increasing energy production is the goal of a wind turbine controller to reduce the cost of producing wind energy. Conventional wind turbines use feedback methods to optimize these goals, reacting to wind disturbances after they have already impacted the wind turbine. Lidar sensors offer a means to provide additional inputs to a wind turbine controller, enabling new techniques to improve control methods, allowing a controller to actuate a wind turbine in anticipation of an incoming wind disturbance. This paper will look at the development of lidar-enhanced controls and how they have been used for various turbine load reductions with pitch actuation, as well as increased energy production with improved yaw control. Ongoing work will also be discussed to show that combining pitch and torque control using feedforward nonlinear model predictive control can lead to both reduced loads and increased energy production. Future work is also proposed on extending individual wind turbine controls to the wind plant level and determining how lidars can be used for control methods to further lower the cost of wind energy by minimizing wake impacts in a wind farm.

17 WIND ENERGY↗

Cosmic shear results from the deep lens survey. II. Full cosmological parameter constraints from tomography

We present a tomographic cosmic shear study from the Deep Lens Survey (DLS), which, providing a limiting magnitude r{sub lim}∼27 (5σ), is designed as a precursor Large Synoptic Survey Telescope (LSST) survey with an emphasis on depth. Using five tomographic redshift bins, we study their auto- and cross-correlations to constrain cosmological parameters. We use a luminosity-dependent nonlinear model to account for the astrophysical systematics originating from intrinsic alignments of galaxy shapes. We find that the cosmological leverage of the DLS is among the highest among existing >10 deg{sup 2} cosmic shear surveys. Combining the DLS tomography with the 9 yr results of the Wilkinson Microwave Anisotropy Probe (WMAP9) gives Ω{sub m}=0.293{sub −0.014}{sup +0.012}, σ{sub 8}=0.833{sub −0.018}{sup +0.011}, H{sub 0}=68.6{sub −1.2}{sup +1.4} km s{sup −1} Mpc{sup −1}, and Ω{sub b}=0.0475±0.0012 for ΛCDM, reducing the uncertainties of the WMAP9-only constraints by ∼50%. When we do not assume flatness for ΛCDM, we obtain the curvature constraint Ω{sub k}=−0.010{sub −0.015}{sup +0.013} from the DLS+WMAP9 combination, which, however, is not well constrained when WMAP9 is used alone. The dark energy equation-of-state parameter w is tightly constrained when baryonic acoustic oscillation (BAO) data are added, yielding w=−1.02{sub −0.09}{sup +0.10} with the DLS+WMAP9+BAO joint probe. The addition of supernova constraints further tightens the parameter to w=−1.03±0.03. Our joint constraints are fully consistent with the final Planck results and also with the predictions of a ΛCDM universe.

79 ASTRONOMY AND ASTROPHYSICS↗