Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Simulation”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

LIGHT BRIDGE IN A DEVELOPING ACTIVE REGION. II. NUMERICAL SIMULATION OF FLUX EMERGENCE AND LIGHT BRIDGE FORMATION

Light bridges, the bright structure dividing umbrae in sunspot regions, show various activity events. In Paper I, we reported on an analysis of multi-wavelength observations of a light bridge in a developing active region (AR) and concluded that the activity events are caused by magnetic reconnection driven by magnetconvective evolution. The aim of this second paper is to investigate the detailed magnetic and velocity structures and the formation mechanism of light bridges. For this purpose, we analyze numerical simulation data from a radiative magnetohydrodynamics model of an emerging AR. We find that a weakly magnetized plasma upflow in the near-surface layers of the convection zone is entrained between the emerging magnetic bundles that appear as pores at the solar surface. This convective upflow continuously transports horizontal fields to the surface layer and creates a light bridge structure. Due to the magnetic shear between the horizontal fields of the bridge and the vertical fields of the ambient pores, an elongated cusp-shaped current layer is formed above the bridge, which may be favorable for magnetic reconnection. The striking correspondence between the observational results of Paper I and the numerical results of this paper provides a consistent physical picture of light bridges. The dynamic activity phenomena occur as a natural result of the bridge formation and its convective nature, which has much in common with those of umbral dots and penumbral filaments.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

Numerical Simulations of the XR-5 Hall Thruster for Life Assessment at Different Operating Conditions

NASA's Jet Propulsion Laboratory has been investigating the applicability of Aerojet Rocketdyne's XR-5 thruster, a 4.5 kW class Hall thruster, for deep-space missions. Major considerations for qualifying the XR-5 for deep-space missions are demonstration of a wide throttling envelope and a usable life capability in excess of 10,000 h. Numerical simulations with the 2-D axisymmetric code Hall2De are employed to inform the qualification process by assessing erosion rates at the thruster surfaces in a wide range of throttling conditions without the need for conducting costly endurance testing. In previous work at JPL by Jorns et al., the anomalous collision frequency distribution for 11 different throttling conditions of the XR-5 spanning 0.3-4.5 kW were identified based on probe measurements of the electron temperature in the near plume region. In this paper, we provide estimates for the erosion rates at the channel walls and pole covers for the same 11 conditions. Uncertainties in the plasma measurements and in the anomalous collision frequency distribution are addressed by determining upper and lower bounds of the erosion rates. Results suggest that erosion of the walls only occurs in the last 5% of the acceleration channel and the rate of such erosion decreases as the geometry of the thruster changes in time due to magnetic shielding. A quasi-zero-erosion state is eventually achieved in all the examined throttling conditions. Examination of the results for pole surface erosion and estimated cathode life indicates that the XR-5 propellant throughput capability will exceed 700 kg, which provides 50% margin over the usable throughput capability of 466 kg as already demonstrated in wear testing.

Lopez Ortega, Alejandro↗
DOE OSTI2018

Analysis of Plane-Parallel Electron Beam Propagation in Different Media by Numerical Simulation Methods

Simulation by the Monte Carlo method is widely used to calculate the character of ionizing radiation interaction with substance. A wide variety of programs based on the given method allows users to choose the most suitable package for solving computational problems. In turn, it is important to know exactly restrictions of numerical systems to avoid gross errors. Results of estimation of the feasibility of application of the program PCLab (Computer Laboratory, version 9.9) for numerical simulation of the electron energy distribution absorbed in beryllium, aluminum, gold, and water for industrial, research, and clinical beams are presented. The data obtained using programs ITS and Geant4 being the most popular software packages for solving the given problems and the program PCLab are presented in the graphic form. A comparison and an analysis of the results obtained demonstrate the feasibility of application of the program PCLab for simulation of the absorbed energy distribution and dose of electrons in various materials for energies in the range 1–20 MeV.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Numerical Simulation of the Current Dependence of Emission Spectra of High-Power Pulsed Lasers Based on Separate-Confinement Double Heterostructures

