Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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 307 records · Page 17

ML-Based Pebble Power Reconstruction for Pebble Bed Reactor Analysis

Pebble power reconstruction has been explored to complement the conventional homogenized modeling approach in pebble bed reactor (PBR) analysis, as detailed heterogeneous geometry calculations are computationally expensive. The random distribution of pebble fuels within the core challenges the application of conventional pin power reconstruction methods. To address this, we introduce a machine learning approach based on the transformer model, composed of encoder and decoder layers, to estimate the flux and power form functions for reconstructing individual pebble neutron fluxes and powers. The homogeneous neutron flux distribution within each spectral zone (SZ) is obtained from finite element solutions of global diffusion or transport calculations. Verification tests demonstrate that the trained transformer model accurately predicts power form functions over a range of conditions, including variations in pebble enrichment, location, type, SZ size, and burnup. In particular, verification using a three-dimensional PBR benchmark with burned pebbles shows good agreement in heterogeneous pebble power distributions between Griffin and Serpent. These results highlight the potential of applying conventional pin power reconstruction approaches to PBR cores with randomly distributed pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Theory of plasma contactors in ground-based experiments and low Earth orbit

Previous theoretical work on plasma contactors as current collectors has fallen into two categories: collisionless double layer theory (describing space charge limited contactor clouds) and collisional quasineutral theory. Ground based experiments at low current are well explained by double layer theory, but this theory does not scale well to power generation by electrodynamic tethers in space, since very high anode potentials are needed to draw a substantial ambient electron current across the magnetic field in the absence of collisions (or effective collisions due to turbulence). Isotropic quasineutral models of contactor clouds, extending over a region where the effective collision frequency upsilon sub e exceeds the electron cyclotron frequency omega sub ce, have low anode potentials, but would collect very little ambient electron current, much less than the emitted ion current. A new model is presented, for an anisotropic contactor cloud oriented along the magnetic field, with upsilon sub e less than omega sub ce. The electron motion along the magnetic field is nearly collisionless, forming double layers in that direction, while across the magnetic field the electrons diffuse collisionally and the potential profile is determined by quasineutrality. Using a simplified expression for upsilon sub e due to ion acoustic turbulence, an analytic solution has been found for this model, which should be applicable to current collection in space. The anode potential is low and the collected ambient electron current can be several times the emitted ion current.

Gerver, M. J.↗

Numerical Simulation for Droplet Combustion Using Lagrangian Hydrodynamics

A predictive model of spray combustion must incorporate models for the wide variety of physical environments in a practical combustor. In regions where droplets are closely spaced, combustion resembles a diffusion flame; where they are well separated, an envelope or wake flame results. The relative velocity field between the fuel droplets and oxidizer in influences boundary layer development about the droplet, recirculating flow patterns, and droplet shape and stability. A model must encompass these interacting temporal and spatial effects as well as complicated combustor boundaries. The objective of the current work is to develop the triangular gridding method for describing the individual and collective properties of vaporizing and burning fuel droplets.

Fritts, M. J.↗

A simple biosphere model (SiB) for use within general circulation models

A simple realistic biosphere model for calculating the transfer of energy, mass and momentum between the atmosphere and the vegetated surface of the earth has been developed for use in atmospheric general circulation models. The vegetation in each terrestrial model grid is represented by an upper level, representing the perennial canopy of trees and shrubs, and a lower level, representing the annual cover of grasses and other heraceous species. The vegetation morphology and the physical and physiological properties of the vegetation layers determine such properties as: the reflection, transmission, absorption and emission of direct and diffuse radiation; the infiltration, drainage, and storage of the residual rainfall in the soil; and the control over the stomatal functioning. The model, with prescribed vegetation parameters and soil interactive soil moisture, can be used for prediction of the atmospheric circulation and precipitaion fields for short periods of up to a few weeks.

Sellers, P. J.↗

Thermal Pollution Mathematical Model: Verification of One-Dimensional Numerical Model at Lake Keowee - Volume 2

A one dimensional model for studying the thermal dynamics of cooling lakes was developed and verified. The model is essentially a set of partial differential equations which are solved by finite difference methods. The model includes the effects of variation of area with depth, surface heating due to solar radiation absorbed at the upper layer, and internal heating due to the transmission of solar radiation to the sub-surface layers. The exchange of mechanical energy between the lake and the atmosphere is included through the coupling of thermal diffusivity and wind speed. The effects of discharge and intake by power plants are also included. The numerical model was calibrated by applying it to Cayuga Lake. The model was then verified through a long term simulation using Lake Keowee data base. The comparison between measured and predicted vertical temperature profiles for the nine years is good. The physical limnology of Lake Keowee is presented through a set of graphical representations of the measured data base.

Lee, S. S.↗

Parametric studies with an atmospheric diffusion model that assesses toxic fuel hazards due to the ground clouds generated by rocket launches

