Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution species”

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 289 records · Page 16

Variable Lidar Ratios for Marine and Dusty Marine Using MODIS AOD Constrained Retrievals and GOCART Model Simulations and their Impact on Level 2 CALIPSO Data Products in Version 5

The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard CALIPSO provided global measurements of attenuated backscatter profiles of various tropospheric aerosol species from June 2006 through June 2023. Extinction profiles are retrieved from these backscatter measurements by assuming a value for the extinction-to-backscatter or lidar ratio (LR). Up until version V4.51 of the data products, only a constant global value for each of the species has been used, spatially and temporally. However, it has been well known from various ground-based measurements that LRs of each aerosol species can potentially vary spatially as well seasonally. For the forthcoming Version 5 of the CALIPSO products, the CALIPSO team will implement variable lidar ratios at 532 nm for the species classified as “marine” and “dusty marine”. Here we present the essential elements of this scheme, which involves estimating lidar ratios by constraining Fernald solutions for the CALIOP particulate backscatter profiles using aerosol optical depths (AOD) from collocated MODIS retrievals. To account for the lack of measurements in various regions and seasons (e.g., Arctic winters), we leverage collocated model sea-salt volume fractions (SSVF) simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model and use an empirical relationship between the SSVF and LR to build hybrid climatological maps for the entire globe, i.e., using the modelled values when reliable values of the constrained LRs are not available. Initial evaluation of the resulting extinction profiles and corresponding AODs, both globally and in specific regions are presented using a limited amount of data. The primary improvement in AOD occurs in the coastal areas, like the Arabian Sea and Bay of Bengal. The estimated marine LR values in these regions are significantly higher than the constant value of 23 sr used in previous data releases, likely due to mixing with offshore pollution which results in higher AODs. Initial validation results using the CALIPSO Ocean Derived Column Optical Depth (ODCOD) retrievals and ground-based AERONET AOD measurements are presented.

CALIPSO↗

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗

Computation of iodine species concentrations in water

During an evaluation of the use of iodine as a water disinfectant and the development of methods for measuring various iodine species in water onboard Space Freedom, it became necessary to compute the concentration of the various species based on equilibrium principles alone. Of particular concern was the case when various amounts of iodine, iodide, strong acid, and strong base are added to water. Such solutions can be used to evaluate the performance of various monitoring methods being considered. The authors of this paper present an overview of aqueous iodine chemistry, a set of nonlinear equations which can be used to model the above case, and a computer program for solving this system of equations using the Newton-Raphson method. The program was validated by comparing results over a range of concentrations and pH values with those previously presented by Gottardi for a given pH. Use of this program indicated that there are multiple roots to many cases and selecting an appropriate initial guess is important. Comparison of program results with laboratory results for the case when only iodine is added to water indicates the program gives high pH values for the iodine concentrations normally used for water disinfection. Extending the model to include the effects of iodate formation results in the computer pH values being closer to those observed, but the model with iodate does not agree well for the case in which base is added in addition to iodine to raise the pH. Potential explanations include failure to obtain equilibrium conditions in the lab, inaccuracies in published values for the equilibrium constants, and inadequate model of iodine chemistry and/or the lack of adequate analytical methods for measuring the various iodine species in water.

Schultz, John R.↗

Numerical simulation of hypersonic flow over an aeroassist flight experiment vehicle

Solutions of wind-tunnel and entry-flight flow around the vehicle are obtained from the Navier-Stokes equations coupled with the chemical species continuity equations if needed. The time-iterative method employs several techniques: shock fitting, chemistry-split ADI and an algebraic grid in conformal spherical-polar space. Sensitivities of the results to numerical parameters and to frozen, equilibrium and finite rate reactions are investigated in the forebody computation. Quantitative results are obtained for the shock layer and the near wake for the entire vehicle corresponding to both ground test and flight conditions. Complex flow characteristics are analyzed on the basis of the complete flowfield over the aerobrake and simplified afterbodies. The method is stable and cost effective, and has yielded shock locations and wall pressure distributions which are in good agreement with wind-tunnel data.

Li, C. P.↗

Numerical Simulation of Lean Blowout of Alternative Fuels in 7-element Lean Direct Injector

