Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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 757 records · Page 42

Models for the atmospheres of Jupiter and Saturn

Numerical models of the atmospheres of Jupiter and Saturn have been created for use in the development of design criteria for space vehicles intended to investigate these plants. These models provide a basis for the assessment of entry heating, the transmission and emission of radiation, and structural and aerodynamic interactions. The model atmospheres are based on data and analyses published in the literature through 1970. The composition is duscussed, and the relative amounts of each molecular species are tabulated for nominal and extreme cases. The structure of the atmospheric regions is discussed, and the principal features of one nominal and two limiting models for each planet are indicated. The models are illustrated by pressure-temperature profiles which also display the major cloud layers.

Divine, N.↗

Helium-burning flashes on an accreting neutron star - A model for X-ray burst sources

Detailed numerical models of X-ray bursts resulting from thermonuclear flashes near the surface of an accreting neutron star have been constructed. The models assume a nonrotating nonmagnetized spherically accreting neutron star of 1.4 solar masses, radius 6.6 km, core temperature in the range from 250 million to 570 million K, and accretion rate in the range (0.3-3) x 10 to the 17th power g/s. Under many conditions the helium-burning shell undergoes thermonuclear flashes that result in the emission of X-ray bursts, the gross properties of which are remarkably similar to those of most observed X-ray burst sources. Neutron stars with moderately low core temperatures, low accretion rates, and weak magnetic fields are most likely to produce X-ray bursts.

Joss, P. C.↗

A model of active prominences

A one-dimensional numerical model of active loop prominences is investigated. The model includes the effects of gravity, the geometry of the magnetic field, and conduction losses to the chromosphere. Calculations indicate that, as originally proposed by Goldsmith, the thermal instability mechanism is, by itself, sufficient to account for the appearance of bright H-alpha knots in postflare loops. Under certain conditions, initial perturbations in the loop temperature and density profiles of small, but finite, amplitude (approximately 5%) and large size scale (greater than or approximately equal to 10 to the 9th cm) can grow into condensations with temperature and density differences of over an order of magnitude and size scales of less than 10 to the 8th cm. In agreement with observations, the conditions that must be satisfied are such that loop prominence systems are likely to occur only in large flares and such that knots preferentially form at the tops of loops. The velocities, densities, and lifetimes calculated for the loop material are also in agreement with observations. It is concluded that in order for H-alpha knots to occur, heating of some form must continue into the late cooling phase of a flare loop, and that this heating is more intense near the loop base than near the apex.

Antiochos, S. K.↗

Modeling a variable thickness sea ice cover

A numerical model simulating a variable thickness sea ice cover over a seasonal cycle is presented. The model includes a fixed depth mixed-layer formulation with open water heat absorption and lateral melting terms, and a mechanical distribution function consistent with the physics of the ridging process. The equibrium simulation results in realistic geographical ice thickness variations of April ice along the Canadian Archipelago which exceed 7 m, and thicknesses of about 2 m along the Alaskan North Slope. Ice velocity fields were realistic in shape but 25% larger than the net ice station drift over a year; sensitivity simulations indicated a reduced average annual ice export rate of 0.04 Sv as compared to 0.09 Sv for the equilibrium simulation.

Hibler, W. D., III↗

Characterization and modeling for dilute and dense sprays

The present numerical modeling of dilute and dense spray combustion solves the governing gas-phase equations in Eulerian coordinates by means of a time-marching multiple pressure correction procedure that is based on the operator-splitting technique. A sequence of validating cases involving (1) nonevaporating, (2) evaporating, (3) burning, (4) dilute, and (5) dense spray cases is undertaken; it is found that in the combusting dilute spray case, the present procedure correctly predicts the general features of flows and yields qualitative agreement with experimental data. In the other cases, reasonably good agreement with experimental results is obtained.

Chen, C. P.↗

Thunderstorm and Lightning Development: Modeling and Observations

