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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Effects of high doping levels silicon solar cell performance

The significance of the heavy doping effects (HDE) on the open-circuit voltage of silicon solar cells is assessed. Voltage calculations based on diffusion theory are modified to include the first order features of the HDE. Comparisions of the open-circuit voltage measured for cells of various base resistivities are made with those calculated using the diffusion model with and without the HDE. Results indicate that the observed variation of voltage with base resistivity is predicted by these effects. A maximum efficiency of 19% (AM0) and a voltage of 0.7 volts are calculated for 0.1 omega-cm cells assuming an optimum diffused layer impurity profile.

Godlewski, M. P.↗
NASA NTRS1990

Computation of the unsteady facilitated transport of oxygen in hemoglobin

The transport of a reacting permeant diffusing through a thin membrane is extended to more realistic dissociation models. A new nonlinear analysis of the reaction-diffusion equations, using implicit finite-difference methods and direct block solvers, is used to study the limits of linearized and equilibrium theories. Computed curves of molecular oxygen permeating through hemoglobin solution are used to illustrate higher-order reaction models, the effect of concentration boundary layers at the membrane interfaces, and the transient buildup of oxygen flux.

Davis, Sanford↗
NASA NTRS1980

Turbulence Modeling for Unsteady Transonic Flows

Conditionally sampled, ensemble-averaged velocity measurements, made with a laser velocimeter, were taken in the flowfield over the rear half of an 18% thick circular arc airfoil at zero incidence tested at M = 0.76 and at a Reynolds number based on chord of 11 x 10(exp 6). Data for one cycle of periodic unsteady flow having a reduced frequency f of 0.49 are analyzed. A series of compression waves, which develop in the early stages of the cycle, strengthen and coalesce into a strong shock wave that moves toward the airfoil leading edge. A thick shear layer forms downstream of the shock wave. The kinetic energy and shear stresses increase dramatically, reach a maximum when dissipation and diffusion of the turbulence exceed production, and then decrease substantially. The response lime of the turbulence to the changes brought about by the shock-wave passage upstream depends on the shock-wave strength and position in the boundary layer. The cycle completes itself when the shock wave passes the midchord, weakens, and the shear layer collapses. Remarkably good comparisons are found with computations that employ the time-dependent Reynolds averaged form of the Navier-Stokes equations using an algebraic eddy viscosity model, developed for steady flows.

Marvin, J. G.↗
NASA NTRS2019

Improvements to an Explicit Algebraic Stress Model for Turbulent Jet Mixing Predictions

Modifications to key coefficients in a k − E based explicit algebraic stress model (EASM) are examined with the objective of improving the prediction of turbulent jet flows. The pressure strain coefficient, C2 and the turbulent diffusion coefficients, σk and σE were investigated. For a series of benchmark subsonic jets at heated and unheated conditions, lowering C2 from the default value of 0.36 to 0.10 resulted in a significant improvement in the jet mixing, when compared to experimental data. Changing σk and σE from default values of 1.00 and 1.4489, respectively, to 0.50 and 0.7244, respectively, improved the initial mixing rate, while reducing the farfield mixing rate and the peak turbulent kinetic energy along the centerline. A high-speed mixing layer was also investigated for performance of baseline and modified EASM coefficients, with similar results as for the jet cases. A flat plate boundary layer was briefly examined to determine the effects of changing the coefficients on the turbulent skin friction coefficient. The change to the pressure strain coefficient, C2 = 0.10 is recommended for future EASM calculation of jets flow; however, it is also recommended that the diffusion coefficients remain at their default values.

Turbulent Mixing↗
NASA NTRS1980

Turbulence Modeling for Unsteady Transonic Flows

Conditionally sampled, ensemble-averaged velocity measurements, made with a laser velocimeter, were taken in the flowfield over the rear half of an 18% thick circular arc airfoil at zero incidence tested at M = 0.76 and of a Reynolds number based on chord of 11 x 10(exp 6). Data for one cycle of periodic unsteady flow having a reduced frequency bar-f of 0.49 are analyzed. A series of compression waves, which develop in the early stages of the cycle, strengthen and coalesce into a strong shock wave that moves toward the airfoil leading edge. A thick shear layer forms downstream of the shock wave. The kinetic energy and shear stresses increase dramatically, reach a maximum when dissipation and diffusion of the turbulence exceed production, and then decrease substantially. The response time of the turbulence to the changes brought about by the shock-wave passage upstream depends on the shock-wave strength and position in the boundary layer. The cycle completes itself when the shock wave passes the midchord, weakens, and the shear layer collapses. Remarkably good comparisons are found with computations that employ the time-dependent Reynolds averaged form of the Navier-Stokes equations using an algebraic eddy viscosity model, developed for steady flows.

