Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak transitions”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗
NASA NTRS1996

AB Initio Characterization of MgCCH, MgCCH(+), and MgC2, and Pathways to their Formation in the Interstellar Medium

A study of Mg-bearing compounds has been performed in order to determine molecular properties which are critical for planning new astronomical searches and laboratory studies. The primary focus of the work is on MgCCH, MgCCH(+), and the isomers of MgC2. Only MgCCH has been identified in laboratory studies. Additional calculations have been carried out on MgH, MgNC, MgCN, and their cations in an effort to evaluate pathways to the formation of MgCCH and MgCCH(+) in the InterStellar Medium (ISM) or in circumstellar envelopes. Correlated ab initio methods and correlation-consistent basis sets have been employed. Properties including structures, rotational constants, dipole moments, and harmonic frequencies are reported. A transition state between linear MgCC and cyclic MgC2 has been characterized and was found to yield a minimal barrier (approx. 0.5 kcal/mole), indicating easy interconversion to the cyclic form. Direct reactions in the ISM between Mg or Mg(+) and HCCH are precluded by energetic considerations, but a number of ion- molecule or neutral-neutral exchange reactions between CCH and various Mg-containing species offer plausible pathways to MgCCH or MgCCH(+). Weakly bound MgH may react with CCH to form MgCCH, but MgH has not been detected. Both MgNC and MgCN have been observed, but reactions with CCH are slightly endothermic by 1-3 kcal/mole. Although MgH(+), MgNC(+), and MgCN(+) have not been detected, their reactions with CCH to form MgCCH(+) are all exothermic. With only a small barrier separating linear MgCC and cyclic MgC2, the dissociative recombination of MgCCH(+) with an electron is expected to yield cyclic MgC2, and regenerate Mg and CCH. New astronomical searches for MgCCH, MgCCH(+), cyclic MgC2, MgNC(+), and MgCN(+) will provide further insight into organo-magnesium astrochemistry.

Woon, David E.↗
NASA NTRS2020

Propionitrile in the Two Lowest Excited Vibrational States in the Laboratory and on Titan

The rotational spectrum of propionitrile (ethyl cyanide, C2H5CN) in the two lowest excited vibrational states, ν(sub 13) = 1 and ν(sub 21) = 1, was studied in detail. The compiled experimental data set consisted of published line lists, lines measured from published spectra up to 645 GHz, and new measurements at frequencies from 8 to 331 GHz. The two states are subject to considerable mutual interaction but application of a- and b-axis Coriolis coupling model to the measured A-symmetry internal rotation substate transitions was found to be successful. A total of over 1900 transitions were fitted for J values up to 76 and Ka up to 24, resulting in determined vibrational separation between the two states ΔE = 6:2070246(14) cm(exp -1), and Coriolis coefficients |ζ(sup a)(sub 13,21)| = 0.378, |ζ(sup b) (sub 13,21)| = 0.127, in good agreement with quantum chemistry calculations. The new line catalog was used to model ALMA spectra of Titan’s atmosphere from 2016, fitting many weak lines not previously modeled. By including the vibrationally excited lines of C2H5CN and a new partition function, we derive a disk-averaged volume mixing ratio estimate in 2016 of 7.8 ± 0.2 ppb assumed constant above 300 km.

Analysis of rotational spectrum↗
NASA NTRS2006

An Ultrasensitive Hot-Electron Bolometer for Low-Background SMM Applications

We are developing a hot-electron superconducting transition-edge sensor (TES) that is capable of counting THz photons and operates at T = 0.3K. The main driver for this work is moderate resolution spectroscopy (R approx. 1000) on the future space telescopes with cryogenically cooled (approx. 4 K) mirrors. The detectors for these telescopes must be background-limited with a noise equivalent power (NEP) approx. 10(exp -19)-10(exp -20) W/Hz(sup 1/2) over the range v = 0.3-10 THz. Above about 1 THz, the background photon arrival rate is expected to be approx. 10-100/s), and photon counting detectors may be preferable to an integrating type. We fabricated superconducting Ti nanosensors with a volume of approx. 3x10(exp -3) cubic microns on planar substrate and have measured the thermal conductance G to the thermal bath. A very low G = 4x10(exp -14) W/K, measured at 0.3 K, is due to the weak electron-phonon coupling in the material and the thermal isolation provided by superconducting Nb contacts. This low G corresponds to NEP(0.3K) = 3x10(exp -19) W/Hz(sup 1/2). This Hot-Electron Direct Detector (HEDD) is expected to have a sufficient energy resolution for detecting individual photons with v > 0.3 THz at 0.3 K. With the sensor time constant of a few microseconds, the dynamic range is approx. 50 dB.

