Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinearity”

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 649 records · Page 36

Recent Applications of Higher-Order Spectral Analysis to Nonlinear Aeroelastic Phenomena

Recent applications of higher-order spectral (HOS) methods to nonlinear aeroelastic phenomena are presented. Applications include the analysis of data from a simulated nonlinear pitch and plunge apparatus and from F-18 flight flutter tests. A MATLAB model of the Texas A&MUniversity s Nonlinear Aeroelastic Testbed Apparatus (NATA) is used to generate aeroelastic transients at various conditions including limit cycle oscillations (LCO). The Gaussian or non-Gaussian nature of the transients is investigated, related to HOS methods, and used to identify levels of increasing nonlinear aeroelastic response. Royal Australian Air Force (RAAF) F/A-18 flight flutter test data is presented and analyzed. The data includes high-quality measurements of forced responses and LCO phenomena. Standard power spectral density (PSD) techniques and HOS methods are applied to the data and presented. The goal of this research is to develop methods that can identify the onset of nonlinear aeroelastic phenomena, such as LCO, during flutter testing.

Silva, Walter A.↗

The Effect of Basis Selection on Static and Random Acoustic Response Prediction Using a Nonlinear Modal Simulation

An investigation of the effect of basis selection on geometric nonlinear response prediction using a reduced-order nonlinear modal simulation is presented. The accuracy is dictated by the selection of the basis used to determine the nonlinear modal stiffness. This study considers a suite of available bases including bending modes only, bending and membrane modes, coupled bending and companion modes, and uncoupled bending and companion modes. The nonlinear modal simulation presented is broadly applicable and is demonstrated for nonlinear quasi-static and random acoustic response of flat beam and plate structures with isotropic material properties. Reduced-order analysis predictions are compared with those made using a numerical simulation in physical degrees-of-freedom to quantify the error associated with the selected modal bases. Bending and membrane responses are separately presented to help differentiate the bases.

Rizzi, Stephen A.↗

Nonlinear Reduced Order Random Response Analysis of Structures with Shallow Curvature

The goal of this investigation is to further develop nonlinear modal numerical simulation methods for application to geometrically nonlinear response of structures with shallow curvature under random loadings. For reduced order analysis, the modal basis selection must be capable of reflecting the coupling in both the linear and nonlinear stiffness. For the symmetric shallow arch under consideration, four categories of modal basis functions are defined. Those having symmetric transverse displacements (ST modes) can be designated as transverse dominated (ST-T) modes and in-plane dominated (ST-I) modes. Those having anti-symmetric transverse displacements (AT modes) can similarly be designated as transverse dominated (AT-T) modes and in-plane dominated (AT-I) modes. The response of an aluminum arch under a uniformly distributed transverse random loading is investigated. Results from nonlinear modal simulations made using various modal bases are compared with those obtained from a numerical simulation in physical degrees-of-freedom. While inclusion of ST-T modes is important for all response regimes, it is found that the ST-I modes become increasingly important in the nonlinear response regime, and that AT-T and AT-I modes are critical in the autoparametric regime.

Przekop, Adam↗

On the Response of a Nonlinear Structure to High Kurtosis Non-Gaussian Random Loadings

This paper is a follow-on to recent work by the authors in which the response and high-cycle fatigue of a nonlinear structure subject to non-Gaussian loadings was found to vary markedly depending on the nature of the loading. There it was found that a non-Gaussian loading having a steady rate of short-duration, high-excursion peaks produced essentially the same response as would have been incurred by a Gaussian loading. In contrast, a non-Gaussian loading having the same kurtosis, but with bursts of high-excursion peaks was found to elicit a much greater response. This work is meant to answer the question of when consideration of a loading probability distribution other than Gaussian is important. The approach entailed nonlinear numerical simulation of a beam structure under Gaussian and non-Gaussian random excitations. Whether the structure responded in a Gaussian or non-Gaussian manner was determined by adherence to, or violations of, the Central Limit Theorem. Over a practical range of damping, it was found that the linear response to a non-Gaussian loading was Gaussian when the period of the system impulse response is much greater than the rate of peaks in the loading. Lower damping reduced the kurtosis, but only when the linear response was non-Gaussian. In the nonlinear regime, the response was found to be non-Gaussian for all loadings. The effect of a spring-hardening type of nonlinearity was found to limit extreme values and thereby lower the kurtosis relative to the linear response regime. In this case, lower damping gave rise to greater nonlinearity, resulting in lower kurtosis than a higher level of damping.