This research presents the result of numerically simulating 7-element swirl-venturi Lean Direct Injector (SV-LDI) lean blowout (LBO) experiments conducted at NASA Glenn Research Center in May of 2019. After simulating a cold flow case to confirm the pressure drop agrees well with the experiment, additional cases with two different fuels (an average jet fuel and a Gevo alcohol-to-jet fuel) from the National Jet Fuels Combustion Program (NJFCP) were computed to numerically determine the LBO condition. The procedure to approach the LBO follows the method used in the experiment where the air mass flow rate is gradually increased while the fuel supply is maintained. Transient history of global heat release rate as a function of air flow rate is presented, as well as temperature contours at different conditions to give a visual representation of the flame state. The Open National Combustion Code (OpenNCC) used in this research adopted reduced HyChem (Hybrid Chemistry) models along with k-LES turbulence model and a Lagrangian spray model that takes into account droplet internal temperature distribution affected by the shear force on the droplet surface. The transport equations and chemical reaction terms are integrated together to enhance conservation of chemical species that are especially important in the near LBO conditions. After showing the computed range of LBO agrees well with the experimental measurements, time averaged solutions of both fuels at their initial condition and at their limiting condition just before LBO are compared in detail to facilitate the understanding of LBO mechanism.

combustion↗

Capitol Reef Ecological Conservation: Mapping Vegetation Functional Groups to Inform Invasive Vegetation Management, Ecological Conservation and Restoration in Capitol Reef National Park

Invasive exotic plant (IEP) species have been found within the park boundaries of Capitol Reef National Park (CARE) in Utah. Currently, remotely sensed datasets such as the Rangeland Analysis Platform (RAP) from the United States Department of Agriculture (USDA) have been used to investigate IEP species within the park, but validation of the national RAP program is necessary for informing decisions at a local scale. CARE seeks a remote monitoring solution that can precisely target managerial efforts within the park’s challenging terrain and hard-to-reach locations. To fulfill this objective, we harnessed Landsat 8 Operational Land Imager (OLI) imagery and leveraged Random Forest (RF) modeling to generate classification maps characterizing vegetation functional groups for 2013 and 2022 within the park. Subsequently, the Land Change Modeler (LCM) in Idrisi TerrSet facilitated the production of a predicted classification map for 2033. The team also devised an annual grass probability map to accentuate areas impacted by exotic grasses. A comparative assessment between the RF classification map and the RAP map for 2022 revealed an overall agreement of 47.41%, with disparities primarily arising from differences in bare soil and shrub areas. Significantly, the 2022 RF-generated classification map showcased an impressive overall accuracy of 92.17%. In short, the probability map, the land cover change detection spanning 2013 to 2022, and the forecasting of observed trends into the future aids in the evaluation of invasive plant impacts and facilitation of CARE’s preparedness for potential ecological disturbances. Notably, in comparison to the RAP, the RF classification method generates functional group maps that are more representative of the study area.

Vanchy Li↗

Graphite solubility and co-vesiculation in basalt-like melts at one-ATM

The identity and source of the vapor phase that caused lunar lava-fountaining and vesiculation in lunar basalts continues to be of interest because of its implications for the composition and state of the lunar interior and because of its implications for lunar resources. In light of the apparent near-absence of H2O on the Moon, it has been suggested that the vapor phase may be CO2-CO. This premise is supported by the presence of carbon on the surface of volcanic glass beads. However, although the rapid exsolution of CO2 from a melt during decompression may be consistent with firefountaining, it fails to provide a satisfying explanation for vesiculation in mare basalt where exsolution of the gas phase would more reasonably be related to cooling/crystallization at low pressure rather than decompression from high pressure. Also, geochemical trends in lunar volcanic glasses suggest that their source has an oxygen fugacity more reducing than the iron-wustite buffer, an oxygen fugacity that is inconsistent with presence of dissolved CO2-CO at depth. The results of experiments in which a vesicular 'basalt' is produced from a melt equilibrated with graphite and pure CO gas at one atmosphere pressure are reported. The vesiculation is apparently related to exsolution of CO or a CO species during cooling of the melt or growth of quench crystals. Additionally, particulate carbon dispersed through the quenched sample suggests that elemental carbon is either in solution in the melt prior to quenching or tends to go into suspension perhaps as colloid-like particles. These two observations may provide insight into the nature of fire-fountaining and vesiculation on the Moon.

Colson, R. O.↗

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗

On upstream blocking in a viscous diffusive stratified flow

