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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Modeling of the Influence of a Damaged Thermally Grown Oxide (TGO) Layer in an Environmental Barrier Coating System

Environmental barrier coatings (EBC) are an enabling technology for the successful application of ceramic matrix composites (CMCs) in air-breathing gas turbine engines. EBCs are susceptible to a number of failure modes including oxidation/delamination, recessions, chemical attack and dissolution, thermomechanical degradation, erosion, and foreign object impact damage in a combustion environment. Spallation of environmental barrier coating (EBC), induced by a thermally grown oxide (TGO) layer, is a key EBC failure mode. The TGO layer, resulting from steam oxidation, grows either from a silicon bond coat layer (if present) or from a silicon carbide (SiC) based substrate itself. The TGO layer evolves (i.e., thickness increases with time) as water vapor and oxygen gradually diffuse through the EBC, and the EBC spalls off once the TGO layer reaches some critical thickness. The critical thickness of the TGO layer for failure is in the range of 20-30 microns, but it can vary due to exposure temperature, microstructure etc. Current work at NASA Glenn Research Center, under the Revolutionary Tools and Methods (RTM) subproject is aimed at addressing associated failure modes in EBC systems and developing robust analysis tools to aid in the design/analysis of these systems. The objective of the current work is to conduct a sensitivity study to examine the influence of uniformly and non-uniformly grown oxide layers on the associated driving forces leading to spallation of the EBC when subjected to isothermal loading. The effect of damage in the TGO layer on the resulting stress states is also assessed both in uniform and non-uniform TGO layers.

Mital, Subodh K.↗

NASA NTRS1981

A model of temporal variations in ozone density in the Martian atmosphere

Temporal variations of the ozone density profile in the Martian atmosphere at high latitudes are calculated for the course of a Martian year, taking into account seasonal and diurnal variations in temperature, water vapor and solar radiation. Calculations are based on a model including 35 neutral photochemical reactions, and vertical eddy diffusion using a time step of 12 min for the region from the surface to 240 km altitude. Results are found to be in better agreement with Mariner 9 observations of the time and magnitude of the seasonal maximum than previous model calculations. The diurnal variation is predicted to be small near the solstices, with the nighttime ozone density greater than the daytime and the magnitude of the difference dependent on season. Opposite temporal variations are predicted for ozone densities above and below about 25 km, and an ozone density maximum at 35-40 km is obtained. It is suggested that the effects of an aerosol layer may not be important in enhancing predicted ozone concentration, and may even decrease it.

Shimazaki, T.↗

NASA NTRS1999

Global Mapping of Underwater UV Irradiances and DNA-Weighted Exposures using TOMS and SeaWiFS Data Products

The global stratospheric ozone-layer depletion results In an increase in biologically harmful ultraviolet (UV) radiation reaching the surface and penetrating to ecologically significant depths in natural waters. Such an increase can be estimated on a global scale by combining satellite estimates of UV irradiance at the ocean surface from the Total Ozone Mapping Spectrometer (TOMS) satellite instrument with the SeaWIFS satellite ocean-color measurements in the visible spectral region. In this paper we propose a model of seawater optical properties in the UV spectral region based on the Case I water model in the visible range. The inputs to the model are standard monthly SeaWiFS products: chlorophyll concentration and the diffuse attenuation coefficient at 490nm. Penetration of solar UV radiation to different depths in open ocean waters is calculated using the RT (radiative transfer) quasi-single scattering approximation (QSSA). The accuracy of the QSSA approximation in the water is tested using more accurate codes. The sensitivity study of the underwater UV irradiance to atmospheric and oceanic optical properties have shown that the main environmental parameters controlling the absolute levels of the UVB (280-320nm) and DNA-weighted irradiance underwater are: solar-zenith angle, cloud transmittance, water optical properties, and total ozone. Weekly maps of underwater UV irradiance and DNA-weighted exposure are calculated using monthly-mean SeaWiFS chlorophyll and diffuse attenuation coefficient products, daily SeaWiFS cloud fraction data, and the TOMS-derived surface UV irradiance daily maps. The final products include global maps of weekly-average UVB irradiance and DNA-weighted daily exposures at 3m and 10m, and depths where the UVB irradiance and DNA-weighted dose rate at local noon are equal to 10% of their surface values.