bolometers↗
NASA NTRS1994

Imaging spectroscopy of the solar CO lines at 4.67 microns

We analyze spatially and temporally resolved spectra of the fundamental vibration-rotation transitions of carbon monoxide (CO) in the solar spectrum at 4.67 micrometers. Our observations imply that, in the quiet Sun, spatial variations in CO intensity are largely dynamical in nature, reinforcing the suggestion that dynamical effects play a key role in the formation of the dark CO cores. Time sequences of resolved spectra exhibit mainly 3 minute power in line-core intensity but mainly a 5 minute period in Doppler shift. The weak 7-6 R68 line shows normal Evershed flow in the penumbra of a sunspot; we find evidence for the onset of inverse Evershed flow in the strong 3-2 R14 line. Spectra at the limb indicate that 3-2 R14 emission extends approximately 360 km beyond the continuum limb.

Uitenbroek, Han↗
NASA NTRS2010

Ground Based Observation of Isotopic Oxygen in the Martian Atmosphere Using Infrared Heterodyne Spectroscopy

Infrared heterodyne spectra of isotopic CO2 in the Martian atmosphere were obtained using the Goddard Heterodyne Instrument for Planetary Wind and Composition, HIPWAC, which was interfaced with the 3-meter telescope at the NASA Infrared Telescope Facility- Spectra were colle cted at a resolution of lambda/delta lambda=10(exp 7). Absorption fea tures of the CO2 isotopologues have been identified from which isotop ic ratios of oxygen have been determined. The isotopic ratios O-17/O -16 and O-18/O-16 in the Martian atmosphere can be related to Martian atmospheric evolution and can be compared to isotopic ratios of oxyg en in the Earth's atmosphere. Isotopic carbon and oxygen are importa nt constraints on any theory for the erosion of the Martian primordia l atmosphere and the interaction between the atmosphere and surface o r subsurface chemical reservoirs. This investigation explored the pr esent abundance of the stable isotopes of oxygen in Mars' atmospheric carbon dioxide by measuring rovibrational line absorption in isotop ic species of CO2 using groundbased infrared heterodyne spectroscopy in the vicinity of the 9.6 micron and 10.6 micron CO2 lasing bands. T he target transitions during this observation were O-18 C-12 O-16 as well as O-178 C-12 O-16 and O-16 C-113 O-16 at higher resolving power of lambda/delta lambda=10(exp 7) and with high signal-to-noise ratio (longer integration time) in order to fully characterize the absorpt ion line profiles. The fully-resolved lineshape of both the strong n ormal-isotope and the weak isotopic CO2 lines were measured simultane ously in a single spectrum.

Smith, R. L.↗
NASA NTRS2018

Performance of an X-ray Microcalorimeter with a 240 μm Absorber and a 50 μm TES Bilayer

Superconducting transition-edge sensor (TES) microcalorimeters are being developed for a variety of potential astrophysics missions, including Athena. The X-ray integral field unit instrument on this mission requires close-packed pixels on a 0.25 mm pitch, and high quantum efficiency between 0.2 and 12 keV. In this work, we describe a new approach with 50 μm square TESs consisting of a Mo/Au bilayer, deposited on silicon nitride membranes to provide a weak thermal conductance to a ~50 mK heat bath. Larger TESs usually have additional normal metal stripes on top of the bilayer to reduce the noise. However, we have found that excellent spectral performance can be achieved without the need for any normal metal stripes on top of the TES. A spectral performance of 1.58+/-0.12 eV at 5.9 keV has been achieved, the best resolution seen in any of our devices with this pixel size.

Miniussi, Antoine R.↗
NASA NTRS1982

Skylab observations of H I Lyman-alpha

The H I Lyman-alpha profiles obtained with the Skylab SO82-B spectrograph at spectral and spatial resolutions of 60 mA and 2 arcsec x 60 arcsec are presented. Profile variations are illustrated from center to limb, at different heights above the limb, and over a cell and network boundary. The Lyman-alpha core exhibits weak or no limb darkening, whereas the line wings show weak limb darkening. Above the limb, the Lyman-alpha core integrated intensity and line width remain in essence constant up to plus 12 arcsec and may result from resonance scattering of chromospheric radiation by coronal neutral hydrogen or from spicules. On the disk, intensities in the line core are seen as correlating well with the H-alpha network; however, the Lyman-alpha wing intensities are found to correlate better with the transition region lines. The average of profiles acquired over the cell and network region is found to be similar to profiles obtained in the cell.

