Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulated flux procedure”

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 37 records · Page 2

Application of a flux-split algorithm to chemically relaxing, hypervelocity blunt-body flows

Viscous, nonequilibrium, hypervelocity flow fields over two axisymmetric configurations are numerically simulated using a factored, implicit, flux-split algorithm. The governing gas-dynamic and species-continuity equations for laminar flow are presented. The gas-dynamics/nonequilibrium-chemistry coupling procedure is developed as part of the solution procedure and is described in detail. Numerical solutions are presented for hypervelocity flows over a hemisphere and over an axisymmetric aeroassisted orbital transfer vehicle using three different chemistry models. The gas models considered are those for an ideal gas, for a frozen gas, and for chemically relaxing air consisting of five species. The calculated results are compared with existing numerical solutions in the literature along the stagnation line of the hemisphere. The effects of free-stream Reynolds number on the nonequilibrium flow field are discussed.

Balakrishnan, A.↗

Development and application of computational aerothermodynamics flowfield computer codes

Multiple nozzle plume flow field is computed with a 3-D, Navier-Stokes solver. Numerical simulation is performed with a flux-split, two-factor, time asymptotic viscous flow solver of Ying and Steger. The two factor splitting provides a stable 3-D solution procedure under ideal-gas assumptions. An ad-hoc acceleration procedure that shows promise in improving the convergence rate by a factor of three for steady state problems is utilized. Computed solutions to generic problems at various altitude and flight conditions show flow field complexity and three-dimensional effects due to multiple nozzle jet interactions. Viscous, ideal gas solutions for the symmetric nozzle are compared with other numerical solutions.

Venkatapathy, Ethiraj↗

Cloud forcing and the net longwave surface radiation

A method is suggested which determines the monthly mean longwave radiation at the surface over the oceans using currently available satellite data sets. An independent estimate of the frequency of occurrence of clouds or cloud top heights is not required. A climatological relationship between the LWCRF (longwave cloud radiation forcing) at the top of the atmosphere and surface may be substituted. This ratio is obtained from a GMC simulation but in principle could be obtained observationally, although global coverage is doubtful. A flow diagram shows the procedure used to obtain maps for the monthly mean net upward longwave flux at the surface.

HARSHVARDHAN↗

Cloud Power Spectra-Dependence on Solar Zenith Angle and Wavelength, Implications for Cloud Optical Property Retrievals

Scale breaks (spatial scales at which power-law exponent changes occur) observed in Landsat radiances have proven to be useful indicators of radiative interactions, and have aided the development of improved techniques in the remote sensing of clouds. This work extends previous theoretical studies to absorbing wavelengths by using both Landsat Thematic Mapper (TM) observations and Monte Carlo (MC) simulations to infer the systematic dependencies of power spectral shape on cloud characteristics, illumination conditions, and wavelength. We show that MC simulations operating on a simple fractal model of horizontally inhomogeneous clouds produce power spectra that qualitatively resemble observed spectra. We also show that the decrease in the spectra power-law exponent seen at intermediate scales (referred to as "roughening") as the Sun becomes more oblique is more pronounced at absorbing wavelengths. An automated procedure designed to detect the small scale break location is unable to find systematic differences between TM Band 4 and Band 7, despite the fact that MC simulations point to systematic differences in horizontal fluxes. The effect of these qualitative characteristics of the spatial spectra on the retrieval of cloud optical properties is examined by comparing power spectra of nadir radiances with power spectra of optical properties retrieved using either traditional Independent Pixel Approximation approaches or modifications based on normalized radiance indices and the inverse Non-local Independent Pixel Approximation. Assuming that the actual cloud properties follow perfect scaling behavior at all scales, we show the improvement of the proposed retrieval modifications.

Oreopoulos, Lazaros↗

Quadrupole mass spectrometer for Space Shuttle applications - Flight capabilities and ground calibration

The AFGL Quadrupole Ion/Neutral Mass Spectrometer (QINMS), its flight capabilities, and a few examples of QINMS data from the STS-4 flight are discussed in addition to the beam facility and calibration procedures used at the Los Alamos National Laboratory facility. A 5-eV atomic oxygen beam is produced using a continuous optical discharge technique, and the resulting pumped beam has sufficient high flux and velocity to simulate orbital conditions. The calibrated instrument will be used in the NASA Johnson Space Center Evaluation of Oxygen Interactions with Materials experiment to measure the density of ambient oxygen and to identify the products of oxygen-surface reactions.

Hunton, D. E.↗

A multidomain global modeling technique for analysis of Space Shuttle main engine