We use results of a numerical model of convective clouds with a lightning parameterization to examine thunderstorm development and the location of lightning initiation in a Florida storm. These results are compared with radar observations and interferometric measurements of VHF emissions from both intra-cloud (IC) and cloud-to-ground (CG) lightning flashes. Observed rate and mean duration of flashes are analyzed through model results, and we comment altitudes of VHF sources given by the ONERA 3D interferometric system.

Solomon, Robert↗

On the effect of acoustic coupling on random and harmonic plate vibrations

The effect of acoustic coupling on random and harmonic plate vibrations is studied using two numerical models. In the coupled model, the plate response is obtained by integration of the nonlinear plate equation coupled with the nonlinear Euler equations for the surrounding acoustic fluid. In the uncoupled model, the nonlinear plate equation with an equivalent linear viscous damping term is integrated to obtain the response of the plate subject to the same excitation field. For a low-level, narrow-band excitation, the two models predict the same plate response spectra. As the excitation level is increased, the response power spectrum predicted by the uncoupled model becomes broader and more shifted towards the high frequencies than that obtained by the coupled model. In addition, the difference in response between the coupled and uncoupled models at high frequencies becomes larger. When a high intensity harmonic excitation is used, causing a nonlinear plate response, both models predict the same frequency content of the response. However, the level of the harmonics and subharmonics are higher for the uncoupled model. Comparisons to earlier experimental and numerical results show that acoustic coupling has a significant effect on the plate response at high excitation levels. Its absence in previous models may explain the discrepancy between predicted and measured responses.

Frendi, A.↗

Effect of acoustic coupling on random and harmonic plate vibrations

The effect of acoustic coupling on random and harmonic plate vibrations is studied using two numerical models. In the coupled model, the plate response is obtained by integration of the nonlinear plate equation coupled with the nonlinear Euler equations for the surrounding acoustic fluid. In the uncoupled model, the nonlinear plate equation with an equivalent linear viscous damping term is integrated to obtain the response of the plate subject to the same excitation field. For a low-level, narrow-band excitation, the two models predict the same plate response spectra. As the excitation level is increased, the response power spectrum predicted by the uncoupled model becomes broader and more shifted towards the high frequencies than that obtained by the coupled model. In addition, the difference in response between the coupled and uncoupled models at high frequencies becomes larger. When a high intensity harmonic excitation is used, causing a nonlinear plate response, both models predict the same frequency content of the response. However, the level of the harmonics and subharmonics are higher for the uncoupled model. Comparisons to earlier experimental and numerical results show that acoustic coupling has a significant effect on the plate response at high excitation levels. Its absence in previous models may explain the discrepancy between predicted and measured responses.

Frendi, Abdelkader↗

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗

NASA National Combustion Code Simulations

A systematic effort is in progress to further validate the National Combustion Code (NCC) that has been developed at NASA Glenn Research Center (GRC) for comprehensive modeling and simulation of aerospace combustion systems. The validation efforts include numerical simulation of the gas-phase combustor experiments conducted at the Center for Turbulence Research (CTR), Stanford University, followed by comparison and evaluation of the computed results with the experimental data. Presently, at GRC, a numerical model of the experimental gaseous combustor is built to simulate the experimental model. The constructed numerical geometry includes the flow development sections for air annulus and fuel pipe, 24 channel air and fuel swirlers, hub, combustor, and tail pipe. Furthermore, a three-dimensional multi-block, multi-grid grid (1.6 million grid points, 3-levels of multi-grid) is generated. Computational simulation of the gaseous combustor flow field operating on methane fuel has started. The computational domain includes the whole flow regime starting from the fuel pipe and the air annulus, through the 12 air and 12 fuel channels, in the combustion region and through the tail pipe.

Iannetti, Anthony↗

Thermal Property Estimation of Fibrous Insulation: Heat Transfer Modeling and the Continuous Genetic Algorithm

