Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOMOGENIZATION METHODS”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

SURFACE EXOSURE TIMESCALES OF APOLLO CORE SAMPLE 73002 SPACE WEATHERED GRAINS

Introduction: Space weathering causes the surface regolith on airless bodies like the Moon to be morphologically, microstructurally, and chemically altered due to micrometeoroid bombardment and solar wind irradiation[1]. These processes produce a multitude of microstructural and chemical changes in individual soil grains that accumulate as grains are exposed on the surface over time. One microstructural feature produced by exposure to the solar wind is the development of ion-damaged rims on regolith grains via H+ and He+ ion irradiation. Also present are solar energetic particle (SEP) tracks, which are nanoscale lineations of ionization damage within grain interiors formed by high energy solar flare ions(primarily Fe group nuclei) which penetrate millimeters below the surface [2]. Recent work has confirmed that the thickness of solar wind-damaged amorphous rims on anorthite grains, nanocrystalline rims on olivine grains, and their respective SEP track densities are correlated with their surface exposure ages [3].Core sample 73002, recently released under the Apollo Next Generation Sample Analysis (ANGSA) Program, has provided an opportunity to study material collected from the light mantle formation during Apollo 17. The light mantle is thought to have been deposited via a landslide originating from the neighboring South Massif [4]. Spectral profiles and ferromagnetic resonance measurements of bulk soils sampled at cm-intervals from 73002 indicate the existence of mature regolith extending ~8 cm below the surface, suggesting the presence of an in-situ reworking zone in the core[5,6].Here we present the distribution of SEP track densities, solar wind damaged rim widths, and the corresponding surface exposure ages of grains residing in the proposed in-situ reworking zone of core 73002. Methods: Bulk samples of regolith from the first eight intervals and every following fourth interval down the core (dissection Pass 2) were delivered to Purdue University as <45 μm size fractions. The first eight0.5 cm intervals, representing the top 4cmof regolith, were individually dry sieved to a<20 μm size fraction, and grains were prepared by ultramicrotomy for analysis in the scanning transmission electron microscope (STEM). Bright field (BF) and dark field (DF) STEM images of solar flare tracks and solar wind damaged rims were acquired on a JEOL 2500SE TEM, equipped with a 60 mm2ultra-thin window silicon drift energy dispersive X-ray (EDX) spectrometer at NASA John-son Space Center. Grain compositions were determined by EDX compositional spectrum imaging. SEP track densities and amorphous rim thicknesses were measured on BF and DF images. Results: The fraction of grains with discernable space weathering features decreased with depth, with 100% in the first interval exhibiting space weathering features as opposed to ~60% in Interval 8.BF and DF STEM images show splash-melt and vapor deposited rims on outermost grain margins with embedded Fe-bearing nanoparticles ranging in size up to ~10 nm in diameter. Images also show solar-wind damaged rims below the vapor deposits, and SEP tracks present in grain interiors. The rim thicknesses and SEP track densities were determined using the methods of [3] for 54grains spanning the top 4 cm. ~80% of the grains analyzed were anorthite and ~20% were olivine. The track production rate from[3] was used to estimate the surface exposure times of the grains in this study. The majority of exposure times are in the 1 -5 MY range for all intervals. The lowest value was 3.8x 10¬5years in Interval 8and the highest value was1.1x 107years in Interval 6. Grains with the highest solar wind damage rim thicknesses and SEP track densities are found in Interval 6.The highest rim thickness and track density observed for anorthiteis126.5 nm and4.7 x 1011tracks/cm2, while for olivine they are 98.6 nm and 1.83 x 1011tracks/cm2, respectively. Discussion: Intervals from the top 4 cm of 73002 show similar SEP track density and solar wind damaged rim thickness distributions. Spectral and ferromagnetic resonance measurements of 73002 indicate an in-situ reworking zone up to~8 cm below the surface[5,6].The depth analyzed in this work is within the core’s uppermost in-situ reworking zone and may serve as an explanation to the homogenous distribution of exposure ages. The decreasing abundance of space weathered features with depth is consistent with regolith mixing models showing an exponential decay in regolith maturity with depth over time[7].A maximum surface exposure age determined via SEP track density is approximately 10million years. This timescale is consistent with the lower estimate of surface exposure ages of the light mantle determined in previous studies, which can range from 10s to over 100 Ma [4].It also falls within the core regolith reworking development timescale of ~17 Ma [5,6]. Additional analysis will be performed for a further seven intervals in 73002.