Parametric studies were made with a multilayer atmospheric diffusion model to place quantitative limits on the uncertainty of predicting ground-level toxic rocket-fuel concentrations. Exhaust distributions in the ground cloud, cloud stabilized geometry, atmospheric coefficients, the effects of exhaust plume afterburning of carbon monoxide CO, assumed surface mixing-layer division in the model, and model sensitivity to different meteorological regimes were studied. Large-scale differences in ground-level predictions are quantitatively described. Cloud alongwind growth for several meteorological conditions is shown to be in error because of incorrect application of previous diffusion theory. In addition, rocket-plume calculations indicate that almost all of the rocket-motor carbon monoxide is afterburned to carbon dioxide CO2, thus reducing toxic hazards due to CO. The afterburning is also shown to have a significant effect on cloud stabilization height and on ground-level concentrations of exhaust products.

Stewart, R. B.↗

Similarity theory of the buoyantly interactive planetary boundary layer with entrainment

A similarity model is developed for the vertical profiles of turbulent flow variables in an entraining turbulent boundary layer of arbitrary buoyant stability. In the general formulation the vertical profiles, internal rotation of the velocity vector, discontinuities or jumps at a capping inversion and bulk aerodynamic coefficients of the boundary layer are given by solutions to a system of ordinary differential equations in the similarity variable. To close the system, a formulation for buoyantly interactive eddy diffusivity in the boundary layer is introduced which recovers Monin-Obukhov similarity near the surface and incorporates a hypothesis accounting for the observed variation of mixing length throughout the boundary layer. The model is tested in simplified versions which depend only on roughness, surface buoyancy, and Coriolis effects by comparison with planetary-boundary-layer wind- and temperature-profile observations, measurements of flat-plate boundary layers in a thermally stratified wind tunnel and observations of profiles of terms in the turbulent kinetic-energy budget of convective planetary boundary layers. On balance, the simplified model reproduced the trend of these various observations and experiments reasonably well, suggesting that the full similarity formulation be pursued further.

Hoffert, M. I.↗

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗

Numerical simulation of a supersonic reacting mixing layer

In order to arrive at physical models that can adequately describe supersonic combustion, and develop accurate and efficient numerical techniques for the solution of such models' governing equations, a computer program has been developed for the study of reacting flows which considers the multicomponent diffusion and convection of important chemical species, as well as their finite state reaction and the interaction of the fluid mechanics and the chemistry that occurred. The code employs a hybrid Chebyshev pseudospectral technique for integration of the models' resulting governing equations; the program is here used to study a spatially developing and reacting mixing layer.

Drummond, J. Philip↗

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

Characterizing pulse reverse-pulse uranium recovery in molten LiCl-KCl with intermittent electrochemical impedance spectroscopy

In this work, EIS was applied to characterize the temporal recovery of uranium from UCl 4 onto a tungsten electrode in molten LiCl-KCl salt between sets of pulse reverse-pulse (PRP) deposition steps. Impedance measurements were fitted to an equivalent circuit (EC) model of the electrode interface using Simplex optimization to calculate the double layer capacitance, charge transfer resistance, Warburg impedance, and CPE-restricted diffusion. Diffusion of uranium towards the electrode surface dominated the surface interactions for the first 80 cycles, followed by a transition to kinetic control for the remaining 520 cycles. A decrease in the Warburg impedance showed depletion of uranium from the bulk melt, which was confirmed to be a 52.5% reduction through ICP-OES. Reduction in the phasance from CPE-restricted diffusion suggested the formation of pores in the uranium deposits, which was confirmed through SEM characterization.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Analysis of epitaxial drift field N on P silicon solar cells

The performance of epitaxial drift field silicon solar cell structures having a variety of impurity profiles was calculated. These structures consist of a uniformly doped P-type substrate layer, and a P-type epitaxial drift field layer with a variety of field strengths. Several N-layer structures were modeled. A four layer solar cell model was used to calculate efficiency, open circuit voltage and short circuit current. The effect on performance of layer thickness, doping level, and diffusion length was determined. The results show that peak initial efficiency of 18.1% occurs for a drift field thickness of about 30 micron with the doping rising from 10 to the 17th power atoms/cu cm at the edge of the depletion region to 10 to the 18th power atoms/cu cm in the substrate. Stronger drift fields (narrow field regions) allowed very high performance (17% efficiency) even after irradiation to 3x10 to the 14th power 1 MeV electrons/sq cm.

Baraona, C. R.↗

Study of diamond film growth and properties

The objective was to study diamond film growth and its properties in order to enhance the laser damage threshold of substrate materials. Calculations were performed to evaluate laser induced thermal stress parameter, R(sub T) of diamond. It is found that diamond has several orders of magnitude higher in value for R(sub T) compared to other materials. Thus, the laser induced damage threshold (LIDT) of diamond is much higher. Diamond films were grown using a microwave plasma enhanced chemical vapor deposition (MPECVD) system at various conditions of gas composition, pressure, temperature, and substrate materials. A 0.5 percent CH4 in H2 at 20 torr were ideal conditions for growing of high quality diamond films on substrates maintained at 900 C. The diamond films were polycrystalline which were characterized by scanning electron microscopy (SEM) and Raman scattering spectroscopy. The top surface of the growing film is always rough due to the facets of polycrystalline film while the back surface of the film replicates the substrate surface. An analytical model based on two dimensional periodic heat flow was developed to calculate the effective in-plane (face parallel) diffusivity of a two layer system. The effective diffusivity of diamond/silicon samples was measured using a laser pulse technique. The thermal conductivity of the films was measured to be 13.5 W/cm K, which is better than that of a type Ia natural diamond. Laser induced damage experiments were performed on bare Si substrates, diamond film coated Si, and diamond film windows. Significant improvements in the LIDT were obtained for diamond film coated Si compared to the bare Si.

