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 811 records · Page 45

Constraints on Mars sampling based on models of basaltic flow surfaces and interiors

Recent field observation and numerical modelling of the pattern and origin of vesicle zones and joints in terrestrial basaltic flows has resulted in increased understanding of the processes which affect flow surface morphology. This work has documented the ubiquitous occurrence of three vertical zones in basalt flows: (1) an upper vesicular zone; (2) a middle vesicle-free zone; and (3) a lower vesicular zone. The upper vesicular zone is generally about one-half of the total flow thickness. Computer modeling of the development of these zones confirms that vesicle zonation is a result of the nucleation, growth and rise of bubbles in solidifying lava and can be expected to occur in all basaltic flows. Degradation of basaltic flows, therefore, will produce vesicular blocks until the erosional level reaches the central vesicle-free zone. In addition, observation of terrestrial basaltic flows has shown that most thin (less than 10 m thick) flows have a regular pattern of orthogonal joints in vertical section in which the spacing of joints increases with depth beneath the flow surface. Using these studies we have performed a preliminary analysis of the Viking lander sites.

Aubele, J. C.↗

Modeling and computation of turbulent free-surface jets

A 2D numerical model for a turbulent surface jet issuing into a finite-depth channel is presented. The model uses an unstructured, finite element procedure. The calculations confirm the experimental observation that scale similarity may be apparent in the mean velocity profile even in the presence of a large encroaching recirculation cell.

Miner, E. W.↗

Modeling Cyclic Phase Change and Energy Storage in Solar Heat Receivers

Numerical results pertaining to cyclic melting and freezing of an encapsulated phase change material (PCM), integrated into a solar heat receiver, have been reported. The cyclic nature of the present melt/freeze problem is relevant to latent heat thermal energy storage (LHTES) systems used to power solar Brayton engines in microgravity environments. Specifically, a physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) project was developed and results compared with available experimental data. Multi-conjugate effects such as the convective fluid flow of a low-Prandtl-number fluid, coupled with thermal conduction in the phase change material, containment tube and working fluid conduit were accounted for in the model. A single-band thermal radiation model was also included to quantify reradiative energy exchange inside the receiver and losses through the aperture. The eutectic LiF-CaF2 was used as the phase change material (PCM) and a mixture of He/Xe was used as the working fluid coolant. A modified version of the computer code HOTTube was used to generate results for comparisons with GTD data for both the subcooled and two-phase regimes. While qualitative trends were in close agreement for the balanced orbit modes, excellent quantitative agreement was observed for steady-state modes.

Hall, Carsie A., III↗

Magnetograph response to canopy-type fields

The response of longitudinal-field magnetographs to magnetic fields which are semi-infinite or confined to a horizontal layer is discussed with respect to the interpretation of solar diffuse fields, observed towards the limb, in terms of magnetic canopy models. Numerical results are presented for several reference solar models and typical 'calibration' curves are shown for the CI 9111 A, Fe I 8688 A, and Ca II 8542 A lines in magnetostatic atmospheres derived from a mean model. A procedure is developed for determining the base heights of magnetic canopies from observations with an uncertainty not exceeding the order of a pressure scale height. Until definitive information regarding atmospheric structure inside flux tubes can be developed from theory or observation, reliable field strengths cannot be derived from the data.

Jones, H. P.↗

Shock waves and shock tubes; Proceedings of the Fifteenth International Symposium, Berkeley, CA, July 28-August 2, 1985

A detailed survey is presented of shock tube experiments, theoretical developments, and applications being carried out worldwide. The discussions explore shock tube physics and the related chemical, physical and biological science and technology. Extensive attention is devoted to shock wave phenomena in dusty gases and other multiphase and heterogeneous systems, including chemically reactive mixtures. Consideration is given to techniques for measuring, visualizing and theoretically modeling flowfield, shock wave and rarefaction wave characteristics. Numerical modeling is explored in terms of the application of computational fluid dynamics techniques to describing flowfields in shock tubes. Shock interactions and propagation, in both solids, fluids, gases and mixed media are investigated, along with the behavior of shocks in condensed matter. Finally, chemical reactions that are initiated as the result of passage of a shock wave are discussed, together with methods of controlling the evolution of laminar separated flows at concave corners on advanced reentry vehicles.