Vasilkov, Alexander↗

DOE OSTI2021

Mechanistic study of pH effect on organic solvent nanofiltration using carboxylated covalent organic framework as a modeling and experimental platform

Covalent organic framework (COF) membranes for organic solvent nanofiltration (OSN) have attracted growing interest due to their high stability, well-defined pore size, and controllable pore functionality. In aqueous nanofiltration processes, solution pH has a significant effect on membrane permeability, selectivity, and fouling. However, our understanding of pH effects on OSN processes is still limited. In this work, we present experimental and reactive force field molecular dynamics (ReaxFF MD) simulation results on the effect of pH on OSN using a well-characterized carboxylated COF (C-COF) membrane. Filtration experimental and modeling results both demonstrate that pH changes can affect solvent permeance and solute rejection dramatically. Specifically, methanol permeance decreases with the increase of NaOH and HCl concentration, respectively. When Alcian Blue dye was used as the solute in methanol, the selectivity increased from 23% to 98% when changing the feed pH from 2.2 to 10.1, respectively. While ReaxFF MD results show that C-COF solvated pore size is only marginally affected by the pH of solution, methanol diffusion is affected by the size of methanol-solvated Na + and Cl - ions, which increase with increasing Na + and Cl - concentrations. Solution pH change also affected solute solvation layer thicknesses and C-COF pore charges, contributing to the significant AB dye selectivity increase with pH increase. To summarize, we have shown that dramatic changes in C-COF selectivity are not anomalous in OSN, but are instead controllable by varying feed solution pH to modify solvated solute radii and surface charge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Flue Gas CO 2 Capture via Electrochemically Mediated Amine Regeneration: Desorption Unit Design and Analysis

We analyzed the performance of an electrochemical desorption unit to be implemented in an electrochemically mediated amine regeneration (EMAR) process, which separates CO 2 and regenerates the amine upon modulation of a targeted species that complexes with the amine molecules. Using a simplified one-dimensional boundary layer model, we combined the kinetic expressions of the faradaic reactions with mass transfer effects to analyze the polarization behavior of the flow cell. Mass transfer limitations become dominant at high current densities, leading to increases in energy consumption. An understanding of scaling laws with key parameters, such as the Sherwood number, is required to achieve high power operation with minimal mass transfer losses. Because the EMAR chemistry is sensitive to temperature, heat distribution in the unit is also examined. Our results indicate that the desorption performance increases with addition of a high-temperature heat source. Additionally, incorporation of membrane crossover suggests energy losses because of diffusion and migration of active species. Furthermore, these findings inform on scaling laws for at-scale desorber design and provide a model for scaling up the EMAR system for larger scale operations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Transport and Energetics of Carbon Dioxide in Ionic Liquids at Aqueous Interfaces

A major hurdle in utilizing carbon dioxide (CO 2 ) lies in separating it from industrial flue gas mixtures and finding suitable storage methods that enable its application in various industries. To address this issue, we utilized a combination of molecular dynamics simulations and experiments to investigate the behavior of CO 2 in common room-temperature ionic liquids (RTIL) when in contact with aqueous interfaces. Here, our investigation of RTILs, [EMIM][TFSI] and [OMIM][TFSI], and their interaction with a pure water layer mimics the environment of a previously developed ultrathin enzymatic liquid membrane for CO 2 separation. We analyzed diffusion constants and viscosity, which reveals that CO 2 molecules exhibit faster mobility within the selected ILs compared to what would be predicted solely based on the viscosity of the liquids using the standard Einstein–Stokes relation. Moreover, we calculated the free energy of translocation for various species across the aqueous–IL interface, including CO 2 and HCO 3 – . Free energy profiles demonstrate that CO 2 exhibits a more favorable partitioning behavior in the RTILs compared to that in pure water, while a significant barrier hinders the movement of HCO 3 – from the aqueous layer. Experimental measurement of the CO 2 transport in the RTILs corroborates the model. These findings strongly suggest that hydrophobic RTILs could serve as a promising option for selectively transporting CO 2 from aqueous media and concentrating it as a preliminary step toward storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Modeling the aqueous interface of amorphous TiO 2 using deep potential molecular dynamics

