Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Condensation polyimides

Polyimides belong to a class of polymers known as polyheterocyclics. Unlike most other high temperature polymers, polyimides can be prepared from a variety of inexpensive monomers by several synthetic routes. The glass transition and crystalline melt temperature, thermooxidative stability, toughness, dielectric constant, coefficient of thermal expansion, chemical stability, mechanical performance, etc. of polyimides can be controlled within certain boundaries. This versatility has permitted the development of various forms of polyimides. These include adhesives, composite matrices, coatings, films, moldings, fibers, foams and membranes. Polyimides are synthesized through both condensation (step-polymerization) and addition (chain growth polymerization) routes. The precursor materials used in addition polyimides or imide oligomers are prepared by condensation method. High molecular weight polyimide made via polycondensation or step-growth polymerization is studied. The various synthetic routes to condensation polyimides, structure/property relationships of condensation polyimides and composite properties of condensation polyimides are all studied. The focus is on the synthesis and chemical structure/property relationships of polyimides with particular emphasis on materials for composite application.

Hergenrother, P. M.↗
NASA NTRS1991

Effect of Zn additions on precipitation during aging of alloy 8090

Precipitation events have been observed by TEM under two different aging conditions in three 'stretched' alloys, whose compositions are encompassed by the 8090 composition window but contain Zn additions of up to 1.07 wt pct. DSC was also used to obtain deeper insight of the precipitation-event effects obtainable through Zn content variation; it was thereby revealed that Zn is incorporated into the delta-prime phase, perhaps stabilizing it. Coarse, Zn-containing precipitates can form on the boundaries and within the interiors of the grains, when the Zn content reaches the presently investigated maximum of 1.07 wt pct.

Kilmer, R. J.↗
NASA NTRS1992

A quasi-static magnetospheric convection model in two-dimensions

A self-consistent 2D model of subsonic plasma sheet convection is presented. Specifically, time sequences of static equilibrium solutions for the 2D magnetospheric magnetic field are constructed consistent with adiabatic convection. This model self-consistently includes a dipole field and a reasonable accounting for the effects of inner magnetospheric shielding. Starting from a relaxed magnetospheric equilibrium, the earthward convection of plasma sheet flux tubes results in the stretching of inner plasma sheet field lines, the development of a local minimum in the equatorial magnetic field Be in the near-earth plasma sheet, and an increasing lobe magnetic field. This evolution in time occurs independent of the specific magnetopause or far-tail boundary conditions, provided the plasma sheet is not at marginal interchange stability. This behavior results solely from the convection of flux tubes of increasing entropy into the near-earth plasma sheet. These results are discussed in the general context of earthward convection in earth's plasma sheet and in the specific context of magnetospheric substorms.

Erickson, G. M.↗
NASA NTRS1993

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗
NASA NTRS2005

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗
NASA NTRS2010

A Simple Technique for Creating Regional Composites of Sea Surface Temperature from MODIS for Use in Operational Mesoscale NWP