Rizzi, Stephen A.↗

Nonlinear Oscillators in Space Physics

We discuss dynamical systems that produce an oscillation without an external time dependent source. Numerical results are presented for nonlinear oscillators in the Em1h's atmosphere, foremost the quasi-biennial oscillation (QBOl. These fluid dynamical oscillators, like the solar dynamo, have in common that one of the variables in a governing equation is strongly nonlinear and that the nonlinearity, to first order, has particular form. of 3rd or odd power. It is shown that this form of nonlinearity can produce the fundamental li'equency of the internal oscillation. which has a period that is favored by the dynamical condition of the fluid. The fundamental frequency maintains the oscillation, with no energy input to the system at that particular frequency. Nonlinearities of 2nd or even power could not maintain the oscillation.

Lester,Daniel↗

The Nonlinear Magnetosphere: Expressions in MHD and in Kinetic Models

Like most plasma systems, the magnetosphere of the Earth is governed by nonlinear dynamic evolution equations. The impact of nonlinearities ranges from large scales, where overall dynamics features are exhibiting nonlinear behavior, to small scale, kinetic, processes, where nonlinear behavior governs, among others, energy conversion and dissipation. In this talk we present a select set of examples of such behavior, with a specific emphasis on how nonlinear effects manifest themselves in MHD and in kinetic models of magnetospheric plasma dynamics.

Hesse, Michael↗

Nonlinear Reduced-Order Simulation Using An Experimentally Guided Modal Basis

A procedure is developed for using nonlinear experimental response data to guide the modal basis selection in a nonlinear reduced-order simulation. The procedure entails using nonlinear acceleration response data to first identify proper orthogonal modes. Special consideration is given to cases in which some of the desired response data is unavailable. Bases consisting of linear normal modes are then selected to best represent the experimentally determined transverse proper orthogonal modes and either experimentally determined inplane proper orthogonal modes or the special case of numerically computed in-plane companions. The bases are subsequently used in nonlinear modal reduction and dynamic response simulations. The experimental data used in this work is simulated to allow some practical considerations, such as the availability of in-plane response data and non-idealized test conditions, to be explored. Comparisons of the nonlinear reduced-order simulations are made with the surrogate experimental data to demonstrate the effectiveness of the approach.

Rizzi, Stephen A.↗

Nonlinear Modeling by Assembling Piecewise Linear Models

To preserve nonlinearity of a full order system over a parameters range of interest, we propose a simple modeling approach by assembling a set of piecewise local solutions, including the first-order Taylor series terms expanded about some sampling states. The work by Rewienski and White inspired our use of piecewise linear local solutions. The assembly of these local approximations is accomplished by assigning nonlinear weights, through radial basis functions in this study. The efficacy of the proposed procedure is validated for a two-dimensional airfoil moving at different Mach numbers and pitching motions, under which the flow exhibits prominent nonlinear behaviors. All results confirm that our nonlinear model is accurate and stable for predicting not only aerodynamic forces but also detailed flowfields. Moreover, the model is robustness-accurate for inputs considerably different from the base trajectory in form and magnitude. This modeling preserves nonlinearity of the problems considered in a rather simple and accurate manner.

Yao, Weigang↗

A Nonlinear Modal Aeroelastic Solver for FUN3D

A nonlinear structural solver has been implemented internally within the NASA FUN3D computational fluid dynamics code, allowing for some new aeroelastic capabilities. Using a modal representation of the structure, a set of differential or differential-algebraic equations are derived for general thin structures with geometric nonlinearities. ODEPACK and LAPACK routines are linked with FUN3D, and the nonlinear equations are solved at each CFD time step. The existing predictor-corrector method is retained, whereby the structural solution is updated after mesh deformation. The nonlinear solver is validated using a test case for a flexible aeroshell at transonic, supersonic, and hypersonic flow conditions. Agreement with linear theory is seen for the static aeroelastic solutions at relatively low dynamic pressures, but structural nonlinearities limit deformation amplitudes at high dynamic pressures. No flutter was found at any of the tested trajectory points, though LCO may be possible in the transonic regime.

Goldman, Benjamin D.↗

Curved Displacement Transfer Functions for Geometric Nonlinear Large Deformation Structure Shape Predictions

