Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 1,297 records · Page 72

Numerical Simulation of Ion Thruster Optics

A three-dimensional simulation code (ffx) designed to analyze ion thruster optics is described. It is an extension of an earlier code and includes special features like the ability to model a wide range of grid geometries, cusp details, and mis-aligned aperture pairs to name a few. However, the principle reason for advancing the code was in the study of ion optics erosion. Ground based testing of ion thruster optics, essential to the understanding of the processes of grid erosion, can be time consuming and costly. Simulation codes that can accurately predict grid lifetimes and the physical mechanisms of grid erosion can be of great utility in the development of future ion thruster optics designed for more ambitious applications. Results of simulations are presented that describe wear profiles for several standard and nonstandard aperture geometries, such as those grid sets with square- or slotted-hole layout patterns. The goal of this paper will be to introduce the methods employed in the ffx code and to briefly demonstrate their use.

Rawlin, Vincent K.↗

Large eddy simulation of turbulent channel flow: ILLIAC 4 calculation

The three-dimensional time dependent equations of motion were numerically integrated for fully-developed turbulent channel flow. A large scale flow field was obtained directly from the solution of these equations, and small scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC 4 computer. The computed flow patterns show that the wall layer consists of coherent structures of low speed and high speed streaks alternating in the spanwise direction. These structures were absent in the regions away from the wall. Hot spots, small localized regions of very large turbulent shear stress, were frequently observed. The profiles of the pressure velocity-gradient correlations show a significant transfer of energy from the normal to the spanwise component of turbulent kinetic energy in the immediate neighborhood of the wall ('the splatting effect').

Kim, J.↗

Unsteady separated boundary layer in a transonic diffuser flow with self-excited oscillations

A numerical investigation of two-dimensional unsteady boundary layer in a transonic diffuser flow with self-excited oscillations and strong flow separation by solving the compressible, Reynolds-averaged, thin-layer Navier-Stokes equations with two-equations turbulence model is described. Three different meshes with constant streamwise mesh distribution and varying vertical mesh distribution were used. Results obtained indicate that a refinement of mesh studied here has minimal effect on the mean boundary layer flow but significantly increases the amplitude of oscillation of all flow variables. Comparisons of unsteady wall pressure, velocity profile, terminal shock, and separation pocket among computations and with experiment are presented.

Hsieh, T.↗

Numerical Solutions for the CAWAPI Configuration on Unstructured Grids at NASA LaRC, United States

In support of the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) with its goal of improving the Technology Readiness Level of flow solvers by comparing results with measured F-16XL-1 flight data, NASA Langley employed the TetrUSS unstructured grid solver, USM3D, to obtain solutions for all seven flight conditions of interest. A newly available solver version that incorporates a number of turbulence models, including the two-equation linear and non-linear k- , was used in this study. As a first test, a choice was made to utilize only a single grid resolution with the solver for the simulation of the different flight conditions. Comparisons are presented with three turbulence models in USM3D, flight data for surface pressure, boundary-layer profiles, and skin-friction distribution, as well as limited predictions from other solvers. A result of these comparisons is that the USM3D solver can be used in an engineering environment to predict vortex-flow physics on a complex configuration at flight Reynolds numbers with a two-equation linear k- turbulence model.

Lamar, John E.↗

Control theory based airfoil design using the Euler equations

This paper describes the implementation of optimization techniques based on control theory for airfoil design. In our previous work it was shown that control theory could be employed to devise effective optimization procedures for two-dimensional profiles by using the potential flow equation with either a conformal mapping or a general coordinate system. The goal of our present work is to extend the development to treat the Euler equations in two-dimensions by procedures that can readily be generalized to treat complex shapes in three-dimensions. Therefore, we have developed methods which can address airfoil design through either an analytic mapping or an arbitrary grid perturbation method applied to a finite volume discretization of the Euler equations. Here the control law serves to provide computationally inexpensive gradient information to a standard numerical optimization method. Results are presented for both the inverse problem and drag minimization problem.

