Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data intensive computing”

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 703 records · Page 39

Tracker: Image-Processing and Object-Tracking System Developed

Tracker is an object-tracking and image-processing program designed and developed at the NASA Lewis Research Center to help with the analysis of images generated by microgravity combustion and fluid physics experiments. Experiments are often recorded on film or videotape for analysis later. Tracker automates the process of examining each frame of the recorded experiment, performing image-processing operations to bring out the desired detail, and recording the positions of the objects of interest. It can load sequences of images from disk files or acquire images (via a frame grabber) from film transports, videotape, laser disks, or a live camera. Tracker controls the image source to automatically advance to the next frame. It can employ a large array of image-processing operations to enhance the detail of the acquired images and can analyze an arbitrarily large number of objects simultaneously. Several different tracking algorithms are available, including conventional threshold and correlation-based techniques, and more esoteric procedures such as "snake" tracking and automated recognition of character data in the image. The Tracker software was written to be operated by researchers, thus every attempt was made to make the software as user friendly and self-explanatory as possible. Tracker is used by most of the microgravity combustion and fluid physics experiments performed by Lewis, and by visiting researchers. This includes experiments performed on the space shuttles, Mir, sounding rockets, zero-g research airplanes, drop towers, and ground-based laboratories. This software automates the analysis of the flame or liquid s physical parameters such as position, velocity, acceleration, size, shape, intensity characteristics, color, and centroid, as well as a number of other measurements. It can perform these operations on multiple objects simultaneously. Another key feature of Tracker is that it performs optical character recognition (OCR). This feature is useful in extracting numerical instrumentation data that are embedded in images. All the results are saved in files for further data reduction and graphing. There are currently three Tracking Systems (workstations) operating near the laboratories and offices of Lewis Microgravity Science Division researchers. These systems are used independently by students, scientists, and university-based principal investigators. The researchers bring their tapes or films to the workstation and perform the tracking analysis. The resultant data files generated by the tracking process can then be analyzed on the spot, although most of the time researchers prefer to transfer them via the network to their offices for further analysis or plotting. In addition, many researchers have installed Tracker on computers in their office for desktop analysis of digital image sequences, which can be digitized by the Tracking System or some other means. Tracker has not only provided a capability to efficiently and automatically analyze large volumes of data, saving many hours of tedious work, but has also provided new capabilities to extract valuable information and phenomena that was heretofore undetected and unexploited.

Klimek, Robert B.↗

Normalizing a Relativistic Model of X-Ray Reflection Definition of the Reflection Fraction and Its Implementation in relxill

Aims. The only relativistic reflection model that implements a parameter relating the intensity incident on an accretion disk to the observed intensity is relxill. The parameter used in earlier versions of this model, referred to as the reflection strength, is unsatisfactory; it has been superseded by a parameter that provides insight into the accretion geometry, namely the reflection fraction. The reflection fraction is defined as the ratio of the coronal intensity illuminating the disk to the coronal intensity that reaches the observer. Methods. The relxill model combines a general relativistic ray-tracing code and a photoionization code to compute the component of radiation reflected from an accretion that is illuminated by an external source. The reflection fraction is a particularly important parameter for relativistic models with well-defined geometry, such as the lamp post model, which is a focus of this paper. Results. Relativistic spectra are compared for three inclinations and for four values of the key parameter of the lamp post model,namely the height above the black hole of the illuminating, on-axis point source. In all cases, the strongest reflection is produced for low source heights and high spin. A low-spin black hole is shown to be incapable of producing enhanced relativistic reflection. Results for the relxill model are compared to those obtained with other models and a Monte Carlo simulation. Conclusions. Fitting data by using the relxill model and the recently implemented reflection fraction, the geometry of a system can be constrained. The reflection-fraction is independent of system parameters such as inclination and black hole spin. The reflection-fraction parameter was implemented with the name reflec_frac all flavours of the relxill model, and the non-relativistic reflection model xillver, in v0.4a (18 January 2016).

accretion↗

A stochastic precipitation disaggregation scheme for GCM applications

