Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discontinuity Factors”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Controls on Variations of Surface Energy, Water, and Carbon Budgets within Large-Scale Amazon Basin

A key research focus of the LBA Research Program is understanding the space-time variations in interlinked surface energy, water, and carbon budgets, the controls on these variations, and the implications of these controls on the carbon sequestering capacity of the large scale forest-pasture system that dominates the Amaz6nia landscape. Quantification of these variations and controls are investigated by a combination of in situ measurements, remotely sensed measurements from space, and a realistically forced hydrometeorological model coupled to a carbon assimilation model, capable of simulating details within the surface energy and water budgets along with the principle processes of photosynthesis and respiration. Herein we describe the results of an investigation concerning the space-time controls of carbon sources and sinks distributed over the large scale Amazon basin. The results are derived from a carbon-water-energy budget retrieval system for the large scale Amazon basin, which uses a coupled carbon assimilation-hydrometeorological model as an integrating system, forced by both in situ meteorological measurements and remotely sensed radiation and precipitation fluxes obtained from a combination of GOES, SSM/I, TOMS, and TRh4M satellite measurements. Results include validation of (a) retrieved surface radiation and precipitation fluxes based on 30-min averaged surface measurements taken at Ji-Parani in Rondania and Manaus in Amazonas, and (b) modeled sensible, latent, and C02 fluxes based on tower measurements taken at Reserva Jaru, Manaus and Fazenda Nossa Senhora. The space-time controls on carbon sequestration are partitioned into sets of factors classified by: (1) above canopy meteorology, (2) incoming surface radiation, (3) precipitation interception, and (4) indigenous stomatal processes varied over the different land covers of pristine rainforest, partially, and fully logged rainforests, and pasture lands. These are the principle meteorological, thermodynamical, hydrological, and biophysical control paths which perturb net carbon fluxes and sequestration, produce time-space switching of carbon sources and sinks, undergo modulation through atmospheric boundary layer feedbacks, and respond to any discontinuous intervention on the landscape itself such as produced by human intervention in converting rainforest to pasture or conducting selective/clearcut logging operations. The results demonstrate how relative carbon sequestration capacity of the Amazonian ecosystem responds to these controls, and how interpretation of space-time heterogeneities in carbon sequestration depends on a fairly exact quantification of the interacting non-linear properties of photosynthesis in response to incoming solar flux, air-canopy temperatures, and leaf water interception -- and soil respiration in response to upper layer soil temperature and water content. The results also show how the interpretation of the control processes is highly sensitive to the scales at which the surface fluxes are analyzed.

Smith, Eric A.↗
DOE OSTI2022

Subcell limiting strategies for discontinuous Galerkin spectral element methods

Here, we present a general family of subcell limiting strategies to construct robust high-order accurate nodal discontinuous Galerkin (DG) schemes. The main strategy is to construct compatible low order finite volume (FV) type discretizations that allow for convex blending with the high-order variant with the goal of guaranteeing additional properties, such as bounds on physical quantities and/or guaranteed entropy dissipation. For an implementation of this main strategy, four main ingredients are identified that may be combined in a flexible manner: (i) a nodal high-order DG method on Legendre–Gauss–Lobatto nodes, (ii) a compatible robust subcell FV scheme, (iii) a convex combination strategy for the two schemes, which can be element-wise or subcell-wise, and (iv) a strategy to compute the convex blending factors, which can be either based on heuristic troubled-cell indicators, or using ideas from flux-corrected transport methods. By carefully designing the metric terms of the subcell FV method, the resulting methods can be used on unstructured curvilinear meshes, are locally conservative, can handle strong shocks efficiently while directly guaranteeing physical bounds on quantities such as density, pressure or entropy. We further show that it is possible to choose the four ingredients to recover existing methods such as a provably entropy dissipative subcell shock-capturing approach or a sparse invariant domain preserving approach. We test the versatility of the presented strategies and mix and match the four ingredients to solve challenging simulation setups, such as the KPP problem (a hyperbolic conservation law with non-convex flux function), turbulent and hypersonic Euler simulations, and MHD problems featuring shocks and turbulence.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1967

Theoretical studies on the relationship between the thermionic work functions of refractory-metal intermetallic compounds and their electronic and crystal structures Final summary report

