Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,621 records · Page 90

Nonlinear axisymmetric liquid currents in spherical annuli

A numerical analysis of non-linear axisymmetric viscous flows in spherical annuli of different gap sizes is presented. Only inner sphere was supposed to rotate at a constant angular velocity. The streamlines, lines of constant angular velocity, kinetic energy spectra, and spectra of velocity components are obtained. A total kinetic energy and torque needed to rotate the inner sphere are calculated as functions of Re for different gap sizes. In small-gap annulus nonuniqueness of steady solutions of Navier-Stokes equations is established and regions of different flow regime existences are found. Numerical solutions in a wide-gap annulus and experimental results are used in conclusions about flow stability in the considered range of Re. The comparison of experimental and numerical results shows close qualitative and quantitative agreement.

Astafyeva, N. M.↗

Nonlinear Aerodynamics and the Design of Wing Tips

The analysis and design of wing tips for fixed wing and rotary wing aircraft still remains part art, part science. Although the design of airfoil sections and basic planform geometry is well developed, the tip regions require more detailed consideration. This is important because of the strong impact of wing tip flow on wing drag; although the tip region constitutes a small portion of the wing, its effect on the drag can be significant. The induced drag of a wing is, for a given lift and speed, inversely proportional to the square of the wing span. Concepts are proposed as a means of reducing drag. Modern computational methods provide a tool for studying these issues in greater detail. The purpose of the current research program is to improve the understanding of the fundamental issues involved in the design of wing tips and to develop the range of computational and experimental tools needed for further study of these ideas.

Kroo, Ilan↗

Development of Renormalization Group Analysis of Turbulence

The renormalization group (RG) procedure for nonlinear, dissipative systems is now quite standard, and its applications to the problem of hydrodynamic turbulence are becoming well known. In summary, the RG method isolates self similar behavior and provides a systematic procedure to describe scale invariant dynamics in terms of large scale variables only. The parameterization of the small scales in a self consistent manner has important implications for sub-grid modeling. This paper develops the homogeneous, isotropic turbulence and addresses the meaning and consequence of epsilon-expansion. The theory is then extended to include a weak mean flow and application of the RG method to a sequence of models is shown to converge to the Navier-Stokes equations.

L M Smith↗

Thermoviscoplastic analysis of fibrous periodic composites using triangular subvolumes

The nonlinear viscoplastic behavior of fibrous periodic composites is analyzed by discretizing the unit cell into triangular subvolumes. A set of these subvolumes can be configured by the analyst to construct a representation for the unit cell of a periodic composite. In each step of the loading history, the total strain increment at any point is governed by an integral equation which applies to the entire composite. A Fourier series approximation allows the incremental stresses and strains to be determined within a unit cell of the periodic lattice. The nonlinearity arising from the viscoplastic behavior of the constituent materials comprising the composite is treated as fictitious body force in the governing integral equation. Specific numerical examples showing the stress distributions in the unit cell of a fibrous tungsten/copper metal matrix composite under viscoplastic loading conditions are given. The stress distribution resulting in the unit cell when the composite material is subjected to an overall transverse stress loading history perpendicular to the fibers is found to be highly heterogeneous, and typical homogenization techniques based on treating the stress and strain distributions within the constituent phases as homogeneous result in large errors under inelastic loading conditions.

Walker, Kevin P.↗

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗

Influence of thickness and camber on the aeroelastic stability of supersonic throughflow fans: An engineering approach