Bershader, D.↗

Impulsive acceleration and scatter-free transport of about 1 MeV per nucleon ions in (He-3)-rich solar particle events

Impulsive solar energetic particle (SEP) events are studied to: (1) describe a distinct class of SEP ion events observed in interplanetary space, and (2) test models of focused transport through detailed comparisons of numerical model prediction with the data. An attempt will also be made to describe the transport and scattering properties of the interplanetary medium during the times these events are observed and to derive source injection profiles in these events. ISEE 3 and Helios 1 magnetic field and plasma data are used to locate the approximate coronal connection points of the spacecraft to organize the particle anisotropy data and to constrain some free parameters in the modeling of flare events.

Mason, G. M.↗

Vibration transmission through rolling element bearings. I - Bearing stiffness formulation. II - System studies

How vibratory motion can be transmitted from the rotating shaft to the casing and other connecting structures in rotating mechanical equipment is addressed here by developing a new mathematical model of precision rolling element bearings. A new grating stiffness matrix is proposed in order to demonstrate a coupling between the shaft bending motion and the flexural motion of the casing plate. It is shown that the translational bearing stiffness coefficients currently used in rotor dynamic models are a small subset of the proposed matrix. The theory is validated by examples, and the proposed bearing formulation is then extended to demonstrate its superiority over existing models in vibration transmission analyses. It is shown that the model can easily be incorporated into analytical or numerical models typically used for dynamic analyses.

Lim, T. C.↗

An observational and theoretical study of highly supercooled altocumulus

The microphysics of two altocumulus clouds sampled at a temperature of about -30 C through the NCAR King Air aircraft and coincident lidar has been characterized. The clouds are structurally similar to stratocumulus and have extensive cloudtop entrainment, a capping temperature inversion, and a dry layer above. Microphysical and radiative properties of both clouds have been modeled numerically. When entrainment effects are incorporated in the model, calculations of droplet concentration and mean diameter profiles compare favorably with the measurements. Radiative transfer calculations indicate that radiative cooling causes sufficient negative buoyancy in cloudtop parcels for producing convective instability and reproducing the observed downdraft velocities. It is shown that entrainment of warmer, drier air near a cloudtop counteracts the radiatively induced negative buoyancy in the downdrafts.

Heymsfield, Andrew J.↗

Particle-gas dynamics in the protoplanetary nebula

In the past year we made significant progress in improving our fundamental understanding of the physics of particle-gas dynamics in the protoplanetary nebula. Having brought our code to a state of fairly robust functionality, we devoted significant effort to optimizing it for running long cases. We optimized the code for vectorization to the extent that it now runs eight times faster than before. The following subject areas are covered: physical improvements to the model; numerical results; Reynolds averaging of fluid equations; and modeling of turbulence and viscosity.

Cuzzi, Jeffrey N.↗

The cosmic gamma-ray background from Type Ia supernovae

We present an improved calculation of the cumulative gamma-ray spectrum of Type Ia supernovae during the history of the universe. We follow Clayton & Ward (1975) in using a few Friedmann models and two simple histories of the average galaxian nucleosynthesis rate, but we improve their calculation by modeling the gamma-ray scattering in detailed numerical models of SN Ia's. The results confirm that near 1 MeV the SN Ia background may dominate, and that it is potentially observable, with high scientific importance. A very accurate measurement of the cosmic background spectrum between 0.1 and 1.0 MeV may reveal the turn-on time and the evolution of the rate of Type Ia supernova nucleosynthesis in the universe.

The, Lih-Sin↗

Analysis of supersonic flows using k-epsilon model and the RPLUS code: Progress towards high speed combustor analysis

