Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear Analysis”

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 343 records · Page 19

Nonlinear dynamic analysis of flexible multibody systems

Two approaches are developed to analyze the dynamic behavior of flexible multibody systems. In the first approach each body is modeled with a modal methodology in a local non-inertial frame of reference, whereas in the second approach, each body is modeled with a finite element methodology in the inertial frame. In both cases, the interaction among the various elastic bodies is represented by constraint equations. The two approaches were compared for accuracy and efficiency: the first approach is preferable when the nonlinearities are not too strong but it becomes cumbersome and expensive to use when many modes must be used. The second approach is more general and easier to implement but could result in high computation costs for a large system. The constraints should be enforced in a time derivative fashion for better accuracy and stability.

Bauchau, Olivier A.↗

Nonlinear Dynamic Analysis of Deploying Flexible Space Booms

A fundamental investigation of the planar deployment and lock-up of two flexible boom type appendages which have attached tip masses and are connected to a central rigid body through a rotational spring is presented. Nondimensional parameters are identified and it is shown that, in general, the solution depends only on two mass ratios and one nondimensional stiffness parameter. Results are presented for boom tip deflections, deployment time and root moments at lock-up. A threshold value of the nondimensional stiffness parameter is identified beyond which boom deflections become large. Also, a thorough examination of the effect of nonlinear terms in the equations governing the deployment phase is performed. Nonlinear terms in the deployment equations due to kinematics and structural deformation are required to predict more accurately boom deflections, but retention of an inconsistent set of nonlinear terms leads to erroneous results. In particular, retaining nonlinear kinematic terms while neglecting nonlinear structural terms can produce inaccurate results even below the threshold stiffness value.

Mcgowan, P. E.↗

Identification and quantification analysis of nonlinear dynamics properties of combustion instability in a diesel engine

The cycling combustion instabilities in a diesel engine have been analyzed based on chaos theory. The objective was to investigate the dynamical characteristics of combustion in diesel engine. In this study, experiments were performed under the entire operating range of a diesel engine (the engine speed was changed from 600 to 1400 rpm and the engine load rate was from 0% to 100%), and acquired real-time series of in-cylinder combustion pressure using a piezoelectric transducer installed on the cylinder head. Several methods were applied to identify and quantitatively analyze the combustion process complexity in the diesel engine including delay-coordinate embedding, recurrence plot (RP), Recurrence Quantification Analysis, correlation dimension (CD), and the largest Lyapunov exponent (LLE) estimation. The results show that the combustion process exhibits some determinism. If LLE is positive, then the combustion system has a fractal dimension and CD is no more than 1.6 and within the diesel engine operating range. We have concluded that the combustion system of diesel engine is a low-dimensional chaotic system and the maximum values of CD and LLE occur at the lowest engine speed and load. This means that combustion system is more complex and sensitive to initial conditions and that poor combustion quality leads to the decrease of fuel economy and the increase of exhaust emissions.

42 ENGINEERING↗

Mixed models and reduced/selective integration displacement models for nonlinear shell analysis

One of the objectives of the present investigation is to identify classes of equivalent mixed models and reduced/selective integration displacement models for curved shell structures. A second objective is concerned with the identification of the spurious modes exhibited by various mixed models and their equivalent reduced integration displacement models. The merits of using mixed models are also evaluated. The mixed elements developed in the current investigation differ from those presented by Noor and Hartley (1977) by the fact that the stress resultants are discontinuous at interelement boundaries and, therefore, are eliminated on the element level. The high accuracy and effectiveness of the elements developed is demonstrated by means of numerical examples, taking into account three geometrically nonlinear problems of shallow shells.

Noor, A. K.↗

Analysis of nonlinear optical properties in donor–acceptor materials

Time-dependent density functional theory has been used to calculate nonlinear optical (NLO) properties, including the first and second hyperpolarizabilities as well as the two-photon absorption cross-section, for the donor-acceptor molecules p-nitroaniline and dimethylamino nitrostilbene, and for respective materials attached to a gold dimer. The CAMB3LYP, B3LYP, PBE0, and PBE exchange-correlation functionals all had fair but variable performance when compared to higher-level theory and to experiment. The CAMB3LYP functional had the best performance on these compounds of the functionals tested. However, our comprehensive analysis has shown that quantitative prediction of hyperpolarizabilities is still a challenge, hampered by inadequate functionals, basis sets, and solvation models, requiring further experimental characterization. Attachment of the Au{sub 2}S group to molecules already known for their relatively large NLO properties was found to further enhance the response. While our calculations show a modest enhancement for the first hyperpolarizability, the enhancement of the second hyperpolarizability is predicted to be more than an order of magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear modal analysis of penetrative convection

