Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Yield Prediction”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Critical velocity experiments in space

Published data from active space experiments designed to demonstrate the Alfven critical-velocity effect are compiled in graphs and compared with the predictions of numerical simulations. It is found that the discrepancies in the ionization yields obtained in shaped-charge releases of alkali metals are related to the macroscopic limits of time and energy in such releases. It is argued that the total ionization yield is an inadequate measure of the critical-velocity effect, and a new criterion based on eta, the efficiency of energy transfer from the recently ionized neutrals to a heated electron population, is proposed: the effect would be verified if eta values of 10 percent or greater were observed.

Torbert, R. B.↗
NASA NTRS1993

Evaporation kinetics of Mg2SiO4 crystals and melts from molecular dynamics simulations

Computer simulations based on the molecular dynamics (MD) technique were used to study the mechanisms and kinetics of free evaporation from crystalline and molten forsterite (i.e., Mg2SiO4) on an atomic level. The interatomic potential employed for these simulations reproduces the energetics of bonding in forsterite and in gas-phase MgO and SiO2 reasonably accurately. Results of the simulation include predicted evaporation rates, diffusion rates, and reaction mechanisms for Mg2SiO4(s or l) yields 2Mg(g) + 20(g) + SiO2(g).

Kubicki, J. D.↗
NASA NTRS1993

The dynamics of hydroponic crops for simulation studies of the CELSS initial reference configurations

During the past several years, the NASA Program in Controlled Ecological Life Support Systems (CELSS) has continued apace with crop research and logistic, technological, and scientific strides. These include the CELSS Test Facility planned for the space station and its prototype Engineering Development Unit, soon to be active at Ames Research Center (as well as the advanced crop growth research chamber at Ames); the large environmental growth chambers and the planned human test bed facility at Johnson Space Center; the NSCORT at Purdue with new candidate crops and diverse research into the CELSS components; the gas exchange data for soy, potatoes, and wheat from Kennedy Space Center (KSC); and the high-precision gas exchange data for wheat from Utah State University (USU). All these developments, taken together, speak to the need for crop modeling as a means to connect the findings of the crop physiologists with the engineers designing the system. A need also exists for crop modeling to analyze and predict the gas exchange data from the various locations to maximize the scientific yield from the experiments. One fruitful approach employs what has been called the 'energy cascade'. Useful as a basis for CELSS crop growth experimental design, the energy cascade as a generic modeling approach for CELSS crops is a featured accomplishment in this report. The energy cascade is a major tool for linking CELSS crop experiments to the system design. The energy cascade presented here can help collaborations between modelers and crop experimenters to develop the most fruitful experiments for pushing the limits of crop productivity. Furthermore, crop models using the energy cascade provide a natural means to compare, feature for feature, the crop growth components between different CELSS experiments, for example, at Utah State University and Kennedy Space Center.

Volk, Tyler↗
NASA NTRS1994

Chemical transitions for interstellar C2 and CN in cloud envelopes

Observations were made of absorption from CH, C2, and CN toward moderately reddened stars in Sco, OB2, Ceo OB3, and Taurus/Auriga. For these directions, most of the reddening is associated with a single cloud complex, for example, the rho Ophiuchus molecular cloud, and as a result, the observations probe moderately dense material. When combined with avaliable data for nearby directions, the survey provides the basis for a comprehensive analysis of the chemistry for these species. The chemical transitions affecting C2 and CN in cloud envelopes were analyzed. The depth into a cloud at which a transition takes place was characterized by tau(sub uv), the grain optical depth at 1000 A. One transition at tau(sub uv) approx. = 2, which arises from, the conversion of C(+) into CO, affects the chemistries for both molecules because of the key role this ion plays. A second one involving production terms in the CN chemistry occurs at tau(sub uv) of approx. = 3; neutral reactions which C2 and CH is more important at larger values for tau(sub uv). The transition from photodissociation to chemical destruction takes place at tau(sub uv) approx. = 4.5 for C2 and CN. The observational data for stars in Sco OB2, Cep OB3, and Taurus/Auriga were studied with chemical rate equations containing the most important production and destruction mechanisms. Because the sample of stars in Sco OB2 includes sight lines with A(sub v) ranging from 1-4 mag, sight lines dominated by photochemistry could be analyzed separately from those controlled by gas-phase destruction. The analysis yielded values for two poorly known rate constants for reactions involved in the production of CN; the reactions are C2 + N yields CN + C and C(+) + NH yields all products. The other directions were analyzed with the inferred values. The predicted column densities for C2 and CN agree with the observed values to better than 50%, and in most instances 20%. When combining the estimates for density and temperature derived from chemical modeling and molecular excitation for a specific cloud, such as the rho Ophiuchus molecular cloud, the portion of the cloud envelope probed by C2 and CN absorption was found to be in pressure equilibrium.