The Group VIA-Re binary systems are characterized by the presence of several intermediate phases. The u phase is present in all of the systems, but the a-Mn phases occur in only the Mo-Re and W-Re systems. The occurrence of these intermediate phases appears to be governed largely by electronic factors and is predictable on the basis of electron/atom ratio. The thermionic work functions of the intermediate phase materials were experimentally determined. function and the screening of the outermost electrons in the Group VIA metals was found. A correlation between the work As might be expected for purely metallic substances, the work functions were relatively high, in the range 4. 5 to 5. 0 eV, and in each case were bracketed by the work functions of the pure constituents. A discontinuous change in the average work function with composition is proposed, based on a thermodynamical treatment of the electron vapor.

REFRACTORY METAL↗
NASA NTRS1993

New developments in the method of space-time conservation element and solution element: Applications to the Euler and Navier-Stokes equations

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods--i.e., finite difference, finite volume, finite element, and spectral methods. It is conceptually simple and designed to avoid several key limitations to the above traditional methods. An explicit model scheme for solving a simple 1-D unsteady convection-diffusion equation is constructed and used to illuminate major differences between the current method and those mentioned above. Unexpectedly, its amplification factors for the pure convection and pure diffusion cases are identical to those of the Leapfrog and the DuFort-Frankel schemes, respectively. Also, this explicit scheme and its Navier-Stokes extension have the unusual property that their stabilities are limited only by the CFL condition. Moreover, despite the fact that it does not use any flux-limiter or slope-limiter, the Navier-Stokes solver is capable of generating highly accurate shock tube solutions with shock discontinuities being resolved within one mesh interval. An accurate Euler solver also is constructed through another extension. It has many unusual properties, e.g., numerical diffusion at all mesh points can be controlled by a set of local parameters.

Chang, Sin-Chung↗
NASA NTRS1995

The tip of the red giant branch as a distance indicator for resolved galaxies. 2: Computer simulations

Based on both empirical data for the nearby galaxies, and on computer simulations, we show that measuring the position of the tip of the first-ascent red-giant branch provides a means of obtaining the distances to nearby galaxies with a precision and accuracy comparable to using Cepheids and/or RR Lyrae variables. We present an analysis of synthetic I vs (V-I) color magnitude diagrams of Population 2 systems to investigate the use of the observed discontinuity in the I-band luminosity function as a primary distance indicator. In the simulations we quantify the effects (1) signal to noise, (2) crowding, (3) population size, and (4) non-giant-branch-star contamination, on the method adopted for detecting the discontinuity,, measuring its luminosity, and estimating its uncertainity. We discuss sources of systematic error in the context of observable parameters, such as the signal-to-noise ratio and/or surface brightness. The simulations are then scaled to observed color-magnitude diagrams. It is concluded, that from the ground the tip of the red-giant-branch method can be sucessfully used to determine distances accurate to +/- 10% for galaxies out to 3 Mpc (mu approximately 27.5 mag); and from space a factor of four further in distance (mu approximately 30.6 mag) can be reached using HST. This method can be applied whereever a metal-poor population (-2.0 less than Z less than -0.7) of red-giant stars is detected (whose age is in the range 7-17 Gyr), whether that population resides in the halo of a spiral galaxy, the extended outer disk of a dwarf irregular, or in the outer periphery of an elliptical galaxy.

Madore, Barry F.↗
NASA NTRS1989

Structural development of laminar flow control aircraft chordwise wing joint designs

For laminar flow to be achieved, any protuberances on the surface must be small enough to avoid transition to turbulent flow. However, the surface must have joints between the structural components to allow assembly or replacement of damaged parts, although large continuous surfaces can be utilized to minimize the number the number of joints. Aircraft structural joints usually have many countersunk bolts or rivets on the outer surface. To maintain no mismatch on outer surfaces, it is desirable to attach the components from the inner surface. It is also desirable for the panels to be interchangeable, without the need for shims at the joint, to avoid surface discontinuities that could cause turbulence. Fabricating components while pressing their outer surfaces against an accurate mold helps to ensure surface smoothness and continuity at joints. These items were considered in evaluating the advantages and disadvantages of the joint design concepts. After evaluating six design concepts, two of the leading candidates were fabricated and tested using many small test panels. One joint concept was also built and tested using large panels. The small and large test panel deflections for the leading candidate designs at load factors up to +1.5 g's were well within the step and waviness requirements for avoiding transition.The small panels were designed and tested for compression and tension at -65 F, at ambient conditions, and at 160 F. The small panel results for the three-rib and the sliding-joint concepts indicated that they were both acceptable. The three-rib concept, with tapered splice plates, was considered to be the most practical. A modified three-rib joint that combined the best attributes of previous candidates was designed, developed, and tested. This improved joint met all of the structural strength, surface smoothness, and waviness criteria for laminar flow control (LFC). The design eliminated all disadvantages of the initial three-rib concept except for unavoidable eccentricity, which was reduced and reacted satisfactorily by the rib supports. It should also result in a relatively simple low-cost installation, and makes it easy to replace any panels damaged in the field.

