SEARCH · Engineering Papers

Results for “Boundary Layer Stability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Interface of an uncoupled boundary layer algorithm with an inviscid core flow algorithm for unsteady supersonic engine inlets

An uncoupled boundary layer algorithm was combined with an inviscid core flow algorithm to model flows within supersonic engine inlets. The inviscid flow algorithm that was used was the LArge Perturbation INlet Code (LAPIN). The boundary layer and inviscid core flow algorithms were formulated in different manners. The boundary layer algorithm was two dimensional and solved in nonconservation form, while the core flow algorithm was one dimensional and solved in conservation form. In order to interface the two codes, the following modifications were important. The coordinate system was set up to maintain the parabolic nature of the boundary layer algorithm while approaching the one dimensional core flow solution far from a wall. The pressure gradient used in the boundary layer equation was calculated using the core flow values and the boundary layer equations, so the boundary layer solution smoothly approached the core flow values far from the wall. Flaring was used for the advection terms perpendicular to the core flow to maintain the stability of the algorithm. With these modifications, the combined viscous/inviscid algorithm matched well with experimental observations of pressure distributions with a supersonic inlet.

Darling, Douglas↗

NASA NTRS1992

Interface of an uncoupled boundary layer algorithm with an inviscid core flow algorithm for unsteady supersonic engine inlets

An uncoupled boundary layer algorithm was combined with an inviscid core flow algorithm to model flows within supersonic engine inlets. The inviscid flow algorithm that was used was the LArge Perturbation INlet Code (LAPIN). The boundary layer and inviscid core flow algorithms were formulated in different manners. The boundary layer algorithm was two dimensional and solved in nonconservation form, while the core flow algorithm was one dimensional and solved in conservation form. In order to interface the two codes, the following modifications were important. The coordinate system was set up to maintain the parabolic nature of the boundary layer algorithm while approaching the one dimensional core flow solution far from a wall. The pressure gradient used in the boundary layer equation was calculated using the core flow values and the boundary layer equations, so the boundary layer solution smoothly approached the core flow values far from the wall. Flaring was used for the advection terms perpendicular to the core flow to maintain the stability of the algorithm. With these modifications, the combined viscous/inviscid algorithm matched well experimental observations of pressure distributions with a supersonic inlet.

Darling, Douglas↗

NASA NTRS1992

Convection of tin in a Bridgman system. I - Flow characterization by effective diffusivity measurements

An electrochemical titration method was used to investigate the dynamic states in a cylindrical layer of convecting tin. The liquid tin was contained in a cell, with curved boundaries made of quartz and flat boundaries made of a solid state electrolyte - yttria-stabilized zirconia (YSZ). The electrolyte acted as a window through which a trace amount of oxygen could be pumped in or out by the application of a constant voltage. The concentration at the YSZ interface was monitored by operating the electrochemical cell in the galvanic mode. Experimentally determined effective diffusivities of oxygen were compared with the molecular diffusivity. Dynamic states in the convective flow were thus inferred. Temperature measurements were simultaneously made in order to identify the onset of oscillations from a steady convective regime. The experiments were conducted for two different aspect ratios for various imposed temperature gradients and two different orientations with respect to gravity. Transcritical states were identified and comparison to two-dimensional numerical models were made.

Sears, B.↗

NASA NTRS1993

Multigrid methods for flow transition in three-dimensional boundary layers with surface roughness

The efficient multilevel adaptive method has been successfully applied to perform direct numerical simulations (DNS) of flow transition in 3-D channels and 3-D boundary layers with 2-D and 3-D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semi-coarsening multigrid method associated with line distributive relaxation scheme, and an improved outflow boundary-condition treatment, which needs only a very short buffer domain to damp all order-one wave reflections, are developed. These approaches make the multigrid DNS code very accurate and efficient. This allows us not only to be able to do spatial DNS for the 3-D channel and flat plate at low computational costs, but also to do spatial DNS for transition in the 3-D boundary layer with 3-D single and multiple roughness elements, which would have extremely high computational costs with conventional methods. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments. The contribution of isolated and distributed roughness to transition is analyzed.

Liu, Chaoqun↗

DOE OSTI2024

Engineering 3D Floquet Codes by Rewinding