Federman, S. R.↗
NASA NTRS1995

The radial component of the heliospheric magnetic field: Ulysses observations

The radial field component, B(sub R), has been monitored continuously since the Ulysses spacecraft left the ecliptic plane in February 1992 travelling toward the southern pole of the Sun. In order to separate spatial from temporal changes, the Ulysses measurements from 0 to 80 heliographic latitude were compared with in-ecliptic measurements of B(sub R) being made simultaneously by IMP-8. The data revealed essentially the same field strengths and time variations at both locations. The conclusion was drawn that there was no significant latitude gradient in B(sub R) and that the stronger polar cap coronal magnetic fields were being transported equatorward to yield a uniform field in the solar wind. The results contrasted with the predictions of the various source surface models which ignore magnetic stresses within 2.5 solar radii. Since the maximum south latitude was attained in September, 1994, Ulysses has traveled northward toward an ecliptic crossing in March 1995 and onward into the north solar hemisphere. The recent results will be presented and compared with those obtained in the southern hemisphere.

Smith, E. J.↗
NASA NTRS1994

Molecular Dynamics of a Water-Lipid Bilayer Interface

We present results of molecular dynamics simulations of a glycerol 1-monooleate bilayer in water. The total length of analyzed trajectories is 5ns. The calculated width of the bilayer agrees well with the experimentally measured value. The interior of the membrane is in a highly disordered fluid state. Atomic density profile, orientational and conformational distribution functions, and order parameters indicate that disorder increases toward the center of the bilayer. Analysis of out-of-plane thermal fluctuations of the bilayer surfaces occurring at the time scale of the present calculations reveals that the distribution of modes agrees with predictions of the capillary wave model. Fluctuations of both bilayer surfaces are uncorrelated, yielding Gaussian distribution of instantaneous widths of the membrane. Fluctuations of the width produce transient thinning defects in the bilayer which occasionally span almost half of the membrane. The leading mechanism of these fluctuations is the orientational and conformational motion of head groups rather than vertical motion of the whole molecules. Water considerably penetrates the head group region of the bilayer but not its hydrocarbon core. The total net excess dipole moment of the interfacial water points toward the aqueous phase, but the water polarization profile is non-monotonic. Both water and head groups significantly contribute to the surface potential across the interface. The calculated sign of the surface potential is in agreement with that from experimental measurements, but the value is markedly overestimated. The structural and electrical properties of the water-bilayer system are discussed in relation to membrane functions, in particular transport of ions and nonelectrolytes across membranes.

Wilson, Michael A.↗
NASA NTRS2004

Discrete Spring Model for Predicting Delamination Growth in Z-Fiber Reinforced DCB Specimens

Beam theory analysis was applied to predict delamination growth in DCB specimens reinforced in the thickness direction with pultruded pins, known as Z-fibers. The specimen arms were modeled as cantilever beams supported by discrete springs, which were included to represent the pins. A bi-linear, irreversible damage law was used to represent Z-fiber damage, the parameters of which were obtained from previous experiments. Closed-form solutions were developed for specimen compliance and displacements corresponding to Z-fiber row locations. A solution strategy was formulated to predict delamination growth, in which the parent laminate mode I fracture toughness was used as the criterion for delamination growth. The solution procedure was coded into FORTRAN 90, giving a dedicated software tool for performing the delamination prediction. Comparison of analysis results with previous analysis and experiment showed good agreement, yielding an initial verification for the analytical procedure.

Ratcliffe, James G.↗
NASA NTRS2004

Virtual Sensors: Using Data Mining to Efficiently Estimate Spectra