An engineering approach was used to include the nonlinear effects of thickness and camber in an analytical aeroelastic analysis of cascades in supersonic acial flow (supersonic leading-edge locus). A hybrid code using Lighthill's nonlinear piston theory and Lanes's linear potential theory was developed to include these nonlinear effects. Lighthill's theory was used to calculate the unsteady pressures on the noninterference surface regions of the airfoils in cascade. Lane's theory was used to calculate the unsteady pressures on the remaining interference surface regions. Two airfoil profiles was investigated (a supersonic throughflow fan design and a NACA 66-206 airfoil with a sharp leading edge). Results show that compared with predictions of Lane's potential theory for flat plates, the inclusion of thickness (with or without camber) may increase or decrease the aeroelastic stability, depending on the airfoil geometry and operating conditions. When thickness effects are included in the aeroelastic analysis, inclusion of camber will influence the predicted stability in proportion to the magnitude of the added camber. The critical interblade phase angle, depending on the airfoil profile and operating conditions, may also be influenced by thickness and camber. Compared with predictions of Lane's linear potential theory, the inclusion of thickness and camber decreased the aerodynamic stifness and increased the aerodynamic damping at Mach 2 and 2.95 for a cascade of supersonic throughflow fan airfoils oscillating 180 degrees out of phase at a reduced frequency of 0.1.

Ramsey, John K.↗

Influence of thickness and camber on the aeroelastic stability of supersonic throughflow fans

An engineering approach was used to include the nonlinear effects of thickness and camber in an analytical aeroelastic analysis of casecades in supersonic acial flow (supersonic leading-edge locus). A hybrid code using Lighthill's nonlinear piston theory and Lane's linear potential theory was developed to include these nonlinear effects. Lighthill's theory was used to calculate the unsteady pressures on the noninterference surface regions of the airfoils in cascade. Lane's theory was used to calculate the unsteady pressures on the remaining interference surface regions. Two airfoil profiles were investigated (a supersonic throughflow fan design and a NACA 66-206 airfoil with a sharp leading edge). Results show that compared with predictions of Lane's potential theory for flat plates, the inclusion of thickness (with or without camber) may increase or decrease the aeroelastic stability, depending on the airfoil geometry and operating conditions. When thickness effects are included in the aeroelastic analysis, inclusion of camber will influence the predicted stability in proportion to the magnitude of the added camber. The critical interblade phase angle, depending on the airfoil profile and operating conditions, may also be influenced by thickness and camber. Compared with predictions of Lane's linear potential theory, the inclusion of thickness and camber decreased the aerodynamic stiffness and increased the aerodynamic damping at Mach 2 and 2.95 for a cascade of supersonic throughflow fan airfoils oscillating 180 degrees out of phase at a reduced frequency of 0.1.

Ramsey, John K.↗

Dynamical Approach Study of Spurious Numerics in Nonlinear Computations

The last two decades have been an era when computation is ahead of analysis and when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems and non-traditional fields. Ensuring a higher level of confidence in the predictability and reliability (PAR) of these numerical simulations could play a major role in furthering the design, understanding, affordability and safety of our next generation air and space transportation systems, and systems for planetary and atmospheric sciences, and in understanding the evolution and origin of life. The need to guarantee PAR becomes acute when computations offer the ONLY way of solving these types of data limited problems. Employing theory from nonlinear dynamical systems, some building blocks to ensure a higher level of confidence in PAR of numerical simulations have been revealed by the author and world expert collaborators in relevant fields. Five building blocks with supporting numerical examples were discussed. The next step is to utilize knowledge gained by including nonlinear dynamics, bifurcation and chaos theories as an integral part of the numerical process. The third step is to design integrated criteria for reliable and accurate algorithms that cater to the different multiscale nonlinear physics. This includes but is not limited to the construction of appropriate adaptive spatial and temporal discretizations that are suitable for the underlying governing equations. In addition, a multiresolution wavelets approach for adaptive numerical dissipation/filter controls for high speed turbulence, acoustics and combustion simulations will be sought. These steps are corner stones for guarding against spurious numerical solutions that are solutions of the discretized counterparts but are not solutions of the underlying governing equations.

Yee, H. C.↗

Nonlinear resonances in a class of multi-degree-of-freedom systems