Thermal properties of high-temperature fibrous insulation materials were estimated from transient thermal tests using an inverse heat transfer technique. Transient temperature data from an experimental set up was collected and a simple, one-dimensional numerical model was constructed to replicate the temperatures within the test assembly. The Continuous Genetic Algorithm optimization technique in conjunction with a numerical thermal model and transient test data was used to estimate coefficients of a functional representation of the thermal property. The thermal properties, i.e., thermal conductivity and specific heat, of an alumina insulation felt were estimated over the temperature range of 300 K to 1700 K at various constant static pressures in nitrogen gas and compared with published data. The resulting thermal property estimates were within 10% of published values over the entire temperature range at various pressures. The methodology, application, and results are presented.

Frye, Elora↗

A detailed model for simulation of catchment scale subsurface hydrologic processes

A catchment scale numerical model is developed based on the three-dimensional transient Richards equation describing fluid flow in variably saturated porous media. The model is designed to take advantage of digital elevation data bases and of information extracted from these data bases by topographic analysis. The practical application of the model is demonstrated in simulations of a small subcatchment of the Konza Prairie reserve near Manhattan, Kansas. In a preliminary investigation of computational issues related to model resolution, we obtain satisfactory numerical results using large aspect ratios, suggesting that horizontal grid dimensions may not be unreasonably constrained by the typically much smaller vertical length scale of a catchment and by vertical discretization requirements. Additional tests are needed to examine the effects of numerical constraints and parameter heterogeneity in determining acceptable grid aspect ratios. In other simulations we attempt to match the observed streamflow response of the catchment, and we point out the small contribution of the streamflow component to the overall water balance of the catchment.

Paniconi, Claudio↗

A technique for measuring vertically and horizontally polarized microwave brightness temperatures using electronic polarization-basis rotation

This technique for electronically rotating the polarization basis of an orthogonal-linear polarization radiometer is based on the measurement of the first three feedhorn Stokes parameters, along with the subsequent transformation of this measured Stokes vector into a rotated coordinate frame. The technique requires an accurate measurement of the cross-correlation between the two orthogonal feedhorn modes, for which an innovative polarized calibration load was developed. The experimental portion of this investigation consisted of a proof of concept demonstration of the technique of electronic polarization basis rotation (EPBR) using a ground based 90-GHz dual orthogonal-linear polarization radiometer. Practical calibration algorithms for ground-, aircraft-, and space-based instruments were identified and tested. The theoretical effort consisted of radiative transfer modeling using the planar-stratified numerical model described in Gasiewski and Staelin (1990).

Gasiewski, Albin J.↗

Dynamic fragmentation in impacts - Hydrocode simulation of laboratory impacts

The dynamic fragmentation in impacts into solids is examined using a physical model for the formation and growth of cracks in rocks. The physical model is then inserted into a numerical model (hydrocode) of stress wave propagation and interaction, from which the outcome of a given impact event can be computed. The hydrocode model predicts fragment sizes due to impact in terms of shock waves propagating in a homogeneous elastic medium containing a distribution of crack nucleation centers known as Weibull flaws.

Melosh, H. J.↗

Forward Monte Carlo Computations of Polarized Microwave Radiation

Microwave radiative transfer computations continue to acquire greater importance as the emphasis in remote sensing shifts towards the understanding of microphysical properties of clouds and with these to better understand the non linear relation between rainfall rates and satellite-observed radiance. A first step toward realistic radiative simulations has been the introduction of techniques capable of treating 3-dimensional geometry being generated by ever more sophisticated cloud resolving models. To date, a series of numerical codes have been developed to treat spherical and randomly oriented axisymmetric particles. Backward and backward-forward Monte Carlo methods are, indeed, efficient in this field. These methods, however, cannot deal properly with oriented particles, which seem to play an important role in polarization signatures over stratiform precipitation. Moreover, beyond the polarization channel, the next generation of fully polarimetric radiometers challenges us to better understand the behavior of the last two Stokes parameters as well. In order to solve the vector radiative transfer equation, one-dimensional numerical models have been developed, These codes, unfortunately, consider the atmosphere as horizontally homogeneous with horizontally infinite plane parallel layers. The next development step for microwave radiative transfer codes must be fully polarized 3-D methods. Recently a 3-D polarized radiative transfer model based on the discrete ordinate method was presented. A forward MC code was developed that treats oriented nonspherical hydrometeors, but only for plane-parallel situations.