Detecting clouds within a satellite image is essential for retrieving surface geophysical parameters, such as albedo and temperature, from optical and thermal imagery because the retrieval methods tend to be valid for clear skies only. Thus, routine satellite data processing requires reliable automated cloud detection algorithms that are applicable to many surface types. Unfortunately, cloud detection over snow and ice is difficult due to the lack of spectral contrast between clouds and snow. Snow and clouds are both highly reflective in the visible wavelen,ats and often show little contrast in the thermal Infrared. However, at 1.6 microns, the spectral signatures of snow and clouds differ enough to allow improved snow/ice/cloud discrimination. The recent Terra and Aqua Moderate Resolution Imaging Spectro-Radiometer (MODIS) sensors have a channel (channel 6) at 1.6 microns. Presently the most comprehensive, long-term information on surface albedo and temperature over snow- and ice-covered surfaces comes from the Advanced Very High Resolution Radiometer ( AVHRR) sensor that has been providing imagery since July 1981. The earlier AVHRR sensors (e.g. AVHRR/2) did not however have a channel designed for discriminating clouds from snow, such as the 1.6 micron channel available on the more recent AVHRR/3 or the MODIS sensors. In the absence of the 1.6 micron channel, the AVHRR Polar Pathfinder (APP) product performs cloud detection using a combination of time-series analysis and multispectral threshold tests based on the satellite's measuring channels to produce a cloud mask. The method has been found to work reasonably well over sea ice, but not so well over the ice sheets. Thus, improving the cloud mask in the APP dataset would be extremely helpful toward increasing the accuracy of the albedo and temperature retrievals, as well as extending the time-series of albedo and temperature retrievals from the more recent sensors to the historical ones. In this work, we use data mining methods to construct a model of MODIS channel 6 as a function of other channels that are common to both MODIS and AVHRR. The idea is to use the model to generate the equivalent of MODIS channel 6 for AVHRR as a function of the AVHRR equivalents to MODIS channels. We call this a Virtual Sensor because it predicts unmeasured spectra. The goal is to use this virtual channel 6. to yield a cloud mask superior to what is currently used in APP . Our results show that several data mining methods such as multilayer perceptrons (MLPs), ensemble methods (e.g., bagging), and kernel methods (e.g., support vector machines) generate channel 6 for unseen MODIS images with high accuracy. Because the true channel 6 is not available for AVHRR images, we qualitatively assess the virtual channel 6 for several AVHRR images.

Srivastava, Ashok↗
NASA NTRS2004

Discrete Spring Model for Predicting Delamination Growth in Z-Fiber Reinforced DCB Specimens

Beam theory analysis was applied to predict delamination growth in Double Cantilever Beam (DCB) specimens reinforced in the thickness direction with pultruded pins, known as Z-fibers. The specimen arms were modeled as cantilever beams supported by discrete springs, which were included to represent the pins. A bi-linear, irreversible damage law was used to represent Z-fiber damage, the parameters of which were obtained from previous experiments. Closed-form solutions were developed for specimen compliance and displacements corresponding to Z-fiber row locations. A solution strategy was formulated to predict delamination growth, in which the parent laminate mode I critical strain energy release rate was used as the criterion for delamination growth. The solution procedure was coded into FORTRAN 90, giving a dedicated software tool for performing the delamination prediction. Comparison of analysis results with previous analysis and experiment showed good agreement, yielding an initial verification for the analytical procedure.

Ratcliffe, James G.↗
NASA NTRS2011

Crowd Sourcing for Challenging Technical Problems and Business Model