It is pointed out that thermal convection in many astrophysical and geophysical settings occurs in an unstable layer bounded above and below by regions which are stably stratified. The convective motions may extend a substantial distance into the adjacent stable zones. If the motions have little direct effect upon the mean stratification of the stable zone, then they are usually referred to as convective overshooting. The primary objective of the present investigation is to study the dynamics of overshooting at very large Rayleigh numbers, mainly with stellar applications in mind. Numerically this is only feasible if severe simplifications are made in the description of what are likely to be turbulent motions. The approach employed in the investigation utilizes nonlinear modal equations in which the vertical and temporal structure of the convection is described accurately at the expense of the horizontal structure. A summary of the salient properties of penetrative convection is provided on the basis of the results of the conducted studies.

Zahn, J.-P.↗

Superconvergence analysis for nonlinear parabolic equation with EQ{sub 1}{sup rot} nonconforming finite element

EQ{sub 1}{sup rot} nonconforming finite element method (FEM) applied to a class of nonlinear parabolic equation is discussed. First, by use of two typical characters of this element (one is that the associated FE interpolation operator is identical to its traditional Ritz projection operator; the other is that the consistency error is of order O(h{sup 2}), one order higher than the interpolation error, when the exact solution of the problems belongs to H{sup 3}(Ω)), the supercloseness of order O(h{sup 2}) in broken H{sup 1}-norm for semi-discrete scheme is obtained without the boundedness of numerical solution in L{sup ∞}-norm through splitting the numerical solution into several parts. Moreover, we get the desired result with requirement of u,u{sub t}∈H{sup 3}(Ω) only. Secondly, a linearized Crank–Nicolson fully discrete scheme is proposed and the superclose property of order O(h{sup 2}+τ{sup 2}) is derived by constructing a suitable auxiliary problem. Finally, a numerical example is carried out to confirm our theoretical analysis. Here, h is the subdivision parameter and τ is the time step.

97 MATHEMATICS AND COMPUTING↗

Nonlinear Aeroelastic Analysis of the HIAD TPS Coupon in the NASA 8' High Temperature Tunnel: Theory and Experiment

The purpose of this work is to develop a set of theoretical and experimental techniques to characterize the aeroelasticity of the thermal protection system (TPS) on the NASA Hypersonic Inflatable Aerodynamic Decelerator (HIAD). A square TPS coupon experiences trailing edge oscillatory behavior during experimental testing in the 8' High Temperature Tunnel (HTT), which may indicate the presence of aeroelastic flutter. Several theoretical aeroelastic models have been developed, each corresponding to a different experimental test configuration. Von Karman large deflection theory is used for the plate-like components of the TPS, along with piston theory for the aerodynamics. The constraints between the individual TPS layers and the presence of a unidirectional foundation at the back of the coupon are included by developing the necessary energy expressions and using the Rayleigh Ritz method to derive the nonlinear equations of motion. Free vibrations and limit cycle oscillations are computed and the frequencies and amplitudes are compared with accelerometer and photogrammetry data from the experiments.

Benjamin D. Goldman↗

Direct numerical simulation of electrokinetic instability and transition to chaotic motion

