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 901 records · Page 50

Structure of a swirl-stabilized spray flame by imaging, laser Doppler velocimetry, and phase Doppler anemometry

Data are presented which describe the mean structure of a steady, swirl-stabilized, kerosene spray flame in the near-injector region of a research furnace. The data presented include ensemble-averaged results of schlieren, luminosity, and extinction imaging, measurement of the gas phase velocity field by laser Doppler velocimetry, and characterization of the condensed phase velocity by phase Doppler anemometry. The results of these studies define six key regions in the flame: the dense spray region; the rich, two-phase, fuel jet; the main air jet; the internal product recirculation zone; the external product recirculation zone; and the gaseous diffusion flame zone. The first five of these regions form a conical mixing layer which prepares the air and fuel for combustion. The air and fuel jets comprise the central portion of this mixing layer and are bounded on either side by the hot product gases of the internal and external recirculation zones. Entrainment of these product gases into the air/fuel streams provides the energy required to evaporate the fuel spray and initiate combustion. Intermittency of the internal recirculation and spray jet flows accounts for unexpected behavior observed in the aerodynamics of the two phases. The data reported herein are part of the database being accumulated on this spray flame for the purpose of detailed comparison with numerical modeling.

Edwards, C. F.↗

Structure of Propagating and Attached Hydrocarbon Flames

Direct numerical simulations with C3-chemistry and radiative heat-loss models have been performed to reveal the internal structure of propagating and attached flames in an axisymmetric fuel jet of methane, ethane, ethylene, acetylene, or propane in air under normal and zero gravity. Observations of the flames were also made at the NASA Glenn 2.2-Second Drop Tower. In computations, the fuel issued into quasi-quiescent air for a fixed mixing time before it was ignited along the centerline at stoichiometry. The edge of the flame propagated through a flammable layer at the laminar flame speed of the stoichiometric fuel-air mixture independent of gravity. For all cases, a peak reactivity spot, i.e., reaction kernel, was formed in the flame base, thereby holding a trailing diffusion flame. The location of the reaction kernel in the attached flames depended inversely on the reactivity. The reaction-kernel correlations between the reactivity and the velocity were developed further using variables related to local Damkahler and Peclet numbers.

Takahashi, Fumiaki↗

Revealing a Pathway for Low-Temperature Recrystallization in Germanium

Thermally activated annealing in semiconductors faces inherent limitations, such as dopant diffusion. Here, a nonthermal pathway is demonstrated for a complete structural restoration in predamaged germanium via ionization-induced recovery. By combining experiments and modeling, this study reveals that the energy transfer of only 2.4 keV nm −1 from incident ions to target electrons can effectively annihilate pre-existing defects and restore the original crystalline structure at room temperature. Moreover, it is revealed that the irradiation-induced crystalline-to-amorphous (c/a) transformation in Ge is reversible, a phenomenon previously considered unattainable without additional thermal energy imposed during irradiation. For partially damaged Ge, the overall damage fraction decreases exponentially with increasing fluence. Surprisingly, the recovery process in preamorphized Ge starts with defect recovery outside the amorphous layer and a shrinkage of the amorphous thickness. After this initial stage, the remaining damage decreases slowly with increasing fluence, but full restoration of the pristine state is not achieved. These differences in recovery are interpreted in the framework of structural differences in the initial defective layers that affect recovery kinetics. This study provides new insights on reversing the c/a transformation in Ge using highly-ionizing irradiation and has broad implications across materials science, radiation damage mitigation, and fabrication of Ge-based devices.

athermal recovery↗

Semiconductor photoelectrochemistry

Semiconductor photoelectrochemical reactions are investigated. A model of the charge transport processes in the semiconductor, based on semiconductor device theory, is presented. It incorporates the nonlinear processes characterizing the diffusion and reaction of charge carriers in the semiconductor. The model is used to study conditions limiting useful energy conversion, specifically the saturation of current flow due to high light intensity. Numerical results describing charge distributions in the semiconductor and its effects on the electrolyte are obtained. Experimental results include: an estimate rate at which a semiconductor photoelectrode is capable of converting electromagnetic energy into chemical energy; the effect of cell temperature on the efficiency; a method for determining the point of zero zeta potential for macroscopic semiconductor samples; a technique using platinized titanium dioxide powders and ultraviolet radiation to produce chlorine, bromine, and iodine from solutions containing their respective ions; the photoelectrochemical properties of a class of layered compounds called transition metal thiophosphates; and a technique used to produce high conversion efficiency from laser radiation to chemical energy.