Fischler, J. E.↗
NASA NTRS1993

Inelastic deformation of metal matrix composites

A theoretical model capable of predicting the thermomechanical response of continuously reinforced metal matrix composite laminates subjected to multiaxial loading was developed. A micromechanical model is used in conjunction with nonlinear lamination theory to determine inelastic laminae response. Matrix viscoplasticity, residual stresses, and damage to the fiber/matrix interfacial zone are explicitly included in the model. The representative cell of the micromechanical model is considered to be in a state of generalized plane strain, enabling a quasi two-dimensional analysis to be performed. Constant strain finite elements are formulated with elastic-viscoplastic constitutive equations. Interfacial debonding is incorporated into the model through interface elements based on the interfacial debonding theory originally presented by Needleman, and modified by Tvergaard. Nonlinear interfacial constitutive equations relate interfacial tractions to displacement discontinuities at the interface. Theoretical predictions are compared with the results of an experimental program conducted on silicon carbide/titanium (SiC/Ti) unidirectional, (O4), and angle-ply, (+34)(sub s), tubular specimens. Multiaxial loading included increments of axial tension, compression, torque, and internal pressure. Loadings were chosen in an effort to distinguish inelastic deformation due to damage from matrix plasticity and separate time-dependent effects from time-independent effects. Results show that fiber/matrix debonding is nonuniform throughout the composite and is a major factor in the effective response. Also, significant creep behavior occurs at relatively low applied stress levels at room temperature.

Lissenden, C. J.↗
DOE OSTI2022

Fatigue-crack propagation behavior in a high-carbon chromium SUJ2 bearing steel: Role of microstructure

High-carbon chromium martensite steels are commonly selected for bearing components in the power generation, automotive and aerospace industries where fatigue failure is a major concern. Accordingly, it is of importance to elucidate the structure-property relationships governing the fatigue properties of such bearing steels. Here, the role of microstructure in influencing the fatigue-crack propagation behavior of a high-carbon chromium SUJ2 bearing steel is examined with emphasis on three martensitic structures: (i) a fine-grained (~10 μm) structure with intragranular (but no intergranular) carbides, (ii) a coarser-grained (~25 μm) structure with discontinuous grain-boundary carbides primarily near triple junctions, and (iii) a coarse-grained (~31 μm) structure with continuous grain-boundary carbides. Although growth rates were fairly similar above ~10 -8 m/cycle, significant differences in behavior were observed at lower, near-threshold growth rates; specifically at a load ratio of 0.1, the ΔKth fatigue threshold was increased from 3.8 MPa m ½ in the fine-grained structure to 5.9 to 6.3 MPa m ½ in the two coarser-grained structures. The dominant factor governing such near-threshold behavior was found to be crack-tip shielding from roughness-induced crack closure, which was most pronounced in the coarser-grained structures; additionally, shielding from crack deflection and branching was apparent in these structures as the crack propagated through the network of intersecting martensite lamellae within the coarse grains. Compared to high-strength stainless steels and low-strength carbon steels, the fine-grained SUJ2 steel displays an exceptional combination of high strength and fatigue-crack growth resistance and is considered to be a superior structural alloy for bearing applications.

36 MATERIALS SCIENCE↗
NASA NTRS1950

Hydrodynamic Properties of Planing Surfaces and Flying Boats

The study of the hydrodynamic properties of planing bottom of flying boats and seaplane floats is at the present time based exclusively on the curves of towing tests conducted in tanks. In order to provide a rational basis for the test procedure in tanks and practical design data, a theoretical study must be made of the flow at the step and relations derived that show not only qualitatively but quantitatively the inter-relations of the various factors involved. The general solution of the problem of the development of hydrodynamic forces during the motion of the seaplane float or flying boat is very difficult for it is necessary to give a three-dimensional solution, which does not always permit reducing the analysis to the form of workable computation formulas. On the other had, the problem is complicated by the fact that the object of the analysis is concerned with two fluid mediums, namely, air and water, which have a surface of density discontinuity between them. The theoretical and experimental investigations on the hydrodynamics of a ship cannot be completely carried over to the design of floats and flying-boat hulls, because of the difference in the shape of the contour lines of the bodies, and, because of the entirely different flow conditions from the hydrodynamic viewpoint.

