Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space 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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Development of Scalable Coastal and Offshore Macroalgal Farming

The US DOE/ARPA-E MARINER program funded a four-year project to determine an optimal way to grow kelps in large, nearshore and offshore arrays for the eventual purpose of biofuel production with the goal of keeping the cost below $\$$80 USD per dry metric ton of kelp. This project specifically looked at how Saccharina latissima can be grown in the Gulf of Alaska to reach that goal. There were three major aspects of the research: 1) optimizing nursery production and seeding lines for outplanting, 2) designing an economical, modular outplanting structure and 3) developing methods to efficiently harvest the product. Farm designs were based on catenary structures and the use of spreader bars with variable spacing of grow lines and line types. The spacing of the grow lines makes a difference in the yield. Grow line spacing of ≥1.5m showed about a 50% increase in kg/m. There was no statistical difference in the growth of Saccharina latissima whether in the middle or the outside of the array, but the line type and perhaps line thickness can make a difference in yield. Sagging caused by the weight of the mature fronds resulted in lower growth at depth. Various harvesting approaches for mature kelps were tested by collaborating farmers. One promising innovation is the use of large bags for holding the freshly harvested fronds. Although the weight of the fronds on the growlines causes the lines to sink, the bags packed with the harvested fronds float, allowing for easy loading to the transfer vessel. Another advance in harvesting is a specially built harvest vessel, the Harvest Buddy, allowing a more mechanized and faster way to harvest. A significant aspect of this project was the TTO/T2M. Several different outreach activities were performed by the Alaska Fisheries Development Foundation and GreenWave. A techno-economic assessment (TEA) using our data has pointed to solutions to reach the $\$$80 goal. A second phase of this project involved the co-culture of two different species of kelps. Interest in farming kelps has grown beyond using kelp for food, feed or biofuels. There is considerable interest in generating biomass from seaweed for use in bioplastics and other products that would substitute for petroleum-derived products. For these uses to be viable, large amounts of biomass are needed. Very large kelp farms can be expensive to build and maintain, leading to the need to optimize the biomass per unit area. Although close spacing of growlines can lead to poor growth, a viable approach may be to grow two species of kelps together: one that hangs down and one that is buoyant, growing up. This system would increase the spacing in three dimensions. In Alaska, Saccharina latissima is commonly grown hanging down from longlines. One of the buoyant Alaskan kelps is Nereocystis luetkeana. Because there are commercial uses for wild-harvested Nereocystis in Alaska, we undertook a preliminary trial in Kodiak, Alaska that grew both Saccharina and Nereocystis in the same longline array. Closely spaced lines were seeded the first week of February 2023 and set at 3m below the surface. The arrays were harvested in late June 2023. Despite having 45% fewer grow-lines, the total yield of the Nereocystis on the combined arrays was statistically similar to the Nereocystis only arrays. Total yields were greatest on the combined arrays, followed by the Nereocystis only and Saccharina only arrays. These results may have significance for large scale macroalgal production.

09 BIOMASS FUELS↗
DOE OSTI2020

Beam characterization at NSRL for radiobiological experiments—phase 1

An experimental campaign was carried out at the NASA Space Radiation Laboratory to perform an additional, independent dosimetric characterization of the beams of protons, helium and carbon ions for radiobiological experiments. The campaign was undertaken by the request and with the support from the National Cancer Institute, U.S. In this initial phase, the goals were to obtain a first assessment of the dosimetric reproducibility of the beam control system, including analysis of spatial homogeneity and evaluation of ion beam contamination. They should facilitate the design of further experimental campaigns for beam characterization for radiobiological experiments. Measurements included reference dosimetry with comparison of in-house and external ionization chambers and electrometers, lateral-dose profile measurements in air, depth-dose profile in a water tank, evaluation of water equivalent thickness of a HDPE binary range shifter and estimation of impurities of the investigated helium-ion beam. The experiments and results are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Exact Landau Level Description of Geometry and Interaction in a Flatband