In the surface hydrologic parameterization of general circulation models (GCMs), it is commonly assumed that the precipitation processes are homogeneous over a GCM grid square and that the precipitation intensity is uniformly distributed. Based on evidence that the spatial distribution of precipitation within a GCM grid square is crucial for the land surface hydrology parameterization, a few researchers have explored the impacts of assuming that the precipitation is exponentially distributed. This paper explores the suitability of the aforementioned assumptions. First, a statistical analysis is conducted of historical precipitation data for three GCM grids in different regions of the United States. The analysis suggests that neither the uniform nor the exponential distribution assumption may be suitable at the GCM grid scale and, that instead, the spatial variability in precipitation is characterized by statistical patterns that are inhomogeneous. These patterns vary from grid to grid and are induced by the interaction between atmospheric conditions and various land surface characteristics, such as topographical features, surface properties, etc. Within the same grid square, however, the statistical patterns are generally constant from year to year. Based on this analysis, a computationally viable (i.e., usable with GCMs) stochastic precipitation disaggregation scheme that utilizes these stable statistical patterns is proposed. The method was used to generate spatially distributed hourly rainfall for a summer season in the southwestern region of the continental United States. Analysis of the results shows that the methodology preserves the seasonal characteristics of spatial variability in precipitation that is observed in the long-term historical data.

Gao, Xiaogang↗

Convective flow effects on protein crystal growth

A high-resolution microscopic interferometric setup for the monitoring of protein morphologies has been developed. Growth or dissolution of a crystal can be resolved with a long-term depth resolution of 200 A and a lateral resolution of 2 microns. This capability of simultaneously monitoring the interfacial displacement with high local depth resolution has yielded several novel results. We have found with lysozyme that (1) the normal growth rate is oscillatory, and (2) depending on the impurity content of the solution, the growth step density is either greater or lower at the periphery of a facet than in its center. The repartitioning of Na plus and Cl minus ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed, to interpret the large body of data in unified way. The results strongly suggest that (1) the ion to lysozyne ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter electron microscopy results appear to confirm this finding, which could have far-reaching consequences for x-ray diffraction studies. A computational model for diffusive-convective transport in protein crystallization has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies and our kinetics data for the growth of lysozyme. The results show that even in the small cell employed, protein concentration nonuniformities and gravity-driven solutal convection can be significant. The calculated convection velocities are of the same order to magnitude as those found in earlier experiments. As expected, convective transport, i.e., at Og, lysozyme crystal growth remains kinetically limited. The salt distribution in the crystal is predicted to be non-uniform at both 1g and 0g, as a consequence of protein depletion in the solution. Static and dynamic light scattering studies in undersaturated and supersaturated solutions have been performed. Diffusivities in undersaturated solutions, were found to vary with lysozyme concentrations. Depending on the salt concentration, the diffusivities either increase or decrease. Interestingly, the corresponding static scattering intensities behave oppositely, Our current analysis indicates that these changes are inconsistent with aggregation in undersaturated solutions. However, the data are compatible with concentration-dependent changes of the interactions between protein and salt.

Rosenberger, Franz↗

A Statistical Study of the Global Structure of the Ring Current

In this paper we derive the average configuration of the ring current as a function of the state of the magnetosphere as indicated by the Dst index. We sort magnetic field data from the Combined Release and Radiation Effects Satellite (CRRES) by spatial location and by the Dst index in order to produce magnetic field maps. From these maps we calculate local current systems by taking the curl of the magnetic field. NN7e find both the westward (outer) and the eastward (inner) components of the ring current. We find that the ring current intensity varies linearly with D.St as expected, and that the ring current is asymmetric for all Dst values. The azimuthal peak of the ring current is located in the afternoon sector for quiet conditions, and near midnight for disturbed conditions. The ring current also moves closer to the Earth during disturbed conditions. We are able to recreate the Dst index by integrating the magnetic perturbations caused by the ring current. We find that we needed to apply a 20 nT offset to Dst, and assume a perfectly conducting Earth to obtain an optimal agreement between the computed and the observed Dst. We interpret the 20 nT offset as the magnetic field generated by the quiet time ring current used as baseline in computing Dst.

Jorgensen, A. M.↗

Particulate Emission and Optical Measurements in n-Heptane Low-Swirl Non-Premixed Flames at Elevated Pressures

