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 613 records · Page 34

Water cycle and its management for plant habitats at reduced pressures

Experimental and mathematical models were developed for describing and testing temperature and humidity parameters for plant production in bioregenerative life support systems. A factor was included for analyzing systems operating at low (10-101.3 kPa) pressure to reduce gas leakage and structural mass (e.g., inflatable greenhouses for space application). The expected close relationship between temperature and relative humidity was observed, along with the importance of heat exchanger coil temperature and air circulation rate. The presence of plants in closed habitats results in increased water flux through the system. Changes in pressure affect gas diffusion rates and surface boundary layers, and change convective transfer capabilities and water evaporation rates. A consistent observation from studies with plants at reduced pressures is increased evapotranspiration rates, even at constant vapor pressure deficits. This suggests that plant water status is a critical factor for managing low-pressure production systems. The approach suggested should help space mission planners design artificial environments in closed habitats.

NASA Center KSC↗

Numerical simulation of transonic flows in diffusers

Numerical simulations were made of two-dimensional transonic flows in diffusers, including flow separation induced by a shock or adverse pressure gradient. The mass-averaged, time-dependent, compressible Navier-Stokes equations, simplified by the thin-layer approximation, were solved using MacCormack's hybrid method. The eddy-viscosity formulation was described by the Wilcox-Rubesin's two-equation, k-omega model. Detailed comparison of the computed results with measurements showed good agreement in all cases, including one with massive separation induced by a strong shock. The computation correctly predicted the details of a distinct lambda shock pattern, closely duplicating the configuration observed experimentally in spark-schlieren photographs.

Liou, M.-S.↗

On the two-dimensional modelling of horizontal CVD reactors and its limitations

In order to clarify under what conditions two-dimensional theoretical models can realistically account for the actual growth rate distribution in horizontal CVD reactors, a systematic numerical study was performed, and the results were compared with published experimental data for the MOCVD of gallium arsenide. Three models of increasing complexity were studied. It is shown that two-dimensional models can produce realistic predictions only for reactors with large width-to-height (aspect) ratios, that are operated at subcritical Rayleigh number,i.e., under dominant forced flow conditions. But even then, depending on the carrier gas, thermal (Soret) diffusion must be included in the model. Furthermore, velocity corrections for finite aspect ratios must be made, and buoyancy effects can be significant in the entrance region. On the other hand, the results are relatively insensitive to axial diffusion, temperature dependence of transport properties (versus averaged properties), and an assumed Poiseuille velocity profile (versus a rigorous solution of the velocity field). Yet, boundary layer approximations are prone to give unrealistic results at the low Reynolds numbers typical for CVD operations.

Ouazzani, Jalil↗

Modelling of the dynamics of HgCdTe growth by the vertical Bridgman method

The analysis of Kim and Brown (1989), which described the pseudo-steady state of the HgCdTe melt/crystal system grown by the vertical Bridgman method, is extended to the study of the transients during the crystal growth and, in particular, at the beginning of the growth. The time-dependent numerical simulations developed show the development of a diffusion layer in the melt and the direct relationship between this layer and the intensity and the pattern of convection in the melt. It is shown that the spread of the solute layer occurs on a diffusive time scale, with the initial transient to the steady-state growth rate lasting several hours for the slow growth rates that are typical of the HgCdTe growth. The simulations demonstrate all the features observed experimentally by Szofran et al. (1984).

Kim, Do H.↗

The effect of adding roughness and thickness to a transonic axial compressor rotor

The performance deterioration of a high speed axial compressor rotor due to surface roughness and airfoil thickness variations is reported. A 0.025 mm (0.001 in.) thick rough coating with a surface finish of 2.54-3.18 RMS microns (100-125 RMS microinches) is applied to the pressure and suction surface of the rotor blades. Coating both surfaces increases the leading edge thickness by 10% at the hub and 20% at the tip. Application of this coating results in a loss in efficiency of 6 points and a 9% reduction in the pressure ratio across the rotor at an operating condition near the design point. To separate the effect of thickness and roughness, a smooth coating of equal thickness is also applied to the blade. The smooth coating surface finish is 0.254-0.508 RMS microns (10-20 RMS microinches), compared to the bare metal blade surface finish of 0.508 RMS microns (20 RMS microinches). The smooth coating results in approximately half of the performance deterioration found from the rough coating. Both coatings are then applied to different portions of the blade surface to determine which portions of the airfoil are most sensitive to thickness/roughness variations. Aerodynamic performance measurements are presented for a number of coating configurations at 60%, 80%, and 100% of design speed. The results indicate that thickness/roughness over the first 10% of blade chord accounts for virtually all of the observed performance degradation for the smooth coating, compared to about 70% of the observed performance degradation for the rough coating. The performance deterioration is investigated in more detail at design speed using laser anemometer measurements as well as predictions generated by a quasi-3D Navier-Stokes flow solver which includes a surface roughness model. Measurements and analysis are performed on the baseline blade and the full-coverage smooth and rough coatings. The results indicate that coating the blade causes a thickening of the blade boundary layers. The interaction between the rotor passage shock and the thickened suction surface boundary layer then results in an increase in blockage which reduces the diffusion level in the rear half of the blade passage, thus reducing the aerodynamic performance of the rotor.