The problem of interest here is the analysis of chemically reacting flow inside of supersonic ramjet combustors. Two key parameters need to be determined: mixing/combustion efficiency and kinetic energy efficiency. In order to determine these parameters, several potential loss mechanisms must be modeled correctly, including mixing, shear, turbulence, vorticity, shock waves, heat transfer, fuel injector drag, poor wall pressure integral, and chemical dissociation. This document discusses k-epsilon turbulence models for use in the RPLUS code for numerical modeling of combustor flow fields.

Lee, J.↗

The impact of land-surface wetness heterogeneity on mesoscale heat fluxes

Vertical heat fluxes associated with mesoscale circulations generated by land-surface wetness discontinuities are often stronger than turbulent fluxes, especially in the upper part of the atmospheric planetary boundary layer. As a result, they contribute significantly to the subgrid-scale fluxes in large-scale atmospheric models. Yet they are not considered in these models. To provide some insights into the possible parameterization of these fluxes in large-scale models, a state-of-the-art mesoscale numerical model was used to investigate the relationships between mesoscale heat fluxes and atmospheric and land-surface characteristics that play a key role in the generation of mesoscale circulations. The distribution of land-surface wetness, the wavenumber and the wavelength of the land-surface discontinuities, and the large-scale wind speed have a significant impact on the mesoscale heat fluxes. Empirical functions were derived to characterize the relationships between mesoscale heat fluxes and the spatial distribution of land-surface wetness. The strongest mesoscale heat fluxes were obtained for a wavelength of forcing corresponding approximately to the local Rossby deformation radius. The mesoscale heat fluxes are weakened by large-scale background winds but remain significant even with moderate winds.

Chen, Fei↗

Target & Propagation Models for the FINDER Radar

Finding persons still alive in piles of rubble following an earthquake, a severe storm, or other disaster is a difficult problem. JPL is currently developing a victim detection radar called FINDER (Finding Individuals in Emergency and Response). The subject of this paper is directed toward development of propagation & target models needed for simulation & testing of such a system. These models are both physical (real rubble piles) and numerical. Early results from the numerical modeling phase show spatial and temporal spreading characteristics when signals are passed through a randomly mixed rubble pile.

propagation↗

Splashing Droplets

Current data on droplet breakup is scarce for the sizes and velocities typical of practical applications such as in spray combustion processes and coating processes. While much more representative of practical applications, the small spatial scales and rapid time-scales prevent detailed measurement of the internal fluid dynamics and liquid property gradients produced by impinging upon surfaces. Realized through the extended spatial and temporal scales afforded by a microgravity environment, an improved understanding of drop breakup dynamics is sought to understand and ultimately control the impingement dynamics of droplets upon surfaces in practical situations. The primary objective of this research will be to mark the onset of different 'splashing modes' and to determine their temperature, pressure and angle dependence for impinging droplets representative of practical fluids. In addition, we are modeling the evolution of droplets that do not initially splash but rather undergo a 'fingering' evolution observed on the spreading fluid front and the transformation of these fingers into splashed products. An example of our experimental data is presented below. These images are of Isopar V impacting a mirror-polished surface. They were acquired using a high-speed camera at 1000 frames per second. They show the spreading of a single droplet after impact and ensuing finger instabilities. Normal gravity experimental data such as this will guide low gravity measurements in the 2.2 second drop tower and KC-135 aircraft as available. Presently we are in the process of comparing the experimental data of droplet shape evolution to numerical models, which can also capture the internal fluid dynamics and liquid property gradients such as produced by impingement upon a heated surface. To-date isothermal numerical data has been modeled using direct numerical simulations of representative splashing droplets. The data obtained so far indicates that the present model describes well the droplet wall interactions to a point in time just before splash. Additional information is included in the original extended abstract.

VanderWal, Randall L.↗

Advanced Turbulence Modeling Concepts

