Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 541 records · Page 30

Numerical Study of the Transient Vaporization of an Oxygen Droplet at Sub- and Super-Critical Conditions

Unsteady vaporization of a droplet in a high-pressure quiescent environment has been studied. With spherical symmetry and constant pressure, the process is diffusion controlled. At low pressures droplet heating is significant for most of the droplet lifetime, and an unsteady analysis of the gas phase is required. Then, the vaporization of a liquid oxygen droplet in gaseous hydrogen at moderate and high pressures is considered. At super-critical pressures, the surface temperature reaches the computed critical mixture value where a model for super-critical combustion is needed. The details of the diffusion layer of dissolved hydrogen do not significantly affect the results. Various methods are explored to compute the gas-phase density. Under the model's assumptions, homogeneous nucleation is not likely to occur.

Delplanque, J.-P.↗

Thermal effects in rapid directional solidification - Linear theory

We study the morphological instability of the planar solid/liquid interface for a unidirectionally-solidified dilute binary mixture. We use a model developed by Boettinger et al. (1985, 1986), Aziz (1982), and Jackson et al. (1980), which allows for nonequilibrium effects on the interface through velocity-dependent segregation and attachment kinetics. Two types of instabilities are found in the linear stability analysis: (1) a cellular instability, and (2) an oscillatory instability driven by disequilibrium effects. Merchant and Davis (1990) characterized these instabilities subject to the frozen-temperature approximation (FTA). The present work relaxes the FTA by including the effects of latent heat and the full temperature distribution. Thermal effects slightly postpone the onset of the cellular instability but dramatically postpone the onset of the oscillatory instability; however, the absolute-stability conditions, at which at high speed the cellular and oscillatory instabilities are suppressed, remain unchanged from the FTA.

Huntley, D. A.↗

Modelling early failures in Space Station Freedom

A major problem encountered in planning for Space Station Freedom is the amount of maintenance that will be required. To predict the failure rates of components and systems aboard Space Station Freedom, the logical approach is to use data obtained from previously flown spacecraft. In order to determine the mechanisms that are driving the failures, models can be proposed, and then checked to see if they adequately fit the observed failure data obtained from a large variety of satellites. For this particular study, failure data and truncation times were available for satellites launched between 1976 and 1984; no data past 1984 was available. The study was limited to electrical subsystems and assemblies, which were studied to determine if they followed a model resulting from a mixture of exponential distributions.

Navard, Sharon E.↗

Progress Toward an Efficient and General CFD Tool for Propulsion Design/Analysis

The simulation of propulsive flows inherently involves chemical activity. Recent years have seen substantial strides made in the development of numerical schemes for reacting flowfields, in particular those involving finite-rate chemistry. However, finite-rate calculations are computationally intensive and require knowledge of the actual kinetics, which are not always known with sufficient accuracy. Alternatively, flow simulations based on the assumption of local chemical equilibrium are capable of obtaining physically reasonable results at far less computational cost. The present study summarizes the development of efficient numerical techniques for the simulation of flows in local chemical equilibrium, whereby a 'Black Box' chemical equilibrium solver is coupled to the usual gasdynamic equations. The generalization of the methods enables the modelling of any arbitrary mixture of thermally perfect gases, including air, combustion mixtures and plasmas. As demonstration of the potential of the methodologies, several solutions, involving reacting and perfect gas flows, will be presented. Included is a preliminary simulation of the SSME startup transient. Future enhancements to the proposed techniques will be discussed, including more efficient finite-rate and hybrid (partial equilibrium) schemes. The algorithms that have been developed and are being optimized provide for an efficient and general tool for the design and analysis of propulsion systems.

Cox, C. F.↗

Did the Pre-RNA World Rest Upon DNA Molecules?