A series of ground combustion tests conducted at elevated pressures demonstrates in-chamber optical measurements of soot emissions, combined with conventional gaseous and particle emissions sampling techniques. Utilizing a micro-radial-entry counter-swirl (MRX) burner, we successfully stabilized an n-heptane low-swirl, non-premixed spray flame across a wide range of equivalence ratios at pressures reaching 250 psia. Our optical diagnostics feature a temporally gated diffused-backlight illumination extinction imaging (DBI-EI) system, equipped with a 470-nm strobe high-intensity LED array as the illumination source, along with a line-of-sight telecentric camera. By accessing the flame tube, the pulsed DBI-EI provides illuminating optical extinction images, enabling experimental estimation of soot volume fractions within the reaction zone. This imaging technique will allow for real-time observation of soot formation dynamics occurring above the air-blast fuel injector. To further enhance our analysis, we employed simultaneous gaseous and particulate extractive diagnostics, interpreting in-flame particulate information through CO2 sampling, exhaust particle sizing, and particle counting. The processed dataset reveals a promising application of DBI-EI for propulsion emission analysis in high-pressure environments, where optical beam-steering effects can significantly impact results. This report is essential to bridging the divide between in-flame optical diagnostics and exhaust extractive measurements. By achieving an understanding of non-premixed combustion emission characteristics, we attempt to provide invaluable data for computer code validation, advancing the field of combustion research and supporting the development of cleaner, more efficient engines.

high-pressure flame↗

Particulate Emission and Optical Measurements in n-Heptane Low-Swirl Non-Premixed Flames at Elevated Pressures

A series of ground combustion tests conducted at elevated pressures demonstrates in-chamber optical measurements of soot emissions, combined with conventional gaseous and particle emissions sampling techniques. Utilizing a micro-radial-entry counter-swirl (MRX) burner, we successfully stabilized an n-heptane low-swirl, non-premixed spray flame across a wide range of equivalence ratios at pressures reaching 250 psia. Our optical diagnostics feature a temporally gated diffused-backlight illumination extinction imaging (DBI-EI) system, equipped with a 470-nm strobe high-intensity LED array as the illumination source, along with a line-of-sight telecentric camera. By accessing the flame tube, the pulsed DBI-EI provides illuminating optical extinction images, enabling experimental estimation of soot volume fractions within the reaction zone. This imaging technique will allow for real-time observation of soot formation dynamics occurring above the air-blast fuel injector. To further enhance our analysis, we employed simultaneous gaseous and particulate extractive diagnostics, interpreting in-flame particulate information through CO2 sampling, exhaust particle sizing, and particle counting. The processed dataset reveals a promising application of DBI-EI for propulsion emission analysis in high-pressure environments, where optical beam-steering effects can significantly impact results. This report is essential to bridging the divide between in-flame optical diagnostics and exhaust extractive measurements. By achieving an understanding of non-premixed combustion emission characteristics, we attempt to provide invaluable data for computer code validation, advancing the field of combustion research and supporting the development of cleaner, more efficient engines.

high-pressure flame↗

Time-Averaged Velocity, Temperature and Density Surveys of Supersonic Free Jets

A spectrally resolved molecular Rayleigh scattering technique was used to simultaneously measure axial component of velocity U, static temperature T, and density p in unheated free jets at Mach numbers M = 0.6,0.95, 1.4 and 1.8. The latter two conditions were achieved using contoured convergent-divergent nozzles. A narrow line-width continuous wave laser was passed through the jet plumes and molecular scattered light from a small region on the beam was collected and analyzed using a Fabry-Perot interferometer. The optical spectrum analysis air density at the probe volume was determined by monitoring the intensity variation of the scattered light using photo-multiplier tubes. The Fabry-Perot interferometer was operated in the imaging mode, whereby the fringe formed at the image plane was captured by a cooled CCD camera. Special attention was given to remove dust particles from the plume and to provide adequate vibration isolation to the optical components. The velocity profiles from various operating conditions were compared with that measured by a Pitot tube. An excellent comparison within 5m's demonstrated the maturity of the technique. Temperature was measured least accurately, within 10K, while density was measured within 1% uncertainty. The survey data consisted of centerline variations and radial profiles of time-averaged U, T and p. The static temperature and density values were used to determine static pressure variations inside the jet. The data provided a comparative study of jet growth rates with increasing Mach number. The current work is part of a data-base development project for Computational Fluid Dynamics and Aeroacoustics codes that endeavor to predict noise characteristics of high speed jets. A limited amount of far field noise spectra from the same jets are also presented. Finally, a direct experimental validation was obtained for the Crocco-Busemann equation which is commonly used to predict temperature and density profiles from known velocity profiles. Data presented in this paper are available in ASCII format upon request.