Sokolov, N. A.↗
NASA NTRS1977

Geology of the Imbrium Basin Apennine Mountains and relation to the Apollo 15 landing site

Detailed mapping of the Imbrium Basin Apennine Mountains has emphasized several important factors in the geologic evolution of this region. Pre-Imbrian topography was composed of two types of landforms: (1) the nearby Serenitatis, South Imbrium, Aestuum, and Vaporum basins which produced widely spaced regional lows and a sequence of basin deposits; (2) pre-Imbrian craters that reworked existing basin deposits and produced local topographic lows. The sequence of events during the Imbrium impact can be summarized as follows: intensive structural uplift along the Apennine Front occurred prior to the emplacement of ballistic ejecta. Ejecta became more discontinuous toward the backslope and produced radial texture by sedimentary erosional and depositional processes. Post-ejecta-emplacement downfaulting exposed pre-Imbrian deposits along the front and caused severe seismic disturbances that reworked many of the previous units and accentuated pre-basin structures. Tectonic mapping and orbital geophysical data support the idea that little structural modification has occurred in the Apennines since their formation 3.9 b.y. ago.

Spudis, P.↗
NASA NTRS1992

Navier-Stokes simulations of unsteady transonic flow phenomena

Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.

Atwood, C. A.↗
NASA NTRS1991

Observations of plasma dynamics in the coma of P/Halley by the Giotto Ion Mass Spectrometer

Observations in the coma of P/Halley by the Giotto Ion Mass Spectrometer (IMS) are reported. The High Energy Range Spectrometer (HERS) of the IMS obtained measurements of protons and alpha particles from the far upstream region to the near ionopause region and of ions from mass 12 to 32 at distances of about 250,000 to 40,000 km from the nucleus. Plasma parameters from the High Intensity Spectrometer (HIS) of the IMS obtained between 150,000 to 5000 km from the nucleus are also discussed. The distribution functions of water group ions (water group will be used to refer to ions of 16 to 18 m/q, where m is in AMU and q is in unit charges) are observed to be spherically symmetric in velocity space, indicating strong pitch angle scattering. The discontinuity known as the magnetic pile-up boundary (MPB) is apparent only in proton, alpha, and magnetometer data, indicating that it is a tangential discontinuity of solar wind origin. HERS observations show no significant change in the properties of the heavy ions across the MPB. A comparison of the observations to an MHD model is made. The plasma flow directions at all distances greater than 30,000 km from the nucleus are in agreement with MHD calculations. However, despite the agreement in flow direction, within 200,000 km of the nucleus the magnitude of the velocity is lower than predicted by the MHD model and the density is much larger (a factor of 4). Within 30,000 km of the nucleus there are large theoretical differences between the MHD model flow calculations for the plane containing the magnetic field and for the plane perpendicular to the magnetic field. The observations agreed much better with the pattern calculated for the plane perpendicular to the magnetic field. The data obtained by the High Energy Range Spectrometer (HERS) of the IMS that are published herein were provided to the International Halley Watch archive.

Goldstein, B. E.↗
NASA NTRS1993

Compression strength failure mechanisms in unidirectional composite laminates containing a hole

Experiments on graphite-epoxy laminated plates containing unloaded small holes show that these laminates are notch insensitive. That is, the uniaxial strength of these laminates with small holes exceeds the strength predicted by a point stress criterion using the stress concentration factor for the in-plane stress field. Laminates containing large holes exhibit notch sensitive behavior and consequently their strength is reasonably well predicted by the stress concentration effect. This hole size effect is manifested both in tension and in compression. Apparently, some mechanism must cause in-plane stress relief for laminates containing small holes. The purpose of this research was to study the influence of geometric nonlinearity on the micromechanical response of a filamentary composite material in the presence of a strain gradient caused by a discontinuity such as a hole. A mathematical model was developed at the micromechanical level to investigate this geometrically nonlinear effect.

Johnson, Eric R.↗
DOE OSTI2025

Measurement of correlated charge noise in superconducting qubits at an underground facility