Battaglia, A.↗

Spectral Measurements of PMCs from SBUV/2 Instruments

The SBUV/2 (Solar Backscattered Ultraviolet, model 2) instrument is designed to monitor ozone stratospheric profile and total column ozone using measurements of the Earth's backscattered ultraviolet albedo. We have previously demonstrated that the normal radiance measurements from SBUV/2 instruments, which sample 12 discrete wavelengths between 252 and 340 nm during each scan, can be used to identify polar mesospheric clouds (PMCs). Some SBUV/2 instruments also periodically view the earth in continuous scan mode, covering the wavelength range 160-400 nm with 0.15 nm sampling. Analysis of these data show PMC occurrence rates similar to the normal discrete scan results, although the observation technique reduces the number of daily measurements by a factor of six. PMC observed by SBUV/2 instruments show a monotonic variation in the residual spectral albedo over the wavelength range 250 300 nm, with maximum enhancements of 10 15% at 250 nm. This result is consistent with microphysical model predictions from Jensen [1989. A numerical model of polar mesospheric cloud formation and evolution, Ph. D. Thesis, University of Colorado]. We find no evidence for a systematic localized increase in PMC residual albedo for wavelengths near 260 nm, in contrast to the recently reported results from the MSX UVISI instrument [Carbary J.F., et al., 2004. Evidence for bimodal particle distribution from the spectra of polar mesospheric clouds. Geophysics Research. Letters 31, L13108]. This result is observed for three different SBUV/2 instruments in both Northern and Southern Hemisphere data over a 13-year span. Our Mie scattering calculations show that the location and magnitude of the 260 nm hump feature is dependent upon the specific scattering angles appropriate to the MSX measurements. Although it explains the MSX spectrum, the bimodal size distribution proposed by Carbary et al. (2004), cannot explain the lack of scattering angle dependence of the SBUV/2 spectral shapes. The spectral signature of the SBUV/2 continuous scan PMC data is thus inconsistent with the bimodal particle size distribution suggested by Carbary et al. (2004).

DeLand, Matthew T.↗

TPSAS-NF1676L-17931-DND

Many geophysical variables including temperature, clouds, and precipitation exhibit robust diurnal cycles in response to the daily cycle of solar insolation. Due to the fundamental nature of this variability, it is critical that weather and climate models accurately represent the diurnal cycle. Numerical model, however, have difficulty reproducing this observed diurnal cycle behavior, which leads to systematic errors in model representation of earth energy budget terms: including TOA and surface radiation, precipitation, and surface latent and sensible heat fluxes. This study evaluates the regional diurnal cycle in the Tropics within reanalysis models (ERA-Interim and MERRA) and quantifies systematic errors in the simulated TOA flux and precipitation due to biases in the diurnal cycle representation. The focus of this study is to quantify the importance of diurnal cycle simulation to systematic bias in reanalysis data set climatologies.

Patrick C Taylor↗

Computational Analysis of a Laminar Jet in a Cross-Flow

A jet in cross-flow (JIFC) consists of a jet exhausting at a large angle into a freestream flow. It is a flow field which is relevant to a wide variety of technologies and applications. Despite the nearly 65 years of JIFC research there are few results available for laminar hypersonic flows, a combination which will be encountered by re-entry and high altitude vehicles over some portion of their flight path. This research consists of developing a numerical model to investigate the interaction of a normal sonic jet exhausting into a hypersonic cross-flow. The model was validated by comparing experimental measurements with corresponding numerical results generated by the model.

Jeffries, Damon K.↗