A new physical model is proposed for describing the observed current dependence of the emission spectra of high-power pulse lasers based on double separate-confinement heterostructures. The model takes into account the presence of isotype n–n heterojunctions at the interface between the waveguide and the active region. A numerical simulation of the lasing processes was performed. The resulting analytical current dependence of the 2D concentration of holes in the quantum well made it possible to evaluate for the first time the current dependence of the short-wavelength limit of the lasing spectrum. A good agreement was reached between the calculated and experimental characteristics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1993

Numerical simulations of a pulsed detonation wave augmentation device

We present here the concept of a hybrid engine for Single Stage To Orbit (SSTO) air-breathing hypersonic vehicle. This concept relies on the use of pulsed detonation waves, both for thrust generation and mixing/combustion augmentation. We describe the principles behind the engine concept, which we call the Pulsed Detonation Wave Augmentor (PDWA). We demonstrate the principles of operation for two possible configurations through numerical simulations. We also attempt a first approximation to engine design, and propose various applications.

Cambier, Jean-Luc↗
NASA NTRS1992

Numerical simulation of plasma transport driven by the Io torus

The Rice convection model (RCM) has been modified to a form suitable for Jupiter (RCM-J) to study plasma interchange motion in and near the Io plasma torus. The net result of the interchange is that flux tubes, heavily loaded with torus plasma, are transported outward, to be replaced by tubes containing little low-energy (less than 1 keV) plasma. The process is numerically simulated in terms of time evolution from an initial torus that is longitudinally asymmetric and with gradually decreasing density outward from Io's orbit. In the simulations, the nonlinear stage of the instability characteristically exhibits outreaching fingers of heavily-loaded flux tubes that lengthen at an accelerating rate. The principal finding is that the primary geometrical form of outward transport of torus plasma in Jupiter's magnetosphere is through long, outward-moving fingers of plasma. In the simulations, the fingers mainly form in the active sector of the Io torus (the heavier side of the asymmetric torus), and they are spaced longitudinally roughly 20 deg apart.

Yang, Y. S.↗
NASA NTRS1978

Carbonaceous materials subjected to extreme heating - A comparison of numerical simulation and experiments

The ablation of carbonaceous materials in a hydrogen-helium stream has been simulated using a charring materials ablator computer code. These results are compared with the first ablation data to be obtained from the Ames-NASA Giant Planet Pilot Facility. Test stream diagnostics and ablation effects on convective and radiative heat transfer are discussed since these parameters constitute important input data to the numerical simulation. Graphite ablation was predicted to within 10 to 20%, and carbon-phenolic somewhat less accurately.

Davy, W. C.↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part II: The Effect of Roughness on the Performance of a High-Pressure Turbine Blade

Abstract Turbine blades generally present surface roughness introduced in the manufacturing process or caused by in-service degradation, which can have a significant impact on aero-thermal performance. A better understanding of the fundamental physical mechanisms arising from the interaction between the roughness and the turbine flow at engine-relevant conditions can provide insights for the design of blades with improved efficiency and longer operational life. To this end, a high-fidelity numerical framework combining a well-validated solver for direct numerical simulation and a second-order accurate immersed boundary method is employed to predict roughness-induced aero-thermal effects on an LS89 high-pressure turbine (HPT) blade at engine-relevant conditions. Different amplitudes and distributions of surface roughness are investigated and a reference smooth-blade simulation under the same flow conditions is conducted for comparison. Roughness of increasing amplitude progressively shifts the blade suction side boundary layer transition upstream, producing larger values of the turbulent kinetic energy and higher total wake losses. The on-surface data-capturing capabilities of the numerical framework provide direct measurements of the heat flux and the skin friction coefficient, hence offering quantitative information between the surface topology and engineering-relevant performance parameters. This work may provide a benchmark for future numerical studies of turbomachinery flows with roughness.

Engineering↗
NASA NTRS2000

Direct Numerical Simulation of Automobile Cavity Tones

