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 901 records · Page 50

Accurate calculation and instability of supersonic wake flows

This study is concerned with the computation and linear stability of a class of laminar compressible wake flows. The emphasis is on correct basic flow profiles that satisfy the steady equations of motion, and to this end the unperturbed state is obtained through numerical integration of the compressible boundary-layer equations. The linear stability of the flow is examined via the Rayleigh equation that describes evolution of inviscid disturbances. Analytical results are given for short- and long-wavelength disturbances and some numerical results of the general eigenvalue problem are also reported.

Papageorgiou, Demetrius T.↗

Effects of streamwise variations in noise levels and spectra on supersonic boundary-layer transition

Transition data for sharp cones in two quiet wind tunnels at Mach numbers 3.5 and 5.0 have been correlated in terms of noise parameters with data from several conventional wind tunnels and from the flight data for the AEDC transition cone. The noise parameters were developed to account for the large axial variations of the rms stream noise and the high frequency noise spectra that occurred in the quiet tunnels for some test conditions. The correlation results indicated transition in the quiet tunnels was dominated by the local stream noise that was incident on the cone boundary layer upstream of the neutral stability point. The correlation results also suggested that the energy in high frequency components of the quiet tunnel noise spectra had significant adverse effects on transition when the noise was incident on the boundary layer both upstream and downstream of the neutral stability point.

Chen, F.-J.↗

Study of Stratus-Lowering Marine-Fog Events Observed During C-FOG

Two stratus lowering marine fog events observed on 28 September and 04 October 2018 during the Coastal Fog (C-Fog) field campaign that took place offshore of Eastern Canada during 1 September to 6 October 2018 are described. In-situ, profiling and remote sensing observations were made at selected land sites in eastern Newfoundland (NL) and Nova Scotia (NS) as well as aboard the research vessel (R/V) Hugh R. Sharp that cruised in adjoining coastal waters. Synoptic scale analysis showed that both fog episodes were an outcome of the interaction between synoptic-scale surface level low-pressure systems and a contiguous high-pressure system. At the same time, back trajectories revealed that the bulk of the fog layer is formed due to differential advection. Here, the diameter of the fog droplets at the surface gradually decreased from the centre of the fog layer to its leading/trailing edges. The bimodal fog droplet diameter distribution with peaks at 5-10 µm and 20-25 µm provided clues on droplet collision and coalescence processes. The observed difference between microphysical variables and droplet distribution between the two fog events and within same fog layer may have been governed by atmospheric boundary layer conditions (e.g., humidity conditions and turbulence) prevailed in the fog layer. Overall, it is concluded that the life cycle of observed stratus-lowering coastal-fog episodes is dependent on synoptic conditions as well as atmospheric boundary layer characteristics such as stability, cloud top cooling and entrainment.

54 ENVIRONMENTAL SCIENCES↗

Evolving a Puncture Black Hole with Fixed Mesh Refinement

We present a detailed study of the effects of mesh refinement boundaries on the convergence and stability of simulations of black hole spacetimes. We find no technical problems. In our applications of this technique to the evolution of puncture initial data, we demonstrate that it is possible to simulaneously maintain second order convergence near the puncture and extend the outer boundary beyond 100M, thereby approaching the asymptotically flat region in which boundary condition problems are less difficult.

Imbiriba, Breno↗

Turbulence and heat stratification

With applications to meteorology in view, the method of small oscillations is used to investigate the stability of a plane laminar boundary layer in which the density depends in a special way on the distance from the plate. It is found that for a constant Froude number the flow becomes more stable as the Richardson number increases until, at a critical value of the Richardson number that depends on the Froude number, the flow is stable for all disturbance wave lengths and Reynolds numbers. The calculated results are verified by experiment.

FLOW VELOCITY - DISTRIBUTION↗

Transition In Supersonic Flows With Corners

