Search NASA⌕ Search

SEARCH · Search NASA

Results for “edge instability”

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 415 records · Page 23

Centrifugal Compressor Surge Margin Improved With Diffuser Hub Surface Air Injection

Aerodynamic stability is an important parameter in the design of compressors for aircraft gas turbine engines. Compression system instabilities can cause compressor surge, which may lead to the loss of an aircraft. As a result, engine designers include a margin of safety between the operating line of the engine and the stability limit line of the compressor. The margin of safety is typically referred to as "surge margin." Achieving the highest possible level of surge margin while meeting design point performance objectives is the goal of the compressor designer. However, performance goals often must be compromised in order to achieve adequate levels of surge margin. Techniques to improve surge margin will permit more aggressive compressor designs. Centrifugal compressor surge margin improvement was demonstrated at the NASA Glenn Research Center by injecting air into the vaned diffuser of a 4:1-pressure-ratio centrifugal compressor. Tests were performed using injector nozzles located on the diffuser hub surface of a vane-island diffuser in the vaneless region between the impeller trailing edge and the diffuser-vane leading edge. The nozzle flow path and discharge shape were designed to produce an air stream that remained tangent to the hub surface as it traveled into the diffuser passage. Injector nozzles were located near the leading edge of 23 of the 24 diffuser vanes. One passage did not contain an injector so that instrumentation located in that passage would be preserved. Several orientations of the injected stream relative to the diffuser vane leading edge were tested over a range of injected flow rates. Only steady flow (nonpulsed) air injection was tested. At 100 percent of the design speed, a 15-percent improvement in the baseline surge margin was achieved with a nozzle orientation that produced a jet that was bisected by the diffuser vane leading edge. Other orientations also improved the baseline surge margin. Tests were conducted at speeds below the design speed, and similar results were obtained. In most cases, the greatest improvement in surge margin occurred at fairly low levels of injected flow rate. Externally supplied injection air was used in these experiments. However, the injected flow rates that provided the greatest benefit could be produced using injection air that is recirculating between the diffuser discharge and nozzles located in the diffuser vaneless region. Future experiments will evaluate the effectiveness of recirculating air injection.

Skoch, Gary J.↗

Influence of a Backward-Facing Step on Swept-Wing Boundary-Layer Transition

Experimental measurements were performed on a swept flat-plate model with an airfoil leading edge and imposed chordwise pressure gradient to determine the effects of a backward-facing step on transition in a low-speed stationary crossflow-dominated boundary layer. Detailed hot-wire measurements were performed for three step heights ranging from 36 to 49% of the boundary-layer thickness at the step and corresponding to subcritical, nearly critical, and critical cases. In general, the step had a small localized effect on the growth of the stationary crossflow vortex, whereas the unsteady disturbance amplitudes increased with increasing step height. Intermittent spikes in instantaneous velocity began to appear for the two larger step heights. A physical explanation was provided for the mechanism leading to transition and the sudden movement in the transition front due to the critical steps. The large localized velocity spikes, which ultimately led to an intermittent breakdown of the boundary layer, were the result of nonlinear interactions of the different types of unsteady instabilities with each other and with the stationary crossflow vortices. Thus, the unsteady disturbances played the most important role in transition, but the stationary crossflow vortices also had a significant role via the modulation and the increased amplitude of the unsteady disturbances.

Eppink, Jenna L.↗

Ionospheric Convection in the Postnoon Auroral Oval: SuperDARN and Polar UVI Observations