Marvin, J. G.↗
DOE OSTI2023

Modelling and analysis of salt-convection effect on oxide reduction process for uranium oxides using smoothed particle hydrodynamics

Oxygen ion transport via molten salt continuum is a key kinetic feature of electrochemical reduction process for uranium oxides in molten salt. During the electrolysis, the evolution of oxygen ion concentration field in molten salt continuum appears from porous metallic layer to open bulk salt. Despite the important role affecting process throughput, its effect coupled with controlled salt movement was not adequately investigated. To capture the reactive oxygen ion transport phenomena, a multidimensional reduction model based on smoothed particle hydrodynamics (SPH) was combined with the advection-diffusion model. A pelletized oxide fuel bed capturing interconnected salt continuum is considered as a reactant material geometry, which allows forced convective delivery of molten salt through the reactant materials with a relatively low pressure head. Here, in this study, reduction time of a pelletized-fuel bed was quantified with respect to key design parameters of the reactant material configuration. The simulation results demonstrated the advantage of salt convection which effectively removes oxygen ion products from the reactant materials and consequently accelerates the reduction process. Moreover, the reduction-time correlation of a single pelletized fuel was developed as a function of relevant dimensionless parameters and was extended to the application for a pelletized-fuel bed. Finally, several insights and implications in the use of a pelletized oxide fuel were presented. The developed SPH-based framework enabled multi-physics analysis of governing dynamics in appreciating and screening candidate reactant material configurations in a cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Ion specific, thin film confinement effects on conductivity in polymerized ionic liquids.

Acrylate-based polymerized ionic liquids (PILs) with ammonium (Am) or imidazolium (Im) cations and tetrafluoroborate (BF4) or bis(trifluoromethanesulfonyl)imide (TFSI) anions were synthesized and spin coated onto gold interdigitated electrodes on silica to investigate nanoconfinement effects on ion transport. The film thickness ranged from 23 to 313 nm. A significant reduction of the in-plane conductivity was observed in some PIL thin films with thickness below 100 nm. Specifically, Am BF4 PIL showed the largest conductivity drop (ca. 50% difference between a 22 nm and a 261 nm film) while Im TFSI PIL showed almost no change under confinement. The difference in conductivity drop is discussed in terms of (i) differences in interfacial layer thickness by fitting a two-layer conductivity model and (ii) potential changes in glass transition temperature (T-g) under confinement, which were estimated using the Vogel-Fulcher-Tammann fits of bulk conductivities. Decreasing film thickness also caused the dielectric loss peaks to shift to lower frequency, indicating that the ion diffusion process slowed under confinement.

Zhao, Qiujie↗
DOE OSTI2024

SEI Formation and Lithium-Ion Electrodeposition Dynamics in Lithium Metal Batteries via First-Principles Kinetic Monte Carlo Modeling

The stabilization and enhanced performance of lithium metal batteries (LMBs) depend on the formation and evolution of the Solid Electrolyte Interphase (SEI) layer as a critical component for regulating the Li metal electrodeposition processes. This study employs a first-principles kinetic Monte Carlo (kMC) model to simulate the SEI formation and Li + electrodeposition processes on a lithium metal anode, integrating both the electrochemical electrolyte reduction reactions and the diffusion events giving place to the SEI aggregation processes during battery charge and discharge processes. The model replicates the competitive interactions between organic and inorganic SEI components, emphasizing the influence of the cycling regime. Results indicate that grain boundaries within the SEI facilitate faster lithium-ion transport compared to crystalline regions, crucial for improving the performance and stability of LMBs. The findings underscore the importance of dynamic SEI modeling for further development of next-generation high-energy-density batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract 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, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗
NASA NTRS1989

A model for including thermal conduction in molecular dynamics simulations

A technique is introduced for including thermal conduction in molecular dynamics simulations for solids. A model is developed to allow energy flow between the computational cell and the bulk of the solid when periodic boundary conditions cannot be used. Thermal conduction is achieved by scaling the velocities of atoms in a transitional boundary layer. The scaling factor is obtained from the thermal diffusivity, and the results show good agreement with the solution for a continuous medium at long times. The effects of different temperature and size of the system, and of variations in strength parameter, atomic mass, and thermal diffusivity were investigated. In all cases, no significant change in simulation results has been found.

Wu, Yue↗
NASA NTRS1993