A new type of instability—electrokinetic instability—and an unusual transition to chaotic motion near a charge-selective surface (semiselective electric membrane, electrode, or system of micro-/nanochannels) was studied by the numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. A special finite-difference method was used for the space discretization along with a semi-implicit 31/3 -step Runge-Kutta scheme for the integration in time. Two kinds of initial conditions were considered: (a) white-noise initial conditions to mimic “room disturbances” and subsequent natural evolution of the solution, and (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. Our weakly nonlinear analysis and numerical integration of the entire system predict possibility of both kinds of bifurcations at the critical point, supercritical and subcritical, depending on the system parameters. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution, two-dimensional steady electroconvective vortices (stationary point in a proper phase space), unsteady vortices aperiodically changing their parameters (homoclinic contour), periodic motion (limit cycle), and chaotic motion. The transition to chaotic motion does not include Hopf bifurcation. The numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the experimentally observed absence of regular (near-sinusoidal) oscillations for the overlimiting regimes. There is a qualitative agreement of the experimental and the theoretical values of the threshold of instability, the dominant size of the observed coherent structures, and the experimental and theoretical volt–current characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development of Advanced Numerical Tools for the Nonlinear Seismic Analysis of Nuclear Systems

In the framework of Probabilistic Risk Assessment (PRA) of nuclear facilities, a rigorous dynamic characterization of all structural and non-structural components is of paramount concern. The common features of nuclear facilities lies in the combination of three characteristics: 1. low aspect ratio; 2. large mass; and 3. significant embedment in soil. These characteristics make the interaction between the structural system and the surrounding soil a key aspect to consider for seismic risk evaluations. The methodologies and numerical tools currently in use for simulating the seismic response of a nuclear facility with significant Soil Structure Interaction (SSI) date back to the 1970's and rely on a number of simplifying idealizations, mainly related to the characterization of the incident seismic waveforms, the assumption of linear-elastic behavior of the superstructures, and idealization of equivalent linear behavior of the supporting soils. Given this background and taking advantage of the substantial progress in computational mechanics in the last decades, the ongoing research aims at developing advanced numerical tools for the nonlinear, time history simulation of nuclear facilities. This translates into the availability of a robust nonlinear time-domain simulation tool that can accurately model the nonlinearities of coupled soil-structure systems, and whose development includes robust verification and validation. The work described herein is supported by the U.S. Department of Energy (DOE) and focuses on delivering the next-generation, high performance simulation tool in support of DOE needs and mission. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of Nonlinear Shrinkage for the Bound Metal Deposition Manufacturing using Multi-scale Approach

We consider problem of nonlinear shrinkage of the metal part during bound metal deposition manufacturing on the ground and in zero-G. To analyze this problem we developed multi-scale physics-based approach that spans atomistic dynamics at the scale of nanoseconds and the full part shrinkage at the time scale of hours. Using this approach we estimated the key parameters of the problem including grain boundary width, coefficient of surface diffusion, initial redistribution of particles during debinding stage, micro-structure evolution from round particles to densely packed grains and corresponding change of the total and chemical free energy, and sintering stress. The introduced method was used to predict shrinkage at the level of two particles, filament cross-section, sub-model, and the whole green, brown, and metal parts. To further improve accuracy and reliability of the shrinkage predictions we propose concept of intelligent additive manufacturing of metal powders in space that combines the strengths of both physics-based and data-driven methods of analysis of AM.

bound metal deposition↗

Nonlinear pattern analysis of ventricular premature beats by mutual information

The frequency of ventricular premature beats (VPBs) has been related to the risk of mortality. However, little is known about the temporal pattern of occurrence of VPBs and its relationship to autonomic activity. Hence, we applied a general correlation measure, mutual information, to quantify how VPBs are generated over time. We also used mutual information to determine the correlation between VPB production and heart rate in order to evaluate effects of autonomic activity on VPB production. We examined twenty subjects with more than 3000 VPBs/day and simulated random time series of VPB occurrence. We found that mutual information values could be used to characterize quantitatively the temporal patterns of VPB generation. Our data suggest that VPB production is not random and VPBs generated with a higher value of mutual information may be more greatly affected by autonomic activity.

NASA Discipline Regulatory Physiology↗

Solution and sensitivity analysis of nonlinear equations using a hypercomplex-variable Newton-Raphson method

Here, the classical Newton-Raphson (NR) method for solving nonlinear equations is enhanced in two ways through the use of hypercomplex variables and algebra. In particular, i) the Jacobian is computed in a highly accurate and automated way, and ii) the derivative of the solution to the nonlinear equations is computed with respect to any parameter contained within the system of equations. These advances provide two significant enhancements in that it is straightforward to provide an accurate Jacobian and to construct a reduced order model (ROM) of arbitrary order with respect to any parameter of the system. The ROM can then be used to approximate the solution for other parameter values without requiring additional solutions of the nonlinear equations. Several case studies are presented including 1D and 2D academic examples with fully functioning Python code provided. Additionally, a case of study of the catenary of an elastic cable subject to its own weight and a vertical point load. Derivatives up to 10th order were computed with respect to material, loading, and geometrical parameters. The derivatives were used to generate reduced order models of the cable deformation and reaction forces at its ends with respect to multiple input parameters. Results show that from a single hypercomplex evaluation of the cable under a single vertical point load, it is possible to generate an accurate reduced order model capable of predicting the cable deformation with 1.5 times the load in the opposite direction and with 3.5 times the load in the same direction without resolving the system of equations.