Super Dual Auroral Radar Network (SuperDARN) observations, ultraviolet imaging from the Polar satellite (UVI), and particle precipitation data from DMSP satellites have been used to investigate the electrodynamics of the postnoon auroral oval in the Northern hemisphere. We show that: (1) For negative IMF By, the convection reversal (CR) was co-located with the maximum of auroral luminosity, but during positive IMF By the convection reversal was poleward of the auroral oval up to several degrees in latitude; (2) Postnoon auroral oval was associated with a large-scale upward field-aligned current (FAC) of the order of 6x10(exp -7). A m(exp -2) in magnitude (the FAC was inferred from the SuperDARN and UVI data). For negative IMF By, maximum of the auroral intensity coincides in latitude with the maximum of the upward field-aligned current. However, for positive IMF By. the maximum of the upward FAC was shifted to the poleward edge of the auroral oval; (3) In response to the IMF By turning from positive to negative, the maximum of the auroral luminosity did not change its position noticeably, but the position of the convection reversal changed considerably from 80-81 degs to about 76 degs MLAT, and the maximum of FAC moved from 77-78 degs to about 76 degs MLAT. Thus, after IMF By turns negative, both the FAC maximum and CR tend to coincide with the auroral maximum; (4) The IMF Bz positive deflection was followed by a decrease in both field-aligned current intensity and auroral luminosity. However, the decrease in the auroral luminosity lags behind the FAC decrease by about 12 min. Firstly, these observations allow us to suggest that the IMF By-related electric field can penetrate into the closed magnetosphere and produce convection and FAC changes in the region of the postnoon auroral oval. Secondly, we suggest that the interchange instability is a promising mechanism for the postnoon auroras.

Kozlovsky, A.↗

Sloshing of the Magnetized Cool Gas in the Cores of Galaxy Clusters

X-ray observations of many clusters of galaxies reveal the presence of edges in surface brightness and temperature, known as "cold fronts". In relaxed clusters with cool cores, these edges have been interpreted as evidence for the "sloshing" of the core gas in the cluster's gravitational potential. The smoothness of these edges has been interpreted as evidence for the stabilizing effect of magnetic fields "draped" around the front surfaces. To check this hypothesis, we perform high-resolution magnetohydrodynamics simulations of magnetized gas sloshing in galaxy clusters initiated by encounters with subclusters. We go beyond previous works on the simulation of cold fronts in a magnetized intracluster medium by simulating their formation in realistic, idealized mergers with high resolution ((Delta)x approx. 2 kpc). Our simulations sample a parameter space of plausible initial magnetic field strengths and field configurations. In the simulations, we observe strong velocity shears associated with the cold fronts amplifying the magnetic field along the cold front surfaces, increasing the magnetic field strength in these layers by up to an order of magnitude, and boosting the magnetic pressure up to near-equipartition with thermal pressure in some cases. In these layers, the magnetic field becomes strong enough to stabilize the cold fronts against Kelvin-Helmholtz instabilities, resulting in sharp, smooth fronts as those seen in observations of real clusters. These magnetic fields also result in strong suppression of mixing of high and low-entropy gas in the cluster, seen in our simulations of mergers in the absence of a magnetic field. As a result, the heating of the core due to sloshing is very modest and is unable to stave off a cooling catastrophe.

ZuHone, J. A.↗

Excitation of high-frequency electromagnetic waves by energetic electrons with a loss cone distribution in a field-aligned potential drop

The electron cyclotron maser instability (CMI) driven by momentum space anisotropy (df/dp (sub perpendicular) greater than 0) has been invoked to explain many aspects, such as the modes of propagation, harmonic emissions, and the source characteristics of the auroral kilometric radiation (AKR). Recent satellite observations of AKR sources indicate that the source regions are often imbedded within the auroral acceleration region characterized by the presence of a field-aligned potential drop. In this paper we investigate the excitation of the fundamental extraordinary mode radiation due to the accelerated electrons. The momentum space distribution of these energetic electrons is modeled by a realistic upward loss cone as modified by the presence of a parallel potential drop below the observation point. On the basis of linear growth rate calculations we present the emission characteristics, such as the frequency spectrum and the emission angular distribution as functions of the plasma parameters. We will discuss the implication of our results on the generation of the AKR from the edges of the auroral density cavities.

Fung, Shing F.↗

Transition in a Supersonic Boundary-Layer Due to Roughness and Acoustic Disturbances