It is proposed to investigate the effect of sudden change in the slope of the surface on the stability and transition onset point in supersonic and hypersonic boundary layers. Since we can not use the linear stability theory in flows where there are discontinuities we have to solve the full Navier-Stokes equations to determine the amplification of the disturbances across the oblique shock which exists at the corner. We will use the Navier-Stokes code which was developed by Dr. Harold Atkins in the aerodynamic and acoustic methods branch. The code is developed using the higher order ENO scheme and is currently set up for two-dimensional flows. In the first year of our investigation, we will investigate the evolution of small amplitude two-dimensional disturbances and determine the N-factors for the most amplified disturbances. We will perform the analysis for the parameter similar to X-31 model and flight experiment conditions. From these computation we will infer what is the maximum amplification rate possible in quite environments without any tripping devices. In the second part of the investigation, we will perform nonlinear computations to determine what is the minimum amplitudes necessary to cause the transition at a designed location. In the third part, we have to determine what kind of trip devices necessary to excite the required disturbances. The principal investigator of this project is Dr. P. Balakumar, Associate Professor of Department of Aerospace Engineering, Old Dominion University. We will actively collaborate with research scientists in the Aero-Thermodynamics Branch at NASA/Langley Research Center. A half-time graduate student will also participate in the project.

Balakumar, P.↗

The Kelvin–Helmholtz Instability at the Boundary of Relativistic Magnetized Jets

We study the linear stability of a planar interface separating two fluids in relative motion, focusing on conditions appropriate for the boundaries of relativistic jets. The jet is magnetically dominated, whereas the ambient wind is gas-pressure-dominated. We derive the most general form of the dispersion relation and provide an analytical approximation of its solution for an ambient sound speed much smaller than the jet Alfvén speed v A , as appropriate for realistic systems. The stability properties are chiefly determined by the angle ψ between the wavevector and the jet magnetic field. For ψ = π/2, magnetic tension plays no role, and our solution resembles the one of a gas-pressure-dominated jet. Here, only sub-Alfvénic jets are unstable ($0\lt {M}_{e}\equiv (v/{v}_{{\rm{A}}})\cos \theta \lt 1$ , where v is the shear velocity and θ the angle between the velocity and the wavevector). For ψ = 0, the free energy in the velocity shear needs to overcome the magnetic tension, and only super-Alfvénic jets are unstable ($1\lt {M}_{e}\lt \sqrt{{(1+{{\rm{\Gamma }}}_{w}^{2})/[1+({v}_{{\rm{A}}}/c)}^{2}{{\rm{\Gamma }}}_{w}^{2}]}$, with Γ w the wind adiabatic index). Our results have important implications for the propagation and emission of relativistic magnetized jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal Growth in Hypersonic Boundary Layers

The linear form of the parabolized linear stability equations is used in a variational approach to extend the previous body of results for the optimal, nonmodal disturbance growth in boundary-layer flows. This paper investigates the optimal growth characteristics in the hypersonic Mach number regime without any high-enthalpy effects. The influence of wall cooling is studied, with particular emphasis on the role of the initial disturbance location and the value of the spanwise wave number that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary-layer equations, mean flow solutions based on the full Navier-Stokes equations are used in select cases to help account for the viscous- inviscid interaction near the leading edge of the plate and for the weak shock wave emanating from that region. Using the full Navier-Stokes mean flow is shown to result in further reduction with Mach number in the magnitude of optimal growth relative to the predictions based on the self-similar approximation to the base flow.

Navier Stokes equations↗

Stationary disturbances in three-dimensional boundary layers over concave surfaces

A two-dimensional boundary layer on a concave surface is known to be susceptible to centrifugal instability which manifests itself in the form of stationary streamwise counter-rotating vortices commonly known as Goertler vortices. In this paper, the problem of the stability of a three-dimensional boundary layer on a concave surface is considered. Linear stability equations, including streamline and surface curvature effects, are solved for an infinitely swept wing. The results indicate that the Goetler vortex structure begins to transform into co-rotating vortices when the crossflow Reynolds number is increased. This transition is almost complete when the crossflow Reynolds number is in excess of 45. It is shown that the centrifugal effects destabilize the crossflow vortices. Some comparisons with available experimental results are made.