The effect of diffusion of specie upon the flow about a transverse flat plate moving horizontally in a viscous stratified medium is considered. Asymptotic expansions are used to define a parameter regime where a viscous-diffusive-buoyancy balance is dominant. The solution, expressed in terms of an inverse Fourier transform, is numerically integrated. The results show that, as in the non-diffusive problem, a region of closed streamlines exists ahead of the body. However, unlike the case where diffusion is neglected, the density field within this recirculating region is uniquely determined and found to be statically stable. It is also found that varying the relative amount of diffusion affects not only the density distribution, but the velocity profile as well, indicating a strong coupling between the vorticity and specie equation.

Koop, C. G.↗

Moving wall, continuous flow electronphoresis apparatus

This invention relates generally to electrophoresis devices and more particularly to a moving wall, continuous flow device in which an electrophoresis chamber is angularly positionable with respect to the direction of moving belt walls. A frame with an electrophoresis chamber is rotatably supported between two synchronously driven belt walls. This allows the chamber to be angularly positionable with respect to the direction of belt travel, which compensates for electroosmotic flow within the electrophoresis chamber. Injection of a buffer solution via an opening and a homogenous sample stream via another opening is performed at the end of a chamber, and collection of buffer and the fractionated species particles is done by a conventional collection array at an opposite end of the chamber. Belts are driven at a rate which exactly matches the flow of buffer and sample through the chamber, which entrains the buffer to behave as a rigid electrophoretic medium, eliminating flow distortions (Poiseuille effect). Additionally, belt material for each belt is stored at one end of the device and is taken up by drive wheels at an opposite end. The novelty of this invention particularly lies in the electrophoresis chamber being angularly positionable between two moving belt walls in order to compensate for electroosmotic flow. Additionally, new belt material is continuously exposed within the chamber, minimizing flow distortion due to contamination of the belt material by the sample.

Rhodes, Percy H.↗

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

Toward a Micro Gas Chromatograph/Mass Spectrometer (GC/MS) System

Miniature mass filters (e.g., quadrupoles, ion traps) have been the subject of several miniaturization efforts. A project is currently in progress at JPL to develop a miniaturized Gas Chromatograph/Mass Spectrometer (GC/MS) system, incorporating and/or developing miniature system components including turbomolecular pumps, scroll type roughing pump, quadrupole mass filter, gas chromatograph, precision power supply and other electronic components. The preponderance of the system elements will be fabricated using microelectromechanical systems (MEMS) techniques. The quadrupole mass filter will be fabricated using an X-ray lithography technique producing high precision, 5x5 arrays of quadrupoles with pole lengths of about 3 mm and a total volume of 27 cubic mm. The miniature scroll pump will also be fabricated using X-ray lithography producing arrays of scroll stages about 3 mm in diameter. The target detection range for the mass spectrometer is 1 to 300 atomic mass units (AMU) with are solution of 0.5 AMU. This resolution will allow isotopic characterization for geochronology, atmospheric studies and other science efforts dependant on the understanding of isotope ratios of chemical species. This paper will discuss the design approach, the current state-of-the art regarding the system components and the progress toward development of key elements. The full system is anticipated to be small enough in mass, volume and power consumption to allow in situ chemical analysis on highly miniaturized science craft for geochronology, atmospheric characterization and detection of life experiments applicable to outer planet roadmap missions.

Wiberg, D. V.↗

Role of the area postrema in three putative measures of motion sickness in the rat

After thermal cauterization of the area postrema in rats, the absence of conditioned taste aversion of sucrose paired with lithium chloride (0.15M, 3.3 ml/kg) was used as a pharmacologic/behavioral index of area postrema damage. In a subsequent experiment the effects of area postrema lesions on three measures proposed as species-relevant measures of motion sickness were studied, using off-vertical rotation at 150 deg/s for either 30 or 90 min. Lesions of area postrema did not alter postrotational suppression of drinking or amount of defecation during motion. The initial acquisition of conditioned taste aversion to a novel cider vinegar solution paired with motion was not affected by lesioning of the area postrema, but these taste aversions extinguished more slowly in lesioned rats than in sham-operates or intact controls. Results are discussed in terms of proposed humoral factors which may induce motion sickness and in light of recent data on the role of the area postrema in similar measures in species possessing the complete emetic reflex.

Richard L Sutton↗

Numerical simulation of a supersonic reacting mixing layer

In order to arrive at physical models that can adequately describe supersonic combustion, and develop accurate and efficient numerical techniques for the solution of such models' governing equations, a computer program has been developed for the study of reacting flows which considers the multicomponent diffusion and convection of important chemical species, as well as their finite state reaction and the interaction of the fluid mechanics and the chemistry that occurred. The code employs a hybrid Chebyshev pseudospectral technique for integration of the models' resulting governing equations; the program is here used to study a spatially developing and reacting mixing layer.

