Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

On some recent advances in multidisciplinary analysis of hypersonic vehicles

This paper presents pertinent details of the development and application of an integrated, multidisciplinary finite element analysis tool for the effective modeling and simulation of aerospace vehicles including hypersonic ones. Recent advances in this connection include more efficient CFD solution techniques and also aeroservoelastic stability analysis methods that employ unsteady aerodynamic forces for computation of flutter and divergence speeds. An accelerated Euler solution procedure based on the Aitken acceleration technique has recently been implemented that effects considerable improvement in solution efficiency for the same level of solution accuracy, as evidenced by example problems presented in the paper. A summary of finite element numerical formulations for the various individual disciplines as well as for the unified aeroservoelastic analysis and some relevant numerical examples are also presented in the paper.

Gupta, K. K.↗
NASA NTRS2012

An Experimental and Numerical Study of Roughness-Induced Instabilities in a Mach 3.5 Boundary Layer

Progress on a joint experimental and numerical study of laminar-to-turbulent transition induced by an isolated roughness element in a high-speed laminar boundary layer is reported in this paper. The numerical analysis suggests that transition is driven by the instability of high- and low-speed streaks embedded in the wake of the isolated roughness element. In addition, spatial stability analysis revealed that the wake flow supports multiple modes (even and odd) of convective instabilities that experience strong enough growth to cause transition. The experimental measurements, which included hot-wire and pitot-probe surveys, confirmed the existence of embedded high- and low-speed streaks in the roughness wake. Furthermore, the measurements indicate the presence of both even and odd modes of instability, although their relative magnitude depends on the specifics of the roughness geometry and flow conditions (e.g., the value of Re(sub kk) or k/delta. For the two test cases considered in the measurements (Re(sub kk) values of 462 and 319), the even mode and the odd mode were respectively dominant and appear to play a primary role in the transition process.

Kegerise, Michael A.↗
NASA NTRS1986

The effect of circumferential aerodynamic detuning on coupled bending-torsion unstalled supersonic flutter

A mathematical model developed to predict the enhanced coupled bending-torsion unstalled supersonic flutter stability due to alternate circumferential spacing aerodynamic detuning of a turbomachine rotor. The translational and torsional unsteady aerodynamic coefficients are developed in terms of influence coefficients, with the coupled bending-torsion stability analysis developed by considering the coupled equations of this aerodynamic detuning on coupled bending-torsion unstalled supersonic flutter as well as the verification of the modeling are then demonstrated by considering an unstable 12 bladed rotor, with Verdon's uniformly spaced Cascade B flow geometry as a baseline. However, with the elastic axis and center of gravity at 60 percent of the chord, this type of aerodynamic detuning has a minimal effect on stability. For both uniform and nonuniform circumferentially space rotors, a single degree of freedom torsion mode analysis was shown to be appropriate for values of the bending-torsion natural frequency ratio lower than 0.6 and higher 1.2. When the elastic axis and center of gravity are not coincident, the effect of detuning on cascade stability was found to be very sensitive to the location of the center of gravity with respect to the elastic axis. In addition, it was determined that when the center of gravity was forward of an elastic axis located at midchord, a single degree of freedom torsion model did not accurately predict cascade stability.

Hoyniak, D.↗
NASA NTRS1973

BUCLAP2: A computer program for instability analysis of laminated long plates subjected to combined inplane loads

This program description document describes the program structure and the design details of a CDC 6600 FORTRAN 4 digital computer program, BUCLAP2, which uses minimum energy principles to do an elastic stability analysis of curved and flat laminated rectangular long plates subjected to combined inplane normal and shear loads. Given the geometry, the material properties, and slected boundary conditions for the plate element, the program calculates the minimum buckling load for various wave lengths. The two parallel ends of the program calculates the minimum buckling load for various wave lengths. The two parallel ends of the long plate must be simply supported and arbitrary elastic boundary conditions may be imposed along either one or both external longitudinal sides. For guide to program use, see

Halstead, D. W.↗
NASA NTRS1974

Conical isogrid adapter structural test results. Appendix C: Cylindrical isogrid structures

Analysis techniques and data for evaluating the load carrying capabilities of cylindrical isogrid structures subjected to compression and body bending loads are presented. Structural tests were conducted on a 120-inch diameter by 37-inch long cylindrical adapter. Photographs and drawings of the flanged isogrid adapter structure are included. The subjects considered are: (1) summary of the test data resulting from structural analysis, (2) techniques for predicting local and general instability, (3) a method for predicting extended local instability, (4) the effects of basic parameter and constant variations in extended local stability analysis, and (5) evaluation of effective skin width and compression buckling of skin panels on the basis of test data.

Source record↗
NASA NTRS1986

The predicted effect of aerodynamic detuning on coupled bending-torsion unstalled supersonic flutter

A mathematical model is developed to predict the enhanced coupled bending-torsion unstalled supersonic flutter stability due to alternate circumferential spacing aerodynamic detuning of a turbomachine rotor. The translational and torsional unsteady aerodynamic coefficients are developed in terms of influence coefficients, with the coupled bending-torsion stability analysis developed by considering the coupled equations of motion together with the unsteady aerodynamic loading. The effect of this aerodynamic detuning on coupled bending-torsion unstalled supersonic flutter as well as the verification of the modeling are then demonstrated by considering an unstable 12 bladed rotor, with Verdon's uniformly spaced Cascade B flow geometry as a baseline. However, with the elastic axis and center of gravity at 60 percent of the chord, this type of aerodynamic detuning has a minimal effect on stability. For both uniform and nonuniform circumferentially space rotors, a single degree of freedom torsion mode analysis was shown to be appropriate for values of the bending-torsion natural frequency ratio lower than 0.6 and higher 1.2. When the elastic axis and center of gravity are not coincident, the effect of detuning on cascade stability was found to be very sensitive to the location of the center of gravity with respect to the elastic axis. In addition, it was determined that when the center of gravity was forward of an elastic axis located at midchord, a single degree of freedom torsion model did not accurately predict cascade stability.

Hoyniak, D.↗
NASA NTRS1992

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗
DOE OSTI2024

Conceptual Design of a Tension Leg Platform With 22.3 MW Vertical Axis Turbine

Here, this paper presents the conceptual design of a tension leg platform (TLP) for the ARCUS “towerless” vertical-axis wind turbine (VAWT). VAWTs are ideal for floating offshore sites and have several advantages over horizontal-axis wind turbines (HAWT) including reduced top mass, lower center of gravity, increased energy capture, and in turn lower cost. The towerless ARCUS VAWT drives these advantages further through increased structural efficiency and by enabling more optimized TLP designs with simplified installation procedures. For hull sizing, we have studied three turbine sizes with corresponding power ratings of 5.1 MW, 10.4 MW and 22.3 MW. The largest turbine was identified as having the greatest potential to reduce the levelized cost of energy (LCOE) and is the reference size used for the further detailed design process. The conceptual design of the VAWT TLP has been awarded with an ABS Approval in Principle Certificate. This paper contains brief analysis results and design findings for a TLP designed to house a VAWT, including the following topics: • Applicable Design Codes • Metocean Conditions • ARCUS Turbine Loads • Design Load Cases and Requirements - Pre-service TLP Stability - In-place TLP Global Performance • Platform Configurations, Hull Structure Scantling Design, Weight and CG Estimation, and General Arrangement Drawings • Hull Ballast Plan for both Pre-service and In-place Conditions • Pre-service Quayside Integration, Transportation and Wet Tow Stability Analysis • Global Performance Analysis for Motions and Tendon tensions • Summary of cost components and system levelized cost of energy

17 WIND ENERGY↗
NASA NTRS1972

Hydrodynamic stability of viscoelastic fluids - Importance of fluid model, overstability, and form of disturbance.

A linearized stability analysis has been applied to a fluid flowing in a gravity field between horizontal planes in Couette flow under conditions such that the temperature of the bottom plane exceeds that of the top. It is shown that, under conditions likely to be encountered with polymer solutions, oscillatory instabilities will not be controlling. Criteria are offered for ascertaining when an analysis based upon a second-order fluid model may be expected to yield physically meaningful results. It is also shown that, for the fluid model considered, critical conditions for stability are not changed when disturbances which vary in the flow direction are substituted for those which are a function of the coordinate transverse to the flow.

Mcintire, L. V.↗
NASA NTRS2017

Non-Linear Analysis of the NASA Super Pressure Balloons: Whole Flight Simulations

Tensys have a long-established background in the shape generation and load analysis of architectural stressed membrane structures. Founded upon their inTENS finite element analysis suite, these activities have broadened to encompass 'lighter than air' structures such as aerostats, airships, hybrid air vehicles and stratospheric balloons. Tensys have acted as consultants to the NASA Super Pressure Balloon (SPB) Program since 2004. Previous papers have focussed upon the application of inTENS to the overall structural and stability analysis of pumpkin type balloons as used in the SPB Program. Particular emphasis has been placed upon the ability to study both stress and stability at all stages of a flight. As the program has developed, increasing modelling fidelity has been introduced as design refinement has moved emphasis from the overall shape to the performance of individual details. Examples include the introduction of a contact capability into inTENS to better represent the separate cap layer of film. Localised investigations have considered the consequence of local geometric anomalies introduced during fabrication and the effects of debonding of the PBO tendon within its sleeve. Analysis to date has used a material model for the polyethylene shell film developed by Dr Rand of Winzen Engineering, supported by a program of fine resolution material tests at the Balloon Research and Development Laboratory (BRDL) at the GSFC Wallops Island Facility. Based upon work first reported by Schapery, this model has been incorporated into so-called 'snapshot' analyses within inTENS. For a given elapsed time point and temperature distribution, non-linear material properties are iteratively updated for the current stress state in each individual element until changes in those material properties are insignificant. This process has to be iterative as the film stresses and material properties are interdependent. Attempts to derive a time-stepping incremental viscoelastic capability using the Schapery Rand model encountered problems when dealing with the varying temperatures associated with balloon deployment, pressurisation and diurnal behaviour. An alternative approach has been developed by Pellegrino et al at Caltech, again with support from the NASA Balloon Program Office. This is a large strain non-linear viscoelastic model that includes film out-of-plane mechanical and thermal effects. Working from a unified base, this model can be utilised for time-stepping analyses in both a modulus or compliance mode, or a combination of both. This new capability is intended to enhance the current NASA balloon design process. This paper presents the implementation of the Caltech model into a specialist finite element suite, inTENS and its application to whole flight simulations.

Wakefield, David↗
DOE OSTI2024

Velocity Dependence of Rate‐And‐State Friction in Granular Fault Gouge and Implications for Slow‐Slip Events

Abstract The standard rate‐and‐state friction (RSF) has extensively captured frictional behaviors, but it fails to explain the velocity dependence of frictional stability transition and widespread slow‐slip events (SSEs) in experiments and nature adequately. An alternative microphysical Chen‐Niemeijer‐Spiers (CNS) model can well describe the velocity dependence of frictional behaviors of granular gouges. Using the original CNS model, standard RSF parameters can be quantified microphysically. However, some micro‐parameters are not easy to estimate quantitatively, making it difficult to extrapolate to natural and experimental conditions. Here, we simplify the microphysically‐derived RSF parameters including direct effect a , evolution effect b , and critical slip distance D c , as well as equivalent values ( a eq , b eq , and D eq ). The simplified friction parameters directly illustrate their velocity dependence, namely the essentially constant a , a eq , and D c , negatively velocity‐dependent b and b eq , as well as varying D eq for different laws. They are roughly consistent with experimental results in various fault gouges. A modified CNS model is further derived from the original CNS model, establishing a direct link between the standard RSF and CNS models. The modified CNS model exhibits virtually identical frictional behaviors to the original CNS, but differs from the standard RSF at large velocity perturbations. Moreover, the linearized stability analysis indicates that the critical stiffness for the modified CNS model is velocity‐dependent. Compared with the standard RSF, the modified CNS model not only explains the velocity dependence of frictional stability transition, but also exhibits a more gradual transition for SSEs with a broader range of stiffness ratios.

Geochemistry & Geophysics↗
NASA NTRS1996

Oscillatory Thermocapillary Convection

Stability analysis of thermocapillary convection in rectangular cavities is performed using direct numerical simulation. The influence of the Reynolds number (Re), the fluid Prandtl number (Pr) and cavity aspect ratio (Ar) on the motion is investigated. Neutral stability curves for transition to time-dependent convection are delineated in the Re-Ar plane for fluids with Pr=1.0, 4.4, 6.78, and 10. Several interesting features of these diagrams are discussed. One important conclusion is that Ar(epsilon tau) increases as Pr decreases. Thus large values of both Ar and Re are necessary to induce thermocapillary oscillations for small Pr fluids such as liquid metals and semiconductor melts. Energy analysis is also performed for the oscillatory flow in the neighborhood of critical points in order to gain insight into the mechanisms leading to instability.

Xu, Jieyong↗
DOE OSTI2023

Insights into ENSO Diversity from an Intermediate Coupled Model. Part I: Uniqueness and Sensitivity of the ENSO Mode

Abstract The basic dynamics of the spatiotemporal diversity for El Niño–Southern Oscillation (ENSO) has been the subject of extensive research and, while several hypotheses have been proposed, remains elusive. One promising line of studies suggests that the observed eastern Pacific (EP) and central Pacific (CP) ENSO may originate from two coexisting leading ENSO modes. We show that the coexistence of unstable EP-like and CP-like modes in these studies arises from contaminated linear stability analysis due to unnoticed numerical scheme caveats. In this two-part study, we further investigate the dynamics of ENSO diversity within a Cane–Zebiak-type model. We first revisit the linear stability issue to demonstrate that only one ENSO-like linear leading mode exists under realistic climate conditions. This single leading ENSO mode can be linked to either a coupled recharge-oscillator (RO) mode favored by the thermocline feedback or a wave-oscillator (WO) mode favored by the zonal advective feedback at the weak air–sea coupling end. Strong competition between the RO and WO modes for their prominence in shaping this ENSO mode into a generalized RO mode makes it sensitive to moderate changes in these two key feedbacks. Modulations of climate conditions yield corresponding modulations in spatial pattern, amplitude, and period associated with this ENSO mode. However, the ENSO behavior undergoing this linear climate condition modulations alone does not seem consistent with the observed ENSO diversity, suggesting the inadequacy of linear dynamics in explaining ENSO diversity. A nonlinear mechanism for ENSO diversity will be proposed and discussed in Part II.

Geng, Licheng↗
DOE OSTI2024

An Experimental Setup for Mechanical Vibration Analysis Using VLC

This study explores the potential applications across various domains, including earthquake detection and warning systems, where the system’s sensitivity to ground vibrations can contribute to early seismic event detection. Additionally, the study paves the way of developing applications of VLC/T in mechanical vibration and stability analysis of engines and platforms, offering insights into structural integrity and performance optimization. These multifaceted applications underscore the adaptability and potential of VLC/T systems in diverse fields, heralding advancements in sensing, communication, and security technologies. To achieve this, in this study, Peak to Average Power Ratio (PAPR) is proposed to represent the impact of mechanical shocks and vibrations generated by several weights dropped onto the platform with which the receiver is fixed. Even though non-contact measurement methodology is preferred for various reasons, the proposed measurement campaign obtains the data in contact form; however, the system and signal model proposed in this study could easily be extended into non-contact form. Considering the fact that the proposed measurement campaign employs off-the-shelf products, it is cost-effective and very scalable.

Yilmaz, Ahmet Mucahit↗
NASA NTRS1975

Design and application of squeeze film dampers for turbomachinery stabilization

The steady-state transient response of the squeeze film damper bearing was investigated. Both the steady-state and transient equations for the hydrodynamic bearing forces are derived; the steady-state equations were used to determine the damper equivalent stiffness and damping coefficients. These coefficients are used to find the damper configuration which will provide the optimum support characteristics based on a stability analysis of the rotor-bearing system. The effects of end seals and cavitated fluid film are included. The transient analysis of rotor-bearing systems was conducted by coupling the damping and rotor equations and integrating forward in time. The effects of unbalance, cavitation, and retainer springs are included. Methods of determining the stability of a rotor-bearing system under the influence of aerodynamic forces and internal shaft friction are discussed.

Gunter, E. J.↗
NASA NTRS1991

Nonlinear stability and control study of highly maneuverable high performance aircraft

The purpose was to develop and apply new nonlinear system methodologies to the stability analysis and adaptive control of high angle of attack (alpha) aircraft such as the F-18. Considerable progress is documented on nonlinear adaptive control and associated model development, identification, and simulation. The analysis considered linear and nonlinear, longitudinal, high alpha aircraft dynamics with varying degrees of approximation dependent on the purpose. In all cases, angle of attack or pitch rate was controlled primarily by a horizontal stabilizer. In most cases studied, a linear adaptive controller provided sufficient stability. However, it has been demonstrated by simulation of a simplified nonlinear model that certain large rapid maneuvers were not readily stabilized by the investigated linear adaptive control, but were controlled instead by means of a nonlinear time-series based adaptive control.

Mohler, R. R.↗