For shape predictions of structures under large geometrically nonlinear deformations, Curved Displacement Transfer Functions were formulated based on a curved displacement, traced by a material point from the undeformed position to deformed position. The embedded beam (depth-wise cross section of a structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution could be described with a piecewise linear or a piecewise nonlinear function. The discretization approach enabled piecewise integrations of the embedded-beam curvature equations to yield the Curved Displacement Transfer Functions, expressed in terms of embedded beam geometrical parameters and surface strains. By entering the surface strain data into the Displacement Transfer Functions, deflections along each embedded beam can be calculated at multiple points for mapping the overall structural deformed shapes. Finite-element linear and nonlinear analyses of a tapered cantilever tubular beam were performed to generate linear and nonlinear surface strains and the associated deflections to be used for validation. The shape prediction accuracies were then determined by comparing the theoretical deflections with the finiteelement- generated deflections. The results show that the newly developed Curved Displacement Transfer Functions are very accurate for shape predictions of structures under large geometrically nonlinear deformations.

displacement theory↗

A Nonlinear Viscoelastic-Viscoplastic Constitutive Model for Epoxy Polymers

The objective of this paper is to develop a nonlinear viscoelastic–viscoplastic constitutive model for epoxy polymers. The classic nonlinear viscoelasticity model is reformulated to yield a closed-form incremental constitutive relation, which relates the stress increments to the viscoelastic strain increments. A viscoplasticity model, which consists of the Drucker–Prager yield function, nonlinear isotropic and kinematic hardening laws, and the Perzyna viscosity function, is subsequently developed. The nonlinear viscoelasticity and the vis-coplasticity models are then implemented in a radial return algorithm. A closed-form algorithmic tangent operator is derived to facilitate model implementation. The nonlinear viscoelasticity model is validated by reproducing the test data for polymethyl methacrylate(PMMA). The present constitutive model’s capabilities are demonstrated through modeling viscoelastic–viscoplastic PMMA loaded at different strain rates. The mechanics of structure genome-based micromechanics approach, embedded with the present constitutive model, is used to homogenize a unidirectional fiber-reinforced composite with a PMMA matrix, subjected to uniaxial and shear loading at different strain rates. The viscoelastic characterization of PMT-F7 epoxy is evaluated by comparing the predicted creep–recovery responses by the present constitutive model with the experimental ones.

Liang Zhang↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Data-Centric Approach to Capture Non-Polynomial Nonlinear Dynamics

We propose an analytical construction of observable functions in the extended dynamic mode decomposition (EDMD) algorithm. EDMD is a numerical method for approximating the spectral properties of the Koopman operator. The choice of observable functions is fundamental for applying EDMD to nonlinear problems arising in systems and control. Existing methods either start from a set of dictionary functions and look for the subset that best fits the underlying nonlinear dynamics or rely on machine learning algorithms to “learn” observable functions. Conversely, in this paper, we start from the dynamical system model and lift it through the Lie derivatives, rendering it into a polynomial form. This proposed transformation into a polynomial form is exact and provides an adequate set of observable functions. The strength of the proposed approach is its applicability to a broader class of nonlinear dynamical systems, particularly those with nonpolynomial functions and compositions thereof. Moreover, it retains the physical interpretability of the underlying dynamical system and can be readily integrated into existing numerical libraries. We demonstrate the proposed approach with an application to electric power systems. The modeled system consists of a single generator connected to an infinite bus, where nonlinear terms include sine and cosine functions. The results demonstrate the effectiveness of the proposed procedure in off-attractor nonlinear dynamics for estimation and prediction; the observable functions obtained from the proposed construction outperform methods that use dictionary functions comprising monomials or radial basis functions.

extended dynamic mode decomposition↗

Computing resonant modes of accelerator cavities by solving nonlinear eigenvalue problems via rational approximation

Here, we present an efficient and reliable algorithm for solving a class of nonlinear eigenvalue problems arising from the modeling of particle accelerator cavities. The eigenvalue nonlinearity in these problems results from the use of waveguides to couple external power sources or to allow certain excited electromagnetic modes to exit the cavity. We use a rational approximation to reduce the nonlinear eigenvalue problem first to a rational eigenvalue problem. We then apply a special linearization procedure to turn the rational eigenvalue problem into a larger linear eigenvalue problem with the same eigenvalues, which can be solved by existing iterative methods. By using a compact scheme to represent both the linearized operator and the eigenvectors to be computed, we obtain a numerical method that only involves solving linear systems of equations of the same dimension as the original nonlinear eigenvalue problem. We refer to this method as a compact rational Krylov (CORK) method. We implemented the CORK method in the Omega3P module of the Advanced Computational Electromagnetic 3D Parallel (ACE3P) simulation suite and validated it by comparing the computed cavity resonant frequencies and damping Q factors of a small model problem to those obtained from a fitting procedure that uses frequency responses computed by another ACE3P module called S3P. We also used the CORK method to compute trapped modes damped in an ideal eight 9-cell SRF cavity cryomodule. This was the first time it was possible to compute these modes directly. The damping Q factors of the computed modes match well with those measured in experiments and the difference in resonant frequencies is within the range introduced by cavity imperfection. Therefore, the CORK method is an extremely valuable tool for computational cavity design.

43 PARTICLE ACCELERATORS↗

New nonlinear optical potassium iodate K[IO{sub 3}] and borates K{sub 3}[B{sub 6}O{sub 10}]Br, KTa[B{sub 4}O{sub 6}(OH){sub 4}](OH){sub 2}{center_dot}1.33H{sub 2}O-Synthesis, structures and relation to the properties

Three optically uniaxial non-centrosymmetric potassium compounds, iodate K[IO{sub 3}] and two borates K{sub 3}[B{sub 6}O{sub 10}]Br, KTa[B{sub 4}O{sub 6}(OH){sub 4}](OH){sub 2}{center_dot}1.33H{sub 2}O have been synthesized and characterized by single-crystal X-ray diffraction. The materials were synthesized through hydrothermal techniques using initial reagents and mineralizers. All the compounds are trigonal-hexagonal: the space groups are R3 (K[IO{sub 3}]), R3m (K[B{sub 6}O{sub 10}]Br) and P-62m (KTa[B{sub 4}O{sub 6}(OH){sub 4}](OH){sub 2}{center_dot}1.33H{sub 2}O). Powder second-harmonic generation (SHG) measurements on crystals, using 1064 nm Nd:YAG laser radiation, indicate the materials are all phase-matchable and have strong second-order nonlinearities. The correlation with the perovskite structure has been found and described for K,Br-borate and K-iodate. Structure-properties relation is discussed and attributed to stereo-active lone-pair on I{sup 5+} and asymmetrical bonds in the compounds. The role of K atoms is pronounced from crystal chemistry point of view, contributing optical nonlinearities. - Graphical abstract: Structure-properties relation is attributed to stereo-active lone-pair on I{sup 5+} and asymmetrical bonds in the compounds with the crucial role of K atoms, K[IO{sub 3}] structure along 3-fold axis with umbrella-like IO{sub 3}-group and with octahedral coordination-analog with perovskite. Highlights: Black-Right-Pointing-Pointer Syntheses in hydrothermal conditions trigonal KIO{sub 3}, K{sub 3}[B{sub 6}O{sub 10}]Br and KTa[B{sub 4}O{sub 6}(OH){sub 4}](OH){sub 2}. Black-Right-Pointing-Pointer Correlation with the perovskite, found for K-iodate and K,Br-borate. Black-Right-Pointing-Pointer Structure-properties relation, attributed to I{sup 5+} lone-pair and asymmetric bonds including K-atoms. Black-Right-Pointing-Pointer The compounds are phase matchable and strongly nonlinear from powder SHG measurements. Black-Right-Pointing-Pointer Averaged second-order nonlinear coefficients Left-Pointing-Angle-Bracket d Right-Pointing-Angle-Bracket are 24, 2.8 and 1.9 pm/V.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear normal modes modal interactions and isolated resonance curves

The objective of the present study is to explore the connection between the nonlinear normal modes of an undamped and unforced nonlinear system and the isolated resonance curves that may appear in the damped response of the forced system. To this end, an energy balance technique is used to predict the amplitude of the harmonic forcing that is necessary to excite a specific nonlinear normal mode. A cantilever beam with a nonlinear spring at its tip serves to illustrate the developments. Furthermore, the practical implications of isolated resonance curves are also discussed by computing the beam response to sine sweep excitations of increasing amplitudes.

42 ENGINEERING↗

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