J A McFadden↗
DOE OSTI2020

Methodology for inferring reactor core power distribution from an optical fiber based gamma thermometer array

Here, a 4-step data analytic methodology has been devised for the purpose of inferring the distribution of power in a reactor core, based on the response of an array of optical fiber based gamma thermometers (OFBGTs). This data analytic methodology is crucial for the development of a system of OFBGTs for the purpose of calibrating local power range monitors in boiling water reactors. Such a system would be an improvement to the present calibration system in boiling water reactors, in terms of safety, efficiency, and permanence. The first step of this methodology is to establish an energy balance method. In this method, one uses MCNP to determine response functions, which allow one to convert from gamma thermometer response to power. The gamma thermometers and the reactor core are segmented, such that each gamma thermometer segment provides an estimate of power for each reactor core segment. The estimates of power for the reactor core segments (hereafter referred to as fuel assembly segments) are calculated based on the measured response of the gamma thermometers. The second step of the methodology is to employ a weighting scheme to combine the estimates of the power of the various fuel assembly segments, based on the response of the various OFGBT segments and the incremental dose rates of the various fuel assembly segments to the OFBGT segments. The third step is to iteratively calculate the estimates of the power of the fuel assembly segments, until a convergence criterion is met, which indicates that the calculation has converged. The fourth and final step of the methodology is to estimate the uncertainty of the power for each fuel assembly segment. The mathematical basis for this step is not the focus of this paper. We have used a 3-D homogeneous reactor model to demonstrate the data analytic methodology; and have found that the data analytic methodology operates as intended.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Silica Supported Organometallic Ir I Complexes Enable Efficient Catalytic Methane Borylation

Catalytic C–H borylation is an attractive method for the conversion of the most abundant hydrocarbon, methane (CH 4 ), to a mild nucleophilic building block. However, existing CH 4 borylation catalysts often suffer from low turnover numbers and conversions, which is hypothesized to result from inactive metal hydride agglomerates. Herein we report that the heterogenization of a bisphosphine molecular precatalyst, [(dmpe)Ir(cod)CH 3 ], onto amorphous silica dramatically enhances its performance, yielding a catalyst that is 12-times more efficient than the current standard for CH 4 borylation. The catalyst affords over 2000 turnovers at 150 °C in 16 h with a selectivity of 91.5% for mono- vs diborylation. Higher catalyst loadings improve yield and selectivity for the monoborylated product (H 3 CBpin) with 82.8% yield and >99% selectivity being achieved with 1255 turnovers. X-ray absorption and dynamic nuclear polarization-enhanced solid-state NMR spectroscopic studies identify the supported precatalyst as an Ir I species, and indicate that upon completion of catalysis, multinuclear Ir polyhydrides are not formed. This is consistent with the hypothesis that immobilization of the organometallic Ir species on a surface prevents bimolecular decomposition pathways. Importantly, immobilization of the homogeneous Ir I fragment onto amorphous silica represents a unique and simple strategy to improve the TON and longevity of a CH 4 borylation catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2013

Conductive Carbon Nanotube Inks for Use with Desktop Inkjet Printing Technology