The Navier Stokes equation is solved computationally by the Dispersion-Relation-Preserving (DRP) scheme for the flow and acoustic fields associated with a laminar boundary layer flow over an automobile door cavity. In this work, the flow Reynolds number is restricted to R(sub delta*) < 3400; the range of Reynolds number for which laminar flow may be maintained. This investigation focuses on two aspects of the problem, namely, the effect of boundary layer thickness on the cavity tone frequency and intensity and the effect of the size of the computation domain on the accuracy of the numerical simulation. It is found that the tone frequency decreases with an increase in boundary layer thickness. When the boundary layer is thicker than a certain critical value, depending on the flow speed, no tone is emitted by the cavity. Computationally, solutions of aeroacoustics problems are known to be sensitive to the size of the computation domain. Numerical experiments indicate that the use of a small domain could result in normal mode type acoustic oscillations in the entire computation domain leading to an increase in tone frequency and intensity. When the computation domain is expanded so that the boundaries are at least one wavelength away from the noise source, the computed tone frequency and intensity are found to be computation domain size independent.

Kurbatskii, Konstantin↗
DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Numerical simulation of the interaction of the plasma sheet with the lobes of the Earth's magnetotail

Codes involving one and two spatial dimensions and three velocity dimensions were used to model the Earth's magnetotail. It was shown that the magnetotail can become inflated as a consequence of low energy plasma convection toward the neutral plane. The computer study exhibits a conversion of both magnetic field energy and of energy supplied by the convection electric field into particle energy. The numerical simulations suggest that much of the magnetotail substorm morphology may be a simple consequence of an increase, followed by a decrease, in the convection electric field, without the requirement of any magnetospheric size scale plasma instability or other disruptive processes. It is also concluded that the presence of the convection electric field and a continuing replenishment of low energy particles in the magnetotail are both necessary for maintenance of the magnetotail.

Swift, D. W.↗
NASA NTRS2001

Detailed numerical simulations of laser cooling processes

We developed a detailed semiclassical numerical code of the forces applied on atoms in optical and magnetic fields to increase the understanding of the different roles that light, atomic collisions, background pressure, and number of particles play in experiments with laser cooled and trapped atoms.

numerical simulations laser cooling processes↗
DOE OSTI2016

The thermal regime of the Campi Flegrei magmatic system reconstructed through 3D numerical simulations

In this paper, we illustrate a quantitative conductive/convective thermal model incorporating a wide range of geophysical, petrological, geological, geochemical and isotopical observations that constrain the thermal evolution and present state of the Campi Flegrei caldera (CFc) magmatic system. The proposed model has been computed on the basis of the current knowledge of: (1) the volcanic and magmatic history of the volcano over the last 44 ka, (2) its underlying crustal structure, and (3) the physical properties of the erupted magmas. 3D numerical simulations of heat conduction and convection within heterogeneous rock/magma materials with evolving heat sources and boundary conditions that simulate magma rise from a deep (≥ 8 km depth) to shallow (2–6 km) reservoirs, magma chamber formation, magma extrusion, caldera collapse, and intra-caldera hydrothermal convection, have been carried out. The evolution of the CFc magmatic system through time has been simulated through different steps related to its changes in terms of depth, location and size of magma reservoirs and their replenishment. The thermal modeling results show that both heat conduction and convection have played an important role in the CFc thermal evolution, although with different timing. Finally, the simulated present heat distribution is in agreement with the measured geothermal profiles (Agip, 1987), reproduces the thermal gradient peaks at the CFc margins in correspondence to the anomalies in surface gradients (Corrado et al., 1998), and suggests temperatures of 700 °C at depth of 4 km in the central portion of the caldera, in agreement with the estimated temperature for the brittle-ductile transition (Hill, 1992).

3D conductive/convective model↗
DOE OSTI2015

Study of coherent structures of turbulence with large wall-normal gradients in thermophysical properties using direct numerical simulation