Jameson, Antony↗

Crystal Growth of ZnSe by Physical Vapor Transport: A Modeling Study

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of Physical Vapor Transport (PVT) growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe arc reported in this paper. The PVI process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The preliminary results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided.

Ramachandran, Narayanan↗

Assessment Study of Small Space Debris Removal by Laser Satellites

Space debris in Earth orbit poses significant danger to satellites, humans in space, and future space exploration activities. In particular, the increasing number of unidentifiable objects, smaller than 10 cm, presents a serious hazard. Numerous technologies have been studied for removing unwanted objects in space. Our approach uses a short wavelength laser stationed in orbit to vaporize these small objects. This paper discusses the power requirements for space debris removal using lasers. A short wavelength laser pumped directly or indirectly by solar energy can scan, identify, position, and illuminate the target, which will then be vaporized or slow down the orbital speed of debris by laser detonation until it re-enters the atmosphere. The laser-induced plasma plume has a dispersive motion of approximately 105 m/sec with a Lambertian profile in the direction of the incoming beam [1-2]. The resulting fast ejecting jet plume of vaporized material should prevent matter recombination and condensation. If it allows any condensation of vaporized material, the size of condensed material will be no more than a nanoscale level [3]. Lasers for this purpose can be indirectly pumped by power from an array of solar cells or directly pumped by the solar spectrum [4]. The energy required for vaporization and ionization of a 10 cm cube (~ 2700 gm) of aluminum is 87,160 kJ. To remove this amount of aluminum in 3 minutes requires a continuous laser beam power of at least 5.38 MW under the consideration of 9% laser absorption by aluminum [5] and 5% laser pumping efficiency. The power needed for pumping 5.38 MW laser is approximately 108 MW, which can be obtained from a large solar array with 40% efficiency solar cells and a minimal area of 450 meters by 450 meters. This solar array would collect approximately 108 MW. The power required for system operation and maneuvering can be obtained by increasing solar panel size. This feasibility assessment covers roughly the power requirement, laser system, and a potential operational scenario.

Choi, Sang H.↗

The role of forced planetary waves in the annual cycle of the zonal mean circulation of the middle atmosphere

A severely truncated semispectral numerical model is used to simulate the annual cycle of the zonally averaged circulation in the middle atmosphere (16-96 km). The model includes only a single zonal harmonic wave component which interacts with the mean flow; the circulation is driven by diabatic heating and by a specified perturbation in the topography of the lower boundary, which is taken to be the 100 mb surface. A comparison of the annual cycle simulated by this model with the results of an analogous two-dimensional model indicates that planetary waves have relatively little influence on the zonal mean temperature profiles and on the solstice mean zonal winds at high latitudes. The primary effects of the forced waves are in decelerating the mean winds at low latitudes in the winter hemisphere to produce a region of weak westerlies, and in generating final warmings at the spring equinoxes.

Holton, J. R.↗

A comparison of radiative transfer models for predicting the microwave emission from soils

Noncoherent and coherent numerical models for predicting emission from soils are compared. Coherent models use the boundary conditions on the electric fields across the layer boundaries to calculate the radiation intensity, and noncoherent models consider radiation intensities directly. Interference may cause different results in the two approaches when coupling between soil layers in coherent models causes greater soil moisture sampling depths. Calculations performed at frequencies of 1.4 and 19.4 GHz show little difference between the models at 19.4 GHz, although differences are apparent at the lower frequency. A definition for an effective emissivity is also given for when a nonuniform temperature profile is present, and measurements made from a tower show good agreement with calculations from the coherent model.

Schmugge, T. J.↗

Evolution of rotationally and tidally distorted low-mass, close binary systems - Implications for the minimum orbital period of cataclysmic variables