Panda, Jayanta↗

Mercury's helium exosphere after Mariner 10's third encounter

From a comparison of the Mariner 10 third encounter UV spectrometer data with intensities generated from a newly constructed model exosphere, a new value of 4.5 x 10 to the -4th power for the fraction of the solar wind He(++) flux to be intercepted and captured by Mercury's magnetosphere, if the observed He atmosphere is maintained by the solar wind, has been derived. If an internal source for He prevails, the corresponding upper bound for the global outgassing rate is estimated to be 4.5 x 10 to the 22nd power/s. These values differ from those given earlier because of the present use of a surface temperature distribution satisfying the heat equation over Mercury's entire surface which employs Mariner-10-determined mean surface thermal characteristics. The mean standoff distance of Mercury's magnetopause averaged over Mercury's orbit is also used. Agreement between the observed and calculated intensities is found to be good. Minor discrepancies on the nightside of the terminator are explicable in terms of differences between actual and computed temperatures and scale height structure changes.

Curtis, S. A.↗

Heat transfer measurements from a NACA 0012 airfoil in flight and in the NASA Lewis icing research tunnel

Local heat transfer coefficients from a smooth and roughened NACA 0012 airfoil were measured using a steady state heat flux method. Heat transfer measurements on the specially constructed 0.533 meter chord airfoil were made both in flight on the NASA Lewis Twin Otter Research Aircraft and in the NASA Lewis Icing Research Tunnel (IRT). Roughness was obtained by the attachment of small, 2 mm diameter, hemispheres of uniform size to the airfoil surface in four distinct patterns. The flight data was taken for the smooth and roughened airfoil at various Reynolds numbers based on chord in the range of 1.24x10(exp 6) to 2.50x10(exp 6) and at various angles of attack up to 4 degrees. During these flight tests the free stream velocity turbulence intensity was found to be very low (less than 0.1 percent). The wind tunnel data was taken in the Reynolds number range of 1.20x10(exp 6) to 4.52x10(exp 6) and at angles of attack from -4 degrees to +8 degrees. The turbulence intensity in the IRT was 0.5 to 0.7 percent with the cloud making spray off. Results for both the flight and tunnel tests are presented as Frossling number based on chord versus position on the airfoil surface for various roughnesses and angle of attack. A table of power law curve fits of Nusselt number as a function of Reynolds number is also provided. The higher level of turbulence in the IRT versus flight had little effect on heat transfer for the lower Reynolds numbers but caused a moderate increase in heat transfer at the higher Reynolds numbers. Turning on the cloud making spray air in the IRT did not alter the heat transfer. Roughness generally increased the heat transfer by locally disturbing the boundary layer flow. Finally, the present data was not only compared with previous airfoil data where applicable, but also with leading edge cylinder and flat plate heat transfer values which are often used to estimate airfoil heat transfer in computer codes.

Poinsatte, Philip E.↗

Application of engineering codes to the assessment of spacecraft charging induced hazards

One of the major hurdles to developing a better understanding of the environmentally induced interactions in spacecraft systems (spacecraft charging) is the lack of simple tools to analyze the behavior of the spacecraft in a sub-storm environment. The available tool, NASA Charging Analyzer Program (NASCAP), is expensive to use and treats only half of the problem - surface effects only. In response to this need, we have developed, under our IRAD Program, a set of engineering tools that operate on personal computers and addresses both the surface and bulk charging aspects of the interactions. These engineering codes utilize circuit analysis techniques and model simplifications to allow rapid computations. The codes were compared to both NASCAP predictions and to bulk charging data with excellent agreement for both phenomena. The models used in our engineering codes are shown, and in the engineering codes, the models are defined in terms of plates or triangles in an unconfined space. The predictions are based upon current balances to each plate. The code has routines to allow interactions between the plates to simulate 'barrier' effects. The agreement for the predicted differential surface voltages for the model is given. In conducting the comparison of engineering codes to NASCAP predictions, it was found that some inconsistencies existed in the NASCAP code that are not obvious to the user. The use of the grid imposes constraints that voltage gradients exist between the grid points. This allows the existence of voltage gradients in conductors and influences surface resistivity predictions. These points are discussed along with the predictions for treating these inconsistencies properly. The model predictions for a variable intensity substorm is also presented.

Stevens, N. J.↗