Forced-convection heat transfer in a heated working fluid at a thermodynamic state near its pseudocritical point is poorly predicted by correlations calibrated with data at subcritical temperatures and pressures. This is suggested to be primarily due to the influence of large wall-normal thermophysical property gradients that develop in proximity of the pseudocritical point on the concentration of coherent turbulence structures near the wall. The physical mechanisms dominating this influence remain poorly understood. In the present study, direct numerical simulation is used to study the development of coherent vortical structures within a turbulent spot under the influence of large wall-normal property gradients. A turbulent spot rather than a fully turbulent boundary layer is used for the study, for the coherent structures of turbulence in a spot tend to be in a more organized state which may allow for more effective identification of cause-and-effect relationships. Large wall-normal gradients in thermophysical properties are created by heating the working fluid which is near the pseudocritical thermodynamic state. It is found that during improved heat transfer, wall-normal gradients in density accelerate the growth of the Kelvin-Helmholtz instability mechanism in the shear layer enveloping low-speed streaks, causing it to roll up into hairpin vortices at a faster rate. It is suggested that this occurs by the baroclinic vorticity generation mechanism which accelerates the streamwise grouping of vorticity during shear layer roll-up. The increased roll-up frequency leads to reduced streamwise spacing between hairpin vortices in wave packets. The density gradients also promote the sinuous instability mode in low-speed streaks. The resulting oscillations in the streaks in the streamwise-spanwise plane lead to locally reduced spanwise spacing between hairpin vortices forming over adjacent low-speed streaks. The reduction in streamwise and spanwise spacing between hairpin vortices causes them to interact more frequently by merging together and by breaking apart into smaller turbulence structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

The Impact of TRMM Data on Mesoscale Numerical Simulation of Super Typhoon Paka

Accurate measurement of the spatial and temporal variations of tropical rainfall around the globe had remained as a critical problem in meteorology until the recent launch of the Tropical Rainfall Measuring Mission (TRMM). TRMM offers a unique opportunity to improve understanding of tropical meteorology and also offers a great opportunity to evaluate the impact of rainfall data on tropical weather forecasts. This study assesses the impact of TRMM Microwave Imager (TMI) derived surface rainfall data on the numerical simulation of Super Typhoon Paka (1997).

Pu, Zhao-Xia↗
NASA NTRS1987

Numerical simulation of three dimensional transonic flows

The three-dimensional flow over a projectile has been computed using an implicit, approximately factored, partially flux-split algorithm. A simple composite grid scheme has been developed in which a single grid is partitioned into a series of smaller grids for applications which require an external large memory device such as the SSD of the CRAY X-MP/48, or multitasking. The accuracy and stability of the composite grid scheme has been tested by numerically simulating the flow over an ellipsoid at angle of attack and comparing the solution with a single grid solution. The flowfield over a projectile at M = 0.96 and 4 deg angle-of-attack has been computed using a fine grid, and compared with experiment.

Sahu, Jubaraj↗
NASA NTRS1993

Development of a computational testbed for numerical simulation of combustion instability

A synergistic hierarchy of analytical and computational fluid dynamic techniques is used to analyze three-dimensional combustion instabilities in liquid rocket engines. A mixed finite difference/spectral procedure is employed to study the effects of a distributed vaporization zone on standing and spinning instability modes within the chamber. Droplet atomization and vaporization are treated by a variety of classical models found in the literature. A multi-zone, linearized analytical solution is used to validate the accuracy of the numerical simulations at small amplitudes for a distributed vaporization region. This comparison indicates excellent amplitude and phase agreement under both stable and unstable operating conditions when amplitudes are small and proper grid resolution is used. As amplitudes get larger, expected nonlinearities are observed. The effect of liquid droplet temperature fluctuations was found to be of critical importance in driving the instabilities of the combustion chamber.

Grenda, Jeffrey↗
DOE OSTI2017

Numerical simulation of a full-loop circulating fluidized bed under different operating conditions

Both experimental and computational studies of the fluidization of high-density polyethylene (HDPE) particles in a small-scale full-loop circulating fluidized bed are conducted. Experimental measurements of pressure drop are taken at different locations along the bed. The solids circulation rate is measured with an advanced Particle Image Velocimetry (PIV) technique. The bed height of the quasi-static region in the standpipe is also measured. Comparative numerical simulations are performed with a Computational Fluid Dynamics solver utilizing a Discrete Element Method (CFD-DEM). This paper reports a detailed and direct comparison between CFD-DEM results and experimental data for realistic gas-solid fluidization in a full-loop circulating fluidized bed system. The comparison reveals good agreement with respect to system component pressure drop and inventory height in the standpipe. In addition, the effect of different drag laws applied within the CFD simulation is examined and compared with experimental results.

01 COAL, LIGNITE, AND PEAT↗