A (1 + 0.4) solar mass close binary system consisting of a compact primary and a red dwarf secondary has been evolved numerically. Such a binary system should effectively model cataclysmic variables for which a minimum orbital period cutoff of about 81 minutes has been observed. The influence of gravitational radiation losses which drive Roche lobe overflow has been studied, and the effects of rotational and tidal distortion have also been incorporated in the calculation. The evolution of the He-3 abundance, which has been suggested as a possible explanation for the upper limit of the apparent orbital period gap exhibited by cataclysmic variables, is also considered. It is found that both the distortional effects and the He-3 chemical profile can play an important role in determining the subsequent evolution of these systems. Specifically, when distortion is included, the theoretical minimum orbital period is increased by about 10 percent, yielding better agreement with observations.

Nelson, L. A.↗

Modelling of three-dimensional shock wave turbulent boundary layer interactions

The accuracy of numerical computations of the three-dimensional sharp fin flow field using the Reynolds-averaged Navier-Stokes equations is examined. Computed flows are compared with a set of benchmark experiments at Mach 3 for two different Reynolds numbers. A detailed comparison is performed between experimental data and separate computational results for the three-dimensional sharp fin for a 10 deg fin angle in order to examine the accuracy of two different turbulence models. The computed surface pressure for the fin at a 20 deg angle is compared with recent experimental data. The computed yaw angle profiles are shown to be sensitive to the turbulence model. In the immediate vicinity of the surface, the calculated results using the Baldwin-Lomax (1978) model are in better agreement with the data than the Jones-Launder (1972) model. Within the inner portion of the boundary layer, both models tend to underpredict the yaw angle.

Knight, D. D.↗

Mixing control in a plane shear layer

An investigation of mixing processes in a plane shear layer by a direct numerical simulation with a high-order numerical technique, the isoparametric spectral element method, is presented. The computational domain includes a region with the splitter plate that allows to start the flow as a laminar (Blasius) boundary layers with different velocities and concentrations, merging in the shear layer and undergoing transition to highly-unsteady structures with intensive mixing. The augmentation of mixing as a function of a downstream location, the intensity and type of vortical structures in the layer is investigated by an introduction of disturbance at the high-velocity inflow. The character of the dynamic mixing is displayed along with time averaged profiles and statistical characteristics.

Korczak, K. Z.↗

A Stereo Imaging Velocimetry Technique for Analyzing Structure of Flame Balls at Low Lewis-Number (SOFBALL) Data

Stereo Imaging Velocimetry (SIV) is a NASA Glenn Research Center (GRC) developed fluid physics technique for measuring threedimensional (3-D) velocities in any optically transparent fluid that can be seeded with tracer particles. SIV provides a means to measure 3-D fluid velocities quantitatively and qualitatively at many points. This technique provides full-field 3-D analysis of any optically clear fluid or gas experiment using standard off-the-shelf CCD cameras to provide accurate and reproducible 3-D velocity profiles for experiments that require 3-D analysis. A flame ball is a steady flame in a premixed combustible atmosphere which, due to the transport properties (low Lewis-number) of the mixture, does not propagate but is instead supplied by diffusive transport of the reactants, forming a premixed flame. This flame geometry presents a unique environment for testing combustion theory. We present our analysis of flame ball phenomena utilizing SIV technology in order to accurately calculate the 3-D position of a flame ball(s) during an experiment, which can be used as a direct comparison of numerical simulations.

McDowell, Mark↗

Evaluating the Impact of Geostationary Sounders in the Context of International Coordination