The isolation of a DNA sequence that catalyzes the ligation of oligodeoxynucleotides via the formation of 3' - 5' phosphodiester linkage significance in selection experiments has been reported. Ball recently used this to discuss the possibility that natural DNA molecules may have formed in the primitive Earth leading to the origin of life. As noted by Ferris and Usher, if metabolic pathways evolved backwards, it could be argued that the biosynthesis of 2-deoxyribose from ribose suggests that RNA came from DNA. As summarized elsewhere, there are several properties of deoxyribose which could be interpreted to support the possibility that DNA-like molecules arose prior to the RNA world. For example, 2-deoxyribose is slightly more soluble than ribose (which may have been an advantage in a drying pool scenario), may have been more reactive under possible prebiotic conditions (it forms a nucleoside approx. 150 times faster than ribose with the alternative base urazole at 25 C), while it decomposes in solution (approximately 2.6 times more slowly than ribose at 100 C). Other advantages of DNA over RNA are that it has one fewer chiral center, has a greater stability at the 8.2 pH value of the current oceans, and does not has the 2'5' and 3'5' ambiguity in polymerizations. Yet, there is strong molecular biological and biochemical evidence that RNA was featured in the biology well before the last common ancestor. The presence of sugar acids, including both ribo- and deoxysugar acids, in the 4.6 Ga old Murchison meteorite suggest that both may have been available in the primitive Earth, derived from the accretion of extraterrestrial sources and/or from endogenous processes involving formaldehyde and its derivatives. However, the abiotic synthesis of deoxyribose, ribose, and other sugars from glyceraldehyde and acetaldehyde under alkaline conditions is inefficient and unespecific. Although sugars are labile compounds, the role of cyanamide or borate minerals in the stabilization of the cyclic forms of ribose and other pentoses has recently been demonstrated. Nonetheless, the assumption either RNA or DNA was the first genetic material needs to be supplemented by laboratory models demonstrating that the prebiotic synthesis of activated beta-D-(deoxy)ribonucleotides and their polymers was feasible. As of today such evidence is lacking, and there is no convincing synthesis of any nucleotide, since all model experiments produce complex mixtures of products in which there is no preferential synthesis of chiral D-nucleotides. This strongly suggests that both DNA and RNA may have been preceded by pairing structures much simpler than extant nucleic acids. It is doubtful that DNA molecules, or indeed other (de0xy)ribofuranoid oligonucleotides formed the basis of these as yet undescribed pre-RNA worlds.

Lazcano, Antonio↗

A Model Study of the Thermal Evolution of Astrophysical Ices

We address the question of the evolution of ices that have been exposed to radiation from stellar sources and cosmic rays. We studied in the laboratory the thermal evolution of a model ice sample: a mixture of water, hydrogen peroxide, dioxygen, and ozone produced by irradiating solid H2O2 with 50 keV H(+) at 17 K. The changes in composition and release of volatiles during warming to 200 K were monitored by infrared spectroscopy, mass spectrometry, and microbalance techniques. We find evidence for voids in the water component from the infrared bands due to dangling H bonds. The absorption from these bands increases during heating and can be observed at temperatures as high as approx. 155 K. More O2 is stored in the radiolyzed film than can be retained by codeposition of O2 and H2O. This O2 remains trapped until approx. 155 K, where it desorbs in an outburst as water ice crystallizes. Warming of the ice also drastically decreases the intrinsic absorbance of O2 by annealing defects in the ice. We also observe loss of O3 in two stages during heating, which correlates with desorption and possibly chemical reactions with radicals stored in the ice, triggered by the temperature increase.

Loeffler, M. J.↗

Method and Apparatus for Predicting Unsteady Pressure and Flow Rate Distribution in a Fluid Network

A method and apparatus for analyzing steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics, external body forces such as gravity and centrifugal force and conjugate heat transfer. In some embodiments, a graphical user interface provides for the interactive development of a fluid network simulation having nodes and branches. In some embodiments, mass, energy, and specific conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. In some embodiments, contained herein are data objects for computing thermodynamic and thermophysical properties for fluids. In some embodiments, the systems of equations describing the fluid network are solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods.

Majumdar, Alok K.↗

Application of Artificial Neural Networks to the Development of Improved Multi-Sensor Retrievals of Near-Surface Air Temperature and Humidity Over Ocean