Flatbands appear in many condensed matter systems, including the two-dimensional electron gas in a high magnetic field, correlated materials, and moiré heterostructures. They are characterized by intrinsic geometric properties such as the Berry curvature and Fubini-Study metric. The influence of the band geometry on electron-electron interaction is difficult to understand analytically because the geometry is in general nonuniform in momentum space. In this work, we study the topological flatband of Chern number C = 1 with a momentum-dependent but positive definite Berry curvature that fluctuates in sync with Fubini-Study metric. We derive an exact correspondence between such ideal flatbands and Landau levels and show that the band geometry fluctuation gives rise to a new type of interaction in the corresponding Landau levels that depends on the center of mass of two particles. We characterize such interactions by generalizing the usual Haldane pseudopotentials. This mapping gives exact zero-energy ground states for short-ranged repulsive generalized pseudopotentials in flatbands, in analogy to fractional quantum Hall systems. Driving the center-of-mass interactions beyond the repulsive regime leads to a dramatic reconstruction of the ground states towards gapless phases. The generalized pseudopotential could be a useful basis for future numerical studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Towards Precision Photometric Type ia Supernova Cosmology With Machine Learning

The revolutionary discovery of dark energy and accelerating cosmic expansion was made with just 42 type Ia supernovae (SNe Ia) in 1999. Since then, large synoptic surveys, e.g., Dark Energy Survey (DES), have observed thousands more SNe Ia and the upcoming Rubin Legacy Survey of Space and Time (LSST) and Roman Space Telescope promise to deliver millions in the next decade. This unprecedented data volume can produce the required precision to unambiguously test concordance cosmology which could represent a monumental shift in our understanding of dark energy and its role in cosmic history. However, extracting a pure SN Ia sample with accurate redshifts for such a large dataset will be a challenge. Specifically, spectroscopic classification will not be possible for the vast majority of discovered objects, and only ~25% will have spectroscopic redshifts. This thesis presents a series of observational and methodological studies designed to address the questions associated with this new era of photometric SN Ia cosmology. First, we present a machine learning method for photometric classification of SNe, Supernova Classification with a COnvolutional Neural Network (SCONE). Photometric classification enables SNe with no spectroscopic information to be confidently categorized, and is a critical component of current and future analysis pipelines. SCONE achieves >99% accuracy distinguishing simulated SNe Ia from non-Ia SNe, and has been integrated into DES, LSST, and Roman analysis pipelines. We also demonstrate the efficacy of SCONE on early-time photometric classification, which will be vital for optimal allocation of spectroscopic resources. We show that SCONE can distinguish between 6 SN types with 75% accuracy on the night of initial discovery, comparable to results in the literature for full-phase SNe. Next, we study current methods for estimating SN Ia redshifts and propose a machine learning alternative that uses SN photometry alone to extract redshift information. Most SNe Ia inherit redshift information from their host galaxy, but the process of matching SNe to the correct host galaxy can be challenging. We systematically analyze the impact of incorrect redshifts from host galaxy mismatch on 5 years of DES SN data, and conclude that improved host matching or redshift estimation methods can reduce our systematic errors by ~10%. In response to this finding, we present a SN photometry-only method for estimating redshifts independent of host galaxy information, Photo-zSNthesis. We show that Photo-zSNthesis redshift estimates are accurate to within 2% across the full redshift range of LSST, a first in the literature. Finally, we focus on the robustness of machine learning (ML) algorithms for real-world and scientific applications. ML models generalize poorly beyond their training set and often experience severe performance degradation when deployed on new data. We demonstrate a general method for improving robustness that achieves new state-of-the-art results on astronomical object classification, wildlife identification, and tumor detection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Filamentation and focusing of electron beams due to interactions with plasma waves

Results of numerical modeling of the interaction of an electron beam propagating across relativistic plasma waves indicate that electron beam filamentation and focusing may occur under certain conditions. The model is based on solving the relativistic equation of motion in three dimensions for the individual electrons in a tenuous Gaussian beam, as they pass through a relativistic plasma wave. Several electron beam and plasma wave parameters were varied, and the results are summarized. One of the results is that the spacing of the electron beam filaments correlates with the wavelength of the plasma wave. The electron beam filaments appear as vertical slabs after the beam exits the plasma waves. The electron beam also compresses to a focus after it exits the plasma, and the focal distance depends on several parameters including the electron beam energy and phase velocity of the relativistic plasma wave. It is suggested that these focusing and filamentation phenomena may be the basis for diagnostics schemes for laser plasma interactions. The parameters used in the model were electron beam energies in the 5–50 keV range and plasma wave properties typical for the beat-wave produced by CO2 lasers, which correspond to the facilities available in our laboratory. The limitations of these results to lower energy density beam and plasma regimes and to higher energy density regimes will be discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Global quantification of the dispersion effect with POLDER satellite data