Crowd sourcing may be defined as the act of outsourcing tasks that are traditionally performed by an employee or contractor to an undefined, generally large group of people or community (a crowd) in the form of an open call. The open call may be issued by an organization wishing to find a solution to a particular problem or complete a task, or by an open innovation service provider on behalf of that organization. In 2008, the Space Life Sciences Directorate (SLSD), with the support of Wyle Integrated Science and Engineering, established and implemented pilot projects in open innovation (crowd sourcing) to determine if these new internet-based platforms could indeed find solutions to difficult technical challenges. These unsolved technical problems were converted to problem statements, also called "Challenges" or "Technical Needs" by the various open innovation service providers, and were then posted externally to seek solutions. In addition, an open call was issued internally to NASA employees Agency wide (10 Field Centers and NASA HQ) using an open innovation service provider crowd sourcing platform to post NASA challenges from each Center for the others to propose solutions). From 2008 to 2010, the SLSD issued 34 challenges, 14 externally and 20 internally. The 14 external problems or challenges were posted through three different vendors: InnoCentive, Yet2.com and TopCoder. The 20 internal challenges were conducted using the InnoCentive crowd sourcing platform designed for internal use by an organization. This platform was customized for NASA use and promoted as NASA@Work. The results were significant. Of the seven InnoCentive external challenges, two full and five partial awards were made in complex technical areas such as predicting solar flares and long-duration food packaging. Similarly, the TopCoder challenge yielded an optimization algorithm for designing a lunar medical kit. The Yet2.com challenges yielded many new industry and academic contacts in bone imaging, microbial detection and even the use of pharmaceuticals for radiation protection. The internal challenges through NASA@Work drew over 6000 participants across all NASA centers. Challenges conducted by each NASA center elicited ideas and solutions from several other NASA centers and demonstrated rapid and efficient participation from employees at multiple centers to contribute to problem solving. Finally, on January 19, 2011, the SLSD conducted a workshop on open collaboration and innovation strategies and best practices through the newly established NASA Human Health and Performance Center (NHHPC). Initial projects will be described leading to a new business model for SLSD.

Davis, Jeffrey R.↗
NASA NTRS2024

Assessing Coulomb Collision Contributions to Photoelectron Deceleration above Earth’s Polar Caps

NASA’s Endurance rocket mission resolved a subtle energy shift in photoelectrons escaping from Earth’s polar caps, from which the first successful measurement of Earth’s ambipolar electrostatic field was made. Preflight models predicted that coulomb collisions with thermal electrons could add an extra energy shift of up to 1 eV. However, no such additional shift was reported in the initial analysis. Here, we identify a previously unreported effect related to the spacecraft wake, which introduced a systematic +0.25eV energy shift that countered and masked the shift from coulomb collisions. After applying an in-flight calibration for this wake effect, we find a total photopeak energy shift consistent with the combined influence of the ambipolar field and the subtler -0.37eV coulomb shift, predicted by the fully kinetic STET model. Further analysis of these recalibrated measurements yields an ambipolar field estimate in close agreement with the initial analysis.

Glyn A Collinson↗
DOE OSTI2024

MicroBooNE investigations on the photon interpretation of the MiniBooNE low energy excess

The MicroBooNE experiment is a liquid argon time projection chamber with 85-ton active volume at Fermilab, operated from 2015 to 2020 to collect neutrino data from Fermilab's Booster Neutrino Beam. One of MicroBooNE's physics goals is to investigate possible explanations of the low-energy excess observed by the MiniBooNE experiment in $\nu_{\mu}\rightarrow \nu_{e}$ neutrino oscillation measurements. MicroBooNE has performed searches to test hypothetical interpretations of the MiniBooNE low-energy excess, including the underestimation of the photon background or instrinic $\nu_{e}$ background. This thesis presents MicroBooNE's searches for two neutral current (NC) single-photon production processes that contribute to the photon background of the MiniBooNE measurement: NC $\Delta$ resonance production followed by $\Delta$ radiative decay: $\Delta \rightarrow N\gamma$, and NC coherent single-photon production. Both searches take advantage of boosted decision trees to yield efficient background rejection, and a high-statistic NC $\pi^0$ measurement to constrain dominant background, and make use of MicroBooNE's first three years of data. The NC $\Delta \rightarrow N\gamma$ measurement yielded a bound on the $\Delta$ radiative decay process at 2.3 times the predicted nominal rate at 90\% confidence level(C.L.), disfavoring a candidate photon interpretation of the MiniBooNE low-energy excess as a factor of 3.18 times the nominal NC Δ radiative decay rate at the 94.8\% C.L. The NC coherent single photon measurement leads to the world's first experimental limit on the cross-section of this process below 1 GeV, of $1.49 \times 10^{-41} \text{cm}^2$ at 90\% C.L., corresponding to 24.0 times the nominal prediction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Searching for Neutrino Tridents in the NOvA Near Detector