Buoncristiani, A. M.↗

Numerical prediction of turbulent flow over airfoil sections with a new nonequilibrium turbulence model

A new nonequilibrium turbulence closure model has been developed for computing wall bounded two-dimensional turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present model has been evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated turbulent flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the new nonequilibrium turbulence model accurately captures the history effects of convection and diffusion on turbulence.

Ahmed, S.↗

The Liquid-Water Oscillation in Modeling Boundary-Layer Cumuli with Third-Order Turbulence Closure Models

A hierarchy of third-order turbulence closure models are used to simulate boundary-layer cumuli in this study. An unrealistically strong liquid-water oscillation (LWO) is found in the fully prognostic model, which predicts all third moments. The LWO propagates from cloud base to cloud top with a speed of 1 m/s. The period of the oscillation is about 1000 s. Liquid-water buoyancy terms in the third-moment equations contribute to the LWO. The LWO mainly affects the vertical profiles of cloud fraction, mean liquid-water mixing ratio and the fluxes of liquid-water potential temperature and total water, but has less impact on the vertical profiles of other second-moments and third-moments. In order to minimize the LWO, a moderate large diffusion coefficient and a large turbulent dissipation at its originating level are needed. However, this approach distorts the vertical distributions of cloud fraction and liquid-water mixing ratio. A better approach is to parameterize liquid-water buoyancy more reasonably. A minimally prognostic model, which diagnoses all third moments except for vertical velocity, is shown to produce better results, compared to a fully prognostic model.

Cheng, A.↗

Mineral-Catalysed Formation of Marine NO and N 2 O on the Anoxic Early Earth

Microbial denitrification converts fixed nitrogen species into gases in extant oceans. However, it is unclear how such transformations occurred within the early nitrogen cycle of the Archean. Here 5we present experimental anoxic surface-catalyzed reduction of nitrite and nitrate via green rust and magnetite combined with diffusion and photochemical modeling. We find that in a Fe 2+ -rich marine environment, Fe minerals could have catalyzed abiotic denitrification reactions leading to the formation of nitric oxide (NO) and nitrous oxide (N 2 O). Nitrate did not exhibit reactivity in the presence of either mineral or aqueous Fe2+, however, both minerals induced rapid nitrite reduction 10to NO and N 2 O. While N 2 O escaped into the atmosphere (63% of nitrite-nitrogen, with green rust as catalyst), NO remained associated with precipitates (7%) serving as a potential shuttle to the benthic ocean.The modeling suggests that marine N 2 O emissions would have sustained 0.8-6 ppb of atmospheric N 2 O without a protective ozone layer. Our findings add detail to the as yet incompletely documented Archean nitrogen cycle, implying a globally distributed process driven 15by chemical kinetics similar to those of modern enzymatically mediated conversions.

Steffen Buessecker↗

A Decentralized Approach for Modeling Organized Convection Based on Thermal Populations on Microgrids

Abstract In this study, a spectral model for convective transport is coupled to a thermal population model on a two‐dimensional horizontal “microgrid,” covering the typical gridbox size of general circulation models. The goal is to explore new ways of representing impacts of spatial organization in cumulus cloud fields. The thermals are considered the smallest building block of convection, with thermal life cycle and movement represented through binomial functions. Thermals interact through two simple rules, reflecting pulsating growth and environmental deformation. Long‐lived thermal clusters thus form on the microgrid, exhibiting scale growth and spacing that represent simple forms of spatial organization and memory. Size distributions of cluster number are diagnosed from the microgrid through an online clustering algorithm, and provided as input to a spectral multiplume eddy‐diffusivity mass flux scheme. This yields a decentralized transport system, in that the thermal clusters acting as independent but interacting nodes that carry information about spatial structure. The main objectives of this study are (a) to seek proof of concept of this approach, and (b) to gain insight into impacts of spatial organization on convective transport. Single‐column model experiments demonstrate satisfactory skill in reproducing two observed cases of continental shallow convection. Metrics expressing self‐organization and spatial organization match well with large‐eddy simulation results. We find that in this coupled system, spatial organization impacts convective transport primarily through the scale break in the size distribution of cluster number. The rooting of saturated plumes in the subcloud mixed layer plays a key role in this process.