The charge environment of superconducting qubits may be studied through the introduction of controlled, quantified amounts of ionizing radiation. We measure space- and time-correlated charge jumps on a four-qubit device, operating 107 meters below the Earth’s surface in a low-radiation, cryogenic facility designed for the characterization of low-threshold particle detectors. The rock overburden of this facility reduces the cosmic ray muon flux by over 99% compared to laboratories at sea level. Combined with 4π coverage of a movable lead shield, this facility enables quantifiable control over the flux of ionizing radiation on the qubit device. Long-time-series charge tomography measurements on these weakly charge-sensitive qubits capture discontinuous jumps in the induced charge on the qubit islands, corresponding to the interaction of ionizing radiation with the qubit substrate. The rate of these charge jumps scales with the flux of ionizing radiation on the qubit package, as characterized by a series of independent measurements on another energy-resolving detector operating simultaneously in the same cryostat with the qubits. Using lead shielding, we achieve a minimum charge jump rate of $0.19^{+0.04}_{-0.03}$ mHz, almost an order of magnitude lower than that measured in surface tests, but a factor of roughly seven higher than expected based on reduction of ambient gammas alone. We operate four qubits for over 22 consecutive hours with zero correlated charge jumps at length scales above three millimeters.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2011

Getting around Antarctica: New High-Resolution Mappings of the Grounded and Freely-Floating Boundaries of the Antarctic Ice Sheet Created for the International Polar Year

Two ice-dynamic transitions of the Antarctic ice sheet - the boundary of grounded ice features and the freely-floating boundary - are mapped at 15-m resolution by participants of the International Polar Year project ASAID using customized software combining Landsat-7 imagery and ICESat/GLAS laser altimetry. The grounded ice boundary is 53 610 km long; 74% abuts to floating ice shelves or outlet glaciers, 19% is adjacent to open or sea-ice covered ocean, and 7% of the boundary ice terminates on land. The freely-floating boundary, called here the hydrostatic line, is the most landward position on ice shelves that expresses the full amplitude of oscillating ocean tides. It extends 27 521 km and is discontinuous. Positional (one-sigma) accuracies of the grounded ice boundary vary an order of magnitude ranging from +/- 52m for the land and open-ocean terminating segments to +/- 502m for the outlet glaciers. The hydrostatic line is less well positioned with errors over 2 km. Elevations along each line are selected from 6 candidate digital elevation models based on their agreement with ICESat elevation values and surface shape inferred from the Landsat imagery. Elevations along the hydrostatic line are converted to ice thicknesses by applying a firn-correction factor and a flotation criterion. BEDMAP-compiled data and other airborne data are compared to the ASAID elevations and ice thicknesses to arrive at quantitative (one-sigma) uncertainties of surface elevations of +/-3.6, +/-9.6, +/-11.4, +/-30 and +/-100m for five ASAID-assigned confidence levels. Over one-half of the surface elevations along the grounded ice boundary and over one-third of the hydrostatic line elevations are ranked in the highest two confidence categories. A comparison between ASAID-calculated ice shelf thicknesses and BEDMAP-compiled data indicate a thin-ice bias of 41.2+/-71.3m for the ASAID ice thicknesses. The relationship between the seaward offset of the hydrostatic line from the grounded ice boundary only weakly matches a prediction based on beam theory. The mapped products along with the customized software to generate them and a variety of intermediate products are available from the National Snow and Ice Data Center.

Bindschadler, R.↗
NASA NTRS2023

Ultrasonic Inspection During Autoclave Cure of Reflowable-Interface Composite Joints

In structural bonds, the interface between adherend and adhesive is nearly two-dimensional making the interface susceptible to minute quantities of contamination, which can cause weak bonds. Regulatory organizations such as the Federal Aviation Administration (FAA) often require redundant load paths in secondary-bonded primary-structures due to uncertainty in bond performance. To address this issue, the NASA Convergent Aeronautics Solutions (CAS): Adhesive Free Bonding of Composites (AERoBOND) project is investigating reformulated aerospace epoxy-matrix resins to enable reflow and diffusion of the resin at the joint interface during a secondary bonding and cure process that can eliminate the material discontinuity at the interface. Implementing in-situ process monitoring enables assessing the bond quality during processing rather than waiting for post-fabrication mechanical testing to be completed. In this paper, an in-situ ultrasonic inspection system is used to monitor the joining of three composite laminates assembled using the AERoBOND technique. For each panel, the amplitude of the wave reflection at the joint was measured throughout the cure cycle. The results indicate the timing and extent of reflow and cure of the epoxy resin at the joint. Factors indicating end of cure were ascertained and a metric was developed to qualitatively predict acceptable mode-II fracture toughness based on the ultrasonic amplitude at the joint line during cure. In addition, since the inspection system scans an area of the joint, high-resolution localized results can be obtained across the joint.

Tyler B Hudson↗
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

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.↗