The transition process induced by the interaction of an isolated roughness with acoustic disturbances in the free stream is numerically investigated for a boundary layer over a flat plate with a blunted leading edge at a free stream Mach number of 3.5. The roughness is assumed to be of Gaussian shape and the acoustic disturbances are introduced as boundary condition at the outer field. The governing equations are solved using the 5'h-~rder accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third- order total-variation-diminishing (TVD) Runge- Kutta scheme for time integration. The steady field induced by the two and three-dimensional roughness is also computed. The flow field induced by two-dimensional roughness exhibits different characteristics depending on the roughness heights. At small roughness heights the flow passes smoothly over the roughness, at moderate heights the flow separates downstream of the roughness and at larger roughness heights the flow separates upstream and downstream of the roughness. Computations also show that disturbances inside the boundary layer is due to the direct interaction of the acoustic waves and isolated roughness plays a minor role in generating instability waves.

Balakumar, P.↗

The clumpy circumstellar medium around young supernova remnants

Each of the youngest supernova remnants known in the Milky Way, Cas A, Kepler's SNR, and Tycho's SNR, shows a different morphological structure caused by different conditions in the progenitor stars and their surroundings. In all three cases, however, the observed shells have a thickness of about 1/4 the radius, a sharp outer edge, and significant brightness irregularities. These features require that the circumstellar medium be highly clumped. To investigate the phenomenon, models of the expansion have been constructed using a one-dimensional spherical hydrodynamic code. As a supernova shock moves down the external density gradient of the star, material behind the shock begins to go into free expansion. Then as surrounding material is encountered a reverse shock moving back into the ejectum will be formed. Until the expansion has swept up about eight times the ejected mass when the situation can be considered as a point explosion in its surroundings, the dynamics are controlled by conditions between the shocks. The region is also where the synchrotron radio emission from relativistic electrons trapped in magnetic fields arises. Initial particles and fields are accelerated and amplified by eddy motion at the interface between the ejected and swept-up material and at the boundaries of clumps. Polarimetry shows that these SNR have a net radial orientation of their magnetic fields apparently from stretching by Rayleigh-Taylor instabilities at the contact surfaces. Without clumps the observed shell is much too narrow and steep on the inside.

Dickel, J. R.↗

Intermittency and Cessation of Vortex Shedding in the Wake of a Thin Flat Plate with a Circular Trailing Edge

The near and very near wake of thin flat plates with both sharp and circular trailing edges (TEs) are investigated with direct numerical simulations (DNSs). The TE is circular in two of the cases (IN & NS) and sharp in one of them (ST). The separating boundary layers are turbulent in all cases. The objectives of this study are twofold. The first is to explore the effect of significantly reducing Re(sub D) (Reynolds number based on circular TE diameter, D) on the flow in the TE region, and the shedding process (Cases IN and NS). The second is to better understand the reasons underlying the findings of an earlier experimental wake investigation (sharp TE) where (1) the center-line values in normal intensity, and the peak in shear stress profiles in the cross-stream direction, were found to first increase in the streamwise direction (x), from that obtained at the TE, before diminishing further downstream, and (2) a broadband peak was observed in centerline cross-stream velocity (v) spectra (indicating quasi-periodicity, possibly due to vortices or wave-like motions). Case ST from the present study showed a near wake instability resulting in spanwise vortices (with a streamwise component). The instability is intermittent and contributes to both the broadband peak in the v spectrum and the initial increases in normal intensity and shear stress (as in the experiment). Case NS, with the lower value of ReD is an "essentially" non-shedding case where the flow in the TE region continually changes direction (upward/downward) because of turbulence. Case IN, with twice the value of ReD as Case NS, also exhibits a swaying motion in the TE region. In addition, vortex shedding is initiated during periods when the flow direction changes rapidly. Shedding in this case is intermittent. It results in a peak in the v spectrum obtained at the centerline (x/D = 1.0).

direct numerical simulations↗

Evolution of an Atmospheric Kármán Vortex Street From High‐Resolution Satellite Winds: Guadalupe Island Case Study