Amorphous titanium dioxide (a-TiO 2 ) is widely used as a coating material in applications such as electrochemistry and self-cleaning surfaces where its interface with water has a central role. However, little is known about the structures of the a-TiO 2 surface and aqueous interface, particularly at the microscopic level. In this work, we construct a model of the a-TiO 2 surface via a cut-melt-and-quench procedure based on molecular dynamics simulations with deep neural network potentials (DPs) trained on density functional theory data. After interfacing the a-TiO 2 surface with water, we investigate the structure and dynamics of the resulting system using a combination of DP-based molecular dynamics (DPMD) and ab initio molecular dynamics (AIMD) simulations. Both AIMD and DPMD simulations reveal that the distribution of water on the a-TiO 2 surface lacks distinct layers normally found at the aqueous interface of crystalline TiO 2 , leading to an ~10 times faster diffusion of water at the interface. Bridging hydroxyls (Ti 2 –O b H) resulting from water dissociation decay several times more slowly than terminal hydroxyls (Ti–O w H) due to fast Ti–O w H 2 → Ti–O w H proton exchange events. These results provide a basis for a detailed understanding of the properties of a-TiO 2 in electrochemical environments. Moreover, the procedure of generating the a-TiO 2 -interface employed here is generally applicable to studying the aqueous interfaces of amorphous metal oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Nanoscale structural correlations in a model cuprate superconductor

Understanding the extent and role of inhomogeneity is a pivotal challenge in the physics of cuprate superconductors. While it is known that structural and electronic inhomogeneity is prevalent in the cuprates, it has proven difficult to disentangle compound-specific features from universally relevant effects. Here, in this study, we combine advanced neutron and x-ray diffuse scattering with numerical modeling to obtain insight into bulk structural correlations in HgBa 2 ⁢ CuO 4+δ . This cuprate exhibits a high optimal transition temperature of nearly 100 K, pristine charge-transport behavior, and a simple average crystal structure without long-range structural instabilities, and is therefore uniquely suited for investigations of intrinsic inhomogeneity. We uncover diffuse reciprocal-space patterns that correspond to prominent nanoscale correlations of atomic displacements perpendicular to the CuO 2 planes. The real-space nature of the correlations is revealed through three-dimensional pair distribution function analysis and complementary numerical refinement. We find that relative displacements of ionic and CuO 2 layers play a crucial role, and that the structural inhomogeneity is not directly caused by the presence of conventional point defects. The observed correlations are therefore intrinsic to HgBa 2 ⁢ CuO 4+δ , and thus likely important for the physics of cuprates more broadly. It is possible that the structural correlations are closely related to the unusual superconducting fluctuations and Mott-localization in these complex oxides. As advances in scattering techniques yield increasingly comprehensive data, the experimental and analysis tools developed here for large volumes of diffuse scattering data can be expected to aid future investigations of a wide range of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1995

Rapid Numerical Simulation of Viscous Axisymmetric Flow Fields

A two-dimensional Navier-Stokes code has been developed for rapid numerical simulation of axisymmetric flow fields, including flow fields with an azimuthal velocity component. The azimuthal-invariant Navier-Stokes equations in a cylindrical coordinate system are mapped to a general body-fitted coordinate system, with the streamwise viscous terms then neglected by applying the thin-layer approximation. Turbulence effects are modeled using an algebraic model, typically the Baldwin-Lomax turbulence model, although a modified Cebeci-Smith model can also be used. The equations are discretized using central finite differences and solved using a multistage Runge-Kutta algorithm with a spatially varying time step and implicit residual smoothing. Results are presented for calculations of supersonic flow over a waisted body-of-revolution, transonic flow through a normal shock wave in a straight circular duct of constant cross sectional area, swirling supersonic (inviscid) flow through a strong shock in a straight radial duct, and swirling subsonic flow in an annular-to-circular diffuser duct. Comparisons between computed and experimental results are in fair to good agreement, demonstrating that the viscous code can be a useful tool for practical engineering design and analysis work.