Inkjet printing is a common commercial process. In addition to the familiar use in printing documents from computers, it is also used in some industrial applications. For example, wire manufacturers are required by law to print the wire type, gauge, and safety information on the exterior of each foot of manufactured wire, and this is typically done with inkjet or laser printers. The goal of this work was the creation of conductive inks that can be applied to a wire or flexible substrates via inkjet printing methods. The use of inkjet printing technology to print conductive inks has been in testing for several years. While researchers have been able to get the printing system to mechanically work, the application of conductive inks on substrates has not consistently produced adequate low resistances in the kilohm range. Conductive materials can be applied using a printer in single or multiple passes onto a substrate including textiles, polymer films, and paper. The conductive materials are composed of electrical conductors such as carbon nanotubes (including functionalized carbon nanotubes and metal-coated carbon nanotubes); graphene, a polycyclic aromatic hydrocarbon (e.g., pentacene and bisperipentacene); metal nanoparticles; inherently conductive polymers (ICP); and combinations thereof. Once the conductive materials are applied, the materials are dried and sintered to form adherent conductive materials on the substrate. For certain formulations, increased conductivity can be achieved by printing on substrates supported by low levels of magnetic field alignment. The adherent conductive materials can be used in applications such as damage detection, dust particle removal, smart coating systems, and flexible electronic circuitry. By applying alternating layers of different electrical conductors to form a layered composite material, a single homogeneous layer can be produced with improved electrical properties. It is believed that patterning alternate layers of different conductors may improve electrical pathways through alignment of the conductors and band gap optimization. One feature of this innovation is that flexible conductive traces could be accomplished with a conductive ink having a surface resistivity of less than 10 ohms/square. Another result was that a composite material comprising a mixture of carbon nanotubes and metallic nanoparticles could be applied by inkjet printing to flexible substrates, and the resulting applied material was one to two orders of magnitude more conductive than a material made by printing inks containing carbon nanotubes alone.

Roberson, Luke↗
DOE OSTI2021

Discussion on molecular dynamics (MD) simulations of the asphalt materials

The application of asphalt materials in pavement engineering has been increasingly widespread and sophisticated over the past several decades. Variations in the properties of asphalt binder during mixing, transportation, and paving can affect the performance of asphalt pavement. However, the asphalt material is a non-homogeneous and complex organic substance, consisting of various molecules with widely various molecular weights, elemental compositions, and structures. This complexity leads to difficulties for researchers to clearly and immediately understand the properties of asphalt materials and their variations. The multi-scale research approach combines macroscopic experimental data and microscopic simulation results from a practical engineering perspective. It helps to improve the understanding of asphalt materials. The molecular dynamics (MD) simulation proposes a corresponding molecular model of asphalt material based on experimental data, and the simulation algorithm is able to derive properties similar to those of real asphalt. Here, this paper provides a comprehensive review of the current studies on MD simulation of asphalt materials, including modeling, properties, and multi-scale analysis. As a key part of the computational simulation, this paper discusses the typical asphalt binder and asphalt-aggregate interface models constructed by different groups, and also presents their differences from real samples and their feasibility based on fundamental properties. After the introduction of molecular models, the extensive work made by researchers based on molecular models is categorically reviewed and discussed. The strengths and weaknesses of MD simulation methods in the study of asphalt materials are also summarized in order to provide the reader with a more comprehensive understanding of the relevant contents and to guide subsequent research.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Lifetime Predictions of a Titanium Silicate Glass with Machined Flaws

A dynamic fatigue study was performed on a Titanium Silicate glass to assess its susceptibility to delayed failure and to compare the results with those of a previous study. Fracture mechanics techniques were used to analyze the results for the purpose of making lifetime predictions. The material strength and lifetime was seen to increase due to the removal of residual stress through grinding and polishing. Influence on time-to-failure is addressed for the case with and without residual stress present. Titanium silicate glass otherwise known as ultra-low expansion (ULE)* glass is a candidate for use in applications requiring low thermal expansion characteristics such as telescope mirrors. The Hubble Space Telescope s primary mirror was manufactured from ULE glass. ULE contains 7.5% titanium dioxide which in combination with silica results in a homogenous glass with a linear expansion coefficient near zero. delayed failure . This previous study was based on a 230/270 grit surface. The grinding and polishing process reduces the surface flaw size and subsurface damage, and relieves residual stress by removing the material with successively smaller grinding media. This results in an increase in strength of the optic during the grinding and polishing sequence. Thus, a second study was undertaken using samples with a surface finish typically achieved for mirror elements, to observe the effects of surface finishing on the time-to-failure predictions. An allowable stress can be calculated for this material based upon modulus of rupture data; however, this does not take into account the problem of delayed failure, most likely due to stress corrosion, which can significantly shorten lifetime. Fortunately, a theory based on fracture mechanics has been developed enabling lifetime predictions to be made for brittle materials susceptible to delayed failure. Knowledge of the factors governing the rate of subcritical flaw growth in a given environment enables the development of relations between lifetime, applied stress and failure probability for the material under study. Dynamic fatigue is one method of obtaining the necessary information to develop these relationships. In this study, the dynamic fatigue method was used to construct a time-to-failure diagram for polished ULE glass.