The proposed NOAA/NASA Geostationary eXtended Observations (GeoXO) program plans to include a hyperspectral infrared (IR) sounder (GXS) on its central satellite, joining international counterparts. Ahead of launch, the NASA Global Modeling and Assimilation Office (GMAO) assessed the potential effectiveness of GXS both as a lone sounder in a GEO orbit and as part of a global ring of such instruments. Using an observing system simulation experiment (OSSE) framework from a global numerical weather prediction (NWP) perspective, the ability of GXS and the global ring to improve weather prediction of thermodynamic variables was assessed both globally and regionally. GXS dominated regional analysis and forecast improvements and contributed significantly to global increases in forecast skill. However, more sustained global improvements on the order of 4 days rely on international partnerships. Over CONUS, the FSOI metric showed the GXS observations provide the strongest radiance impact on the moist energy error norm reduction. Additionally, GXS shows the capability to improve hurricane forecast track errors, resulting in improved forecast warnings. Overall, the persistent atmospheric profile information from GXS over much of the western hemisphere provide an opportunity to improve the representation of weather systems and their forecasts.

Erica L. McGrath-Spangler↗

Evaluating the Impact of Geostationary Sounders in the Context of International Coordination

The proposed NOAA/NASA Geostationary eXtended Observations (GeoXO) program plans to include a hyperspectral infrared (IR) sounder (GXS) on its central satellite, joining international counterparts. Ahead of launch, the NASA Global Modeling and Assimilation Office (GMAO) assessed the potential effectiveness of GXS both as a lone sounder in a GEO orbit and as part of a global ring of such instruments. Using an observing system simulation experiment (OSSE) framework from a global numerical weather prediction (NWP) perspective, the ability of GXS and the global ring to improve weather prediction of thermodynamic variables was assessed both globally and regionally. GXS dominated regional analysis and forecast improvements and contributed significantly to global increases in forecast skill. However, more sustained global improvements on the order of 4 days rely on international partnerships. Over CONUS, the FSOI metric showed the GXS observations provide the strongest radiance impact on the moist energy error norm reduction. Additionally, GXS shows the capability to improve hurricane forecast track errors, resulting in improved forecast warnings. Overall, the persistent atmospheric profile information from GXS over much of the western hemisphere provide an opportunity to improve the representation of weather systems and their forecasts.

Erica McGrath-Spangler↗

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗

Radiation damage of the HEAO C-1 germanium detectors

The effects of radiation damage from proton bombardment of the four HEAO C-1 high purity germanium detectors have been measured and compared to predictions. Because of the presence of numerous gamma-ray lines in the detector background spectra and because of the relatively long exposure time of the HEAO 3 satellite to cosmic-ray and trapped protons, it has been possible to measure both the energy and time dependence of radiation damage. After 100 d in orbit, each of the four detectors has been exposed to approximately 3 x 10 to the 7th protons/sq cm, and the average energy resolution at 1460 keV had degraded from 3.2 keV fwhm to 8.6 keV fwhm. The lines were all broadened to the low energy side although the line profile was different for each of the four detectors. The damage-related contribution to the degradation in energy resolution was found to be linear in energy and proton influence.

Mahoney, W. A.↗

Approximate Turbulent Boundary-Layer Development in Plane Compressible Flow Along Thermally Insulated Surfaces with Application to Supersonic-Tunnel Contour Correction

Numerical solutions of the differential equation obtained from the momentum theorem for the development of a turbulent boundary layer along a thermally insulated surface in two-dimensional and in radial shock-free flow are presented in tabular form for a range of Mach numbers from 0.100 to 10. The solution can be used in a step-wise procedure with any given distribution of favorable pressure gradients and for zero pressure gradients. Solutions are also given for use with moderate adverse pressure gradients. The mean velocity in the boundary layer is approximated by a power-law profile. In view of the stepwise integration methods to be used, the exponent designated the profile shape can be varied along the surface between the integral fraction limits 1/5 and 1/11 through interpolation. Agreement obtained between theoretical and experimental boundary-layer development in a supersonic nozzle at a nominal Mach number of 2 indicates the general validity of the approximations used in the analysis - in particular, the method of extrapolating low-speed skin-friction relations to high Mach number flows. The extrapolation method used assumes that the skin-friction coefficient depend primarily on Reynolds number, provided that the density and the kinematic viscosity are evaluated at surface conditions.

BOUNDARY LAYER - WIND TUNNEL WALLS↗