97 MATHEMATICS AND COMPUTING↗

Shaker-structure interaction modeling and analysis for nonlinear force appropriation testing

Nonlinear force appropriation is an extension of its linear counterpart where sinusoidal excitation is applied to a structure with a modal shaker and phase quadrature is achieved between the excitation and response. While a standard practice in modal testing, modal shaker excitation has the potential to alter the dynamics of the structure under test. Previous studies have been conducted to address several concerns, but this work specifically focuses on a shaker-structure interaction phenomenon which arises during the force appropriation testing of a nonlinear structure. Under pure-tone sinusoidal forcing, a nonlinear structure may respond not only at the fundamental harmonic but also potentially at sub- or superharmonics, or it can even produce aperiodic and chaotic motion in certain cases. Shaker-structure interaction occurs when the response physically pushes back against the shaker attachment, producing non-fundamental harmonic content in the force measured by the load cell, even for pure tone voltage input to the shaker. This work develops a model to replicate these physics and investigates their influence on the response of a nonlinear normal mode of the structure. Experimental evidence is first provided that demonstrates the generation of harmonic content in the measured load cell force during a force appropriation test. This interaction is replicated by developing an electromechanical model of a modal shaker attached to a nonlinear, three-mass dynamical system. Several simulated experiments are conducted both with and without the shaker model in order to identify which effects are specifically due to the presence of the shaker. Finally, the results of these simulations are then compared to the undamped nonlinear normal modes of the structure under test to evaluate the influence of shaker-structure interaction on the identified system’s dynamics.

42 ENGINEERING↗

Parabolized Stability Equations analysis of nonlinear interactions with forced eigenmodes to control subsonic jet instabilities

Nonlinear evolution of disturbances in an axisymmetric, high subsonic, high Reynolds number hot jet with forced eigenmodes is studied using the Parabolized Stability Equations (PSE) approach to understand how modes interact with one another. Both frequency and azimuthal harmonic interactions are analyzed by setting up one or two modes at higher initial amplitudes and various phases. While single mode excitation leads to harmonic growth and jet noise amplification, controlling the evolution of a specific mode has been made possible by forcing two modes (m{sub 1}, n{sub 1}), (m{sub 2}, n{sub 2}), such that the difference in azimuth and in frequency matches the desired “target” mode (m{sub 1} − m{sub 2}, n{sub 1} − n{sub 2}). A careful setup of the initial amplitudes and phases of the forced modes, defined as the “killer” modes, has allowed the minimizing of the initially dominant instability in the near pressure field, as well as its estimated radiated noise with a 15 dB loss. Although an increase of the overall sound pressure has been found in the range of azimuth and frequency analyzed, the present paper reveals the possibility to make the initially dominant instability ineffective acoustically using nonlinear interactions with forced eigenmodes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of nonlinear features associated with radio-frequency sheaths

An analytic study is conducted of nonlinear features that arise when radio-frequency (RF) sheath boundary conditions are solved self-consistently with the electromagnetic fields in a magnetized plasma. An instability associated with RF sheath–plasma interactions is identified and shown to be closely linked with sheath–plasma resonance and the occurrence of multiple solutions. It is proposed that this instability serves as a physical mechanism that enables jumping from one solution branch to another, providing a physical basis for the prediction of an RF sheath hysteresis in the time domain. It is further shown that this instability can lead to non-intuitive features such as spontaneous symmetry breaking in an otherwise completely symmetric problem. In plane-stratified geometry, two analytically solvable examples are considered for the cases of an unbounded half-space and a bounded domain

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