54 ENVIRONMENTAL SCIENCES↗

Phase mediated dynamics of self-propagating Co/Al nanolaminate reactions

The ignition of sputter deposited nanolaminate foils comprised of alternating Co and Al layers results in rapid, self-propagating formation reactions. The propagating waves present after ignition of 150 nm-thick foils are characterized in movie mode dynamic transmission electron microscopy where these are found to have reaction speeds and wave morphology that vary with bilayer thickness. High speed videography reveals different bilayer thickness-wave character relationships in 750 nm- and 7500 nm-thick Co/Al foils. The reaction speed dependencies on bilayer thickness are calculated for each total thickness by treating the effect of radiation loss as a perturbation from an analytical model described by the difference in the heat of reaction measured in calorimetry and the adiabatic heat of product formation. From this model, an effective activation energy, diffusion constant, and flame temperatures are obtained, which allows for an interpretation of the reaction phase variations with laminate design and their effects on propagating wave morphology.

36 MATERIALS SCIENCE↗

Improving the Performance of the OMPS Limb Profiling Instrument in Monitoring Aerosols in the Post-Pinatubo Stratosphere

The OMPS Limb Profiling (LP) instrument was launched on NOAA's Suomi NPP (SNPP) weather satellite in October 2011. The instrument measures limb-scattered radiation in the ultraviolet, visible and NIR wavelengths with high vertical resolution (~2 km) and relatively dense spatial and temporal sampling (several thousand measurements/day covering the entire sunlit globe). These measurements have very high sensitivity to stratospheric aerosols due to long path of the incoming solar photons through the aerosol layer. However, unlike the recently launched ISS/SAGE III that directly measures extinction of solar/lunar radiation by gases and particles along the line of sight (LOS) of the instrument, information content of OMPS LP measurements is complex. It consists of solar radiation singly scattered and attenuated by aerosols, air, and trace gases along the LOS of the instrument, as well as diffuse upwelling radiation from the lower atmosphere that also gets scattered and attenuated into the LOS. To unscramble this complex signal one needs to have good knowledge of aerosol size distribution (ASD) and their complex refractive index. We will describe our recent efforts in constraining the ASD parameters by analyzing in-situ balloon data from Laramie, Wyoming; aerosol microphysical model data from CARMA; and spectral information provided by OMPS LP. We will also discuss our proposed methodology to refine this information by adjusting the microphysical properties of the aerosols to make data from various occultation and scattering instruments internally consistent. This is a generalized version of the technique employed by AERONET in which solar extinction data are combined with scattered radiation measured by the almucantar technique to determine aerosol particle shape, complex refractive index and size distribution. We suggest that such closure studies are essential to develop confidence in space-based data to validate aerosol microphysical models, such as CARMA, and aerosol chemistry- transport models, such as GOCART, and for estimating radiative forcing due to stratospheric aerosols in both volcanic and quiescent time periods.

Bhartia, P. K.↗

Understanding particle size effect on fast-charging behavior of graphite anode using ultra-thin-layer electrodes

Extreme fast charging (≤ 15 min) of lithium-ion batteries is highly desirable to accelerate mass-market adoption of electric vehicles. However, significant capacity fading, as well as safety issues due to the lithium plating caused by the fast charging rate, limit its implementation. In this study, we investigated the fast-charging capability of graphite materials with various particle sizes. To eliminate the Li + ion concentration gradient effect across the thickness of the electrode, ultra-thin-layer graphite electrodes were developed to investigate the "real" fast-charging capability of graphite at the particle level. Electrochemical assessments as well as microscopic characterizations revealed that smaller particles exhibited superior fast-charging performance, featuring enhanced capacity reversibility, faster charging rate, and less lithium plating under the same fast-charging conditions. It is shown that small-particle graphite (mean radius of 3.3 μm) could withstand a 4C charge (to 80 % state-of-charge) without plating, with minimal plating occurring at 6C. Thicker particles exhibited plating at lower C-rates. Since the experimental data could not directly explain whether intra-particle diffusion limitations or interfacial reaction limitations dominated the plating mechanism, the pseudo-2-dimensional model was used to evaluate the most likely plating mechanism. The model suggested that particle-level diffusion is the dominant mechanism contributing to plating at high rates. Finally, this work provides comprehensive insights into the particle size effects on fast-charging capability, offering a better understanding of fast-charging behavior and valuable guidance for designing optimal electrode architecture for high-rate lithium-ion batteries.