This paper describes a simple technique for creating regional, high-resolution, daytime and nighttime composites of sea surface temperature (SST) for use in operational numerical weather prediction (NWP). The composites are based on observations from NASA s Moderate Resolution Imaging Spectroradiometer (MODIS) aboard Aqua and Terra. The data used typically are available nearly in real time, are applicable anywhere on the globe, and are capable of roughly representing the diurnal cycle in SST. The composites resolution is much higher than that of many other standard SST products used for operational NWP, including the low- and high-resolution Real-Time Global (RTG) analyses. The difference in resolution is key because several studies have shown that highly resolved SSTs are important for driving the air sea interactions that shape patterns of static stability, vertical and horizontal wind shear, and divergence in the planetary boundary layer. The MODIS-based composites are compared to in situ observations from buoys and other platforms operated by the National Data Buoy Center (NDBC) off the coasts of New England, the mid-Atlantic, and Florida. Mean differences, mean absolute differences, and root-mean-square differences between the composites and the NDBC observations are all within tenths of a degree of those calculated between RTG analyses and the NDBC observations. This is true whether or not one accounts for the mean offset between the skin temperatures of the MODIS dataset and the bulk temperatures of the NDBC observations and RTG analyses. Near the coast, the MODIS-based composites tend to agree more with NDBC observations than do the RTG analyses. The opposite is true away from the coast. All of these differences in point-wise comparisons among the SST datasets are small compared to the 61.08C accuracy of the NDBC SST sensors. Because skin-temperature variations from land to water so strongly affect the development and life cycle of the sea breeze, this phenomenon was chosen for demonstrating the use of the MODIS-based composite in an NWP model. A simulated sea breeze in the vicinity of New York City and Long Island shows a small, net, but far from universal improvement when MODIS-based composites are used in place of RTG analyses. The timing of the sea breeze s arrival is more accurate at some stations, and the near-surface temperature, wind, and humidity within the breeze are more realistic.

Knievel, Jason C.↗
NASA NTRS1977

Dynamical theory of stability for elastic rods with nonlinear curvature and twist

Considering non-linear terms in the curvature as well as in the twist, the governing boundary value problem for lateral bending of elastic, transverse loaded rods is formulated by means of Hamilton's principle. Using the method of small vibrations, the associated linearized equations of stability are derived, which complete the currently accepted relations. The example of the simplest lateral bending problem illustrates the improved effect of the proposed equations.

Wauer, J.↗
NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗
NASA NTRS2025

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗
NASA NTRS2020

Streak Instability Analysis for BOLT Configuration

The BOLT configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. A numerical study of the boundary layer instabilities over this configuration has been performed with the goal of applying state-of-the-art tools in multidimensional stability analysis that complement the other ongoing research related to the BOLT flight campaign. Specifically, the paper outlines modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane of a 1/3 rd scale BOLT configuration that was tested in the low-disturbance hypersonic wind tunnel at Purdue University and both quiet and conventional tunnels at Texas A&M University (TAMU). The basic state computations agree well with the heat flux measurements in the Purdue quiet tunnel as well as with the off-body velocity field described in a previous DNS study at the University of Minnesota. The encouraging agreement between the predicted frequency spectra of the streak instability in the region of increased boundary layer thickness and the surface pressure spectra measured during the TAMU experiments suggests that these streak instabilities may play an important role during the transition process near the symmetry plane. The modest value of the predicted peak N-factor (N max ≈ 5) for the streak instabilities provides a potential explanation for why transition is not observed under quiet flow conditions, but may occur in a noisy wind tunnel. Similar instability analyses are performed for the flight configuration at two selected flow conditions from the ascent portion of the planned flight trajectory. The predicted values of the peak N-factor suggest that transition onset in the vicinity of the center plane is possible at the ascent time of t = 26.98 seconds (N max ≈ 12.5 at X/L = 0.77), but the streak instabilities weaken sufficiently at the later time of t = 28.88 seconds (N max ≈ 11.5 at X/L ≈ 1.00) and are less likely to cause transition at that time.

Laminar turbulent transition↗
NASA NTRS1993

Role of secondary instability theory and parabolized stability equations in transition modeling

In modeling the laminar-turbulent transition region, the designer depends largely on benchmark data from experiments and/or direct numerical simulations that are usually extremely expensive. An understanding of the evolution of the Reynolds stresses, turbulent kinetic energy, and quantifies in the transport equations like the dissipation and production is essential in the modeling process. The secondary instability theory and the parabolized stability equations method are used to calculate these quantities, which are then compared with corresponding quantities calculated from available direct numerical simulation data for the incompressible boundary-layer flow of laminar-turbulent transition conditions. The potential of the secondary instability theory and the parabolized stability equations approach in predicting these quantities is discussed; results indicate that inexpensive data that are useful for transition modeling in the early stages of the transition region can be provided by these tools.