Tucker, Dennis S.↗
DOE OSTI2021

Experimental and computational analysis of bending fatigue failure in chopped carbon fiber chip reinforced composites

With a better balance among good mechanical performance, high freedom of design, and low material and manufacturing cost, chopped carbon fiber chip reinforced sheet molding compound (SMC) composites show great potential in different engineering applications. Here in this paper, bending fatigue behaviors of SMC composites considering the heterogeneous fiber orientation distributions have been thoroughly investigated utilizing both experimental and computational methods. First, four-point bending fatigue tests are performed with designed SMC composites, and the local modulus is adopted as a metric to represent the local fiber orientation of two opposing sides. Interestingly, SMC composites with and without large discrepancy in local modulus of opposing sides show different fatigue behaviors. Interrupted tests are conducted to explore the bending fatigue failure mechanism, and the damage processes of valid specimens are also closely examined. We find that the fatigue failure of SMC composites under four-point bending is governed by crack propagation instead of crack initiation. Because of this, the heterogeneous local fiber orientations of both sides of the specimen influence fatigue life. The microstructure of the lower side shows a direct influence while that of the upper side also exhibiting influence which becomes more prominent for high cycle fatigue cases. Furthermore, a hybrid micro–macro computational model is proposed to efficiently study the cyclic bending behavior of SMC composites. The region of interest is reconstructed with a modified random sequential absorption algorithm to conserve all the microstructural details including the heterogeneous fiber orientation, while the rest of the regions are modeled as homogenized macro-scale continua. Combined with a framework to capture the progressive fatigue damage under cyclic bending, the bending fatigue behaviors of SMC composites are accurately captured by the hybrid computational model comparing with our experimental analysis.

36 MATERIALS SCIENCE↗
NASA NTRS2022

An Evaluation of the Liquid Cloud Droplet Effective Radius Derived from MODIS, Airborne Remote Sensing, and in Situ Measurements from CAMP2Ex

The cloud drop effective radius (Re) of the drop size distribution derived from passive satellite sensors is a key variable used in climate research. Validation of these satellite products has often taken place under stratiform cloud conditions that favor the assumption of cloud horizontal homogeneity used by the retrieval techniques. However, many studies have noted concerns with respect to significant biases in retrieved Re arising from cloud heterogeneity, for example, in cumulus cloud fields. Here, we examine data collected during the 2019 “Cloud, Aerosol and Monsoon Processes Philippines Experiment” (CAMP2Ex), which, in part, targeted the objective of providing the first detailed evaluation of Re retrieved across multiple platforms and techniques in a cumulus and congestus cloud region. Our evaluation consists of cross-comparisons of Re between the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra satellite, the Research Scanning Polarimeter (RSP) onboard the NASA P-3 aircraft, and in situ measurements from both the NASA P-3 and Learjet aircraft that are all taken in close spatiotemporal proximity to the same cloud fields. A particular advantage of our approach lies in the capability of the RSP to retrieve Re using a bi-spectral MODIS approach and a polarimetric approach, which allows for the evaluation of bi-spectral and polarimetric Re retrievals from an airborne perspective using the same samples. Averaged over all P-3 flight segments examined here for warm clouds, the RSP polarimetric method, the in situ method, and the bias-adjusted MODIS method of Fu et al. (2019) show a comparable median (mean ± standard deviation) for the Re samples of 9.6 (10.2 ± 4.0) µm, 11.0 (13.6 ± 11.3) µm, and 10.4 (10.8 ± 3.8) µm, respectively. These values are far lower than the values of 15.1 (16.2 ± 5.5) µm and 17.2 (17.7 ± 5.7) µm from the bi-spectral retrievals of RSP and MODIS, respectively. Similar results are observed when Re is segregated by cloud-top height and in detailed case studies. The clouds sampled during CAMP2Ex consist of mostly small (mean transect length ∼ 1.4 km) and low clouds (mean cloud-top height ∼ 1 km), which had more numerous small clouds than the trade wind cumuli sampled in past field campaigns such as Rain in Shallow Cumulus over the Ocean (RICO) and the Indian Ocean Experiment (INDOEX). The overestimates of Re from the RSP bi-spectral technique compared with the polarimetric technique increased as cloud size and cloud optical depth decreased. Drizzle, cloud-top bumpiness, and solar zenith angle, however, are not closely correlated with the overestimate of bi-spectral Re. For shallow clouds that dominated the liquid cloud cover for the CAMP2Ex region and period, we show that 3-D radiative transfer and cloud heterogeneity, particularly for the optically thin and small clouds, appear to be the leading cause of the large positive biases in bi-spectral retrievals. Because this bias varies with the underlying structure of the cloud field, caution continues to be warranted in studies that use bi-spectral Re retrievals in cumulus cloud fields.