Albin, Sacharial↗

The JPL/NASA/TAMU nickel-cadmium battery model development status

A discussion of the development of a fundamental cell model is presented in vugraph format. The nickel oxide layer is described in terms of the electronic conductivity of the oxide layer and proton diffusion through the oxide layer. The kinetic and conductivity expressions for the cadmium electrode were improved. The development process yielded performance predictions that are significantly improved.

Timmerman, Paul↗

Turbulent mixing layers in the interstellar medium of galaxies

It is proposed that turbulent mixing layers are common in the ISM of the Milky Way and selected external galaxies, with many layers per kiloparsec along typical lines of sight. All of the diffuse C I 1550-A background emission and a significant fraction of the diffuse H-alpha background at high latitude can be explained by mixing layers cooling at pressure of about 3000/cu cm K. These models also produce the correct ratio of semiforbidden C IV 15500 III 1663-A emission. Only 10 percent of the disk H-alpha is likely to arise from mixing layers. The observed Galactic absorption-line column densities of C IV, N V, Si IV, and O VI are roughly consistent with mixing-layer models with an intermediate temperature of about 10 exp 5.3 K and depleted abundances.

Slavin, Jonathan D.↗

Vertical Transport Processes for Inert and Scavenged Species: TRACE-A Measurements

The TRACE-A mission of the NASA DC-8 aircraft made a large-scale survey of the tropical and subtropical atmosphere in September and October of 1992. Both In-situ measurements of CO (G. Sachsen NASA Langley) and aerosol size (J. Browell group, NASA Langley) provide excellent data sets with which to constrain vertical transport by planetary boundary layer mixing and deep-cloud cumulus convection. Lidar profiles of aerosol-induced scattering and ozone (also by Bremen) are somewhat require more subtle interpretation as tracers, but the vertical information on layering largely compensates for these complexities. The reason this DC-8 dataset is so useful is that very large areas of biomass burning over Africa and South America provide surface sources of appropriate sizes with which to characterize vertical and horizontal motions; the major limitation of our source description is that biomass burning patterns move considerably every few days, and daily burning inventories are a matter of concurrent, intensive research. We use the Penn State / NCAR MM5 model in an assimilation mode on the synoptic and intercontinental scale, and assess the success it shows in vertical transport descriptions. We find that the general level of emissions suggested by the climatological approach (Will. Has, U. of Montana) appears to be approximately correct, possibly a bit low, for this October, 1992, time period. Vertical transport in planetary boundary layer mixing to 5.5 kin was observed and reproduced in our simulations. Furthermore we find evidence that Blackader "transilient" or matrix-transport scheme is needed, but may require some adaptation in our tracer model: CO seems to exhibit very high values at the top of the planetary boundary layer, a process that stretches the eddy-diffusion parameterization. We will report on progress in improving the deep convective transport of carbon monoxide: the Grail scheme as we used it at 100 kin resolution did not transport enough material to the upper troposphere. We expect to be able to attribute this to either parameterization reasons (inadequacy of this parameterization at the large 100km scale) or other reasons. Nevertheless, the qualitative nature of deep transport by clouds shows up well in the simulations. As for scavengable species, the simulations predict tens of micrograms per standard cubic meter of smoke aerosol in the boundary layer. In a straightforward illustration of our simple bulk-mass scavenging parameterization, to one or two micrograms per standard cubic meter of smoke aerosol in the free troposphere just above the source regions: very high concentrations for the free troposphere. We expect to report on comparisons of these predictions to a variety of observations.

Chatfield, Robert B.↗

Blended Wing Body Systems Studies: Boundary Layer Ingestion Inlets With Active Flow Control

A CFD analysis was performed on a Blended Wing Body (BWB) aircraft with advanced, turbofan engines analyzing various inlet configurations atop the aft end of the aircraft. The results are presented showing that the optimal design for best aircraft fuel efficiency would be a configuration with a partially buried engine, short offset diffuser using active flow control, and a D-shaped inlet duct that partially ingests the boundary layer air in flight. The CFD models showed that if active flow control technology can be satisfactorily developed, it might be able to control the inlet flow distortion to the engine fan face and reduce the powerplant performance losses to an acceptable level. The weight and surface area drag benefits of a partially submerged engine shows that it might offset the penalties of ingesting the low energy boundary layer air. The combined airplane performance of such a design might deliver approximately 5.5% better aircraft fuel efficiency over a conventionally designed, pod-mounted engine.

Geiselhart, Karl A.↗