The ZCET program developed at NASA Glenn Research Center is to study hydrogen/air injection concepts for aircraft gas turbine engines that meet conventional gas turbine performance levels and provide low levels of harmful NOx emissions. A CFD study for ZCET program has been successfully carried out. It uses the most recently enhanced National combustion code (NCC) to perform CFD simulations for two configurations of hydrogen fuel injectors (GRC- and Sandia-injector). The results can be used to assist experimental studies to provide quick mixing, low emission and high performance fuel injector designs. The work started with the configuration of the single-hole injector. The computational models were taken from the experimental designs. For example, the GRC single-hole injector consists of one air tube (0.78 inches long and 0.265 inches in diameter) and two hydrogen tubes (0.3 inches long and 0.0226 inches in diameter opposed at 180 degree). The hydrogen tubes are located 0.3 inches upstream from the exit of the air element (the inlet location for the combustor). To do the simulation, the single-hole injector is connected to a combustor model (8.16 inches long and 0.5 inches in diameter). The inlet conditions for air and hydrogen elements are defined according to actual experimental designs. Two crossing jets of hydrogen/air are simulated in detail in the injector. The cold flow, reacting flow, flame temperature, combustor pressure and possible flashback phenomena are studied. Two grid resolutions of the numerical model have been adopted. The first computational grid contains 0.52 million elements, the second one contains over 1.3 million elements. The CFD results have shown only about 5% difference between the two grid resolutions. Therefore, the CFD result obtained from the model of 1.3-million grid resolution can be considered as a grid independent numerical solution. Turbulence models built in NCC are consolidated and well tested. They can handle both coarse and fine grids near the wall. They can model the effect of anisotropy of turbulent stresses and the effect of swirling. The chemical reactions of Magnusson model and ILDM method were both used in this study.

Shih, Tsan-Hsing↗

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

Two general types of numerical models for predicting microwave emission from soils are compared-coherent and noncoherent. In the former, radiation in the soil is treated coherently, and the boundary conditions on the electric fields across the layer boundaries are used to calculate the radiation intensity. In the latter, the radiation is assumed to be noncoherent, and the intensities of the radiation are considered directly. The results of the two approaches may be different because of the effects of interference, which can cause the transmitted intensity at the surface (i.e., emissivity) to be sometimes higher and sometimes lower for the coherent case than for the noncoherent case, depending on the relative phases of reflected fields from the lower layers. This coupling between soil layers in the coherent models leads to greater soil moisture sampling depths observed with this type of model, and is the major difference that is found between the two types of models. In noncoherent models, the emissivity is determined by the dielectric constraint at the air/soil interface. The subsequent differences in the results are functions of both the frequency of the radiation being considered and the steepness of the moisture gradient near the surface. The calculations were performed at frequencies of 1.4 and 19.4 GHz and for two sets of soil profiles. Little difference was observed between the models at 19.4 GHz; and only at the lower frequency were differences apparent because of the greater soil moisture sampling depth at this frequency.

Schmugge, T. J.↗

Time-dependent models of single- and double-layer mantle convection

Time-dependent numerical models of two-dimensional convection which simulate mantle convection are presented. Starting from a fluid initially at rest with a purely conductive temperature profile, it is found that for a Rayleigh number of ten million, the model experiences a transient period of double-layer convection which lasts on the order of hundreds of millions of years when scaled to the earth's mantle. It is suggested that transient periods of double-layer convection in the earth's mantle may have existed in the past, whether or not such a period exists today.

Boss, A. P.↗

Computer modeling of inversion layer MOS solar cells and arrays

A two dimensional numerical model of the inversion layer metal insulator semiconductor (IL/MIS) solar cell is proposed by using the finite element method. The two-dimensional current flow in the device is taken into account in this model. The electrostatic potential distribution, the electron concentration distribution, and the hole concentration distribution for different terminal voltages are simulated. The results of simple calculation are presented. The existing problems for this model are addressed. Future work is proposed. The MIS structures are studied and some of the results are reported.

Ho, Fat Duen↗