Floquet codes are a novel class of quantum error-correcting codes with dynamically generated logical qubits arising from a periodic schedule of noncommuting measurements. We utilize the interpretation of measurements in terms of condensation of topological excitations and the rewinding of measurement sequences to engineer new examples of Floquet codes. In particular, rewinding is advantageous for obtaining a desired set of instantaneous stabilizer groups on both toric and planar layouts. Our first example is a Floquet code with instantaneous stabilizer codes that have the same topological order as the three-dimensional (3D) toric code(s). This Floquet code also exhibits a splitting of the topological order of the 3D toric code under the associated sequence of measurements, i.e., an instantaneous stabilizer group of a single copy of the 3D toric code in one round transforms into an instantaneous stabilizer group of two copies of the 3D toric code up to nonlocal stabilizers in the following round. We further construct boundaries for this 3D code and argue that stacking it with two copies of the 3D subsystem toric code allows for a transversal implementation of the logical non-Clifford controlled-controlled- Z gate. We also show that the coupled-layer construction of the X-cube Floquet code can be modified by a rewinding schedule such that each of the instantaneous stabilizer codes is finite depth equivalent to the X-cube model up to toric codes; the X-cube Floquet code exhibits a splitting of the X-cube model into a copy of the X-cube model and toric codes under the measurement sequence. Our final 3D example is a generalization of the 2D Floquet toric code on the honeycomb lattice to three dimensions, which has instantaneous stabilizer codes with the same topological order as the 3D fermionic toric code. Published by the American Physical Society2024

Physics↗

DOE OSTI2021

Adaptive mesh based combustion simulations of direct fuel injection effects in a supersonic cavity flame-holder

Here, we present high-fidelity reacting simulations of a supersonic cavity flame-holder configuration. The focus of this work is on flame stabilization brought about by varying the location of fuel injection in a cavity stabilized supersonic flow of air. Central to our approach is a compressible multi-species reacting flow solver that uses adaptive-mesh-refinement (AMR), enabling the resolution of flame, shock-waves, boundary-layers, and small-scale structures in the computational domain. Our analysis indicates that fuel injection closer to the ramp at the aft end of the cavity allows for greater mixing and lower peak temperatures compared to fuel injection upstream that is closer to the backward facing step of the cavity. This difference is mainly due to greater turbulent fluctuations generated from the shear-layer towards the cavity ramp, thereby enhancing the mixing of fuel and air. A low frequency oscillatory behaviour in heat-release and pressure was also observed for the upstream injection case while a much higher-frequency phenomena was observed in the near-ramp injection case. By identifying the important physical determinants of the combustion processes, this study illustrates a promising pathway to design and optimize direct fuel injection strategies in supersonic cavity flame-holders that can improve flame stability, combustion efficiency, and reduce emissions.

42 ENGINEERING↗

NASA NTRS1987

The fully nonlinear development of Goertler vortices in growing boundary layers

The fully nonlinear development of small wavelength Goertler vortices in a growing boundary layer is investigated using a combination of asymptotic and numerical methods. The starting point for the analysis is the weakly nonlinear theory of Hall (1982b) who discussed the initial development of small amplitude vortices in a neighborhood of the location where they first become linearly unstable. That development is unusual in the context of nonlinear stability theory in that it is not described by the Stuart-Watson approach. In fact the development is governed by a pair of coupled nonlinear partial differential evolution equations for the vortex flow and the mean flow correction. Here the further development of this interaction is considered for vortices so large that the mean flow correction driven by them is as large as the basic state. Surprisingly it is found that such a nonlinear interaction can still be described by asymptotic means. It is shown that the vortices spread out across the boundary layer and effectively drive the boundary layer. In fact the system obtained by writing down the equations for the fundamental component of the vortex generate a differential equation for the basic state. Thus the mean flow adjusts so as to make these large amplitude vortices locally neutral. Moreover in the region where the vortices exist the mean flow has a square-root profile and the vortex velocity field can be written down in closed form. The upper and lower boundaries of the region of vortex activity are determined by a free-boundary problem involving the boundary layer equations. In general it is found that this region ultimately includes almost all of the original boundary layer and much of the free-stream. In this situation the mean flow has essentially no relationship to the flow which exists in the absence of the vortices.

Hall, Philip↗

NASA NTRS1978

Some unsteady separation problems for slender bodies

The unsteady Kutta-Joukowski condition, dynamic stall on oscillating airfoils, and unsteady shock wave-boundary layer interaction are discussed. Emphasis is placed on developing reliable prediction techniques and suppression of unsteady separation on oscillating control surfaces, wings, and rotating blades to improve aerodynamic stability.

Mccroskey, W. J.↗

NASA NTRS1985

Instability of time-periodic flows

The instabilities of some spatially and/or time-periodic flows are discussed, in particular, flows with curved streamlines which can support Taylor-Gortler vortices are described in detail. The simplest flow where this type of instability can occur is that due to the torsional oscillations of an infinitely long circular cylinder. For more complicated spatially varying time-periodic flows, a similar type of instability can occur and is spatially localized near the most unstable positions. When nonlinear effects are considered it is found that the instability modifies the steady streaming boundary layer induced by the oscillatory motion. It is shown that a rapidly rotating cylinder in a uniform flow is susceptible to a related type of instability; the appropriate stability equations are shown to be identical to those which govern the instability of a boussinesq fluid of Prandtl number unity heated time periodically from below.