Suder, Kenneth L.↗

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗

Heat Transfer and Fluid Mechanics Institute, Meeting, 25th, University of California, Davis, Calif., June 21-23, 1976, Proceedings

Major research advances in heat transfer and fluid dynamics are outlined, with particular reference to relevant energy problems. Of significant importance are such topics as synthetic fuels in combustion, turbulence models, combustion modeling, numerical methods for interacting boundary layers, and light-scattering diagnostics for gases. The discussion covers thermal convection, two-phase flow and boiling heat transfer, turbulent flows, combustion, and aerospace heat transfer problems. Other areas discussed include compressible flows, fluid mechanics and drag, and heat exchangers. Featured topics comprise heat and salt transfer in double-diffusive systems, limits of boiling heat transfer in a liquid-filled enclosure, investigation of buoyancy-induced flow stratification in a cylindrical plenum, and digital algorithms for dynamic analysis of a heat exchanger. Individual items are announced in this issue.

Mckillop, A. A.↗

Trends and Techniques for Space Base Electronics

Simulations of various phosphorus and boron diffusions in SOS were completed and a sputtering system, furnaces, and photolithography related equipment were set up. Double layer metal experiments initially utilized wet chemistry techniques. By incorporating ultrasonic etching of the vias, premetal cleaning a modified buffered HF, phosphorus doped vapox, and extended sintering, yields of 98% were obtained using the standard test pattern. A two dimensional modeling program was written for simulating short channel MOSFETs with nonuniform substrate doping. A key simplifying assumption used is that the majority carriers can be represented by a sheet charge at the silicon dioxide silicon interface. Although the program is incomplete, the two dimensional Poisson equation for the potential distribution was achieved. The status of other Z-D MOSFET simulation programs is summarized.

Trotter, J. D.↗

Flow-induced morphological instabilities due to temporally-modulated stagnation-point flow

The influence of periodically-modulated planar stagnation-point flow on the morphological stability of a directionally-solidifying interface is presently considered with a view to the effect of unsteady nonparallel flows on single-crystal growth. The modeling of the system assumes that the viscous boundary layer thickness is much greater than that of the solute boundary layer, and that the modulation frequency is much smaller than the strength of plane stagnation-point flow. The solidifying interface is either stabilized or destabilized depending on the ratio of the period of modulation to the solute-diffusion time.

Merchant, G. J.↗

Phase-field model for isothermal phase transitions in binary alloys

A new phase field model is described which models isothermal phase transitions between ideal binary alloy solution phases. Equations are developed for the temporal and spatial variation of the phase field, which describes the identity of the phase, and of the composition. An asymptotic analysis, as the gradient energy coefficient of the phase field becomes small, was conducted. From the analysis, it is shown that the model recovers classical sharp interface models of this situation when the interfacial layers are thin, and they show how to relate the parameters appearing in the phase field model to material and growth parameters in real systems. Further, three stages of temporal evolution are identified: the first corresponding to interfacial genesis which occurs very rapidly; the second to interfacial motion controlled by the local energy difference across the interface and diffusion; the last taking place on a long time scale in which curvature effects are important and which correspond to Ostwald ripening. The results of the numerical calculations are presented.

Wheeler, A. A.↗

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗

A procedure for computation of fully stalled flows in two-dimensional passages