Roussel-Dupre, D.↗
NASA NTRS2022

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

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

Boundary Layer Stability↗
NASA NTRS1994

The HST quasar absorption line key project. 4: HST faint-object spectrograph and ground-based observations of the unusual low-redshift broad absorption-line quasi-stellar object PG 0043+039

Hubble Space Telescope Faint Object Spectrograph (HST FOS) observations have shown that the spectrum of the low-redshift (z(sub em) approximately equal to 0.384) QSO PG 0043+039 exhibits weak broad absorption lines (BALs). The BALs were discovered during the course of UV spectrophotometry made for the HST Quasar Absorption Line Key Project. The HST data are analyzed along with ground-based optical and IUE spectrophotometry. The object is found to have a number of atypical properties relative to normal non-BAL QSOs. The observed continuum is atypical in the sense that it is much weaker than that of a normal optically selected QSO at rest wavelengths approximately less than 2200 A. Intrinsic reddening of E(B-V) approximately equal to 0.11 mag by dust similar to that found in the SMC at the redshift of PG 0043+039 conservatively accounts for the observed continuum shape moderately well. These observed characteristics are typical of low-ionization BAL QSOs, but convincing evidence for BALs due to low-ionization transitions of Mg II, Al III, Al II, or C II does not exist. Therefore, this object may be a misaligned BAL QSO having many of the characteristics of low-ionization BAL QSOs with the sight line passing through a putative dusty region, but evidently missing clouds of high enough column density to produce observable low-ionization BALs. If the intrinsic dust-extinction model is correct, the observations suggest that the dust is not confined to the presumably higher density, low-ionization BAL clouds, but that it has drifted to nearby high-ionization BAL regions. We also consider other possible mechanisms for producing the shape of the continuous energy distribution which cannot be ruled out. We compare the Fe II emission in PG 0043+039 with that in another Key Project QSO, NGC 2841-UB 3, which has optical Fe II emission comparable in strength to that in PG 0043+039, but has anomalously weak UV Fe II emission. In addition, from an analysis of UV and optical spectrophotometric data at 5 epochs over approximately 11 yr, there is tentative evidence that PG 0043+039 has varied in brightness by as much as 1.1 mag during this time interval. Two different interpretations of PG 0043+039 and the low-ionization BAL QSOs are considered. Various model scenarios for explaining the weak narrow-line (O III) emission are considered, but there is no definitive explanation.

Turnshek, David A.↗
NASA NTRS2016

Exhibition of Stochastic Resonance in Vestibular Perception

Astronauts experience sensorimotor changes during spaceflight, particularly during G-transitions. Post flight sensorimotor changes include spatial disorientation, along with postural and gait instability that may degrade operational capabilities of the astronauts and endanger the crew. A sensorimotor countermeasure that mitigates these effects would improve crewmember safety and decrease risk. The goal of this research is to investigate the potential use of stochastic vestibular stimulation (SVS) as a technology to improve sensorimotor function. We hypothesize that low levels of SVS will improve sensorimotor perception through the phenomenon of stochastic resonance (SR), when the response of a nonlinear system to a weak input signal is enhanced by the application of a particular nonzero level of noise. This study aims to advance the development of SVS as a potential countermeasure by 1) demonstrating the exhibition of stochastic resonance in vestibular perception, a vital component of sensorimotor function, 2) investigating the repeatability of SR exhibition, and 3) determining the relative contribution of the semicircular canals (SCC) and otolith (OTO) organs to vestibular perceptual SR. A constant current stimulator was used to deliver bilateral bipolar SVS via electrodes placed on each of the mastoid processes, as previously done. Vestibular perceptual motion recognition thresholds were measured using a 6-degree of freedom MOOG platform and a 150 trial 3-down/1-up staircase procedure. In the first test session, we measured vestibular perceptual thresholds in upright roll-tilt at 0.2 Hz (SCC+OTO) with SVS ranging from 0-700 μA. In a second test session a week later, we re-measured roll-tilt thresholds with 0, optimal (from test session 1), and 1500 μA SVS levels. A subset of these subjects, plus naive subjects, participated in two additional test sessions in which we measured thresholds in supine roll-rotation at 0.2 Hz (SCC) and upright y-translation at 1 Hz (OTO) with SVS up to 700 μA. A sinusoidal galvanic vestibular stimulation (GVS) perceptual threshold was also measured on each test day and used to normalize the SVS levels across subjects. In roll-tilt thresholds with SVS, the characteristic SR curve was qualitatively exhibited in 10 of 12 subjects, and the improvement in motion threshold was significant in 6 subjects, indicating that optimal SVS improved passive body motion perception in a way that is consistent with classical SR theory. A probabilistic comparison to numeric simulations further validated these experimental results. On the second test session, 4 out of the 10 SR exhibitors showed repeated improvement with SVS compared to the no SVS condition. Data collection is ongoing for the last two test sessions in which SCC and OTO only perceptual motion recognition thresholds are being measured with SVS. The final results of these test sessions will give insight into whether vestibular perceptual SR can occur when only one type of vestibular sensor is sensing motion or if it is more evident when sensory integration between the SCC and OTO is occurring during the motion. The overall purpose of this research is to further quantify the effects of SVS on various sensorimotor tasks and to gain a more fundamental understanding of how SVS causes SR in the vestibular system. In the context of human space flight, results from this research will help in understanding how SVS may be practically implemented in the future as a component of a comprehensive countermeasure plan for G-transition adaptation.