Vortex streets formed in the stratocumulus‐capped wake of mountainous islands are the atmospheric analogues of the classic Kármán vortex street observed in laboratory flows past bluff bodies. The quantitative analysis of these mesoscale unsteady atmospheric flows has been hampered by the lack of satellite wind retrievals of sufficiently high spatial and temporal resolution. Taking advantage of the cutting‐edge Advanced Baseline Imager, we derived kilometer‐scale cloud‐motion winds at 5‐min frequency for a vortex street in the lee of Guadalupe Island imaged by Geostationary Operational Environmental Satellite‐16. Combined with Moderate Resolution Imaging Spectroradiometer data, the geostationary imagery also provided accurate stereo cloud‐top heights. The time series of geostationary winds, supplemented with snapshots of ocean surface winds from the Advanced Scatterometer, allowed us to capture the wake oscillations and measure vortex shedding dynamics. The retrievals revealed a markedly asymmetric vortex decay, with cyclonic eddies having larger peak vorticities than anticyclonic eddies at the same downstream location. Drawing on the vast knowledge accumulated about laboratory bluff body flows, we argue that the asymmetric island wake arises from the combined effects of Earth's rotation and Guadalupe's nonaxisymmetric shape resembling an inclined flat plate at low angle of attack. However, numerical simulations will need to establish whether or not the selective destabilization of the shallow atmospheric anticyclonic eddies is caused by the same mechanisms that destabilize the deep columnar anticyclones of laboratory flows, such as three‐dimensional vertical perturbations due to centrifugal or elliptical instabilities.

Horváth, Á.↗

Planar Strain-Rate-Free Diffusion Flames: Initiation, Properties, and Extinction

An effectively strain-rate-free diffusion flame constitutes the most vigorous laminar combustion of initially unmixed reactive gases. Such a diffusion flame is characterized by a relatively long residence time and by a relatively large characteristic length scale. If such a flame were also planar, providing high symmetry, it would be particularly suitable for experimental and theoretical investigations of key combustion phenomena, such as multicomponent diffusion, chemical kinetics, and soot inception, growth, and oxidation. Unfortunately, a planar strain-rate-free diffusion flame is highly disrupted in earth-gravity (e.g., in a counterflow-diffusion-flame apparatus) because of the very rapid onset (approx. 100 ms) of gravity-induced instability. Accordingly, a specially dedicated apparatus was designed, fabricated, and initially checked out for the examination of a planar strain-rate-free diffusion flame in microgravity. Such a diffusion flame may be formed within a hollowed-out squat container (initially configured as 25 cm x 25 cm x 9 cm), with isothermal, noncatalytic, impervious walls. At test initiation, a thin metallic sheet (approx. 1 mm in thickness) that separates the internal volume into two equal portions, each of dimensions 25 cm x 25 cm x 4.5 cm, is withdrawn, by uniform translation (approx. 50 cm/s) in its own plane, through a tightly fitting slit in one side wall. Thereupon, diluted fuel vapor (initially confined to one half-volume of the container) gains access to diluted oxygen (initially with the same pressure, density, and temperature as the fuel, but initially confined to the other half-volume). After a brief delay (approx. 10 ms), to permit limited but sufficient-for-flammability diffusional interpenetration of fuel vapor and oxidizer, burning is initiated by discharge of a line igniter, located along that side wall from which the trailing edge of the separator withdraws. The ignition spawns a triple-flame propagation across the 25 cm x 25 cm centerplane. When a diffusion flame is emplaced in the centerplane, any subsequent travel, and change in temperature, of that planar diffusion flame may be tracked, along with the effectively spatially uniform but temporally evolving pressure within the container. Eventually, nearly complete depletion of the stoichiometrically deficient reactant, along with heat loss to the container surfaces, effects extinction. These data afford an opportunity to check theoretical models of diffusion and chemical kinetics under conditions ranging from intense burning to flame out, or, alternatively, to evolve simple empirical representations of these phenomena. Thus, the project sought to utilize microgravity testing to elucidate commonly encountered phenomenology, arising in the commonly-encountered mode of combustion (whether related to heating, manufacturing, boiling, and propulsion, or to uncontrolled, free-burning fire in structures and wildland vegetation), of those commonly utilized fuels usually categorized as gaseous fuels (such as hydrogen, natural gas, and propane, which are gaseous under atmospheric conditions).

Fendell, Francis↗

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗

Dynamics of Liquids in Edges and Corners (DYLCO): IML-2 Experiment for the BDPU