A procedure is described for computation of incompressible, steady, two-dimensional flows in fully stalled diffusers with plenum exit. The procedure is successful in predicting pressure distributions and patterns to the accuracy of the data. The procedure employs a zonal model; this maintains close connections between the modeling and the physics thereby providing insight into critical aspects of modeling separated flows. The procedure presented is also convenient for computing unstalled flows in passages with turbulent boundary layers for either direct or indirect design problems. Computing times are well within engineering feasibility. The concepts developed can be extended to other classes of separated flows; some of these extensions have already been completed and are referenced.

Woolley, R. L.↗

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗

Modeling Extreme Deformations in Lithium Ion Batteries

A simultaneously coupled modeling approach to study the electrochemical and thermal behavior of lithium-ion batteries under large mechanical deformation has been developed. The thermo-electrochemical pseudo-2D (P2D) battery model is coupled with a mechanical material model. Mechanical, thermal, and electrochemical models are implemented as user-defined sub-routines in the commercial multi-physics code LS-DYNA. The mechanical strain experienced by anode, cathode and separator results in thickness and porosity changes in each layer which in turn influences electrochemical behavior. The evolution of concentration profiles and cell potential are studied under different mechanical loading conditions. Internal short-circuits caused by mechanical deformation and corresponding physical behaviors are also elucidated. We discuss the competing effects of improved transport at higher temperature due to the internal short-circuit versus a drop in the effective ionic conductivity and electrolyte diffusivity due to mechanical deformation.

47 OTHER INSTRUMENTATION↗

Calculations of separation bubbles using boundary-layer-type equations. I. II

Two-dimensional models of separated flows are presented with stream-function boundary equations. A patching procedure is employed to solve for unknowns in viscous-inviscid interactions simultaneously in order to derive strong coupling. The technique features simplified Navier-Stokes equations and is useful for flows around flat plates and more complex geometries such as airfoils and cascades. Sample calculations are provided for separated flows in diffusers and in the neighborhood of trailing edges. Massive bubbles and unclosed bubbles are noted to present difficulties for the models.

Halim, A.↗

Virtual Volumetric Additive Manufacturing (VirtualVAM)

Abstract Tomographic volumetric additive manufacturing (VAM) produces arbitrary 3D geometries by exposure of a rotating volume of photopolymer resin to tomographically‐patterned illumination. This enables high speed, layer‐less printing of parts from a wide range of photopolymers not amenable to layer‐by‐layer processes. Since the entire geometry is produced at once over the course of a few seconds to minutes, molecular diffusion length scales become significant to the printing process. Understanding these molecular reaction and diffusion processes is imperative for advancing VAM to a usable technology. These processes are experimentally very difficult to monitor and measure. Herein, VirtualVAM ‐ a simulation framework for modeling the tomographic VAM process, is developed and experimentally validated. VirtualVAM simulates reaction, diffusion, and heat generation processes over the course of a print with single‐voxel resolution. From a few experimentally‐determined input parameters and a set of images for projection, VirtualVAM is able to generate a large spatio‐temporal data set for any given tomographic VAM print. Using VirtualVAM, a number of experimentally‐unattainable aspects of the VAM process are investigated such as single‐voxel conversion profiles, effect of molecular oxygen, and stopping time determination. VirtualVAM also enables the optimization of exposure patterns to further improve contrast between in‐part and out‐of‐part delivered dose.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mass transport in nanoarray monolithic catalysts: An experimental-theory study

Reducing the mass transfer resistance globally of a catalyst is a key to enhancing the catalytic reaction kinetics and fully utilizing the catalyst activity. Despite the success in tailoring the external mass transfer in the widely studied washcoat monoliths, the internal mass transfer resistance is difficult to be reduced due to the requirement of increasing macroporosity while maintaining high specific surface area and mechanical stability. Therefore, nanostructured array-based monolithic catalysts (nanoarray catalysts) have been developed in the past decade as a promising class of structured catalysts that may complement or substitute washcoat catalysts. This work fundamentally elucidates the enhanced mass transport properties of the nanoarray monolithic catalysts by a combination of experimental measurements and theoretical modeling. Using a low-dimensional model, the relative contributions of resistances were quantified in terms of chemical kinetics, internal and external mass transfers based on a probe model of C2H4 oxidation over the TiO2 supported Pt-based monolithic catalysts. The nanoarray catalysts displayed a lower internal mass transfer resistance than the washcoat counterparts as a result of the high macroporosity and small thickness of nanoarray layers. Finally, the nanoarray configuration provides a new pathway towards designing high-performance monolithic reactors and catalysts with low internal diffusion limitations for various gas phase reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