Improved estimates of near-surface air temperature and air humidity are critical to the development of more accurate turbulent surface heat fluxes over the ocean. Recent progress in retrieving these parameters has been made through the application of artificial neural networks (ANN) and the use of multi-sensor passive microwave observations. Details are provided on the development of an improved retrieval algorithm that applies the nonlinear statistical ANN methodology to a set of observations from the Advanced Microwave Scanning Radiometer (AMSR-E) and the Advanced Microwave Sounding Unit (AMSU-A) that are currently available from the NASA AQUA satellite platform. Statistical inversion techniques require an adequate training dataset to properly capture embedded physical relationships. The development of multiple training datasets containing only in-situ observations, only synthetic observations produced using the Community Radiative Transfer Model (CRTM), or a mixture of each is discussed. An intercomparison of results using each training dataset is provided to highlight the relative advantages and disadvantages of each methodology. Particular emphasis will be placed on the development of retrievals in cloudy versus clear-sky conditions. Near-surface air temperature and humidity retrievals using the multi-sensor ANN algorithms are compared to previous linear and non-linear retrieval schemes.

Roberts, J. Brent↗

T-Matrix Modeling of Linear Depolarization by Morphologically Complex Soot and Soot-Containing Aerosols

We use state-of-the-art public-domain Fortran codes based on the T-matrix method to calculate orientation and ensemble averaged scattering matrix elements for a variety of morphologically complex black carbon (BC) and BC-containing aerosol particles, with a special emphasis on the linear depolarization ratio (LDR). We explain theoretically the quasi-Rayleigh LDR peak at side-scattering angles typical of low-density soot fractals and conclude that the measurement of this feature enables one to evaluate the compactness state of BC clusters and trace the evolution of low-density fluffy fractals into densely packed aggregates. We show that small backscattering LDRs measured with groundbased, airborne, and spaceborne lidars for fresh smoke generally agree with the values predicted theoretically for fluffy BC fractals and densely packed near-spheroidal BC aggregates. To reproduce higher lidar LDRs observed for aged smoke, one needs alternative particle models such as shape mixtures of BC spheroids or cylinders.

atmospheric radiation↗

Validation of the aerosol optical property products derived by the GRASP/Component approach from multi-angular polarimetric observations

The multi-angle polarimetric satellite observations are helpful for improving the retrievals of aerosol parameters. However, practical applications of polarization technology are still limited because of complexity of measurement and interpretation of polarimetric observations. In this study, we analyze the performance of a new component approach developed in the frame of the Generalized Retrieval of Atmosphere and Surface Properties (GRASP) algorithm. In addition to aerosol optical properties including particle size distribution, non-sphericity and index of refraction that are commonly derived from multi-angle radiance and polarization measurements, the GRASP/Component approach also provides some information about aerosol composition. Specifically, in this approach aerosol is modeled as an internal mixture of several components with distinctly different chemical compositions and known refractive indices. The approach is intended not only to provide additional insight on aerosol composition but also to improve retrieval of basic aerosol optical properties. This study presents comprehensive validation and evaluation of Aerosol Optical Depth (AOD), Ångstrom exponent (AE), fine mode AOD (AODF), coarse mode AOD (AODC), and single scattering albedo (SSA) as retrieved by the GRASP/Component approach. The GRASP/Component products include aerosol retrievals using two different aerosol component mixing rules, i.e., Maxwell-Garnett (MG) effective medium approximation and a simple Volume- Weighted averaging (VW). The differences between the results obtained using these two assumptions are also discussed. The obtained results show that the aerosol optical property products of GRASP/Component approach have good agreement with the ground-based AERONET measurements, which is comparable to other PARASOL/ GRASP approaches. Specifically, the AOD retrieved by GRASP/Component approach show high correlation and nearly no bias both over land and ocean, as compared with AERONET. The more detailed aerosol properties such as AE, AODF, AODC and SSA also show one of the best comparisons with AERONET. These improvements can probably be attributed to the use of the additional physical constraints on spectral dependence of the complex refractive index and the reduction of total number of aerosol parameters directly retrieved in the GRASP/Component approach. In addition, the choice of mixing rules had no significant effect on optical retrievals. With the exception of SSA, the results obtained based on the MG mixing rule were found to be slightly better over those obtained using VW mixing rule, especially for bias.

Brent N Holben↗

The CO2-Broadened H2O Continuum in the 100-1500 cm(sup -1) Region: Measurements, Predictions and Empirical Model

