Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Chemistry”

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 397 records · Page 22

X-ray Mapping of Terrestrial and Extraterrestrial Materials Using the Electron Microprobe

Lunar samples returned from the Apollo program motivated development of the Bence-Albee algorithm for the rapid and accurate analysis of lunar materials, and established interlaboratory comparability through its common use. In the analysis of mineral and rock fragments it became necessary to combine micro- and macroscopic analysis by coupling electron-probe microanalysis (EPMA) with automated stage point counting. A coarse grid that included several thousand points was used, and initially wavelength-dispersive (WDS) and later energydispersive (EDS) data were acquired at discrete stage points using approx. 5 sec count times. A approx 50 micrometer beam diameter was used for WDS and up to 500 micrometer beam diameter for EDS analysis. Average analyses of discretely sampled phases were coupled with the point count data to calculate the bulk composition using matrix algebra. Use of a defocused beam resulted in a contribution from multiple phases to each analytical point, and the analytical data were deconvolved relative to end-member phase chemistry on the fly. Impressive agreement was obtained between WDS and EDS measurements as well as comparison with bulk chemistry obtained by other methods. In the 30 years since these methods were developed, significant improvements in EPMA automation and computer processing have taken place. Digital beam control allows routine collection of x-ray maps by EDS, and stage mapping for WDS is conducted continuously at slew speed and incrementally by sampling at discrete points. Digital pulse processing in EDS systems has significantly increased the throughput for EDS mapping, and the ongoing development of Si-drift detector systems promises mapping capabilities rivaling WDS systems. Spectrum imaging allows a data cube of EDS spectra to be acquired and sophisticated processing of the original data is possible using matrix algebra techniques. The study of lunar and meteoritic materials includes the need to conveniently: (1) Characterize the sample at microscopic and macroscopic scales with relatively high sensitivity, (2) Determine the modal abundance of minerals, and (3) Identify and relocate discrete features of interest in terms of size and chemistry. The coupled substitution of cations in minerals can result in significant variation in mineral chemistry, but at similar average Z, leading to poor backscattered-electron (BSE) contrast discrimination of mineralogy. It is necessary to discriminate phase chemistry at both the trace element level and the major element level. To date, the WDS of microprobe systems is preferred for mapping due to high throughput and the ability to obtain the necessary intensity to discriminate phases at both trace and major element concentrations. It is desirable to produce fully quantitative compositional maps of geological materials, which requires the acquisition of k-ratio maps that are background and dead-time corrected, and which have been corrected by phi(delta z> or an equivalent algorithm at each pixel. To date, turnkey systems do not allow the acquisition of k-ratio maps and the rigorous correction in this manner. X-ray maps of a chondrule from the Ourique meteorite, and a comb-layered xenolith from the San Francisco volcanic field, have been analyzed and processed to extract phase information. The Ourique meteorite presents a challenge due to relatively low BSE contrast, and has been studied using spectrum imaging. X-ray maps for Si, Mg, and FeK(alpha) were used to produce RGB images. The xenolith sample contains sector-zoned augite, olivine, plagioclase, and basaltic glass. X-ray maps were processed using Lispix and ImageJ software to produce mineral phase maps. The x-ray maps for Mg, Ca, and Ti were used with traceback to generate binary images that were converted to RGB images. These approaches are successful in discriminating phases, but it is desirable to achieve the methods that were used on lunar samples 30 years ago on current microprobe systems. Curnt research includes x-ray mapping analysis of the Dalgety Downs chondrite by micro x-ray fluorescence and spectrum imaging, in collaboration with Kenny Witherspoon of IXRF Systems and Dale Newbury of NIST.

Carpenter, P.↗

A new technique for low-to-high altitude predictions of ablative hypersonic flowfields

A parabolized Navier-Stokes (PNS) technique for predicting chemically reacting reentry flows over a wide range of altitudes extending from zero to 259 kft is presented. This coupled-chemistry PNS scheme can use either a strongly implicit solution approach or a two-step solution approach and is used to study a hypersonic flow over various configurations over a Mach-number range of 6.89 to 20. It is demonstrated that the two-step solution provides a combination of accuracy, stability, and computational efficiency. Nonablating low-altitude calculations show that the flowfield is in chemical equilibrium only in the near-wall region, while finite-rate chemistry effects are significant over a major portion of the shock layer, and the Blottner-type linearization of the production terms is unstable for low-altitude high-Reynolds-number conditions.

Bhutta, Bilal A.↗

CFD validation of subsonic turbulent planar shear layers

The primary objective of the present study is to assess the limitations and capabilities of RPLUS in predicting the entrainment, mixing, and burning characteristics of a high subsonic planar shear layer; these computations are performed in conjunction with an experiment being conducted at NASA Lewis Research Center. Turbulence is modeled by a two-equation k-epsilon closure modified for compressibility effects. Combustion is modeled by two different models: one in which a finite-rate laminar chemistry model is used for H2-O2 oxidation, and the other based on a composition joint pdf approach to account for the turbulence-chemistry interaction effects. Between upwind and centered differencing schemes, minimal differences are observed for the flow properties of a nonreacting shear layer case studied. In comparison with experimental data, computed growth rates are underpredicted in both reacting and nonreacting cases of a shear layer, but more severely for the reacting flows. However, the qualitative agreement between computation and experiment is reasonably good for the other observable characteristics, and indicative of potential for further improvements.

Lai, H. T.↗

Compositional Effects on Nickel-Base Superalloy Single Crystal Microstructures

Fourteen nickel-base superalloy single crystals containing 0 to 5 wt% chromium (Cr), 0 to 11 wt% cobalt (Co), 6 to 12 wt% molybdenum (Mo), 0 to 4 wt% rhenium (Re), and fixed amounts of aluminum (Al) and tantalum (Ta) were examined to determine the effect of bulk composition on basic microstructural parameters, including gamma' solvus, gamma' volume fraction, volume fraction of topologically close-packed (TCP) phases, phase chemistries, and gamma - gamma'. lattice mismatch. Regression models were developed to describe the influence of bulk alloy composition on the microstructural parameters and were compared to predictions by a commercially available software tool that used computational thermodynamics. Co produced the largest change in gamma' solvus over the wide compositional range used in this study, and Mo produced the largest effect on the gamma lattice parameter and the gamma - gamma' lattice mismatch over its compositional range, although Re had a very potent influence on all microstructural parameters investigated. Changing the Cr, Co, Mo, and Re contents in the bulk alloy had a significant impact on their concentrations in the gamma matrix and, to a smaller extent, in the gamma' phase. The gamma phase chemistries exhibited strong temperature dependencies that were influenced by the gamma and gamma' volume fractions. A computational thermodynamic modeling tool significantly underpredicted gamma' solvus temperatures and grossly overpredicted the amount of TCP phase at 982 C. Furthermore, the predictions by the software tool for the gamma - gamma' lattice mismatch were typically of the wrong sign and magnitude, but predictions could be improved if TCP formation was suspended within the software program. However, the statistical regression models provided excellent estimations of the microstructural parameters based on bulk alloy composition, thereby demonstrating their usefulness.

MacKay, Rebecca A.↗

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

An Improved Retrieval of Tropospheric Nitrogen Dioxide from GOME

We present a retrieval of tropospheric nitrogen dioxide (NO2) columns from the Global Ozone Monitoring Experiment (GOME) satellite instrument that improves in several ways over previous retrievals, especially in the accounting of Rayleigh and cloud scattering. Slant columns, which are directly fitted without low-pass filtering or spectral smoothing, are corrected for an artificial offset likely induced by spectral structure on the diffuser plate of the GOME instrument. The stratospheric column is determined from NO2 columns over the remote Pacific Ocean to minimize contamination from tropospheric NO2. The air mass factor (AMF) used to convert slant columns to vertical columns is calculated from the integral of the relative vertical NO2 distribution from a global 3-D model of tropospheric chemistry driven by assimilated meteorological data (Global Earth Observing System (GEOS)-CHEM), weighted by altitude dependent scattering weights computed with a radiative transfer model (Linearized Discrete Ordinate Radiative Transfer), using local surface albedos determined from GOME observations at NO2 wavelengths. The AMF calculation accounts for cloud scattering using cloud fraction, cloud top pressure, and cloud optical thickness from a cloud retrieval algorithm (GOME Cloud Retrieval Algorithm). Over continental regions with high surface emissions, clouds decrease the AMT by 20- 30% relative to clear sky. GOME is almost twice as sensitive to tropospheric NO2 columns over ocean than over land. Comparison of the retrieved tropospheric NO2 columns for July 1996 with GEOS-CHEM values tests both the retrieval and the nitrogen oxide radical (NOx) emissions inventories used in GEOS-CHEM. Retrieved tropospheric NO2 columns over the United States, where NOx emissions are particularly well known, are within 18% of GEOS-CHEM columns and are strongly spatially correlated (r = 0.78, n = 288, p less than 0.005). Retrieved columns show more NO2 than GEOS-CHEM columns over the Transvaal region of South Africa and industrial regions of the northeast United States and Europe. They are lower over Houston, India, eastern Asia, and the biomass burning region of central Africa, possibly because of biases from absorbing aerosols.

Martin, Randall V.↗

Atmospheric Teleconnection over Eurasia Induced by Aerosol Radiative Forcing During Boreal Spring

The direct effects of aerosols on global and regional climate during boreal spring are investigated based on simulations using the NASA Global Modeling and Assimilation Office (GMAO) finite-volume general circulation model (fvGCM) with Microphyics of clouds in Relaxed Arakawa Schubert Scheme (McRAS). The aerosol loading are prescribed from three-dimensional monthly distribution of tropospheric aerosols viz., sulfate, black carbon, organic carbon, soil dust, and sea salt from output of the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The aerosol extinction coefficient, single scattering albedo, and asymmetric factor are computed as wavelength-dependent radiative forcing in the radiative transfer scheme of the fvGCM, and as a function of the aerosol loading and ambient relative humidity. We find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excites a planetary scale teleconnection pattern in sea level pressure, temperature and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, i.e., South Asia, East Asia, and northern and western Africa. Additionally, atmospheric heating is found in regions with large loading of dust (over Northern Africa, and Middle East), and black carbon (over South-East Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east-west dipole anomaly with strong cooling over the Caspian Sea, and warming over central and northeastern Asia, where aerosol concentration are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations.

Kim, Maeng-Ki↗

Rare Gases Inserted into Biological Building Blocks: A Theoretical Study of Glycine - Rg Compounds (Rg-Xe, Kr, Ar)

Compounds formed by insertion of rare-gas atoms (Xe, Kr, and Ar) into glycine molecule are investigated using accurate ab initio computational methods. Identification of such insertion compounds may open new frontiers in the field of rare-gas chemistry, such as possible existence of biological molecules that include chemically bound rare gas atoms. The most stable glycine-Rg configuration is found to correspond to insertion of Rg atoms into the 0-H bond of glycine. These NH2CH2COORgH compounds are metastable , but separated by sizable potential barriers from the Rg + glycine dissociation products. Preliminary calculations show that NH2CH2COOXeH compound is energetically stable with respect to another (3-body) dissociation channel (NH2CH2COO + Rg + H), while the corresponding Ar species is not stable in this respect. The compound with the inserted Kr is a borderline case, with the 3-body dissociation products being close in energy to the NH2CH2COOKrH minimum.

Chaban, Galina M.↗

Simulation of the Aerosol Size Distribution Using a Neural Network Surrogate for the Modal Aerosol Module (MAM7)

One objective of atmospheric simulations is to quantify the distribution of aerosols and their properties. Accurate parameterizations of the processes governing aerosol mass, particle number, and particle size distribution are important for predicting the Earth’s net radiative balance and aerosol-cloud interactions. The Modal Aerosol Module (MAM7) is a two-moment aerosol model that simulates mass, number, and size distribution of seven modes comprised of internally mixed aerosol species. The two-moment scheme adds significant computational expense but allows for the prediction of varying particle size distribution relative to the bulk method which predicts only total mass. In this work, we developed a neural network surrogate model for MAM7 (MAMnet) to predict the aerosol number concentration in NASA’s Global Earth Observing System (GEOS) without adding prohibitive computational expense. MAMnet, can be driven by output from a single moment, mass-based, aerosol scheme (Goddard Chemistry Aerosol and Radiation model (GOCART)) or from reanalysis products (Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2)). MAMnet was trained using number concentrations from a 5-year GEOS/MAM7 simulation at 1-degree horizontal resolution and using the total mass calculated across modes as inputs, as well as temperature and air density. The model architecture for MAMnet was based on AlexNet, the 2012 winner of the ImageNet Large Scale Visual Recognition Challenge. While some modifications were necessary to accommodate our problem, important aspects of the network were preserved. MAMnet was able to reproduce zonal dynamics and spatial distributions of the aerosol number concentration however predictability in the upper troposphere was poor.

Katherine H Breen↗

Insights into the general circulation of the lower stratosphere from TOMS

Total ozone is controlled by dynamical advection as well as chemistry. At least for day-to-day variations, dynamical processes appear to be in control of total ozone. There also appears to be good evidence that seasonal and secular changes in total ozone are also dynamically controlled. For example, the zonal mean changes in total ozone in the two hemispheres in spring appear to be quite different. The TOMS total ozone data suggest a south polar spring upwelling while the Northern Hemisphere shows a clear downwelling in the same period. Radiative transfer computations support this conclusion. The secular changes in total ozone over the South Pole in spring indicate a change in dynamics rather than chemistry.

Schoeberl, Mark R.↗

A sensitivity analysis of volcanic aerosol dispersion in the stratosphere

A computer sensitivity analysis was performed to determine the uncertainties involved in the calculation of volcanic aerosol dispersion in the stratosphere using a 2 dimensional model. The Fuego volcanic event of 1974 was used. Aerosol dispersion processes that were included are: transport, sedimentation, gas phase sulfur chemistry, and aerosol growth. Calculated uncertainties are established from variations in the stratospheric aerosol layer decay times at 37 latitude for each dispersion process. Model profiles are also compared with lidar measurements. Results of the computer study are quite sensitive (factor of 2) to the assumed volcanic aerosol source function and the large variations in the parameterized transport between 15 and 20 km at subtropical latitudes. Sedimentation effects are uncertain by up to a factor of 1.5 because of the lack of aerosol size distribution data. The aerosol chemistry and growth, assuming that the stated mechanisms are correct, are essentially complete in several months after the eruption and cannot explain the differences between measured and modeled results.

Butler, C. F.↗

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.

Yungster, Shaye↗

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.

Yungster, Shaye↗

Overview of the NCC

A multi-disciplinary design/analysis tool for combustion systems is critical for optimizing the low-emission, high-performance combustor design process. Based on discussions between then NASA Lewis Research Center and the jet engine companies, an industry-government team was formed in early 1995 to develop the National Combustion Code (NCC), which is an integrated system of computer codes for the design and analysis of combustion systems. NCC has advanced features that address the need to meet designer's requirements such as "assured accuracy", "fast turnaround", and "acceptable cost". The NCC development team is comprised of Allison Engine Company (Allison), CFD Research Corporation (CFDRC), GE Aircraft Engines (GEAE), NASA Glenn Research Center (LeRC), and Pratt & Whitney (P&W). The "unstructured mesh" capability and "parallel computing" are fundamental features of NCC from its inception. The NCC system is composed of a set of "elements" which includes grid generator, main flow solver, turbulence module, turbulence and chemistry interaction module, chemistry module, spray module, radiation heat transfer module, data visualization module, and a post-processor for evaluating engine performance parameters. Each element may have contributions from several team members. Such a multi-source multi-element system needs to be integrated in a way that facilitates inter-module data communication, flexibility in module selection, and ease of integration. The development of the NCC beta version was essentially completed in June 1998. Technical details of the NCC elements are given in the Reference List. Elements such as the baseline flow solver, turbulence module, and the chemistry module, have been extensively validated; and their parallel performance on large-scale parallel systems has been evaluated and optimized. However the scalar PDF module and the Spray module, as well as their coupling with the baseline flow solver, were developed in a small-scale distributed computing environment. As a result, the validation of the NCC beta version as a whole was quite limited. Current effort has been focused on the validation of the integrated code and the evaluation/optimization of its overall performance on large-scale parallel systems.

Liu, Nan-Suey↗

Modeling Dust Mineralogical Composition: Sensitivity to Soil Mineralogy Atlases and Their Expected Climate Impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ∼ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ∼ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

Soil dust↗

Variability in Modelled Airborne Dust Mineralogy Derived from Global Soil Composition Uncertainties

Dust aerosols consist on a variety of minerals with different physic and chemical properties. As such, they interact with short and long wave radiation, potentially form clouds, act as nutrients modulating biogeochemical cycles, or influence atmospheric chemistry, differently. Most current state-of-the-art Earth System Models (ESMs) neglect the complexity in dust composition, mainly due to computational constraints, but also to the existing uncertainties in the size resolved composition of parent soils, the resulting distribution of minerals in airborne dust, and the scarcity of observations to constrain them. Within this work, we assess the variability of global dust composition due to uncertainties in the characterization of the parent soil mineralogy. To that end, we consider two available global soil mineralogy atlases, developed by Claquin et al. (1999) –C1999- and Journet et al. (2014) –J2014-, which represent respectively 8 and 12 relevant minerals for climate (namely: illite, smectite/montmorillonite, kaolinite, calcite, gypsum, hematite, quartz, and feldspars in C1999, and those plus chlorite, vermiculite, goethite, and mica in J2014). Thanks to a recently developed feature of the MONARCH atmospheric-chemistry model, we are able to explicitly resolve the minerals’ atmospheric cycle. Therefore, we define two global experiments to assess changes on airborne dust composition attributed to the soil mineralogy assumptions and provide a measure of their variability. We also perform a preliminary evaluation of the global mineralogy results against available observations of mineral fractions in surface dust concentration. Our results will inform the climate modelling community about the potential variability in dust composition, an aspect that will gain relevance as ESMs continue growing in complexity and new processes to better characterize aerosols’ forcing or biogeochemical cycles are added. Further observational constraints, such as those that will derive from the EMIT NASA mission on soil composition or the FRAGMENT experimental campaigns on airborne dust characterization, will be key in the near future to improve our understanding of the impact of dust mineralogy on fundamental climate features.

Dust aerosols↗

Combined LAURA-UPS solution procedure for chemically-reacting flows

A new procedure seeks to combine the thin-layer Navier-Stokes solver LAURA with the parabolized Navier-Stokes solver UPS for the aerothermodynamic solution of chemically-reacting air flowfields. The interface protocol is presented and the method is applied to two slender, blunted shapes. Both axisymmetric and three dimensional solutions are included with surface pressure and heat transfer comparisons between the present method and previously published results. The case of Mach 25 flow over an axisymmetric six degree sphere-cone with a noncatalytic wall is considered to 100 nose radii. A stability bound on the marching step size was observed with this case and is attributed to chemistry effects resulting from the noncatalytic wall boundary condition. A second case with Mach 28 flow over a sphere-cone-cylinder-flare configuration is computed at both two and five degree angles of attack with a fully-catalytic wall. Surface pressures are seen to be within five percent with the present method compared to the baseline LAURA solution and heat transfers are within 10 percent. The effect of grid resolution is investigated and the nonequilibrium results are compared with a perfect gas solution, showing that while the surface pressure is relatively unchanged by the inclusion of reacting chemistry the nonequilibrium heating is 25 percent higher. The procedure demonstrates significant, order of magnitude reductions in solution time and required memory for the three dimensional case over an all thin-layer Navier-Stokes solution.

Wood, William A.↗

Modeling of InP metalorganic chemical vapor deposition

The growth of InP by metalorganic chemical vapor deposition (MOCVD) in a horizontal reactor is being modeled with a commercially available computational fluid dynamics modeling code. The mathematical treatment of the MOCVD process has four primary areas of concern: 1) transport phenomena, 2) chemistry, 3) boundary conditions, and 4) numerical solution methods. The transport processes involved in CVD are described by conservation of total mass, momentum, energy, and atomic species. Momentum conservation is described by a generalized form of the Navier-Stokes equation for a Newtonian fluid and laminar flow. The effect of Soret diffusion on the transport of particular chemical species and on the predicted deposition rate is examined. Both gas-phase and surface chemical reactions are employed in the model. Boundary conditions are specified at the inlet and walls of the reactor for temperature, fluid flow and chemical species. The coupled set of equations described above is solved by a finite difference method over a nonuniform rectilinear grid in both two and three dimensions. The results of the 2-D computational model is presented for gravity levels of zero- and one-g. The predicted growth rates at one-g are compared to measured growth rates on fused silica substrates.

Black, Linda R.↗