This dissertation presents a search for neutrino trident production in the NOvA near detector through the coherent ``dimuon" channel: $\nu_\mu +\hspace{1pt}\text{X} \rightarrow \nu_\mu + \mu^- + \mu^+ +\hspace{1pt}\text{X}$. Trident production is a rare, purely electroweak process with sensitivity to physics beyond the Standard Model. The theoretical background, motivation for studying the process, and previous experimental measurements are reviewed. The analysis uses data collected by the NOvA near detector (ND) from Fermilab's Neutrinos at the Main Injector (NuMI) beam between November 2014 and February 2024, corresponding to an exposure of $25.5\times 10^{20}$ protons on target. The ND is a segmented tracking calorimeter located 800~m from the beam target, receiving neutrinos with a mean energy of 2~GeV. A multi-pass background reduction strategy is implemented, including the development of a novel dimuon-specific tracking technique. Trident candidates are identified using a boost ed decision tree classifier trained on simulated signal and background events. Limited background Monte Carlo statistics necessitate the use of functional fits to sideband data, which are extrapolated to estimate backgrounds in the signal region. The unblinded data contain 9 trident-like events, with an estimated background of 5.66 $\pm$ 5.15 events. This yields a best fit estimate of 3.34 tridents compared to the Standard Model prediction of 4.66. A profiled Feldman-Cousins method is used to determine a 90\% confidence interval of [0,9.1] on the number of signal events, corresponding to an upper limit of 1.95$\times$ the Standard Model prediction. This result represents the lowest energy search for trident events to date, and the first experimental contribution to the process in 27 years.

Bowles, Reed Scott [Indiana U.]↗
NASA NTRS1994

Hot proton anisotropies and cool proton temperatures in the outer magnetosphere

The hot protons of the outer magnetosphere typically exhibit a temperature anistropy such that T(sub perp)/T(sub parallel) greater than 1, where perpendicular and parallel symbols denote directions relative to the background magnetic field. If this anisotrpy is sufficiently large, the electomagneitc proton cyclotron anistropy instability may be excited. This instability is studied using linear Vlasov theory and one-dimensional hybrid simulations for a homogeneous plasma model representative of conditions in the outer magnetosphere with a hot anisotropic proton component (denoted by subscript h) and a cool, initially isotropic proton component (subscript c). Linear theory yields an instability threshold condition on the hot proton temperature anistropy where as the simulations imply an upper bound on T(sub perp h)/T(sub parallel h); both the threshold and the upper bound have similar scaling with the maximum growth rate gamma (sub m), the parallel beta of the hot component, beta(sub parallel h), and the relative density of the hot component n(sub h)/n(sub e). An anlysis of plasma observations from the Los Alamos magnetospheric plasma analyzer (MPA) in geosynchronous orbits finds that the maximum value of the hot proton temperature anisotropy approximately satisfies the predicted scaling with beta(sub parallel h) and nu(sub h)/n(sub e) and yields the proportionality factor that quantifies this upper bound. The simulations are also used to examine the heating of the cool proton cyclotron instability. The simulations yield a scaling for the dimensionless late-time cool proton average temperature T(sub c)/T(sub parallel h) as (n(sub h)/n(sub e))/beta(sub parallel h exp 0.5). Analysis of MPA data shows that the observed values of T(sub c)/T(sub parallel h) have similar scaling and again yield the proportionality factor which quantifies this relationship.

Gary, S. Peter↗
DOE OSTI2022

Small field campaign: aerosol – ice formation closure pilot study. Final Report