25 ENERGY STORAGE↗

Near-contact diffusion and compensation in extrinsic photoconductors

Carrier diffusion from heavily doped contact regions has been experimentally monitored with resistivity measurements of variable thickness Si:P extrinsic photoconductors (n(+) -n-n(+)). Modeling of free carrier diffusion at the interface of heavily doped contacts and highly resistive bulk at low temperatures predicts extended regions (5-30 microns) of excess carriers in high purity materials so that, in thin device structures, free carrier diffusion profiles from each contact will overlap and determine the resistivity of the device. In this work, a decrease in resistivity of four orders of magnitude was observed in a 5-micron thick structure compared to a 10-micron thick device. The resistivity of an ohmic structure in the thin limit is strongly dependent on the bulk and near-contact compensation, and resistivity measurements can be used as a sensitive measure of compensation at interfaces or in the tail of implanted layers.

Olsen, C. S.↗

Cloud Condensation in Titan's Lower Stratosphere

A 1-D condensation model is developed for the purpose of reproducing ice clouds in Titan's lower stratosphere observed by the Composite Infrared Spectrometer (CIRS) onboard Cassini. Hydrogen cyanide (HCN), cyanoacetylene (HC3N), and ethane (C2H6) vapors are treated as chemically inert gas species that flow from an upper boundary at 500 km to a condensation sink near Titan's tropopause (-45 km). Gas vertical profiles are determined from eddy mixing and a downward flux at the upper boundary. The condensation sink is based upon diffusive growth of the cloud particles and is proportional to the degree of supersaturation in the cloud formation regIOn. Observations of the vapor phase abundances above the condensation levels and the locations and properties of the ice clouds provide constraints on the free parameters in the model. Vapor phase abundances are determined from CIRS mid-IR observations, whereas cloud particle sizes, altitudes, and latitudinal distributions are derived from analyses of CIRS far-IR observations of Titan. Specific cloud constraints include: I) mean particle radii of2-3 J.lm inferred from the V6 506 cm- band of HC3N, 2) latitudinal abundance distributions of condensed nitriles, inferred from a composite emission feature that peaks at 160/cm , and 3) a possible hydrocarbon cloud layer at high latitudes, located near an altitude of 60 km, which peaks between 60 and 80 cm l . Nitrile abundances appear to diminish substantially at high northern latitudes over the time period 2005 to 2010 (northern mid winter to early spring). Use of multiple gas species provides a consistency check on the eddy mixing coefficient profile. The flux at the upper boundary is the net column chemical production from the upper atmosphere and provides a constraint on chemical pathways leading to the production of these compounds. Comparison of the differing lifetimes, vapor phase transport, vapor phase loss rate, and particle sedimentation, sheds light on temporal stability of the clouds.

Romani, Paul N.↗

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗

Navier-Stokes analysis of cold scramjet-afterbody flows

The progress of two efforts in coding solutions of Navier-Stokes equations is summarized. The first effort concerns a 3-D space marching parabolized Navier-Stokes (PNS) code being modified to compute the supersonic mixing flow through an internal/external expansion nozzle with multicomponent gases. The 3-D PNS equations, coupled with a set of species continuity equations, are solved using an implicit finite difference scheme. The completed work is summarized and includes code modifications for four chemical species, computing the flow upstream of the upper cowl for a theoretical air mixture, developing an initial plane solution for the inner nozzle region, and computing the flow inside the nozzle for both a N2/O2 mixture and a Freon-12/Ar mixture, and plotting density-pressure contours for the inner nozzle region. The second effort concerns a full Navier-Stokes code. The species continuity equations account for the diffusion of multiple gases. This 3-D explicit afterbody code has the ability to use high order numerical integration schemes such as the 4th order MacCormack, and the Gottlieb-MacCormack schemes. Changes to the work are listed and include, but are not limited to: (1) internal/external flow capability; (2) new treatments of the cowl wall boundary conditions and relaxed computations around the cowl region and cowl tip; (3) the entering of the thermodynamic and transport properties of Freon-12, Ar, O, and N; (4) modification to the Baldwin-Lomax turbulence model to account for turbulent eddies generated by cowl walls inside and external to the nozzle; and (5) adopting a relaxation formula to account for the turbulence in the mixing shear layer.