Hall, P.↗

NASA NTRS1987

Instability of time-periodic flows

The instabilities of some spatially and/or time-periodic flows are discussed, in particular, flows with curved streamlines which can support Taylor-Gortler vortices are described in detail. The simplest flow where this type of instability can occur is that due to the torsional oscillations of an infinitely long circular cylinder. For more complicated spatially varying time-periodic flows, a similar type of instability can occur and is spatially localized near the most unstable positions. When nonlinear effects are considered it is found that the instability modifies the steady streaming boundary layer induced by the oscillatory motion. It is shown that a rapidly rotating cylinder in a uniform flow is susceptible to a related type of instability; the appropriate stability equations are shown to be identical to those which govern the instability of a Boussinesq fluid of Prandtl number unity heated time periodically from below.

Hall, Philip↗

NASA NTRS1989

Upwind algorithm for the parabolized Navier-Stokes equations

A new upwind algorithm based on Roe's scheme has been developed to solve the two-dimensional parabolized Navier-Stokes equations. This method does not require the addition of user-specified smoothing terms for the capture of discontinuities such as shock waves. Thus, the method is easy to use and can be applied without modification to a wide variety of supersonic flowfields. The advantages and disadvantages of this adaptation are discussed in relation to those of the conventional Beam-Warming (1978) scheme in terms of accuracy, stability, computer time and storage requirements, and programming effort. The new algorithm has been validated by applying it to three laminar test cases, including flat-plate boundary-layer flow, hypersonic flow past a 15-deg compression corner, and hypersonic flow into a converging inlet. The computed results compare well with experiment and show a dramatic improvement in the resolution of flowfield details when compared with results obtained using the conventional Beam-Warming algorithm.

Lawrence, Scott L.↗

NASA NTRS1993

Transition prediction and control in subsonic flow over a hump

The influence of a surface roughness element in the form of a two-dimensional hump on the transition location in a two-dimensional subsonic flow with a free-stream Mach number up to 0.8 is evaluated. Linear stability theory, coupled with the N-factor transition criterion, is used in the evaluation. The mean flow over the hump is calculated by solving the interacting boundary-layer equations; the viscous-inviscid coupling is taken into consideration, and the flow is solved within the separation bubble. The effects of hump height, length, location, and shape; unit Reynolds number; free-stream Mach number, continuous suction level; location of a suction strip; continuous cooling level; and location of a heating strip on the transition location are evaluated. The N-factor criterion predictions agree well with the experimental correlation of Fage; in addition, the N-factor criterion is more general and powerful than experimental correlations. The theoretically predicted effects of the hump's parameters and flow conditions on transition location are consistent and in agreement with both wind-tunnel and flight observations.

Masad, Jamal A.↗

NASA NTRS1994

Numerical simulation of two-dimensional spatially-developing mixing layers

Two-dimensional, incompressible, spatially developing mixing layer simulations are performed at Re = 10(exp 2) and 10(exp 4) with two classes of perturbations applied at the inlet boundary; combinations of discrete modes from linear stability theory, and a broad spectrum of modes derived from experimentally measured velocity spectra. The effect of the type and strength of inlet perturbations on vortex dynamics and time-averaged properties are explored. Two-point spatial velocity and autocorrelations are used to estimate the size and lifetime of the resulting coherent structures and to explore possible feedback effects. The computed time-averaged properties such as mean velocity profiles, turbulent statistics, and spread rates show good agreement with experimentally measured values. It is shown that by forcing with a broad spectrum of modes derived from an experimental energy spectrum many experimentally observed phenomena can be reproduced by a 2-D simulation. The strength of the forcing merely affected the length required for the dominant coherent structures to become fully-developed. Thus intensities comparable to those of the background turbulence in many wind tunnel experiments produced the same results, given sufficient simulation length.

Wilson, R. V.↗

NASA NTRS1994

Linear stability of a layered fluid with mobile surface plates