Increased aerosols can modify the shape of the cloud Particle Size Distribution (PSD), thereby influencing the radiative properties of clouds, known as the Dispersion Effect (DE). However, a global, observation-based quantification of its impact on Aerosol-Cloud Interactions (ACI) is lacking, leading to DE being typically ignored in satellite-based estimates of ACI forcing. Here we propose a physics-based method that combines polarimetric satellite data on cloud PSD to achieve global observational quantification of DE’s impact on ACI in liquid-phase stratiform clouds. Globally, DE offsets ACI changes induced by droplet number concentration variation and liquid water path adjustment by 7% and −1.4%, respectively. Furthermore, a parameterization based on the global dataset of PSD shape parameters is developed to improve DE estimation in large-scale models. Both the quantification and parameterization enhance our understanding of DE and facilitate the inclusion of this non-negligible impact of DE on ACI in estimating aerosol climate forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Importance of Engineered and Learned Molecular Representations in Predicting Organic Reactivity, Selectivity, and Chemical Properties

Machine-readable chemical structure representations are foundational in all attempts to harness machine learning for the prediction of reactivities, selectivities, and chemical properties directly from molecular structure. The featurization of discrete chemical structures into a continuous vector space is a critical phase undertaken before model selection, and the development of new ways to quantitatively encode molecules is an active area of research. Here, we highlight the application and suitability of different representations, from expert-guided “engineered” descriptors to automatically “learned” features, in different prediction tasks relevant to organic and organometallic chemistry, where differing amounts of training data are available. These tasks include statistical models of stereo- and enantioselectivity, thermochemistry, and kinetics developed using experimental and quantum chemical data. The use of expert-guided molecular descriptors provides an opportunity to incorporate chemical knowledge, domain expertise, and physical constraints into statistical modeling. In applications to stereoselective organic and organometallic catalysis, where data sets may be relatively small and 3D-geometries and conformations play an important role, mechanistically informed features can be used successfully to obtain predictive statistical models that are also chemically interpretable. We provide an overview of several recent applications of this approach to obtain quantitative models for reactivity and selectivity, where topological descriptors, quantum mechanical calculations of electronic and steric properties, along with conformational ensembles, all feature as essential ingredients of the molecular representations used. Alternatively, more flexible, general-purpose molecular representations such as attributed molecular graphs can be used with machine learning approaches to learn the complex relationship between a structure and prediction target. This approach has the potential to out-perform more traditional representation methods such as “hand-crafted” molecular descriptors, particularly as data set sizes grow. One area where this is particularly relevant is in the use of large sets of quantum mechanical data to train quantitative structure–property relationships. A general approach toward curating useful data sets and training highly accurate graph neural network models is discussed in the context of organic bond dissociation enthalpies, where this strategy outperforms regression using precomputed descriptors. Finally, we describe how graph neural network predictions can be incorporated into mechanistically informed statistical models of chemical reactivity and selectivity. Once trained, this approach avoids the expensive computational overhead associated with quantum mechanical calculations, while maintaining chemical interpretability. We illustrate examples for which fast predictions of bond dissociation enthalpy and of the identities of radicals formed through cleavage of a molecule’s weakest bond are used in simple physical models of site-selectivity and reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗
DOE OSTI2022

Superfluid Weight Bounds from Symmetry and Quantum Geometry in Flat Bands

Flat-band superconductivity has theoretically demonstrated the importance of band topology to correlated phases. In two dimensions, the superfluid weight, which determines the critical temperature through the Berezinksii-Kosterlitz-Thouless criteria, is bounded by the Fubini-Study metric at zero temperature. We show this bound is nonzero within flat bands whose Wannier centers are obstructed from the atoms—even when they have identically zero Berry curvature. Next, we derive general lower bounds for the superfluid weight in terms of momentum space irreps in all 2D space groups, extending the reach of topological quantum chemistry to superconducting states. We find that the bounds can be naturally expressed using the formalism of real space invariants (RSIs) that highlight the separation between electronic and atomic degrees of freedom. Finally, using exact Monte Carlo simulations on a model with perfectly flat bands and strictly local obstructed Wannier functions, we find that an attractive Hubbard interaction results in superconductivity as predicted by the RSI bound beyond mean field. Hence, obstructed bands are distinguished from trivial bands in the presence of interactions by the nonzero lower bound imposed on their superfluid weight.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Impact of particle arrays on phase separation composition patterns