A multidomain global modeling technique is presented. The purpose of this multidomain methodology is to simultaneously analyze the three-dimensional compressible flow field in Space Shuttle main engine components. This technique provides flexibility and modularity for individual domain and also simplifies the grid generation task for each domain. A two-way coupling procedure is adopted for exchanging interface boundary conditions. This algorithm which requires that the adjacent domains share an overlapping grid slab, is fully conservative and implicit. For non-matching grid systems and for non-overlapping interfaces, an interpolation and/or extrapolation procedure is used when transferring the interface boundary conditions at the adjacent domains. To maintain overall conservation, a flux correction method is subsequently used. The coupling procedure was verified by modeling several geometries with features similar to SSME components. Satisfactory results for all test cases were obtained. Finally, this algorithm has been applied to the SSME flow simulation.

Tam, L. T.↗

Simulation Of Unsteady, Viscous, Incompressible Flow

Method for numerical solution of Navier-Stokes equations of viscous, incompressible flow developed based on use of fractional-step procedure. Accurate to second order in both space and time. Attempt made to minimize Poisson-equation difficulties by choosing pressures at centers and volume fluxes across faces of computational cells as dependent variables instead of familiar Cartesian components of velocity. Choice ensures satisfaction of discrete equation of conservation of mass to within round-off errors in any coordinate system and has favorable effects on convergence properties.

Rosenfeld, Moshe↗

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗

On predicting contamination levels of HALOE optics aboard UARS using direct simulation Monte Carlo

A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. Emphasis is placed on simulating a realistic, worst-case set of flowfield and surface conditions and geometric orientations in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. Problems resolving species outgassing and vent flux rates that varied over many orders of magnitude were handled using species weighting factors. Results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface are presented, along with data related to code performance. Using procedures developed in standard contamination analyses, the cumulative level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated to be about 2700A.

Woronowicz, Michael S.↗

3D Simulations of Ion Thruster Accelerator Grid Erosion Accounting for Charge Exchange Ion Space Charge

Accelerator (accel) grid sputtering by ions formed through charge-exchange (CEX) reactions between beam ions and residual neutral gas is a critical life-limiting mechanism for gridded ion thrusters. The three-dimensional ion optics code CEX3D is designed to simulate this grid erosion for a single beamlet, with a particular emphasis on non-axisymmetric features such as the "pits and grooves" erosion commonly observed on the accel grid downstream face in two-grid thrusters. The treatment of CEX ions in the code was recently upgraded with a new particle-in-cell (PIC) module to account for the influence of these ions' space charge on the electrostatic potential downstream of the grids. In order to achieve reasonable computation times while resolving the Debye length near the grids and avoiding gross violations of the Courant-Friedrichs-Lewy (CFL) condition, macroparticle velocities in the PIC calculation are limited through a rescaling procedure that preserves ion trajectories and space charge density. The code accounts for beam divergence, finite momentum transfer in CEX collisions, and radial losses of CEX ions from the beam; these effects are important for determining the CEX ion flux to the accel grid because the calculated potential downstream of the grids can become very flat. The upgraded code has been used to simulate operation of NASA's Evolutionary Xenon Thruster (NEXT) during the 51 kHr Long Duration Test - a selection of results is presented and compared with experimental data.

Chaplin, Vernon H.↗

Windblown sand on Venus - Preliminary results of laboratory simulations

Small particles and winds of sufficient strength to move them have been detected from Venera and Pioneer-Venus data and suggest the existence of aeolian processes on Venus. The Venus wind tunnel (VWT) was fabricated in order to investigate the behavior of windblown particles in a simulated Venusian environment. Preliminary results show that sand-size material is readily entrained at the wind speeds detected on Venus and that saltating grains achieve velocities closely matching those of the wind. Measurements of saltation threshold and particle flux for various particle sizes have been compared with theoretical models which were developed by extrapolation of findings from Martian and terrestrial simulations. Results are in general agreement with theory, although certain discrepancies are apparent which may be attributed to experimental and/or theoretical-modeling procedures. Present findings enable a better understanding of Venusian surface processes and suggest that aeolian processes are important in the geological evolution of Venus.

Greeley, R.↗

An investigation of the energy balance of solar active regions using the ACRIM irradiance data

The correlation between the irradiance (as measured by the Active Cavity Radiometer Irradiance Monitor and the Earth Radiation Budget) as corrected for sunspot flux deficit (which is responsible for most of the variance of the uncorrected signal) and both the 205 nm flux (as measured by Nimbus 7) and a photometric facular index is discussed. The computer program which simulates two-dimensional convection in a compressible, stratified medium is described. Equipment which was acquired to perform high precision, white-light observations of sunspot areas, and procedures were tested. Analysis of observations of large scale convective heat inhomogeneities which were obtained in May 1985 was begun.