Kinetic simulations of the Kelvin-Helmholtz instability at the magnetopause

Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz instability for a magnetopauselike configuration. The unidirectional magnetic field is essentially transverse to the plasma flow velocity, which is the most unstable case according to linear theory and models the flow dynamics in the subsolar region of the magnetopause for northward interplanetary magnetic field. We recover effects analogous to those found in MHD simulations, including a mode cascade to longer wavelengths. The boundary layer consists of coherent structure and is not well described by a diffusive process. Isolated structures on the order of the ion gyroradius are formed which can cross the boundary in either direction. We describe how the time evolution of these structures represents transport across boundary layers, and we consider the possible connection of these entities to flux transfer events and other structure seen in the low-latitude boundary layer at the Earth's magnetopause as well as to flux ropes commonly observed near the ionopause of Venus. We also discuss the relation of the hybrid calculations to previous MHD simulations and to observations.

Thomas, V. A.↗
NASA NTRS1993

Equatorial ground ice on Mars: Steady-state stability

Current Martian equatorial surface temperatures are too warm for water ice to exist at the surface for any appreciable length of time before subliming into the atmosphere. Subsurface temperatures are generally warmer still and, despite the presence of a diffusive barrier of porous regolith material, it has been shown by Smoluchowski, Clifford and Hillel, and Fanale et al. that buried ground ice will also sublime and be lost to the atmosphere in a relatively short time. We investigate the behavior of this subliming subsurface ice and show that it is possible for ice to maintain at a steady-state depth, where sublimation and diffusive loss to the atmosphere is balanced by resupply from beneath by diffusion and recondensation of either a deeper buried ice deposits or ground water. We examine the behavior of equatorial ground ice with a numercial time-marching molecular diffusion model. In our model we allow for diffusion of water vapor through a porous regolith, variations in diffusivity and porosity with ice content, and recondensation of sublimed water vapor. A regolith containing considerable amounts of ice can still be very porous, allowing water vapor to diffuse up from deeper within the ice layer where temperatures are warmer due to the geothermal gradient. This vapor can then recondense nearer to the surface where ice had previously sublimed and been lost to the atmosphere. As a result we find that ice deposits migrate to find a steady-state depth, which represents a balance between diffusive loss to the atmosphere through the overlying porous regolith and diffusive resupply through a porous icy regolith below. This depth depends primarily on the long-term mean surface temperature and the nature of the geothermal gradient, and is independent of the ice-free porosity and the regolith diffusivity. Only the rate of loss of ground ice depends on diffusive properties.

Mellon, Michael T.↗
DOE OSTI2020

Study of turbulence-induced refraction of lower hybrid waves using synthetic scrape-off layer filaments

Here, turbulence-induced refraction effects to lower hybrid (LH) wave propagation and current drive are studied using synthetic scrape-off layer (SOL) blob/filament fields. A synthetic 3D, field-following, blob turbulence model is implemented in the ray-tracing/Fokker-Planck (RTFP) codes GENRAY/CQL3D. In Alcator C-Mod, the blob field is shown to significantly affect LH ray-trajectories, leading to increased on-axis damping and smoother current profiles. This effect depends on the average blob size and amplitude. In addition, the diffusion of ray-trajectories in phase-space caused by turbulence increases the robustness of the RTFP model. A modified $N_{||}$ launch spectrum, acting as a proxy for parametric decay instability (PDI) effects, is included in simulations with the blob model. A synergy between the modified launch spectrum and turbulence-induced refraction results in synthetic hard x-ray profiles that agree with experiment. Lastly, the blob model is used to predict the effect of SOL turbulence on DIII-D high-field side (HFS) LH launch. Assuming low turbulence amplitude in the HFS SOL (~5%), turbulence-induced refraction is predicted to have little effect on current drive efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Design and Performance of McRas in SCMs and GEOS I/II GCMs

The design of a prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) for general circulation models (GCMs) will be discussed. McRAS distinguishes three types of clouds: (1) convective, (2) stratiform, and (3) boundary-layer types. The convective clouds transform and merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds merge into the stratiform clouds instantly. The cloud condensate converts into precipitation following the auto-conversion equations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, as well as diffuse both horizontally and vertically with a fully interactive cloud-microphysics throughout the life-cycle of the cloud, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III and 5-ARN CART datasets has shown that together with the rest of the model physics, McRAS can simulate the observed temperature, humidity, and precipitation without many systematic errors. The time history and time mean incloud water and ice distribution, fractional cloudiness, cloud optical thickness, origin of precipitation in the convective anvil and towers, and the convective updraft and downdraft velocities and mass fluxes all show a realistic behavior. Performance of McRAS in GEOS 11 GCM shows several satisfactory features but some of the remaining deficiencies suggest need for additional research involving convective triggers and inhibitors, provision for continuously detraining updraft, a realistic scheme for cumulus gravity wave drag, and refinements to physical conditions for ascertaining cloud detrainment level.