Knowledge of the behavior of fluids possessing free surfaces is important to many fluid systems, particularly in space, where the normally subtle effects of surface wettability play a more dramatic and often surprising role. DYLCO for the IML-2 mission was proposed as a simple experiment to probe the particular behavior of capillary surfaces in containers of irregular cross section. Temperature control was utilized to vary the fluid-solid contact angle, a questionable thermodynamic parameter of the system, small changes in which can dramatically influence the configuration, stability, and flow of a capillary surface. Container shapes, test fluid, and temperature ranges were selected for observing both local changes in interface curvature as well as a global change in fluid orientation due to a critical wetting phenomenon. The experiment hardware performed beyond what was expected and fluid interfaces could be readily digitized post flight to show the dependence of the interface curvature on temperature. For each of the containers tested surfaces were observed which did not satisfy the classic equations for the prediction of interface shape with constant contact angle boundary condition. This is explained by the presence of contact angle hysteresis arising from expansion and contraction of the liquid during the heating and cooling steps of the test procedure. More importantly, surfaces exceeding the critical surface curvature required for critical wetting were measured, yet no wetting was observed. These findings are indeed curious and pose key questions concerning the role of hysteresis for this critical wetting phenomena. The stability of such surfaces was determined numerically and it is shown that stability is enhance (reduced) when a surface is in its 'advancing' ('receding') state, The analysis shows complete instability as the critical wetting condition is reached. The case of ideal dynamic wetting is addressed analytically in detail with results of significant flow characteristics presented in closed form. The solutions indicate a square root of T dependence of the capillary 'rise' rate which is corroborated by drop tower tests. The analysis clearly shows that infinite time is necessary for surfaces to reorient at the critical wetting transition.

Langbein, D.↗

High speed transition prediction

The main objective of this work period was to develop, acquire and apply state-of-the-art tools for the prediction of transition at high speeds at NASA Ames. Although various stability codes as well as basic state codes were acquired, the development of a new Parabolized Stability Equation (PSE) code was minimal. The time that was initially allocated for development was used on other tasks, in particular for the Leading Edge Suction problem, in acquiring proficiency in various graphics tools, and in applying these tools to evaluate various Navier-Stokes and Euler solutions. The second objective of this work period was to attend the Transition and Turbulence Workshop at NASA Langley in July and August, 1991. A report on the Workshop follows. From July 8, 1991 to August 2, 1991, the author participated in the Transition and Turbulence Workshop at NASA Langley. For purposes of interest here, analysis can be said to consist of solving simplified governing equations by various analytical methods, such as asymptotic methods, or by use of very meager computer resources. From the composition of the various groups at the Workshop, it can be seen that analytical methods are generally more popular in Great Britain than they are in the U.S., possibly due to historical factors and the lack of computer resources. Experimenters at the Workshop were mostly concerned with subsonic flows, and a number of demonstrations were provided, among which were a hot-wire experiment to probe the boundary layer on a rotating disc, a hot-wire rake to map a free shear layer behind a cylinder, and the use of heating strips on a flat plate to control instability waves and consequent transition. A highpoint of the demonstrations was the opportunity to observe the rather noisy 'quiet' supersonic pilot tunnel in operation.

Gasperas, Gediminis↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Phase Transition in Dusty Plasmas: A Microphysical Description

Dust grains immersed in plasma discharges acquire a large negative charge and settle into a dust cloud at the edge of the sheath. In this region, the plasma ions stream toward the electrode at a velocity u approx. cs=(T(sub e)/m(sub i))(exp 1/2). Experimentally at sufficiently high gas pressure P, the random kinetic energy of the grains is damped by gas friction, and the grains are strongly coupled and self-organize into a crystalline configuration. For lower pressures despite the dissipation of grain kinetic energy to gas friction, the dust grains reach a steady-state kinetic temperature T(sub d) which is much larger than the temperature of any other component in the plasma. T(sub d) is so large that the dust acts like a fluid. We have used the dynamically shielded dust (DSD) model to simulate these physical processes. We find that the known experimental features are nicely reproduced in the simulations, and that additional features are revealed. In the figure we plot the variation of T(sub d) as P is continuously varied in a DSD code run. A marked difference is evident between the critical pressure P(sub m) for the melting transition as P is decreased, and the critical pressure P(sub c) for the condensation transition as P is increased. For P(sub m) is less than P is less than P(sub c), mixed phase states are seen. This hysteresis occurs because the instability which triggers melting is different from the instability that heats the dust in the fluid phase and inhibits freezing. At low pressure, the dust is subject to a two-stream instability with the ions. This instability is responsible for the high temperature of the dust at low pressure. The basic physics underlying the melting transition has been elucidated in a series of papers. We are developing a first-principles analytic approach to the melting transition, which embodies the same physics that is present in the DSD code.