Tweedt, Daniel L.↗

NASA NTRS1997

Structure of Microgravity Transitional and Pulsed Jet Diffusion Flames

This paper describes results obtained in a study of pulsed gas jet diffusion flames to better characterize the recently observed vortex/flame interactions in microgravity transitional and turbulent diffusion flames, and to improve the understanding of large-scale structures in corresponding normal-gravity flames. In preparation for a space experiment, tests were conducted in the 5.18-Second Zero-Gravity Facility of the NASA Lewis Research Center. Both unpulsed and pulsed laminar flames were studied and numerical modeling of these flames was carried out for data comparison and model validation. In addition, complementary tests for a series of unpulsed flames were conducted on-board the NASA KC-135 research aircraft. The microgravity transitional and turbulent gas-jet diffusion flames have been observed to be dominated by large-scale disturbances, or structures. These structures first appear intermittently in the flame at Reynolds numbers (based on the cold jet injection properties) of about 2100. With increase in injection Reynolds number, the rate of intermittent disturbances increases until the generation becomes continuous at Reynolds numbers of 3000 and higher. The behavior of these structures depends upon the velocity and temperature characteristics of the jet/flame shear layer. These characteristics are different in normal gravity and microgravity.

Bahadori, M. Yousef↗

NASA NTRS2004

Damped Windows for Aircraft Interior Noise Control

Windows are a significant path for structure-borne and air-borne noise transmission into aircraft. To improve the acoustical performance, damped windows were fabricated using two or three layers of plexiglas with transparent viscoelastic damping material sandwiched between the layers. In this paper, numerical and experimental results are used to evaluate the acoustic benefits of damped windows. Tests were performed in the Structural Acoustic Loads and Transmission Facility at NASA Langley Research Center to measure the transmission loss for diffuse acoustic excitation and radiated sound power for point force excitation. Comparisons between uniform and damped plexiglas windows showed increased transmission loss of 6 dB at the first natural frequency, 6 dB at coincidence, and 4.5 dB over a 50 to 4k Hz range. Radiated sound power was reduced up to 7 dB at the lower natural frequencies and 3.7 dB over a 1000 Hz bandwidth. Numerical models are presented for the prediction of radiated sound power for point force excitation and transmission loss for diffuse acoustic excitation. Radiated sound power and transmission loss predictions are in good agreement with experimental data. A parametric study is presented that evaluates the optimum configuration of the damped plexiglas windows for reducing the radiated sound power.

Buehrle, Ralph D.↗

NASA NTRS1981

Transport properties for a mixture of the ablation products C, C2, and C3

The ablation of carbon-phenolic heat shields upon entry into the atmosphere of one of the outer planets leads to the injection of large amounts of C, C2, and C3 into the shock layer. These species must be included in the calculation of transport properties in the shock layer. The kinetic theory of gases has been used to obtain accurate results for the transport properties of monatomic carbon. The Hulburt-Hirschelder potential, the most accurate general purpose atom-atom potential for states with an attractive minimum, was used to represent such states and repulsive states were represented by fitting quantum mechanical potential energy curves with the exponential repulsive potential. These results were orientation averaged according to the peripheral force model to obtain transport collision integrals for the C-C2 and C2-C2 interaction. Results for C3 were obtained by ignoring the presence of the central carbon atom. The thermal conductivity, viscosity, and diffusion coefficients for pure C, C2, and C3, and for mixtures of these gases, were then calculated from 1000 K - 25,000 K.

Biolsi, L.↗

NASA NTRS1977

Characteristics of laminar gas jet diffusion flames under the influence of elevated gravity