El-Hady, Nabil M.↗
NASA NTRS1986

Computational solid mechanics (finite elements and boundary elements) - Present status and future directions

Computational finite-element and boundary-element methods are reviewed, and their application to the mechanics of solids is discussed. Stability conditions for general FEMs are considered in addition to the use of least-order, stable, invariant, or hybrid/mixed isoparametric elements as alternatives to the displacement-based isoparametric elements. The use of symbolic manipulation, adaptive mesh refinement, transient dynamic response, and boundary-element methods for linear elaslticity and finite-strain problems of inelastic materials are also discussed.

Atluri, S. N.↗
DOE OSTI2021

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part II: Entropy Stability

In the hyperbolic research community, there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability. In the first part of the series, the application of simultaneous approximation terms for linear problems is investigated where the boundary conditions are imposed weakly. By applying this technique, the authors demonstrate that a pure continuous Galerkin scheme is indeed linearly stable if the boundary conditions are imposed in the correct way. In this work, we extend this investigation to the nonlinear case and focus on entropy conservation. Here, by switching to entropy variables, we provide an estimation of the boundary operators also for nonlinear problems, that guarantee conservation. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Numerical experiments on the stability of controlled shear flows

Three-dimensional, nonlinear simulations are presented for the parallel water boundary layer subjected to pressure gradient, suction, or heating control. The focus is on the secondary instability responsible for the K-type transition. Heating control is found to be the least effective in delaying the secondary instability. Flow-field details (including temperature profiles) are presented for both the uncontrolled boundary layer and the heated boundary layer.

Zang, T. A.↗
NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the nonlocal, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, such as those in the vicinity of roughness elements, steps, gaps, or corners. To compute the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 10^(7) to 10^(8) even in relatively simple flows. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear systems. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE. Additionally, inflow-resolvent analysis based on the HLNSE for discrete-roughness-induced disturbances in the nose tip of a blunt cone at Mach 6 demonstrates the importance of including the disturbance amplification along the near vicinity of the roughness element and separation region.

Boundary Layer Stability↗
NASA NTRS2020

Convolutional Neural Network for Transition Modeling Based on Linear Stability Theory

Transition prediction is an important aspect of aerodynamic design because of its impact on skin friction and potential coupling with flow separation characteristics. Traditionally, the modeling of transition has relied on correlation-based empirical formulas based on integral quantities such as the shape factor of the boundary layer. However, in many applications of computational fluid dynamics, the shape factor is not straightforwardly available or not well-defined. We propose using the complete velocity profile along with other quantities (e.g., frequency, Reynolds number) to predict the perturbation amplification factor. While this can be achieved with regression models based on a classical fully connected neural network, such a model can be computationally more demanding. We propose a novel convolutional neural network inspired by the underlying physics as described by the stability equations. Specifically, convolutional layers are first used to extract integral quantities from the velocity profiles, and then fully connected layers are used to map the extracted integral quantities, along with frequency and Reynolds number, to the output (amplification ratio). Numerical tests on classical boundary layers clearly demonstrate the merits of the proposed method. More importantly, we demonstrate that, for Tollmien-Schlichting instabilities in two-dimensional, low-speed boundary layers, the proposed network encodes information in the boundary layer profiles into an integral quantity that is strongly correlated to a well-known, physically defined parameter – the shape factor.

Laminar-turbulent transition↗
NASA NTRS1989

The influence of atmospheric stratification on scatterometer data

The effects of atmospheric stratification and the stability of the atmospheric stratification on the scatterometer data measuring surface winds over the ocean were investigated using the boundary layer model developed by Louis (1979). A variational analysis method is proposed, which allows direct assimilation of scatterometer data. It is shown that the effect of the stability of atmospheric stratification on the wind increment is relatively small. However, it is a systematic effect, and neglecting it would consistently underestimate the winds in stable regions.

Louis, Jean-Francois↗