Collier, F. S., Jr.↗

Guidance Enhancements and Performance Assessments for the Mars Ascent Vehicle Spin-Stabilized Upper Stage Configuration

he objective of the Mars Sample Return (MSR) campaign is to return samples from the surface of Mars to Earth for research. As one element of the MSR campaign, the Mars Ascent Vehicle (MAV) is responsible for transporting the samples from the surface of Mars to a Low-Martian Orbit (LMO) for retrieval. Complete autonomy is required throughout ascent, and orbital insertion is constrained by tight dispersion boundaries. An unguided, spin-stabilized second stage for MAV has been selected over a guided upper-stage to drive mass savings and reduce overall MSR campaign mass risk, at the cost of reduced GNC capability. To address this design change, the MAV GNC team has derived a robust prediction algorithm, building on previous energy management schemes, that solves for a single inertial pointing direction solution for the spin-stabilized 2nd stage burn. Algorithm stability is explored that compared to previous versions of the algorithm. Also, a set of analytical partials was developed to study MAV’s dispersed orbital insertion performance with respect to MAV system uncertainties. These partials were verified through simulation analysis and prove useful for analytical insight into the dynamics of MAV during the 2nd stage maneuver.

Jason M Everett↗

Guidance Enhancements and Performance Assessments for the Mars Ascent Vehicle Spin-Stabilized Upper Stage Configuration

The objective of the Mars Sample Return (MSR) campaign is to return samples from the surface of Mars to Earth for research. As one element of the MSR campaign, the Mars Ascent Vehicle (MAV) is responsible for transporting the samples from the surface of Mars to a Low-Martian Orbit (LMO) for retrieval. Complete autonomy is required throughout ascent, and orbital insertion is constrained by tight dispersion boundaries. An unguided, spin-stabilized second stage for MAV has been selected over a guided upper-stage to drive mass savings and reduce overall MSR campaign mass risk, at the cost of reduced GNC capability. To address this design change, the MAV GNC team has derived a robust prediction algorithm, building on previous energy management schemes, that solves for a single inertial pointing direction solution for the spin-stabilized 2nd stage burn. Algorithm stability is explored that compared to previous versions of the algorithm. Also, a set of analytical partials was developed to study MAV’s dispersed orbital insertion performance with respect to MAV system uncertainties. These partials were verified through simulation analysis and prove useful for analytical insight into the dynamics of MAV during the 2nd stage maneuver.

GNC↗

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

An interaction solution algorithm for viscous aerodynamic performance prediction

A weak-interaction solution algorithm is established for aerodynamic flow field prediction about an isolated airfoil. It requires numerical solution of differential equations governing potential flow, viscous and turbulent boundary layer flow, and the turbulent wake flow downstream of a trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. These in turn alter the viscous flow through imposed pressure gradients. Closure for turbulence is accomplished using a second order model. Numerical evaluations assess factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems as solved using a finite element algorithm.

Baker, A. J.↗

Grain boundary energy control in zinc aluminate nanoceramics

This study investigates the grain boundary energy dependence on segregated dopants in nanocrystalline zinc aluminate ceramics. Atomistic simulations of Σ3 and Σ9 grain boundaries showed that trivalent ions of varying ionic radii [Sc 3+ (74.5 pm), In 3+ (80.0 pm), Y 3+ (90.0 pm), and Nd 3+ (98.3 pm)] have a tendency to segregate to both interfaces, with Y 3+ presenting the highest segregation potentials. The connection between segregation and the reduction of interfacial energies was explored by measuring the grain boundary energy on nanoceramics fabricated via high‐pressure spark plasma sintering (HP‐SPS) using differential scanning calorimetry (DSC). The results revealed that Y 3+ doping at 0.5 mol% reduces the grain boundary energy in zinc aluminate nanoceramics from 1.1–1.3 J/m 2 to 0.6–0.8 J/m 2 ; the range correlates with the observed size dependence of the excess energy, with higher values observed for the smaller grain sizes (∼17 nm). The noted decrease in interfacial energies for doped samples suggests it is indeed possible to alter the stability of zinc aluminate grain boundaries via dopant segregation.