Sud, Yogesh C.↗
NASA NTRS1982

User's manual for the REEDM (Rocket Exhaust Effluent Diffusion Model) computer program

The REEDM computer program predicts concentrations, dosages, and depositions downwind from normal and abnormal launches of rocket vehicles at NASA's Kennedy Space Center. The atmospheric dispersion models, cloud-rise models, and other formulas used in the REEDM model are described mathematically Vehicle and source parameters, other pertinent physical properties of the rocket exhaust cloud, and meteorological layering techniques are presented as well as user's instructions for REEDM. Worked example problems are included.

Bjorklund, J. R.↗
DOE OSTI2023

Particle-in-cell modeling of electron beam generated plasma

Plasmas generated using energetic electron beams are well known for their low electron temperature (T e ) and plasma potential, which makes them attractive for atomic-precision plasma processing applications such as atomic layer etch and deposition. A 2-dimensional particle-in-cell model for an electron beam-generated plasma in argon confined by a constant applied magnetic field is described here in this article. Plasma production primarily occurs in the path of the beam electrons in the center of the chamber. The resulting plasma spreads out in the chamber through non-ambipolar diffusion with a short-circuiting effect allowing unequal electron and ion fluxes to different regions of the bounding conductive chamber walls. The cross-field transport of the electrons (and thus the steady-state characteristics of the plasma) are strongly impacted by the magnetic field. T e is anisotropic in the electron beam region, but low and isotropic away from the plasma production zone. The plasma density increases and the plasma becomes more confined near the region of production when the magnetic field strengthens. The magnetic field reduces both electron physical and energy transport perpendicular to the magnetic field. T e is uniform along the magnetic field lines and slowly decreases perpendicular to it. Electrons are less energetic in the sheath regions where the sheath electric field repels and confines the low-energy electrons from the bulk plasma. Even though electron and ion densities are similar in the bulk plasma due to quasi-neutrality, electron and ion fluxes on the grounded chamber walls are unequal at most locations. Electron confinement by the magnetic field weakens with increasing pressure, and the plasma spread out farther from the electron beam region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fabrication of high-performance gas-diffusion-electrode based membrane-electrode assemblies

This paper demonstrates the fabrication and processing steps required to produce high performance fuel cell membrane electrode assemblies (MEAs) based on spray-coated gas-diffusion electrodes (GDEs). It is demonstrated that coating the catalyst layer with a thin layer of ionomer and then hot pressing the GDEs to the membrane is required to achieve comparable catalyst activity and air performance to catalyst-coated-membrane MEAs. We show that there is a critical amount of ionomer required to achieve maximized performance. Using electron microscopy, we show that the combination of the ionomer overlayer and hot-pressing bonds the catalyst layer to the membrane, increasing the interfacial contact area and quality of this interface. We also find that the ionomer overlayer smooths the surface of the GDE and provides increased contact area between the GDE and the membrane. Additionally, we demonstrate that much less ionomer is required for high-performance than has been previously reported. Through model fitting of electrochemical impedance spectroscopy, we show that this improvement in the catalyst layer - membrane interface reduces the effective catalyst layer resistance, which reduces Ohmic losses and increases catalyst utilization.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Modeling Ni redistribution in the hydrogen electrode of solid oxide cells through Ni(OH) 2 diffusion and Ni-YSZ wettability change

Ni redistribution in the hydrogen electrodes of solid oxide cells is an important degradation mechanism. Its driving force is still under debate. This work focuses on the interplay between Ni(OH) 2 diffusion and Ni-YSZ wettability change. With a Ni(OH) 2 diffusion model developed in our previous work, we further employ three models for the Ni-YSZ contact angle to quantify the Ni-YSZ wettability change. The microstructure evolutions in a reconstructed Ni-YSZ electrode are simulated under selected experimental operating conditions. It is shown that the phenomenological model can capture the driving force of Ni spreading/detachment on YSZ surface and the Ni migration reported in experiments can be qualitatively reproduced through the competition between Ni(OH) 2 diffusion and Ni-YSZ wettability change. It is also found that the initial microstructure of the active layer and the microstructure characteristics in the support layer can strongly affect the distribution of steam partial pressure and overpotential, hence the Ni redistribution, in the active layer, which may explain the inconsistencies in experiments. Further, our results show that there are two missing pieces in current theory of Ni redistribution: a mechanism of fast Ni diffusion under humid atmosphere and the physical reason behind the Ni spreading/detachment.

25 ENERGY STORAGE↗