Galvan-Garza, R. C.↗
NASA NTRS2010

High C/O Chemistry and Weak Thermal Inversion in the Extremely Irradiated Atmosphere of Exoplanet WASP-12b

The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior as opposed to the silicate-dominated composition as found on Earth; the solar C/O is 0.54. Theory, shows that high C/O leads to a diversity of carbon-rich planets that can have very different interiors and atmospheres from those in the solar system. Here we report the detection of C/O greater than or equal to 1 in a planetary atmosphere. The transiting hot Jupiter WASP-12b has a dayside atmosphere depleted in water vapour and enhanced in methane by over two orders of magnitude compared to a solar-abundance chemical equilibrium model at the expected temperatures. The observed concentrations of the prominent molecules CO, CH4, and H2O are consistent with theoretical expectations for an atmosphere with the observed C/O = 1. The C/O ratios are not known for giant planets in the solar system, although they are expected to equal the solar value. If high C/O ratios are common, then extrasolar planets are likely very different in interior composition, and formed very differently, from expectations based on solar composition, potentially explaining the large diversity in observed radii. We also find that the extremely irradiated atmosphere (greater than 2500 K) of WASP-12b lacks a prominent thermal inversion, or a stratosphere, and has very efficient day-night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres.

Madhusudhan, Nikku↗
NASA NTRS2022

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

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

Boundary Layer Stability↗
NASA NTRS2001

Quantitative Studies on the Propagation and Extinction of Near-Limit Premixed Flames under Normal- and Micro-gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion.. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Studies at high strain rates are important in quantifying and understanding the response of vigorously burning flames and determine extinction states. Studies of weakly strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds, S(sup o)(sub u), can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames are stabilized by weak strain rates. Previous studies have shown that near-limit flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In our past studies the laminar flame speeds and extinction strain rates were directly measured at ultra-low strain rates. The laminar flame speeds were measured by having a positively strained planar flame undergoing a transition to a negatively strained Bunsen flame and by measuring the propagation speed during that transition. The extinction strain rates of near-limit flames were measured in mu-g. Results obtained for CH4/air and C3H8/air mixtures are in agreement with those obtained by Maruta et al.

Egolfopoulos, F. N.↗
NASA NTRS2016

Aerospace Applications of Non-Equilibrium Plasma

Nonequilibrium plasma/non-thermal plasma/cold plasmas are being used in a wide range of new applications in aeronautics, active flow control, heat transfer reduction, plasma-assisted ignition and combustion, noise suppression, and power generation. Industrial applications may be found in pollution control, materials surface treatment, and water purification. In order for these plasma processes to become practical, efficient means of ionization are necessary. A primary challenge for these applications is to create a desired non-equilibrium plasma in air by preventing the discharge from transitioning into an arc. Of particular interest is the impact on simulations and experimental data with and without detailed consideration of non-equilibrium effects, and the consequences of neglecting non-equilibrium. This presentation will provide an assessment of the presence and influence of non-equilibrium phenomena for various aerospace needs and applications. Specific examples to be considered will include the forward energy deposition of laser-induced non-equilibrium plasmoids for sonic boom mitigation, weakly ionized flows obtained from pulsed nanosecond discharges for an annular Hall type MHD generator duct for turbojet energy bypass, and fundamental mechanisms affecting the design and operation of novel plasma-assisted reactive systems in dielectric liquids (water purification, in-pipe modification of fuels, etc.).

Nonequilibrium Plasma↗
NASA NTRS1998

Shear-Sensitive Liquid Crystal Coating Method: Surface-Inclination Effects on Shear Vector Measurements