Defect Localization Capabilities of a Global Detection Scheme: Spatial Pattern Recognition Using Full-field Vibration Test Data in Plates

Recently, a conceptually simple approach, based on the notion of defect energy in material space has been developed and extensively studied (from the theoretical and computational standpoints). The present study focuses on its evaluation from the viewpoint of damage localization capabilities in case of two-dimensional plates; i.e., spatial pattern recognition on surfaces. To this end, two different experimental modal test results are utilized; i.e., (1) conventional modal testing using (white noise) excitation and accelerometer-type sensors and (2) pattern recognition using Electronic speckle pattern interferometry (ESPI), a full field method capable of analyzing the mechanical vibration of complex structures. Unlike the conventional modal testing technique (using contacting accelerometers), these emerging ESPI technologies operate in a non-contacting mode, can be used even under hazardous conditions with minimal or no presence of noise and can simultaneously provide measurements for both translations and rotations. Results obtained have clearly demonstrated the robustness and versatility of the global NDE scheme developed. The vectorial character of the indices used, which enabled the extraction of distinct patterns for localizing damages proved very useful. In the context of the targeted pattern recognition paradigm, two algorithms were developed for the interrogation of test measurements; i.e., intensity contour maps for the damaged index, and the associated defect energy vector field plots.

Saleeb, A. F.↗

Level 1 Processing of MODIS Direct Broadcast Data From Terra

In February 2000, an effort was begun to adapt the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 1 production software to process direct broadcast data. Three Level 1 algorithms have been adapted and packaged for release: Level 1A converts raw (level 0) data into Hierarchical Data Format (HDF), unpacking packets into scans; Geolocation computes geographic information for the data points in the Level 1A; and the Level 1B computes geolocated, calibrated radiances from the Level 1A and Geolocation products. One useful aspect of adapting the production software is the ability to incorporate enhancements contributed by the MODIS Science Team. We have therefore tried to limit changes to the software. However, in order to process the data immediately on receipt, we have taken advantage of a branch in the geolocation software that reads orbit and altitude information from the packets themselves, rather than external ancillary files used in standard production. We have also verified that the algorithms can be run with smaller time increments (2.5 minutes) than the five-minute increments used in production. To make the code easier to build and run, we have simplified directories and build scripts. Also, dependencies on a commercial numerics library have been replaced by public domain software. A version of the adapted code has been released for Silicon Graphics machines running lrix. Perhaps owing to its origin in production, the software is rather CPU-intensive. Consequently, a port to Linux is underway, followed by a version to run on PC clusters, with an eventual goal of running in near-real-time (i.e., process a ten-minute pass in ten minutes).

Lynnes, Christopher↗

Cloud Properties Derived From GOES-7 for Spring 1994 ARM Intensive Observing Period Using Version 1.0.0 of ARM Satellite Data Analysis Program

This document describes the initial formulation (Version 1.0.0) of the Atmospheric Radiation Measurement (ARM) program satellite data analysis procedures. Techniques are presented for calibrating geostationary satellite data with Sun synchronous satellite radiances and for converting narrowband radiances to top-of-the-atmosphere fluxes and albedos. A methodology is documented for combining geostationary visible and infrared radiances with surface-based temperature observations to derive cloud amount, optical depth, height, thickness, temperature, and albedo. The analysis is limited to two grids centered over the ARM Southern Great Plains central facility in north-central Oklahoma. Daytime data taken during 5 Apr. - 1 May 1994, were analyzed on the 0.3 deg and 0.5 deg latitude-longitude grids that cover areas of 0.9 deg x 0.9 deg and 10 deg x 14 deg, respectively. Conditions ranging from scattered low cumulus to thin cirrus and thick cumulonimbus occurred during the study period. Detailed comparisons with hourly surface observations indicate that the mean cloudiness is within a few percent of the surface-derived sky cover. Formats of the results are also provided. The data can be accessed through the World Wide Web computer network.

Minnis, Patrick↗

An Investigation of Transonic Resonance in a Mach 2.2 Round Convergent-Divergent Nozzle