We develop a general method of calculating the linear stability of a fluid with homogeneous layers that is heated from below. The method employs a propagator technique to obtain expressions for the fluid velocity, stress, and temperature. The principal advantage of the method is the ease with which solutions are adapted to a wide variety of boundary conditions and fluid properties. We demonstrate the utility of the method using three examples which quantify the effects of (1) rheological layering, (2) mobile plates at the surface, and (3) multiple phase transitions. Each example is presented in the context of Earth's mantle. In the first example, we predict that convection becomes confined to the upper mantle once the viscosity increase between the upper and lower mantle exceeds a factor of 2000, consistent with the nonlinear calculations of Davies (1977). In the second example we find that the heat flux variations in a convecting fluid with variably sized, surface plates can be attributed, in part, to changes in the critical Rayleigh number. The linear stability of a fluid with multiple phase transitions is significantly affects by the locations of the transitions. We find that phase transitions have their largest effect when they are located at the center of the fluid layer and become much less important when they are located near the exterior boundaries.

Buffett, B. A.↗

NASA NTRS1994

Nonlinear stability of Halley comethosheath with transverse plasma motion

Weakly nonlinear Magneto Hydrodynamic (MHD) stability of the Halley cometosheath determined by the balance between the outward ion-neutral drag force and the inward Lorentz force is investigated including the transverse plasma motion as observed in the flanks with the help of the method of multiple scales. The eigenvalues and the eigenfunctions are obtained for the linear problem and the time evolution of the amplitude is obtained using the solvability condition for the solution of the second order problem. The diamagnetic cavity boundary and the adjacent layer of about 100 km thickness is found unstable for the travelling waves of certain wave numbers. Halley ionopause has been observed to have strong ripples with a wavelength of several hundred kilometers. It is found that nonlinear effects have stabilizing effect.

Srivastava, Krishna M.↗

NASA NTRS1986

Stability and natural vibration analysis of laminated plates by using a mixed element based on a refined plate theory

A mixed shear flexible finite element, with relaxed continuity, is developed for the geometrically linear and nonlinear analysis of layered anisotropic plates. The element formulation is based on a refined higher order theory which satisfies the zero transverse shear stress boundary conditions on the top and bottom faces of the plate and requires no shear correction coefficients. The mixed finite element developed herein consists of eleven degrees of freedom per node which include three displacements, two rotations and six moment resultants. The element is evaluated for its accuracy in the analysis of the stability and vibration of anisotropic rectangular plates with different lamination schemes and boundary conditions. The mixed finite element described here for the higher order theory gives very accurate results for buckling loads and natural frequencies.

Putcha, N. S.↗

DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗

NASA NTRS2016

Satellite-Surface Perspectives of Air Quality and Aerosol-Cloud Effects on the Environment: An Overview of 7-SEAS BASELInE

The objectives of 7-SEASBASELInE (Seven SouthEast Asian Studies Biomass-burning Aerosols and Stratocumulus Environment: Lifecycles and Interactions Experiment) campaigns in spring 2013-2015 were to synergize measurements from uniquely distributed ground-based networks (e.g., AERONET (AErosol RObotic NETwork)), MPLNET ( NASA Micro-Pulse Lidar Network)) and sophisticated platforms (e.g.,SMARTLabs (Surface-based Mobile Atmospheric Research and Testbed Laboratories), regional contributing instruments), along with satellite observations retrievals and regional atmospheric transport chemical models to establish a critically needed database, and to advance our understanding of biomass-burning aerosols and trace gases in Southeast Asia (SEA). We present a satellite-surface perspective of 7-SEASBASELInE and highlight scientific findings concerning: (1) regional meteorology of moisture fields conducive to the production and maintenance of low-level stratiform clouds over land; (2) atmospheric composition in a biomass-burning environment, particularly tracers-markers to serve as important indicators for assessing the state and evolution of atmospheric constituents; (3) applications of remote sensing to air quality and impact on radiative energetics, examining the effect of diurnal variability of boundary-layer height on aerosol loading; (4) aerosol hygroscopicity and ground-based cloud radar measurements in aerosol-cloud processes by advanced cloud ensemble models; and (5) implications of air quality, in terms of toxicity of nanoparticles and trace gases, to human health. This volume is the third 7-SEAS special issue (after Atmospheric Research, vol. 122, 2013; and Atmospheric Environment, vol. 78, 2013) and includes 27 papers published, with emphasis on air quality and aerosol-cloud effects on the environment. BASELInE observations of stratiform clouds over SEA are unique, such clouds are embedded in a heavy aerosol-laden environment and feature characteristically greater stability over land than over ocean, with minimal radar surface clutter at a high vertical spatial resolution. To facilitate an improved understanding of regional aerosol-cloud effects, we envision that future BASELInE-like measurement modeling needs fall into two categories: (1) efficient yet critical in-situ profiling of the boundary layer for validating remote-sensing retrievals and for initializing regional transport chemical and cloud ensemble models; and (2) fully utilizing the high observing frequencies of geostationary satellites for resolving the diurnal cycle of the boundary layerheight as it affects the loading of biomass-burning aerosols, air quality and radiative energetics.

AERONET↗