Joyce, Glenn↗

Aeroelastic Airworthiness Assesment of the Adaptive Compliant Trailing Edge Flaps

The Adaptive Compliant Trailing Edge (ACTE) demonstrator is a joint task under the National Aeronautics and Space Administration Environmentally Responsible Aviation Project in partnership with the Air Force Research Laboratory and FlexSys, Inc. (Ann Arbor, Michigan). The project goal is to develop advanced technologies that enable environmentally friendly aircraft, such as adaptive compliant technologies. The ACTE demonstrator flight-test program encompassed replacing the Fowler flaps on the SubsoniC Aircraft Testbed, a modified Gulfstream III (Gulfstream Aerospace, Savannah, Georgia) aircraft, with control surfaces developed by FlexSys. The control surfaces developed by FlexSys are a pair of uniquely-designed unconventional flaps to be used as lifting surfaces during flight-testing to validate their structural effectiveness. The unconventional flaps required a multidisciplinary airworthiness assessment to prove they could withstand the prescribed flight envelope. Several challenges were posed due to the large deflections experienced by the structure, requiring non-linear analysis methods. The aeroelastic assessment necessitated both conventional and extensive testing and analysis methods. A series of ground vibration tests (GVTs) were conducted to provide modal characteristics to validate and update finite element models (FEMs) used for the flutter analyses for a subset of the various flight configurations. Numerous FEMs were developed using data from FlexSys and the ground tests. The flap FEMs were then attached to the aircraft model to generate a combined FEM that could be analyzed for aeroelastic instabilities. The aeroelastic analysis results showed the combined system of aircraft and flaps were predicted to have the required flutter margin to successfully demonstrate the adaptive compliant technology. This paper documents the details of the aeroelastic airworthiness assessment described, including the ground testing and analyses, and subsequent flight-testing performed on the unconventional ACTE flaps.

aeroelasticity↗

A Vortical Dawn Flank Boundary Layer for Near-Radial IMF: Wind Observations on 24 October 2001

We present an example of a boundary layer tailward of the dawn terminator which is entirely populated by rolled-up flow vortices. Observations were made by Wind on 24 October 2001 as the spacecraft moved across the region at the X plane approximately equal to −13 Earth radii. Interplanetary conditions were steady with a near-radial interplanetary magnetic field (IMF). Approximately 15 vortices were observed over the 1.5 hours duration of Wind's crossing, each lasting approximately 5 min. The rolling up is inferred from the presence of a hot tenuous plasma being accelerated to speeds higher than in the adjoining magnetosheath, a circumstance which has been shown to be a reliable signature of this in single-spacecraft observations. A blob of cold dense plasma was entrained in each vortex, at whose leading edge abrupt polarity changes of field and velocity components at current sheets were regularly observed. In the frame of the average boundary layer velocity, the dense blobs were moving predominantly sunward and their scale size along the X plane was approximately 7.4 Earth radii. Inquiring into the generation mechanism of the vortices, we analyze the stability of the boundary layer to sheared flows using compressible magnetohydrodynamic Kelvin-Helmholtz theory with continuous profiles for the physical quantities. We input parameters from (i) the exact theory of magnetosheath flow under aligned solar wind field and flow vectors near the terminator and (ii) the Wind data. It is shown that the configuration is indeed Kelvin-Helmholtz (KH) unstable. This is the first reported example of KH-unstable waves at the magnetopause under a radial IMF.

Magnetosheath↗

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