Baysal, Oktay↗

Navier-Stokes simulation of flow through a highly contoured subsonic diffuser

The flow through a highly offset subsonic diffuser with cross-sectional profiles that varied from rectangular at the duct entrance to circular at the engine face was numerically simulated. A multizonal approach combined with a two-grid topology was used to represent both the internal and external flowfields, and an implicit, approximately-factored, partially flux-split finite-difference algorithm was used to solve the three-dimensional thin-layer Navier-Stokes equations. The computed static pressures along the inlet wall and total pressures on the engine face were compared with experimental data. In addition, the overall flowfield within the duct was examined in detail. Good agreement is shown between experiment and computations, with the limiting factor being the lack of a reliable turbulence model for internal flow problems.

Chyu, W. J.↗

Experimental Investigation of the Pressure Rise Required for the Incipient Separation of Turbulent Boundary Layers in Two-Dimensional Supersonic Flow

An experimental investigation has been made of turbulent boundary-layer separation associated with compression corners, curved surfaces of various radii, and incident shock waves. The purpose of the investigation was to provide design information, and to define significant physical trends, which would aid in the prediction of turbulent separation for various aerodynamic devices, such as compressor blades, flaps, spoilers, and diffusers. A characteristic change in the longitudinal static-pressure distribution (i.e., a change from a curve with one inflection point to a curve with three inflection points) was employed to detect the occurrence of separation. The effects of Reynolds number (10(exp 6) to 10(exp 7) per foot or l.5 x 10(exp 4) to 7.5 x 10(exp 4) based upon boundary-layer thickness) and Mach number (1.6 to 4.2) on the onset of turbulent boundary-layer separation were investigated. The pressure gradient of the boundary-layer flow ahead of the interaction region was essentially zero. The results show a considerable effect of Mach number on the pressure rise for incipient separation for all configurations. For a curved-surface model, the static pressure-rise ratio required to cause separation varied from about 2.5, at a Mach number of 2 to about 16, at a Mach number of 3.5. A substantial effect of Reynolds number on the pressure rise for incipient separation was observed in the upper Mach number range and in the lower Reynolds number range; namely, the pressure rise required for separation decreased with increasing Reynolds number. For low Mach numbers and high Reynolds numbers, there appeared to be no Reynolds number effect. The effects of Mach number and of Reynolds number were similar for all models. Model shape was also found to be an important variable affecting the onset of separation. Large gains were realized in the pressure-rise ratio with no separation when the radius of curvature of the model surface was increased. At a Mach number of 3.4, for instance, the pressure-rise ratio with no separation increased from about 5 to 15 as a result of an increase in the radius of curvature from approximately 0 to 30 boundary-layer thicknesses.

Kuehn, Donald M.↗

Quantitative Analysis of Planar Laser-Induced Fluorescence Measurements in a Hypersonic Boundary Layer

Several quantitative measurements extracted from nitric oxide (NO) planar laser-induced fluorescence (PLIF) data obtained in a hypersonic boundary layer are reported: (a) off-body NO mole fraction; (b) surface heat flux; and (c) near-wall static temperature. The experimental data was obtained at NASA Langley Research Center’s 31 in. Mach 10 air tunnel. NO was seeded into the flow through a spanwise slot on the surface of the 10 degree half-angle wedge model. An ultraviolet planar laser sheet was positioned perpendicular to the wedge surface, downstream of the seeding slot, to excite six fluorescence transitions. A method for extracting the relative NO mole fraction, based on spatial variations of the J”= 0.5 PLIF signal, is presented. Combined with the principle of mass conservation, the absolute NO mole fraction is determined. These measurements were used to assess CFD diffusion modelling, correct previously reported PLIF thermometry results, and develop methods for NO-PLIF heat transfer measurements.

McDougall, Connor C.↗