An analysis is presented of the superharmonic, subharmonic, and combination resonances in a multi-degree-of-freedom system which has cubic nonlinearity and modal viscous damping and is subject to harmonic excitation. It is shown that, in the absence of internal resonances, the steady-state response contains only the modes which are directly excited. It is shown that, in the presence of internal resonances, modes other than those that are directly excited can appear in the response. The strong influence of internal resonances is exhibited in numerical examples involving hinged-clamped beams. It is shown that when a multimode solution exists the lowest mode can dominate the response, even when it is not directly excited.

Sridhar, S.↗

Non-isothermal elastoviscoplastic analysis of planar curved beams

The development of a general mathematical model and solution methodologies, to examine the behavior of thin structural elements such as beams, rings, and arches, subjected to large nonisothermal elastoviscoplastic deformations is presented. Thus, geometric as well as material type nonlinearities of higher order are present in the analysis. For this purpose a complete true abinito rate theory of kinematics and kinetics for thin bodies, without any restriction on the magnitude of the transformation is presented. A previously formulated elasto-thermo-viscoplastic material constitutive law is employed in the analysis. The methodology is demonstrated through three different straight and curved beams problems.

Simitses, G. J.↗

A Bayesian approach to evaluation of soil biogeochemical models

Abstract. To make predictions about the carbon cycling consequences of rising global surface temperatures, Earth system scientists rely on mathematical soil biogeochemical models (SBMs). However, it is not clear which models have better predictive accuracy, and a rigorous quantitative approach for comparing and validating the predictions has yet to be established. In this study, we present a Bayesian approach to SBM comparison that can be incorporated into a statistical model selection framework. We compared the fits of linear and nonlinear SBMs to soil respiration data compiled in a recent meta-analysis of soil warming field experiments. Fit quality was quantified using Bayesian goodness-of-fit metrics, including the widely applicable information criterion (WAIC) and leave-one-out cross validation (LOO). We found that the linear model generally outperformed the nonlinear model at fitting the meta-analysis data set. Both WAIC and LOO computed higher overfitting risk and effective numbers of parameters for the nonlinear model compared to the linear model, conditional on the data set. Goodness of fit for both models generally improved when they were initialized with lower and more realistic steady-state soil organic carbon densities. Still, testing whether linear models offer definitively superior predictive performance over nonlinear models on a global scale will require comparisons with additional site-specific data sets of suitable size and dimensionality. Such comparisons can build upon the approach defined in this study to make more rigorous statistical determinations about model accuracy while leveraging emerging data sets, such as those from long-term ecological research experiments.

54 ENVIRONMENTAL SCIENCES↗

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.↗

Beyond optimization—supervised learning applications in relativistic laser-plasma experiments

We explore the applications of machine learning techniques in relativistic laser-plasma experiments beyond optimization purposes. We predict the beam charge of electrons produced in a laser wakefield accelerator given the laser wavefront change caused by a deformable mirror. Machine learning enables feature analysis beyond merely searching for an optimal beam charge, showing that specific aberrations in the laser wavefront are favored in generating higher beam charges. Supervised learning models allow characterizing the measured data quality as well as recognizing irreproducible data and potential outliers. Furthermore, we also include virtual measurement errors in the experimental data to examine the model robustness under these conditions. This work demonstrates how machine learning methods can benefit data analysis and physics interpretation in a highly nonlinear problem of relativistic laser-plasma interaction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Motion Cueing Algorithm Development: Piloted Performance Testing of the Cueing Algorithms