Transmission spectra of H2O+CO2 mixtures have been recorded, at 296, 325 and 366 K, for various pressures and mixture compositions using two experimental setups. Their analysis enables to retrieve values of the “continuum” absorption by the CO2-broadened H2O line wings between 100 and 1500 cm(exp -1). The results are in good agreement with those, around 1300 cm(exp -1), of the single previous experimental study available. Comparisons are also made with direct predictions based on line-shape correction factors X calculated, almost thirty years ago, using a quasistatic approach and an input H2O-CO2 intermolecular potential. They show that this model quite nicely predicts, with slightly overestimated values, the continuum over a spectral range where it varies by more than three orders of magnitude. An empirical correction is proposed, based on the experimental data, which should be useful for radiative transfer and climate studies in CO2 rich planetary atmospheres.

CO2↗

AL4GAP: Active learning workflow for generating DFT-SCAN accurate machine-learning potentials for combinatorial molten salt mixtures

Machine learning interatomic potentials have emerged as a powerful tool for bypassing the spatiotemporal limitations of ab initio simulations, but major challenges remain in their efficient parameterization. We present AL4GAP, an ensemble active learning software workflow for generating multicomposition Gaussian approximation potentials (GAP) for arbitrary molten salt mixtures. The workflow capabilities include: (1) setting up user-defined combinatorial chemical spaces of charge neutral mixtures of arbitrary molten mixtures spanning 11 cations (Li, Na, K, Rb, Cs, Mg, Ca, Sr, Ba and two heavy species, Nd, and Th) and 4 anions (F, Cl, Br, and I), (2) configurational sampling using low-cost empirical parameterizations, (3) active learning for down-selecting configurational samples for single point density functional theory calculations at the level of Strongly Constrained and Appropriately Normed (SCAN) exchange-correlation functional, and (4) Bayesian optimization for hyperparameter tuning of two-body and many-body GAP models. Here, we apply the AL4GAP workflow to showcase high throughput generation of five independent GAP models for multicomposition binary-mixture melts, each of increasing complexity with respect to charge valency and electronic structure, namely: LiCl–KCl, NaCl–CaCl 2 , KCl–NdCl 3 , CaCl 2 –NdCl 3 , and KCl–ThCl 4 . Our results indicate that GAP models can accurately predict structure for diverse molten salt mixture with density functional theory (DFT)-SCAN accuracy, capturing the intermediate range ordering characteristic of the multivalent cationic melts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Corrected Distortion model for Lagrangian spray simulation of transcritical fuel injection

In this work, we present a detailed implementation and validation of the droplet modeling framework proposed by Dahms and Oefelein (2016) into the engine commercial CFD software CONVERGE using the User Defined Function (UDF) interface. The model accounts for the nonlinear deformation and oscillation experienced by liquid spray droplet injected into high pressure and temperature. Lagrangian spray simulations of Engine Combustion Network (ECN) Spray A are performed. Model validation against standard experimental measurements of liquid velocity, vapor mixture fraction is conducted. To perform more rigorous model validation, new experimental measurements based on Diffused Back Illumination (DBI) are introduced. The new measurements are processed for Projected Liquid Volume (PLV), which offers as close as possible one-to-one model validation for liquid penetration while offering new insights into the spray physics. Comparison with a One-D model based on adiabatic mixing theory by Siebers (1999) and Desantes et al. (2007) are also conducted. Through these model validation exercises, it is shown that the new framework improves liquid-phase penetration predictions, following a tendency for enhanced evaporation, compared to the standard approach for both Reynolds Average Navier Stokes (RANS) and Large Eddy Simulation (LES). Additionally, at the liquid length, maximum mixture fraction values predicted by the new approach are in good agreement those of an adiabatic mixing model. Qualitative analysis of the spray behaviors during the early stage of the injection process reveals that the proposed framework predicts significant increase in droplet evaporation rate with lower drop drag compared to the current standard approach.

42 ENGINEERING↗

Modeling and optimization of steady flow of natural gas and hydrogen mixtures in pipeline networks

Here, we extend the canonical problems of simulation and optimization of steady-state gas flows in pipeline networks with compressors to the transport of mixtures of highly heterogeneous gases injected throughout a network. Our study is motivated by proposed projects to blend hydrogen generated using clean energy into existing natural gas pipeline systems as part of efforts to reduce the reliance of energy systems on fossil fuels. Flow in a pipe is related to endpoint pressures by a basic Weymouth equation model, with an ideal gas equation of state, where the wave speed depends on the hydrogen concentration. At vertices, in addition to mass balance, we also consider mixing of incoming flows of varying hydrogen concentrations. The problems of interest are the heterogeneous gas flow simulation (HGFS), which determines system pressures and flows given fixed boundary conditions and compressor settings, as well as the heterogeneous gas flow optimization (HGFO), which extremizes an objective by determining optimal boundary conditions and compressor settings. We examine conditions for uniqueness of solutions to the HGFS, as well as compare and contrast mixed-integer and continuous nonlinear programming formulations for the HGFO. We develop computational methods to solve both problems, and examine their performance using four test networks of increasing complexity.