The shear-sensitive liquid crystal coating (SSLCC) method is an image-based technique for both visualizing dynamic surface-flow phenomena, such as transition and separation, and for measuring the continuous shear-stress vector distribution acting on an aerodynamic surface. Under proper lighting and viewing conditions (discussed below), the coating changes color in response to an applied aerodynamic shear. This color-change response is continuous and reversible, with a response time of milliseconds, and is a function of both the shear magnitude and the shear vector orientation relative to the observer. The liquid crystal phase of matter is a weakly-ordered, viscous, non-Newtonian fluid state that exists between the nonuniform liquid phase and the ordered solid phase of certain organic compounds. Cholesteric liquid crystal compounds possess a helical molecular arrangement that selectively scatters white light, incident along the helical axis, as a three-dimensional spectrum. This property is linked to the helical pitch length, which is within the range of wavelengths in the visible spectrum. The pitch length, and hence the wavelength of the scattered light, is influenced by shear stress normal to the helical axis. This unique optical property produces a measurable color change in response to an applied shearing force. The full-surface shear stress vector measurement method, developed at NASA-Ames, is schematically illustrated. As with the visualization method, the coated test surface is illuminated from the normal direction with white light and the camera is positioned at an above-plane view angle of approximately 30 deg. Experiments have been initiated at NASA Ames to begin the process of quantifying surface-inclination (surface-curvature) effects on shear vector measurement accuracy. In preliminary experiments, surface-inclination angles theta(sub x), theta(sub y) of 0, +/-5, +/-10, and +/-15 deg were employed. In this arrangement, white-light illumination was positioned normal to the untilted test surface, and the camera above-plane view angle was set at 30 deg relative to the untilted test surface. As can be seen, vector-aligned lambda(sub d) values showed no dependence on theta(sub x) or theta(sub y) for absolute values of these tilt angles is less than or equal to 15 deg. Acquisition and analyses of full-surface color images are presently underway to definitively document the insensitivity limits of the shear vector measurement methodology to surface-slope variations.

Reda, Daniel C.↗
NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗
NASA NTRS1994

The peculiar type II supernova 1993J in M81: Transition to the nebular phase

We present optical spectra of the bright, peculiar Type II supernova 1993J in M81 spanning the first 14 months of its existence, revealing its transition to the nebular phase. Unlike the case in normal Type II supernovae, during the first 2-10 months the H-alpha emission line gradually becomes less prominent relative to other features such as (O I) lambda lambda 6300, 6364 and (Ca II) lambda lambda 7291, 7324, as we had predicted based on early-time (tau less than or approximately equal to 2 months) spectra. The nebular spectrum resembles those of the Type Ib/Ic supernovae 1985F and 1987M, although weak H-alpha emission is easily visible even at late times in SN 1993J. At tau = 8 months a close similarity is found with the spectrum of SN 1987K, the only other Type II supernova known to have undergone such a metamorphosis. The emission lines are considerably broader than those of normal Type II supernovae at comparable phases, consistent with the progenitor having lost a majority of its hydrogen envelope prior to exploding. Consequently, there is now little doubt that Type Ib, and probably Type Ic, supernovae result from core collapse in stripped, massive stars; models of the chemical evolution of galaxies in which these subtypes are ascribed to exploding white dwarfs must be appropriately modified. Although all of the emission lines in spectra of SN 1993J fade roughly exponentially for a considerable time, the fading of H-alpha begins to slow down at tau approximately = 8 months, and in the interval tau = 10-14 months its flux is constant, or even slightly rising in the wings of the line. This behavior, together with the box-like shape and great breadth (full width at half maximum (FWHM) approximately = 17 000 km/s) of the line profile, suggests that the H-alpha emission is being produced by the high-velocity outer layer of hydrogen ejecta interacting with circumstellar gas released by the progenitor prior to its explosion. A similar phenomenon has previously been seen at later phases in several Type II supernovae, most notably SN 1980K. Bumps (FWHM approximately = 1000 km/s, amplitude approximately = 20%) in the H-alpha profile are probably indicative of Rayleigh-Taylor instabilities in the cool gas behind the reverse shock. A very narrow component (unresolved, FWHM less than or approximately equal to 200 km/s) of H-alpha at the symmetric velocity of SN 1993J may instead be produced by a superposed H II region, or perhaps by recombination in a large circumstellar shell or ring that was ionized during the first few hours after outburst. In the near future the spectrum of SN 1993J should become increasingly dominated by broad H-alpha emission.

Filippenko, Alexei V.↗