Cloud microphysics↗
DOE OSTI2026

Achieving high rate performance in hybrid pristine-recycled cathodes using model-informed electrode designs

Direct recycling lithium-ion battery cathodes, a process that retains the engineered oxide structures from end-of-life materials, presents a cost-effective and energy-efficient alternative to other battery recycling methods. However, while direct-recycled cathodes have demonstrated performance comparable to that of pristine materials at low cycling rates, their high-rate performance remains uncertain. Morphology changes in cathode particles, a main mode of degradation, directly impact rate performance by limiting surface kinetics and solid-phase diffusion. If direct recycling processes do not sufficiently restore pristine-like morphologies, the recycled materials may retain structural defects that hinder high-rate performance. The present work uses a physics-based pseudo-2D model to simulate hybrid electrodes with pristine and artificially “aged/recycled” NMC materials to investigate potential impacts of incorporating performance-limited aged cathode materials into cells. The study highlights how differences in transport and kinetic properties can influence rate capabilities in mixed electrodes — particularly in high-loading cells in high-demand applications. However, model results also reveal a possible mitigation strategy via dual-layer electrode architectures with lower-performing materials positioned near the current collector. Simulations of 4.0 mAh cm −2 cells cycled at 4C using a dual-layer architecture provided approximately 5%–30% more capacity in constant-current protocols compared to homogeneously blended electrode architectures with the same loadings and mixed-material compositions. These findings highlight the importance of strategic electrode design in minimizing potential performance losses and facilitating the integration of recycled materials into high-performance batteries, advancing sustainable and cost-effective battery manufacturing.

25 ENERGY STORAGE↗
DOE OSTI2025

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Molecular-Level Design of Heterogeneous Chiral Catalysts

A recent major focus of catalyst research has been on improving the selectivity or atom efficiency of catalytic processes. Perhaps the most subtle and difficult form of selectivity to control is enantioselectivity, a catalyst’s ability to preferentially yield one enantiomer of a chiral product. Since two enantiomers are identical in all respects except for their mirror reflection symmetry, both are formed in non-chiral environments at exactly equal rates; only by introducing some type of chiral modification to the system can the rate of formation of one enantiomer be favored over the other. This has been successfully achieved for homogeneous-phase, organometallic catalysts, but the empirical development of heterogeneous-phase analogs has proven elusive. The design of such heterogeneous-phase catalysts would obviate the need for inefficient purification steps to remove the heavy-metal catalysts. The goal of the proposed research by the Chiral Catalysis Group is to understand the interactions of prochiral reactants on chirally modified surfaces that lead to the preferential formation of one enantiomer over the other. There are two general concepts for designing heterogeneous-phase enantioselective catalysts. The first adopts the strategy used in homogeneous phase, by adding a chiral modifier to the surface. This approach, while conceptually simple has challenges associated with modifying an extended surface. The second concept, which is unique to extended surfaces, is the idea of creating surfaces that are themselves naturally chiral. A fundamental understanding of the interactions that control enantioselectivity will enable the Chiral Catalysis Group to develop the concepts necessary to design novel, heterogenous catalytic materials. Enantioselectivity is also central to catalytic processes that underpin the pharmaceutical and agrochemical industries. With funding from the Department of Energy, the Chiral Catalysis Group has played a leading role in developing a fundamental understanding of chiral catalysis and of the adsorbate structures, interactions and kinetic effects that control heterogenous-phase enantioselectivity. The proposed research program is organized conceptually around three types of enantiospecific interactions: interactions with naturally chiral surfaces; those with modified surfaces categorized as chirally templated surfaces; and one-to-one interactions of prochiral molecules with chiral surface modifiers. Our specific goals are to: (1) determine the structures of chiral surfaces, prochiral reactants, chiral products and chiral modifiers; (2) understand the enantiospecific interactions that occur between prochiral reactants and chiral products with chirally modified surfaces; (3) develop methods for detection of enantiospecific surface phenomena; and (4) elucidate enantioselective reaction pathways and the interactions that control them. The range of expertise within the Chiral Catalysis Group includes molecular-scale theory, surface science experiments, state-of-the-art spectroscopic, kinetic and atomic-scale imaging capabilities carried out in vacuum and under more-realistic liquid- and gas-phase environments complemented by catalytic measurements. By fully integrating such a breadth of approaches and extensive expertise, the Chiral Catalysis Group will, for the first time, be able to understand, rationally design and then control active sites for optimal enantioselectivity. In addition to having a significant impact on the field of chiral catalysis, the proposed research program will further our general understanding of the molecular-scale origins of catalytic selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