Hot-wire and acoustic measurements were taken for a round convergent nozzle and a round convergent-divergent (C-D) nozzle at a jet Mach number of 0.61. The C-D nozzle had a design Mach number of 2.2. Compared to the convergent nozzle jet flow, the Mach 2.2 nozzle jet flow produced excess broadband noise (EBBN). It also produced a transonic resonance tone at 1200 Herz. Computational simulations were performed for both nozzle flows. A steady Reynolds-Averaged Navier-Stokes simulation was performed for the convergent nozzle jet flow. For the Mach 2.2 nozzle flow, a steady RANS simulation, an unsteady RANS (URANS) simulation, and an unsteady Detached Eddy Simulation (DES) were performed. The RANS simulation of the convergent nozzle showed good agreement with the hot-wire velocity and turbulence measurements, though the decay of the potential core was over-predicted. The RANS simulation of the Mach 2.2 nozzle showed poor agreement with the experimental data, and more closely resembled an ideally-expanded jet. The URANS simulation also showed qualitative agreement with the hot-wire data, but predicted a transonic resonance at 1145 Herz. The DES showed good agreement with the hot-wire velocity and turbulence data. The DES also produced a transonic tone at 1135 Herz. The DES solution showed that the destabilization of the shock-induced separation region inside the nozzle produced increased levels of turbulence intensity. This is likely the source of the EBBN.

DES↗

Recommendations for Using Noise Monitors to Estimate Noise Exposure During X-59 Community Tests

A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.

sonic boom↗

The Numerical Analysis of a Turbulent Compressible Jet

A numerical method to simulate high Reynolds number jet flows was formulated and applied to gain a better understanding of the flow physics. Large-eddy simulation was chosen as the most promising approach to model the turbulent structures due to its compromise between accuracy and computational expense. The filtered Navier-Stokes equations were developed including a total energy form of the energy equation. Subgrid scale models for the momentum and energy equations were adapted from compressible forms of Smagorinsky's original model. The effect of using disparate temporal and spatial accuracy in a numerical scheme was discovered through one-dimensional model problems and a new uniformly fourth-order accurate numerical method was developed. Results from two- and three-dimensional validation exercises show that the code accurately reproduces both viscous and inviscid flows. Numerous axisymmetric jet simulations were performed to investigate the effect of grid resolution, numerical scheme, exit boundary conditions and subgrid scale modeling on the solution and the results were used to guide the three-dimensional calculations. Three-dimensional calculations of a Mach 1.4 jet showed that this LES simulation accurately captures the physics of the turbulent flow. The agreement with experimental data was relatively good and is much better than results in the current literature. Turbulent intensities indicate that the turbulent structures at this level of modeling are not isotropic and this information could lend itself to the development of improved subgrid scale models for LES and turbulence models for RANS simulations. A two point correlation technique was used to quantify the turbulent structures. Two point space correlations were used to obtain a measure of the integral length scale, which proved to be approximately 1/2 D(sub j). Two point space-time correlations were used to obtain the convection velocity for the turbulent structures. This velocity ranged from 0.57 to 0.71 U(sub j).

DeBonis, James R.↗

Modeling of the EUV and X-Ray Emission Spectra Induced by the Solar Winds Ions in the Heliosphere

We have carried out investigation of the EUV and X-ray emission spectra induced in interaction between the Solar Wind (SW) and interstellar neutral gas. The spectra of most important SW ions have been computed for the charge-exchange mechanism of X-ray emission using new accurate spectroscopic data from recent laboratory measurements and theoretical calculations. Total spectra have been constructed as a sum of spectra induced in the charge-exchange collisions by individual O(exp q+), C(exp q+), N(exp q+), Ne(exp q+), Mg (exp q+) and Fe(exp q+) ions. Calculations have been performed for X-ray emission from the heliospheric hydrogen and helium gas. X-ray maps of the heliosphere have been computed. The power density of X-ray sources in the heliospheric ecliptic plane is shown for the H gas and for the He gas. Distances from the Sun (0,0) are given in AU. The helium cone is clear seen in the X-ray map of the charge-exchange emission induced by the solar wind. X-ray emission spectra detected by the Chandra X-ray telescope from the "dark" side of Moon has been identified as a X-ray background emission induced by the solar wind from the geocorona. Spectra and intensities of this charge-exchange X-rays have been compared with the heliospheric component of the X-ray background. Observations and modeling of the SW spectra induced from the geocorona indicate a strong presence of emission lines of highly charged oxygen ions. Anisotropy in distribution of heliospheric X-rays has been predicted and calculated for the regions of the fast and slow solar winds.

Kharchenko, Vasili↗