The relative effectiveness in simulating aircraft maneuvers with both current and newly developed motion cueing algorithms was assessed with an eleven-subject piloted performance evaluation conducted on the NASA Langley Visual Motion Simulator (VMS). In addition to the current NASA adaptive algorithm, two new cueing algorithms were evaluated: the optimal algorithm and the nonlinear algorithm. The test maneuvers included a straight-in approach with a rotating wind vector, an offset approach with severe turbulence and an on/off lateral gust that occurs as the aircraft approaches the runway threshold, and a takeoff both with and without engine failure after liftoff. The maneuvers were executed with each cueing algorithm with added visual display delay conditions ranging from zero to 200 msec. Two methods, the quasi-objective NASA Task Load Index (TLX), and power spectral density analysis of pilot control, were used to assess pilot workload. Piloted performance parameters for the approach maneuvers, the vertical velocity upon touchdown and the runway touchdown position, were also analyzed but did not show any noticeable difference among the cueing algorithms. TLX analysis reveals, in most cases, less workload and variation among pilots with the nonlinear algorithm. Control input analysis shows pilot-induced oscillations on a straight-in approach were less prevalent compared to the optimal algorithm. The augmented turbulence cues increased workload on an offset approach that the pilots deemed more realistic compared to the NASA adaptive algorithm. The takeoff with engine failure showed the least roll activity for the nonlinear algorithm, with the least rudder pedal activity for the optimal algorithm.

Houck, Jacob A.↗

Multiscale Static Analysis of Notched and Unnotched Laminates Using the Generalized Method of Cells

The generalized method of cells (GMC) is demonstrated to be a viable micromechanics tool for predicting the deformation and failure response of laminated composites, with and without notches, subjected to tensile and compressive static loading. Given the axial [0], transverse [90], and shear [+45/-45] response of a carbon/epoxy (IM7/977-3) system, the unnotched and notched behavior of three multidirectional layups (Layup 1: [0,45,90,-45](sub 2S), Layup 2: [0,60,0](sub 3S), and Layup 3: [30,60,90,-30, -60](sub 2S)) are predicted under both tensile and compressive static loading. Matrix nonlinearity is modeled in two ways. The first assumes all nonlinearity is due to anisotropic progressive damage of the matrix only, which is modeled, using the multiaxial mixed-mode continuum damage model (MMCDM) within GMC. The second utilizes matrix plasticity coupled with brittle final failure based on the maximum principle strain criteria to account for matrix nonlinearity and failure within the Finite Element Analysis--Micromechanics Analysis Code (FEAMAC) software multiscale framework. Both MMCDM and plasticity models incorporate brittle strain- and stress-based failure criteria for the fiber. Upon satisfaction of these criteria, the fiber properties are immediately reduced to a nominal value. The constitutive response for each constituent (fiber and matrix) is characterized using a combination of vendor data and the axial, transverse, and shear responses of unnotched laminates. Then, the capability of the multiscale methodology is assessed by performing blind predictions of the mentioned notched and unnotched composite laminates response under tensile and compressive loading. Tabulated data along with the detailed results (i.e., stress-strain curves as well as damage evolution states at various ratios of strain to failure) for all laminates are presented.

GMC↗

The Elasto-Plastic Stability of Plates

This article explains results developed from the following research: 'The Stability of Plates and Shells beyond the Elastic Limit.' A significant improvement is found in the derivation of the relations between the stress factors and the strains resulting from the instability of plates and shells. In a strict analysis, the problem reduces to the solution of two simultaneous nonlinear partial differential equations of the fourth order in the deflection and stress function, and in the approximate analysis to a single linear equation of the Bryan type. Solutions are given for the special cases of a rectangular plate buckling into a cylindrical form, and of an arbitrarily shaped plate under uniform compression. These solutions indicate that the accuracy obtained by the approximate method is satisfactory.

Ilyushin, A. A.↗

The dynamic aspects of thermo-elasto-viscoplastic snap-through and creep buckling phenomena

Use of a mathematical model and solution methodology, to examine dynamic buckling and dynamic postbuckling behavior of shallow arches and spherical caps made of a realistic material and undergoing non-isothermal, elasto-viscoplastic deformation was examined. Thus, geometric as well as material type nonlinearities of higher order are included in this analysis. The dynamic stability problem is studied under impulsive loading and suddenly applied loading with loads of constant magnitude and infinite duration. A finite element model was derived directly from the incrementally formulated nonlinear shell equations, by using a tensor-oriented procedure. As an example of the results, the time history of the midspan displacement of a damped shallow circular arch is presented.

Riff, R.↗