DNS and LES of a Shear-Free Mixing Layer

The purpose of this work is twofold. First, given the computational resources available today, it is possible to reach, using DNS, higher Reynolds numbers than in Briggs et al.. In the present study, the microscale Reynolds numbers reached in the low- and high-energy homogeneous regions are, respectively, 32 and 69. The results reported earlier can thus be complemented and their robustness in the presence of increased turbulence studied. The second aim of this work is to perform a detailed and documented LES of the shear-free mixing layer. In that respect, the creation of a DNS database at higher Reynolds number is necessary in order to make meaningful LES assessments. From the point of view of LES, the shear-free mixing-layer is interesting since it allows one to test how traditional LES models perform in the presence of an inhomogeneity without having to deal with difficult numerical issues. Indeed, as argued in Briggs et al., it is possible to use a spectral code to study the shear-free mixing layer and one can thus focus on the accuracy of the modelling while avoiding contamination of the results by commutation errors etc. This paper is organized as follows. First we detail the initialization procedure used in the simulation. Since the flow is not statistically stationary, this initialization procedure has a fairly strong influence on the evolution. Although we will focus here on the shear-free mixing layer, the method proposed in the present work can easily be used for other flows with one inhomogeneous direction. The next section of the article is devoted to the description of the DNS. All the relevant parameters are listed and comparison with the Veeravalli & Warhaft experiment is performed. The section on the LES of the shear-free mixing layer follows. A detailed comparison between the filtered DNS data and the LES predictions is presented. It is shown that simple eddy viscosity models perform very well for the present test case, most probably because the flow seems to be almost isotropic in the small-scale range that is not resolved by the LES.

Knaepen, B.↗
DOE OSTI2024

Atmospheric parameters and chemical abundances within 100 pc: a sample of G, K, and M main-sequence stars

ABSTRACT To date, we have access to enormous inventories of stellar spectra that allow the extraction of atmospheric parameters and chemical abundances essential in stellar studies. However, characterizing such a large amount of data is complex and requires a good understanding of the studied object to ensure reliable and homogeneous results. In this study, we present a methodology to measure homogenously the basic atmospheric parameters and detailed chemical abundances of over 1600 thin disc main-sequence stars in the 100 pc solar neighbourhood, using APOGEE-2 infrared spectra. We employed the code tonalli to determine the atmospheric parameters using a prior on $\log {g}$. The $\log {g}$ prior in tonalli implies an understanding of the treated population and helps to find physically coherent answers. Our atmospheric parameters agree within the typical uncertainties (100 K in $\mathrm{T_{eff}}$, 0.15 dex in $\log {g}$ and [M/H]) with previous estimations of ASPCAP and Gaia DR3. We use our temperatures to determine a new infrared colour–temperature sequence, in good agreement with previous works, that can be used for any main-sequence star. Additionally, we used the bacchus code to determine the abundances of Mg, Al, Si, Ca, and Fe in our sample. The five elements (Mg, Al, Si, Ca, Fe) studied have an abundance distribution centred around slightly subsolar values in agreement with previous results for the solar neighbourhood. The over 1600 main-sequence stars’ atmospheric parameters and chemical abundances presented here are useful in follow-up studies of the solar neighbourhood or as a training set for data-driven methods.

López-Valdivia, Ricardo (ORCID:0000000277950018)↗
NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗
DOE OSTI2021