Drummond, J. Philip↗

Ignition and structure of a laminar diffusion flame in the field of a vortex

The distortion of flames in flows with vortical motion is examined via asymptotic analysis and numerical simulation. The model consists of a constant density, one step, irreversible Arrhenius reaction between initially unmixed species occupying adjacent half-planes which are then allowed to mix and react in the presence of a vortex. The evolution in time of the temperature and mass fraction fields is followed. Emphasis is placed on the ignition time and location as a function of vortex Reynolds number and initial temperature differences of the reacting species. The study brings out the influence of the vortex on the chemical reaction. In all phases, good agreement is observed between asymptotic analysis and the full numerical solution of the model equations.

Macaraeg, Michele G.↗

Testing and Application of Commercial Thermodynamics Software for Solid-Gas Equilibria in Planetary Science Applications

The ability to calculate chemical equilibria from thermochemical constants (e.g., H, S, Cp) or other experimental data is a tool widely employed by scientists. However, phase equilibria for planetary science problems is limited by easy access to libraries of thermochemical data, as well as the software to invert that data. Many academic software packages are focused on specific types of problems (e.g., aqueous solutions or melting and crystallization of magmas). This focus makes them excellent tools for specific uses, but poorer tools for other scenarios, especially those involving elements or chemical species that are less well-studied in that specific scientific field (e.g. NaCl vapor). Amongst the more general thermodynamic software tools, there are closed academic software packages, open-source options, as well as commercial software. We have explored the extent to which commercial chemistry/chemical engineering software may be able to satisfy the need for thermochemical modeling packages that can accommodate the diverse species and conditions of planetary science problems. Here we compare the output of software package HSC Chemistry™ (Metso Outotec) to published vapor-solid phase diagrams for a variety of major and minor elements. HSC Chemistry™ has been used in several recent geochemical studies of terrestrial volcanic systems [2- 4]. However, it has not to our knowledge been benchmarked against peer-reviewed model results for the conditions and compositions relevant to planetary science. We present one such comparison here, the low- pressure environment of the solar nebula, using literature data. In short, we find excellent agreement between HSC Chemistry™ and the well-regarded academic software package CONDOR. These preliminary results suggest that at least some commercial software packages are capable of robust thermochemical calculations for planetary science.

planetary science↗

Computation of Reacting Flows in Combustion Processes

The main objective of this research was to develop an efficient three-dimensional computer code for chemically reacting flows. The main computer code developed is ALLSPD-3D. The ALLSPD-3D computer program is developed for the calculation of three-dimensional, chemically reacting flows with sprays. The ALL-SPD code employs a coupled, strongly implicit solution procedure for turbulent spray combustion flows. A stochastic droplet model and an efficient method for treatment of the spray source terms in the gas-phase equations are used to calculate the evaporating liquid sprays. The chemistry treatment in the code is general enough that an arbitrary number of reaction and species can be defined by the users. Also, it is written in generalized curvilinear coordinates with both multi-block and flexible internal blockage capabilities to handle complex geometries. In addition, for general industrial combustion applications, the code provides both dilution and transpiration cooling capabilities. The ALLSPD algorithm, which employs the preconditioning and eigenvalue rescaling techniques, is capable of providing efficient solution for flows with a wide range of Mach numbers. Although written for three-dimensional flows in general, the code can be used for two-dimensional and axisymmetric flow computations as well. The code is written in such a way that it can be run in various computer platforms (supercomputers, workstations and parallel processors) and the GUI (Graphical User Interface) should provide a user-friendly tool in setting up and running the code.

Keith, Theo G., Jr.↗

Prediction of electron density in the boundary layer on entry vehicles with ablation

A rather general computer code has been developed for the numerical solution of the laminar boundary layer and thin shock layer (stagnation point only) equations for a multi-component gas mixture with finite reaction rates. The purpose of this paper is to indicate the capabilities of this computer program by presenting solutions for pure air flows and flows with carbon ablation. At the body surface the effects of oxidation and sublimation of carbon are taken into account with the following species included in the gas model: 0 2 , N 2 , O, N, NO, NO + , C 1 , C 2 , C 3 , CN, CO, and CO 2 . Results for the boundary layer along a 10° half-angle hyperboloid with ablation of carbon are compared to pure air results. Shock layer results for several wall temperatures and altitudes are presented and compared to pure air results. Other types of problems that can be solved with this computer code are indicated.

F. G. Blottner↗