We examine the symmetry-breaking effect of fixed constellations of particles on the surface-directed spinodal decomposition of binary blends in the presence of particles whose surfaces have a preferential affinity for one of the components. Our phase-field simulations indicate that the phase separation morphology in the presence of particle arrays can be tuned to have a continuous, droplet, lamellar, or hybrid morphology depending on the interparticle spacing, blend composition, and time. In particular, when the interparticle spacing is large compared to the spinodal wavelength, a transient target pattern composed of alternate rings of preferred and non-preferred phases emerges at early times, tending to adopt the symmetry of the particle configuration. We reveal that such target patterns stabilize for certain characteristic length, time, and composition scales characteristic of the pure phase-separating mixture. To illustrate the general range of phenomena exhibited by mixture-particle systems, we simulate the effects of single-particle, multi-particle, and cluster–particle systems having multiple geometrical configurations of the particle characteristic of pattern substrates on phase separation. Our simulations show that tailoring the particle configuration, or substrate pattern configuration, a relative fluid-particle composition should allow the desirable control of the phase separation morphology as in block copolymer materials, but where the scales accessible to this approach of organizing phase-separated fluids usually are significantly larger. Limited experiments confirm the trends observed in our simulations, which should provide some guidance in engineering patterned blend and other mixtures of technological interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Measurement of $d^2 \sigma/d|\vec{q}|d E_{avail}$ and 2p2h contribution using charged current $\nu_\mu$ interactions in the NOvA Near Detector

This Thesis presents the analysis methods and measurement of the $d^2 \sigma/d|\vec{q}|d E_{avail}$-CC inclusive cross section that describes neutrino scattering in a predominantly hydrocarbon medium. Studies that are prerequisite for the analysis including, variable determinations, event selections, efficiency, purity determinations, and systematic uncertainty estimations, are summarized. Tests of the analysis are also detailed, and shown to reproduce the input distributions. The double-differential cross-section measurement provides the foundation for determination of the inclusive rate for the excess induced by 2p2h processes together with nuclear medium effects that are not described by the Fermi gas model. The methodology for the extraction of this rate, defined as the excess observed relative to the rate estimated for known single-nucleon interactions, is described in detail. The following new measurements are reported in this Thesis: (1) The $\nu_\mu-$CC inclusive double-differential cross section as a function of three-momentum transfer and available hadronic energy, for an average $E_\nu$ of 1.8 GeV, is obtained. The double-differential cross section value that obtains when the bin widths are taken to be the dimensions of the entire analysis domain, e.g. $0.2 \le |\vec{q}| \le 2.0 GeV/c$ and $0.0 \le E_{avail} \le 2.0$ GeV, is inclusive cross section is $(6.05 \pm 0.75)\times10^{-39}$ cm$^2$/GeV/GeV/$c$/nucleon. An event excess attributed to 2p2h-MEC processes is observed; most of the rate occurs in a contiguous phase space region $0.3 \le |\vec{q}| \le 1.0$ GeV$/c$ $0.0 \le E_{avail} \le 0.35$ GeV. This rate represents $12.0\pm 6.5\%$ of the observed CC inclusive cross section; its cross section ratio relative to CCQE scattering is estimated to be $44.3 \pm 23.9\%$.

Olson, Travis Grant↗
DOE OSTI2021

Study of structural, transport and magneto-crystalline anisotropy in La 1-x Sr x MnO 3 (0.30 ≤ x ≤ 0.40) perovskite manganites

In this paper we present structural, transport, magnetic and magneto-crystalline anisotropy of La 1–x Sr x MnO 3 (х = 0.30, 0.33, 0.36, and 0.40) synthesized using the solid-state reaction method. The X-ray diffraction pattern of the samples is well indexed to the rhombohedral structure with R 3 - c space group. Lattice parameters and unit cell volume are found to decrease monotonically upon increasing Sr 2+ concentration at La 3+ site. X-ray photoelectron spectroscopy confirms the presence of Mn 4+ and Mn 3+ , the ratio of which (Mn 4+ /Mn 3+ ) increases from 0.44 to 0.64 with increasing Sr 2+ concentration. From the ρ(T) and M(T) data analysis, we find that T MI and M S are larger for x = 0.33 composition, suggesting a weakening of the double-exchange interaction due to competing super-exchange interactions in x = 0.30, 0.36 and 0.40 samples. We calculated MCA constant (K 1 ) by fitting the saturation magnetisation data at a higher magnetic field using the law of approach model, M = M S [1-a/H-b/H 2 ]-χH with K 1 2 = (105/8)bMs 2 and is in the order of 10 5 erg/cm 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗
DOE OSTI2021

Structural, magneto transport and magnetic properties of Ruddlesden–Popper La 2-2x Sr 1+2x Mn 2 O 7 (0.42≤x≤0.52) layered manganites