08 HYDROGEN↗

Application of an Empirical Density Law via Python for Aqueous Plutonium Chloride Systems in MCNP6

Current aqueous plutonium processing models for criticality safety often contain significant bias due to material modeling assumptions. These solutions include plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution. Furthermore, there is no current predictive capability for modeling plutonium metal-water mixtures that is approved for use at Los Alamos National Laboratory (LANL). Recent density measurements for aqueous plutonium chloride systems (PuCl 3 -HCl-H 2 O) now allow for the development and application of a more realistic density law. This work develops a Python-based density law for this ternary solution using an empirical method. This code can be used in conjunction with an MCNP6 input to determine the density and composition of a solution based upon user inputs of plutonium concentration, hydrochloric acid concentration, and temperature. The tool allows users to input plutonium and acid content of a solution in terms of molality, molarity, or concentration, and predicts density within the current data range within 1.4% of experimental data. Current preliminary MCNP6 calculations utilizing this tool have demonstrated a minimum decrease in system reactivity of 5% in comparison to the current modeling conventions. Thus, this tool enables more accurate criticality safety operational limits by better crediting chorine content while still maintaining necessary conservatism.

Density Law↗

Viscosity of rock-ice mixtures and applications to the evolution of icy satellites

Theory and experiments are used to establish lower and upper bounds on the ratio of actual viscosity to pure ice viscosity for a suspension of rock particles in a water ice matrix. A rheological model for rock-ice mixtures is described, establishing bounds for the range of possible viscosity enhancements provided by a suspension of silicate spheres in an ice matrix. A parametrized thermal convection model is described and used to determine a criterion for criticality, defined as the heat flow and/or silicate volume fraction for which the satellite temperature profile intercepts the melting curve of water ice. The consequences of achieving this critical state are examined, and it is shown that under certain circumstances a 'runaway' differentiation can occur in which the silicates settle to form a core and extensive melting of water ice takes place, the latent heat being supplied by the gravitational energy of differentiation. A possible application of these results to Ganymede and Callisto is described.

Friedson, A. J.↗

Tricritical points in multicomponent fluid mixtures

In view of experimental considerations, we give a model-independent argument that the novel tricritical points in multicomponent fluid mixtures, where three phases simultaneously become critical, are points on the boundary of a single two-dimensional surface of critical points. This result is corroborated by the Landau model suggested by Griffiths. The relationship between these tricritical points and the complex 'higher-order' critical points proposed to exist in certain magnetic systems is elucidated.

Hankey, A.↗

State-Specific Kinetic Modeling for Predictions of Radiative Heating in H 2 /He Entry Flows

A 23-species State-to-State (StS) thermochemical model for H 2 /He mixtures is developed and used to simulate 1-D shocks at conditions relevant for ice and gas giant entry flows. To implement this StS model, a literature review of the state-specific excitation and ionization rates of atomic H is first performed. While electron impact rates from various sources are found to be in relatively good agreement, large discrepancies are found in the limited data available on heavy particle impact rates. To investigate the impact of these discrepancies on the predicted shock and radiance profiles, simulations are performed for three different versions of the kinetic model, each with a different set of heavy particle impact rates from the literature. Of these three models, two are found to over-predict and one is found to under-predict the radiance values seen in corresponding experimental data from the NASA Ames Electric Arc Shock Tube (EAST) facility. A sensitivity analysis to the electron and heavy particle impact rates is then performed, from which an improved kinetic model with modified rates is proposed. This modified kinetic model is found to reproduce the total integrated radiance seen in the EAST experiments accurately. This suggests that the associated radiative heating that would be encountered by an entry probe in these environments would be well predicted by this kinetic model. However, some discrepancies are still found i the relative contributions to the total radiance from different spectral features, suggesting that additional improvements to the kinetic model still need to be made.

Outer Planet↗