Petro, L. D.↗

Finite Element Analysis of Poroelastic Composites Undergoing Thermal and Gas Diffusion

A theory for time-dependent thermal and gas diffusion in mechanically time-rate-independent anisotropic poroelastic composites has been developed. This theory advances previous work by the latter two authors by providing for critical transverse shear through a three-dimensional axisymmetric formulation and using it in a new hypothesis for determining the Biot fluid pressure-solid stress coupling factor. The derived governing equations couple material deformation with temperature and internal pore pressure and more strongly couple gas diffusion and heat transfer than the previous theory. Hence the theory accounts for the interactions between conductive heat transfer in the porous body and convective heat carried by the mass flux through the pores. The Bubnov Galerkin finite element method is applied to the governing equations to transform them into a semidiscrete finite element system. A numerical procedure is developed to solve the coupled equations in the space and time domains. The method is used to simulate two high temperature tests involving thermal-chemical decomposition of carbon-phenolic composites. In comparison with measured data, the results are accurate. Moreover unlike previous work, for a single set of poroelastic parameters, they are consistent with two measurements in a restrained thermal growth test.

Salamon, N. J.↗

A redshift survey of IRAS galaxies. III - Reconstruction of the velocity and density fields in N-body model universes

N-body simulations of a 'cold dark matter' universe are presently used to calibrate the accuracy, and assess the limitations, of the procedure previously employed to predict the velocity field within 8000 km/sec of the Local Group through the application of linear gravitational theory to a full-sky, flux-limited sample of IRAS galaxies. The rms difference between the one-dimensional acceleration and velocity of field particles is an increasing function of local density; linear theory can in this way account for all but one-sixth of kinetic energy. A series of artificial IRAS catalogs closely matching the real sample in space density and clustering amplitude is constructed. Velocity correlation functions are used to demonstrate that the predicted velocity fields are in good agreement with the true velocity fields on large scales.

Davis, Marc↗

Developments and applications of dynamic models for large eddy simulation of complex flows

The dynamic modeling procedure for large eddy simulation of turbulent flows is reviewed and recent developments in the theoretical aspects and applications are described. Methods for inclusion of backscatter of energy from small to large scale motions are presented. New formulations of the dynamic procedure are proposed which are optimized based on the subgrid scale flux vector or the energy dissipation rate instead of the subgrid scale stress tensor. Recent results from application of the model to forced isotropic turbulence with an inertial subrange, flow over a backward facing step at Reynolds number of 28000, and flow over a concave curved surface are presented.

Moin, P.↗

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

Experimental study of solar simulator mirror cryocontamination

The background and tasks formulation of the study of Solar Simulator collimation mirror cryocontamination in Large Thermal Vacuum Facility are outlined, research methods and experiment procedures are described, experimental relationships obtained are analyzed and practical recommendations are given. The accepted procedure of thermal vacuum tests as a rule defines the sequence of operations for verifying the spacecraft under test without taking into account measures for preventing Solar Simulator collimation mirror contamination and degradation. On the other hand, evacuation procedures is defined for conditions of achieving the required vacuum in the shortest possible time with using the available evacuation equipment at a regime close to the optimum one. Similarly, cryopanel cooling down cyclogram and test object preparation process are not analyzed from the viewpoint of ways of reducing environmental detrimental effects on thermal vacuum facility contamination-sensitive systems. Solar Simulator mirror contamination and its reflective characteristics change results in degradation of solar flux parameters and reduction of simulator continuous operation time. Methods of consideration of optical effects due to mirror surface contamination are actually missing. The effects themselves are not quite understood and data cited in literature as a rule, were obtained under conditions different from real thermal vacuum facility and therefore should be subjected to additional experimental verification. Only in the last few years contamination effect on optical surfaces degradation has been considered with using empirical relations. Mirror reflective properties degradation leads to the increase of Solar Simulator errors. This ultimately has an adverse effect on S/C ground development, schedule and cost of thermal vacuum tests. Besides, the mirror maintenance in operable state becomes more expensive. The present paper is dedicated to the study of Solar Simulator collimation mirror contamination and to the search of ways for improving the mirror design and thermal vacuum test procedure. On the basis of tests performed, recommendations are devised on reducing chamber inner-optical surfaces cryocontamination and degradation.

Galjaev, V. L.↗

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