Tools for uranium characterization in carbonate samples: case studies of natural U–Pb geochronology reference materials

Laser ablation U–Pb analyses of carbonate (LAcarb) samples has greatly expanded the potential for U–Pb dating to a variety of carbonate-producing settings. Carbonates that were previously considered impossible to date using isotope dilution methods may preserve radiogenic domains that can be dated using spatially resolved laser ablation geochronology techniques. Work is ongoing to identify reference materials and to consider best practices for LAcarb. In this study we apply standard and emerging characterization tool sets on three natural samples with the dual goal of enhancing the study of carbonates and establishing a new set of well-characterized natural reference materials for LAcarb studies. We start with the existing carbonate reference material WC-1 from the Permian Reef Complex of Texas, building on the published description to offer a deeper look at U and associated trace elements. We consider a tufa sample from the Miocene Barstow Formation of the Mojave Block, California, as a possible secondary calcite reference material due to its well-behaved U–Pb systematics. There are currently no natural dolomite standards. We present an unusual dolomite sample with very well-behaved U–Pb systematics from the Miocene of the Turkana Basin of Kenya as a possible dolomite reference material for LAcarb dating. In addition to using X-ray fluorescence (XRF) mapping and spectroscopy to better understand U in these natural samples, we have analyzed multiple aliquots of each of them for 87 Sr/ 86 Sr by thermal ionization mass spectrometry (TIMS). The Sr isotope compositions are analytically homogeneous within petrographically homogeneous regions of all three samples, and thus these materials could be used as Sr isotope standards as well. While not part of the current contribution, this combination could streamline simultaneous LA analyses of 87 Sr/ 86 Sr and U–Pb geochronology.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data-driven selection of stiff chemistry ODE solver in operator-splitting schemes

Most computational fluid dynamics simulations of practical combustion applications employ operator-splitting schemes, where chemistry and transport are separated and integrated with distinct numerical methods. The changes in composition due to chemistry are evaluated by solving ordinary differential equations (ODE) in each cell of the computational domain, which typically dominates the computational cost when detailed chemistry is considered. In this work, a data-driven approach for the selection of chemistry ODE solvers in operator-splitting schemes is presented. Neural networks are used to predict the ODE solvers CPU times and errors for a given thermochemical state. This allows the selection of an optimal ODE solver on a cell-by-cell, timestep-by-timestep basis. The models are trained using a wide set of thermochemical states generated through partially-stirred reactors and flames simulations. The methodology is validated by quantifying the prediction errors, the classification accuracy, and the computational speedup. The model predicts the optimal ODE solver for 70 to 95% of the validation cases and decreases the computional cost by a factor of 3 or more. The generalizability of the methodology to different chemical mechanisms and different fuels is assessed and it is shown that the model’s performance is only slightly degraded and its applicability is significantly enhanced if the inputs to the neural networks are restricted to a small set of thermochemical state variables present in most chemical mechanisms. In conclusion, the models are used in an homogeneous reactor case and a multi-dimensional CFD simulation of a diesel spray at high pressure where a speedup of more than 3 is achieved.

42 ENGINEERING↗
DOE OSTI2020

Reionization history constraints from neural network based predictions of high-redshift quasar continua

ABSTRACT Observations of the early Universe suggest that reionization was complete by z ∼ 6, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Lyα damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Lyα, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening intergalactic medium (IGM). In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Lyα. Our Quasar Spectra from Artificial Neural Network based predictive Regression Algorithm(QSANNdRA) improves the error in this reconstruction compared to the state-of-the-art principal component analysis (PCA) based model in the literature by 14.2 per cent on average, and provides an improvement of 6.1 per cent on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNdRA further achieves an improvement of 22.1 per cent and 16.8 per cent when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at z = 7.0851 and ULAS J1342+0928 at z = 7.5413, respectively. Using our more accurate reconstructions of these two z > 7 quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find $\bar{x}_\mathrm{H\, \small{I}} = 0.25^{+0.05}_{-0.05}$ at z = 7.0851 and $\bar{x}_\mathrm{H\, \small{I}} = 0.60^{+0.11}_{-0.11}$ at z = 7.5413. Our results are consistent with the literature and favour a rapid end to reionization.

Ďurovčíková, Dominika↗