The magneto transport of Ruddlesden–Popper, La 2-2x Sr 1+2x Mn 2 O 7 (0.42≤x≤0.52), 2-dimensional bilayered manganites have been investigated in a broad temperature range. The samples have been synthesized using the solid-state reaction method. Rietveld refinement of the X-ray diffraction data indicates the tetragonal structure formation with I4/mmm space group. The resistivity curves of the samples present a general characteristic of metal-insulator (MI) transition at certain temperature (T MI ). Besides, for samples with x = 0.48, 0.50, and 0.52, at a certain temperature (T CO ) well below T MI the charge ordering is also evident. Furthermore, the samples display a shallow upturn in resistivity below T min due to the Kondo like spin scattering effect, weak localization and electron-phonon interaction. The high temperature semiconducting (T>T MI ) region, the resistivity curve follows the 3D Mott’s variable hopping transport mechanism. The overall suppression of resistivity accounts for a substantial magnetoresistance by applying a magnetic field and the characteristic change in T MI , T CO , and T min are discussed. Temperature-dependent magnetization demonstrates the suppression of ferromagnetism and evident of antiferromagnetic nature with an apparent charge ordering with increasing concentration of Sr 2+ .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
DOE OSTI2022

Insights into Cation Ordering of Double Perovskite Oxides from Machine Learning and Causal Relations

This work investigates origins of cation ordering in double perovskites using first-principles theory computations combined with machine learning (ML) and causal relations. We have considered various oxidation states of A, A', B, and B' from the family of transition metal ions to construct a diverse compositional space. A conventional framework employing traditional ML classification algorithms such as Random Forest (RF) coupled with appropriate features including geometry-driven and key structural modes leads to accurate prediction (~98%) of A-site cation ordering. We have evaluated the accuracy of ML models by employing analyses of decision paths, assignments of probabilistic confidence bound, and finally a direct non-Gaussian acyclic structural equation model to investigate causality. Our study suggests that structural modes are crucial for classifying layered, columnar, and rock-salt ordering. The charge difference between A and A' is the most important feature for predicting clear layered ordering, which in turn depends on the B and B' charge separation. We have also designed mathematical relationships with these features to derive energy differences to form clear layered ordering. Here, the trilinear coupling between tilt, in-phase rotation, and A-site antiferroelectric displacement in the Landau free-energy expansion becomes the necessary condition behind formation of A-site cation ordering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Determining the jet transport coefficient $\hat{q}$ from inclusive hadron suppression measurements using Bayesian parameter estimation

We report a new determination of $\hat{q}$, the jet transport coefficient of the quark-gluon plasma. We use the JETSCAPE framework, which incorporates a novel multistage theoretical approach to in-medium jet evolution and Bayesian inference for parameter extraction. The calculations, based on the Matter and Lbt jet quenching models, are compared to experimental measurements of inclusive hadron suppression in Au + Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and Pb + Pb collisions at the CERN Large Hadron Collider (LHC). The correlation of experimental systematic uncertainties is accounted for in the parameter extraction. The functional dependence of $\hat{q}$ on jet energy or virtuality and medium temperature is based on a perturbative picture of in-medium scattering, with components reflecting the different regimes of applicability of Matter and Lbt. In the multistage approach, the switch between Matter and Lbt is governed by a virtuality scale Q 0 . Comparison of the posterior model predictions to the RHIC and LHC hadron suppression data shows reasonable agreement, with moderate tension in limited regions of phase space. The distribution of $\hat{q}/T^3$ extracted from the posterior distributions exhibits weak dependence on jet momentum and medium temperature T, with 90% credible region (CR) depending on the specific choice of model configuration. The choice of Matter + Lbt, with switching at virtuality Q 0 , has 90% CR of $2<\hat{q}/T^3<4$ for p T jet > 40 GeV=c. The value of Q 0 , determined here for the first time, is in the range 2.0-2.7 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Robust clustering of the local Milky Way stellar kinematic substructures with Gaia eDR3

Understanding local stellar kinematic substructures in the solar neighbourhood helps build a complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. We apply the clustering algorithm HDBSCAN on the Gaia early third data release to identify a list of stable clusters in velocity space and action-angle space by taking into account the measurement uncertainties and studying the stability of the clustering results. We find 1405 (497) stars in 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disc (Sequoia) are identified in velocity space (action-angle space). Although we do not identify any new structures, we find that the HDBSCAN member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries.

79 ASTRONOMY AND ASTROPHYSICS↗