Laminar gas jet diffusion flames of hydrogen, methane, ethane, and propane, stabilized on small cylindrical burners, were studied at elevated gravity in a centrifuge to investigate the role of buoyancy in such flames. Photographic observation of the hydrocarbon flames indicated that length and carbon luminosity decreased with increasing buoyancy. These flames eventually separated from the burner rim and finally extinguished as gravity increased. Hydrogen flames, however, did not separate or extinguish up to about 31 times earth normal gravity. Length data are reported and interpreted with the aid of a simplified model which describes the flame as a heated vertical cylindrical sheet to which oxygen is supplied by a free convection boundary layer.

Altenkirch, R. A.↗

DOE OSTI2020

Graphite Lithiation under Fast Charging Conditions: Atomistic Modeling Insights

The charging of lithium ion batteries in a fast and safe manner is critical for promoting the mass adoption of electric vehicles. Li intercalation in graphite electrodes is known to be one of the bottlenecks during the fast charging process. The mechanism of Li diffusion in highly polarized graphite anode at high current rates remains, however, not well understood. In this work, Density Functional Theory (DFT) calculations are used to gain insights into the Li diffusion process in graphite when it is far from equilibrium under fast charging conditions. The effect of uncompensated charges on Li mobility is determined in the highly polarized regions of the anode close to the interfaces. The extra charge was found to increase the interlayer spacing in the diffusion layer and adjacent channels, increasing the diffusivity and promoting the formation of Li clusters. A concerted diffusion mechanism at the edge of high-concentration Li domains is proposed to enhance the diffusion of Li.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗

DOE OSTI2024

Turbulent drifts of impurity ions as an explanation for anomalous radial transport in the far-SOL of DIII-D

Abstract Successful fusion reactor operation relies on minimal core contamination by impurities, otherwise too much power may be radiated and harm performance. This requires reliable predictions of impurity transport from the scrape-off layer (SOL) into the core, beyond the traditional ‘anomalous’ diffusion approach. We report a set of far-SOL tungsten transport simulations that demonstrate the role of turbulent drifts on radial impurity transport. A turbulent plasma background is simulated using the gyrokinetic SOL code Gkeyll. Tungsten ions are followed within the plasma background using only their drifts. We find that tungsten tends to travel radially outwards with velocities between v r = 300–1200 m s −1 primarily due to polarization drift. We also extract an anomalous radial diffusion coefficient that varies from D r anom = 5–20 m 2 s −1 . These results are compared to and agree with previous interpretive modeling results. We also show how the turbulent polarization drift can transport some tungsten ions from the wall inwards with effective pinch velocities up to 10 000 m s −1 . We conclude that turbulent drifts are a likely explanation for historically anomalous radial impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbomachinery

A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved b a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.

Chima, R. V.↗

NASA NTRS2012

Wind-US Results for the AIAA 1st Propulsion Aerodynamics Workshop

This presentation contains Wind-US results presented at the 1st Propulsion Aerodynamics Workshop. The The workshop was organized by the American Institute of Aeronautics and Astronautics, Air Breathing Propulsion Propulsion Systems Integration Technical Committee with the purpose of assessing the accuracy of computational computational fluid dynamics for air breathing propulsion applications. Attendees included representatives from representatives from government, industry, academia, and commercial software companies. Participants were were encouraged to explore and discuss all aspects of the simulation process including the effects of mesh type and mesh type and refinement, solver numerical schemes, and turbulence modeling. The first set of challenge cases involved computing the thrust and discharge coefficients for a series of convergent convergent nozzles for a range of nozzle pressure ratios between 1.4 and 7.0. These configurations included a included a reference axisymmetric nozzle as well as 15deg , 25deg , and 40deg conical nozzles. Participants were also asked also asked to examine the plume shock structure for two cases where the 25deg conical nozzle was bifurcated by a bifurcated by a solid plate. The final test case was a serpentine inlet diffuser with an outlet to inlet area ratio of 1.52 ratio of 1.52 and an offset of 1.34 times the inlet diameter. Boundary layer profiles, wall static pressure, and total and total pressure at downstream rake locations were examined.

Yoder, Dennis↗