36 MATERIALS SCIENCE↗

The Effects of Acoustic Treatment on Pressure Disturbances From a Supersonic Jet in a Circular Duct

The pressure disturbances generated by an instability wave in the shear layer of a supersonic jet are studied for an axisymmetric jet inside a lined circular duct. For the supersonic jet, locally linear stability analysis with duct wall boundary conditions is used to calculate the eigenvalues and the eigenfunctions at each axial location. These values are used to determine the growth rates and phase velocities of the instability waves and the near field pressure disturbance patterns. The study is confined to the dominant Kelvin-Helmholtz instability mode and to the region just downstream of the nozzle exit where the shear layer is growing but is still small in size compared to the radius of the duct. Numerical results are used to study the effects of changes in the outer flow, growth in the shear layer thickness, wall distance, and wall impedance, and the effects of these changes on non-axisymmetric modes. The primary results indicate that the effects of the duct wall on stability characteristics diminish as the outer flow increases and as the jet azimuthal mode number increases. Also, wall reflections are reduced when using a finite impedance boundary condition at the wall; but in addition, reflections are reduced and growth rates diminished by keeping the imaginary part of the impedance negative when using the negative exponential for the harmonic dependence.

Dahl, Milo D.↗

Spatially Developing Secondary Instabilities and Attachment Line Instability in Supersonic Boundary Layers

This paper reports on progress towards developing a spatial stability code for compressible shear flows with two inhomogeneous directions, such as crossflow dominated swept-wing boundary layers and attachment line flows. Certain unique aspects of formulating a spatial, two-dimensional eigenvalue problem for the secondary instability of finite amplitude crossflow vortices are discussed. A primary test case used for parameter study corresponds to the low-speed, NLF-0415(b) airfoil configuration as tested in the ASU Unsteady Wind Tunnel, wherein a spanwise periodic array of roughness elements was placed near the leading edge in order to excite stationary crossflow modes with a specified fundamental wavelength. The two classes of flow conditions selected for this analysis include those for which the roughness array spacing corresponds to either the naturally dominant crossflow wavelength, or a subcritical wavelength that serves to reduce the growth of the naturally excited dominant crossflow modes. Numerical predictions are compared with the measured database, both as indirect validation for the spatial instability analysis and to provide a basis for comparison with a higher Reynolds number, supersonic swept-wing configuration. Application of the eigenvalue analysis to the supersonic configuration reveals that a broad spectrum of stationary crossflow modes can sustain sufficiently strong secondary instabilities as to potentially cause transition over this configuration. Implications of this finding for transition control in swept wing boundary layers are examined. Finally, extension of the spatial stability analysis to supersonic attachment line flows is also considered.

Li, Fei↗

Effect of 3D Roughness Patch on Instability Amplification in a Supersonic Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of high-speed vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, sinusoidal shape roughness elements on a Mach 3.5 flat plate boundary layer for flow conditions corresponding to the planned conditions of an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at the NASA Langley Research Center. Analysis of convective instabilities in the wake of the roughness patch was reported in a previous paper and the current work extends that analysis to instability amplification across the length of the roughness patch. Quasi-parallel stability analysis of the modified boundary layer flow over the patch indicates two dominant families of unstable disturbances, namely, a group of high frequency modes that amplify in localized regions along the roughness patch and another group of low frequency modes that have smaller peak amplification rates but amplify steadily both above the roughness patch and in the wake region behind it. The results suggest that the amplification factors associated with the high-frequency modes are sufficiently low, at least for the roughness patches considered in this paper, so that these modes are unlikely to have a major influence on the transition process. The amplification of the low-frequency modes within the region of the roughness patch is further quantified via direct numerical simulations. Results confirm the strongly destabilizing influence of the roughness patch on the first mode instabilities, yielding an N-factor increment of ΔN ≈ 3.6 for a roughness patch length of eight wavelengths.

Laminar turbulent transition↗

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