Prediction of atmospheric ice formation from aerosol particles by heterogeneous nucleation represents one of the grand challenges in atmospheric science. Our insufficient predictive understanding of primary ice formation is the reason that climate models typically do not include heterogeneous ice nucleation with subsequent effects on climate uncertainty. Mixed-phase clouds, where supercooled water droplets and ice crystals coexist play globally an important role regulating climate. This is especially the case for the Arctic region that experiences the greatest warming due to climate change compared to other regions in the world. Immersion freezing initiated by ice-nucleating particles (INPs) in supercooled water droplets is recognized as the dominant primary ice formation pathway in mixed-phase cloud regimes. For this reason, it is crucial to evaluate our capability to predict immersion freezing for a given ambient aerosol population. The goal of this project is conducting a field-based pilot study at the U.S. DOE Atmospheric Radiation Measurement (ARM) user facility at Southern Great Plains (SGP) to evaluate our capability to predict the number concentration of aerosol particles that serve as INPs in the immersion freezing mode. Successful prediction of INP number concentrations is also termed “closure”. This field-observational approach represents a first-of-its kind attempt of an aerosol–ice formation closure study (AEROICESTUDY). Very few closure studies related to INPs have been conducted, and to our knowledge, none using robust size-resolved ambient aerosol composition measurements as a starting point. Achievement of aerosol–ice formation closure relies on our ability to characterize the ambient aerosol population with respect to particles size and composition and to determine INP number concentrations for specified freezing temperatures. This requires numerous online and offline instrumentation resulting in this pilot field campaign being a multi-institutional and community-collaborative effort. We chose the ARM SGP megasite for this first aerosol-ice formation closure pilot study due to its significant measurement capabilities available to obtain detailed physical characterization of the local aerosol population including size distribution, mass loading, and chemical composition of non-refractory aerosol particles. The overall objective of this project is to identify ice nucleation parameterizations that produce the most robust predictions of INP numbers and thus are best suited to be included in cloud and climate models. This objective includes the following goals for this field and laboratory-based project: i) What are the crucial aerosol physicochemical property measurements needed to accurately guide ice nucleation representations in models and long-term INP measurements? ii) What level of parameter details needs to be known to achieve aerosol–ice formation closure? iii) What are the leading causes for climate model bias in INP predictions? We found that the advances in our understanding of immersion freezing garnered over the last 20 years allowed us to yield partial and full closures of atmospheric immersion freezing from ambient aerosol particles. When the aerosol population is physicochemically complex and parameterizations for representative INP types are not yet available, we still struggle to accurately predict INP number concentrations. This project clearly demonstrates that with more laboratory and field measurements that are accompanied by particle composition analysis, the necessary datasets to achieve aerosol–ice formation closure for various locations will emerge, thus providing a robust foundation for guiding the representation of INPs in cloud and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Primary and Low-Strain Creep Models for 9Cr Tempered Martensitic Steels Including the Effects of Irradiation Softening and High-Helium Re-Hardening

Primary and low-strain creep represents a very important integrity challenge to large, complex structures, like fusion reactors. Here, we develop a predictive empirical primary creep model for 9Cr tempered martensitic steels (TMS), relating the applied stress (σ) to strain (ε), time (t) and temperature (T). The most accurate model is based on the applied σ normalized by the steel’s T-dependent ultimate tensile stress (σo), σ/σo(T). The model, fit to 17 heats of 9Cr TMS, yielded a σ root mean square error (RMSE) of ≈±11 MPa. Notably, the model also provides robust predictions for all the other TMS, when calibrated only by the fusion candidate Eurofer97 database. The model was extended to explore two possible effects of neutron irradiation, which produces both displacements per atom (dpa) and helium (He in atomic parts per million, appm) damage. These effects, which have not been previously considered, include: (a) softening, as a function of dpa, at T > ≈400–450 °C, in low-He fission environments (<1 He/dpa); and (b) subsequent re-hardening in high-He (≥10 He/dpa) fusion first-wall environments. The irradiation effect models predict (a) accelerated primary creep due to irradiation softening; and (b) fully arrested creep due to high-He re-hardening.

Alam, Md Ershadul (ORCID:0000000345968026)↗
NASA NTRS1992

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗
DOE OSTI2022

Improvement of eukaryotic protein predictions from soil metagenomes

During the last decades, metagenomics has highlighted the diversity of microorganisms from environmental or host-associated samples. Most metagenomics public repositories use annotation pipelines tailored for prokaryotes regardless of the taxonomic origin of contigs. Consequently, eukaryotic contigs with intrinsically different gene features, are not optimally annotated. Using a bioinformatics pipeline, we have filtered 7.9 billion contigs from 6,872 soil metagenomes in the JGI’s IMG/M database to identify eukaryotic contigs. We have re-annotated genes using eukaryote-tailored methods, yielding 8 million eukaryotic proteins and over 300,000 orphan proteins lacking homology in public databases. Comparing the gene predictions we made with initial JGI ones on the same contigs, we confirmed our pipeline improves eukaryotic proteins completeness and contiguity in soil metagenomes. The improved quality of eukaryotic proteins combined with a more comprehensive assignment method yielded more reliable taxonomic annotation. This dataset of eukaryotic soil proteins with improved completeness, quality and taxonomic annotation reliability is of interest for any scientist aiming at studying the composition, biological functions and gene flux in soil communities involving eukaryotes.

54 ENVIRONMENTAL SCIENCES↗