Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data fitting”

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,225 records · Page 68

From Points to Planes: A Workflow for Converting Three‐Dimensional Point Cloud Data Into Discrete Fracture Network Flow and Transport Models

We present the Point cLoud Algorithm for NEtwork Extraction of Discrete Fracture Networks (PLANE-DFN), a point cloud–based algorithm for automatic fracture network extraction designed to support discrete fracture network (DFN) modeling workflows. PLANE-DFN segments three-dimensional fracture planes from raw point cloud data using RANdom SAmple Consensus coupled with statistical outlier removal and density-based clustering to isolate individual fracture features. Each candidate plane is constrained against site-specific structural constraints based on strike and dip. After segmentation, each fracture is converted into a 2-D convex polygon suitable for meshing and simulation. The PLANE-DFN algorithm is validated by comparing geometric and flow and transport data against data from dfnWorks simulations with ensembles of plane-fit networks. We find that the flow and transport in plane-fit networks are comparable to dfnWorks-generated networks when realistic network geometry is maintained. The PLANE-DFN algorithm provides an automated and streamlined workflow to transform point clouds of data into DFN network geometry.

54 ENVIRONMENTAL SCIENCES↗

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗

Influence of living grass Roots and endophytic fungal hyphae on soil hydraulic properties

Soil hydraulic properties are often estimated based on laboratory data or pedotransfer functions dependent on soil physical properties, which often do not consider potential impacts of soil roots or fungal hyphae. Here, we first review current knowledge of how these soil biotic components affect hydraulic properties, then we conducted laboratory experiments to specifically test if the presence of roots and mycorrhizal fungi had a significant effect on the hydraulic properties of two soils with contrasting textures: Flint sand and Hamblen silt loam. Soil cores were seeded with (Panicum virgatum) and grown in a greenhouse over three separate growth periods. The endophytic fungus Serendipita indica was injected as liquid inoculant into designated mycorrhizal cores. Saturated hydraulic conductivity (K sat ) measurements were made with a constant head permeameter, and soil water retention curves were obtained by the evaporation method, supplemented at the dry end for Hamblen silt loam with water activity meter data. Retention curve parameters were obtained by fitting the van Genuchten equation to the resulting measurements. Mean root volume ratios were higher in the mycorrhizal inoculated treatment than in the uninoculated treatment for both soils. For Flint sand, analysis of variance revealed that K sat was reduced by the presence of roots as compared to bare soil. This was likely due to roots clogging soil pores. Results also indicated the presence of roots changed the shape of the water retention curve for Flint sand by increasing water content at saturation and by reducing the slope of the curve. These changes suggested roots created additional porosity and broadened the pore-size distribution. The presence of mycorrhizal fungi accentuated the root effects. The influence of roots and mycorrhizal fungi on hydraulic properties was less obvious for the Hamblen silt loam, as none of the treatments differed from each other at p < 0.05. The results highlight the necessity to consider the impact of root and fungal structures on models of soil hydraulic properties.

59 BASIC BIOLOGICAL SCIENCES↗

Temperature Dependence of Nuclear Quadrupole Resonance and the Observation of Metal–Ligand Covalency in Actinide Complexes: 35 Cl in Cs 2 UO 2 Cl 4

We report a study of the temperature dependence of 35 Cl nuclear quadrupole resonance (NQR) transition energies and spin–lattice relaxation times (T 1 ) for 235 U-depleted dicesium uranyl tetrachloride (Cs 2 UO 2 Cl 4 ) aimed at elucidating electronic interactions between the uranium center and atoms in the equatorial plane of the UO 2 2+ ion. The transition frequency decreases slowly with temperature below 75 K and with a more rapid linear dependence above this temperature. The spin–lattice relaxation time becomes shorter with temperature, and as temperatures increase, the T 1 decrease becomes nearly quadratic. The observed trends are reproduced by a model that assumes phonon-induced fluctuations of the electric field gradient tensor and partial electron delocalization from Cl to U. The fit of the theoretical model to experimental data allows a Debye temperature of 96 K to be estimated. Finally, the generalization of this approach to investigations of covalency in actinide–ligand bonding is examined.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model-Free Approach for Profiling of Polydisperse Soft Matter Using Small Angle Scattering

A strategy for determining the size polydispersity of systems from their small angle coherent scattering is outlined. Here, using the method of moment expansion, we show that the various central moments representing the average particle size, variance of particle size, and skewness of size distribution function (SDF) for polydisperse systems can be extracted from spectral analysis without bias. When the degree of polydispersity is moderate, SDF can be further reconstructed based on the maximum entropy principle. Numerical benchmarking of a model study over a wide range of size nonuniformity demonstrates the validity of this analytical approach for quantifying the size distribution of general soft matter systems in a model-free manner. Furthermore, the efficacy of this method was validated by successfully applying it to the fitting of small-angle neutron scattering data obtained from L64 Pluronic micelles using various form factor models. The numerical and experimental verification underscores the reliability and versatility of this method in accurately characterizing the size distribution of complex soft matter systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exponential vs Gaussian Correlation Functions in the Characterization of Block Copolymer Grain Structure by Depolarized Light Scattering

Block copolymer (BCP) grain structure affects the mechanical, optical, and electrical properties of BCP materials, making the accurate characterization of this grain structure an important goal. In this study, improved BCP grain parameters were obtained by employing an exponentially decaying correlation function within the ellipsoidal grain model, instead of the Gaussian correlation function that was used in previous work. The exponential correlation function provides a better fit to the experimental depolarized light scattering data, which outweighs the disadvantage that it requires numerical integration to obtain the model scattered intensity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TROPOMI NO 2 in the United States: A Detailed Look at the Annual Averages, Weekly Cycles, Effects of Temperature, and Correlation With Surface NO 2 Concentrations

Abstract Observing the spatial heterogeneities of NO 2 air pollution is an important first step in quantifying NO X emissions and exposures. This study investigates the capabilities of the Tropospheric Monitoring Instrument (TROPOMI) in observing the spatial and temporal patterns of NO 2 pollution in the continental United States. The unprecedented sensitivity of the sensor can differentiate the fine‐scale spatial heterogeneities in urban areas, such as emissions related to airport/shipping operations and high traffic, and the relatively small emission sources in rural areas, such as power plants and mining operations. We then examine NO 2 columns by day‐of‐the‐week and find that Saturday and Sunday concentrations are 16% and 24% lower respectively, than during weekdays. We also analyze the correlation of daily maximum 2‐m temperatures and NO 2 column amounts and find that NO 2 is larger on the hottest days (>32°C) as compared to warm days (26°C–32°C), which is in contrast to a general decrease in NO 2 with increasing temperature at moderate temperatures. Finally, we demonstrate that a linear regression fit of 2019 annual TROPOMI NO 2 data to annual surface‐level concentrations yields relatively strong correlation ( R 2 = 0.66). These new developments make TROPOMI NO 2 satellite data advantageous for policymakers and public health officials, who request information at high spatial resolution and short timescales, in order to assess, devise, and evaluate regulations.

54 ENVIRONMENTAL SCIENCES↗

Assessment of Extinction‐, Satellite‐, and Model‐Based Vertical Cloud Condensation Nuclei (CCN) Retrieval Methods Using Airborne CCN Measurements Over the Southern Great Plains

Abstract Accurate estimates of the vertical profile of cloud condensation nuclei (CCN) concentration are crucial to better quantify aerosol‐cloud interactions. We assessed the correlation between the vertical CCN concentrations obtained from extinction‐, satellite‐, and model‐based retrieval methods and airborne CCN concentrations collected at 0.24% supersaturation within the 3, 9, 27, and 81 km regions centered over the U.S. Department of Energy's Atmospheric Radiation Measurement User Facility Southern Great Plains (SGP) site during the spring and summer of 2016. The extinction profiles at a wavelength 355 nm were provided by the ground‐based Raman lidar. Our analysis showed moderate correlation between dry‐corrected extinction and airborne CCN data. We found the retrieved number concentration of CCN (RNCCN) method showed regression best‐fit slopes close to unity and consistent prediction errors for the majority of the data. The Lenhardt et al. (2023, https://doi.org/10.5194/amt‐16‐2037‐2023 ) method showed similar conclusions but only during spring, whereas the Mamouri and Ansmann (2016, https://doi.org/10.5194/acp‐16‐5905‐2016 ) method showed poor correlation. The Shinozuka et al. (2015, https://doi.org/10.5194/acp‐15‐7585‐2015 ) satellite‐based method exhibited reasonable agreement during summer but poor correlation during periods where both high (∼1,400 #/cm 3 ) and low (∼50 #/cm 3 ) airborne CCN concentrations were observed. The Copernicus Atmosphere Monitoring Service reanalysis modeled 3‐D CCN data set showed a moderate to weak positive correlation but performed poorly at high airborne CCN concentrations. Our analysis suggests the extinction‐based RNCCN method performed better than other methods across most observation periods under the diverse meteorological conditions observed at the SGP site.

54 ENVIRONMENTAL SCIENCES↗

Geometrically frustrated interactions drive structural complexity in amorphous calcium carbonate

Amorphous calcium carbonate is an important precursor for biomineralization in marine organisms. Key outstanding problems include understanding the structure of amorphous calcium carbonate and rationalizing its metastability as an amorphous phase. Here we report high-quality atomistic models of amorphous calcium carbonate generated using state-of-the-art interatomic potentials to help guide fits to X-ray total scattering data. Exploiting a recently developed inversion approach, we extract from these models the effective Ca∙∙∙Ca interaction potential governing the structure. This potential contains minima at two competing distances, corresponding to the two different ways that carbonate ions bridge Ca 2+ -ion pairs. We reveal an unexpected mapping to the Lennard-Jones–Gauss model normally studied in the context of computational soft matter. The empirical model parameters for amorphous calcium carbonate take values known to promote structural complexity. We thus show that both the complex structure and its resilience to crystallization are actually encoded in the geometrically frustrated effective interactions between Ca 2+ ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Delocalization error poisons the density-functional many-body expansion

The many-body expansion is a fragment-based approach to large-scale quantum chemistry that partitions a single monolithic calculation into manageable subsystems. This technique is increasingly being used as a basis for fitting classical force fields to electronic structure data, especially for water and aqueous ions, and for machine learning. Here, we show that the many-body expansion based on semilocal density functional theory affords wild oscillations and runaway error accumulation for ion–water interactions, typified by F − (H 2 O) N with N ≳ 15. We attribute these oscillations to self-interaction error in the density-functional approximation. The effect is minor or negligible in small water clusters, explaining why it has not been noticed previously, but grows to catastrophic proportion in clusters that are only moderately larger. This behavior can be counteracted with hybrid functionals but only if the fraction of exact exchange is ≳50%, whereas modern meta-generalized gradient approximations including ωB97X-V, SCAN, and SCAN0 are insufficient to eliminate divergent behavior. Other mitigation strategies including counterpoise correction, density correction (i.e., exchange–correlation functionals evaluated atop Hartree–Fock densities), and dielectric continuum boundary conditions do little to curtail the problematic oscillations. In contrast, energy-based screening to cull unimportant subsystems can successfully forestall divergent behavior. These results suggest that extreme caution is warranted when the many-body expansion is combined with density functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatially resolved x-ray detection with photonic crystal scintillators

We study the self-collimation phenomenon in photonic crystals (PhC) of wide bandgap materials for ultra-fast and high spatial resolution x-ray detection. We work on various heavy inorganic scintillators: BaF 2 , GaN, ZnO, CsI:Tl, NaI:Tl, LYSO, WO 4 compounds, and plastic scintillators. Conventional scintillator detectors do not rely on a direct detection mechanism; hence, they require intricate design and fabrication processes. We offer a PhC design to observe self-collimation phenomena and overcome the ongoing spatial resolution challenges with these types of materials. We investigate the photonic band diagrams and iso-frequency contours. Fourier transforms based on finite-difference time-domain and frequency domain simulations are done for verifying and analyzing the self-collimation with the selected material. Light extraction efficiency at the PhC–air interface, depending on the truncation distance from the excitation point, is measured. Beam divergence values are calculated at 1 mm propagation distance. The vertical field profiles are obtained to observe the confinement. For the spatial resolution analysis, cross-sectional beam profiles have been examined. Furthermore, Gaussian envelopes are fitted to beam profiles for a consistent data analysis, and full-width-at-half-maximum values are considered. As a result, we theoretically prove and demonstrate the spatially resolved x-ray detection at the sub-micrometer level for a wide range of scintillator materials.

36 MATERIALS SCIENCE↗

Mixture Model for Refrigerant Pairs R-32/1234yf, R-32/1234ze(E), R-1234ze(E)/227ea, R-1234yf/152a, and R-125/1234yf

In this work, thermodynamic models based on the corresponding states framework with departure terms are developed for the refrigerant pairs R-32/1234yf, R-32/1234ze(E), R-1234ze(E)/227ea, R-1234yf/152a, and R-125/1234yf. These models are based on new measurements of density, speed of sound, and phase equilibria, combined with the data available in the literature. The model for R-32/1234yf is most comprehensive in its data coverage, with speed of sound deviations within 1%, density deviations within 0.1%, and bubble- and dew-point pressure deviations within 1%. In conclusion, the other mixtures have generally more limited data availability but a similar goodness of fit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic properties of “near-perfect” gas

When applied thermodynamics requires simple representations of thermodynamic quantities, one often finds quantities such as pressure and the thermal energy density expressed as a monomial linear product of density and temperature power laws. This is a simple generalization of the product of density and temperature used for the perfect gas that admits a broader range of thermodynamic behavior into analytic fluid calculations or provides an analytic form that is readily fit to tabulated equation-of-state data for an arbitrary material. This paper reviews the thermodynamic properties of this generalized perfect-gas model, treating it as a class of “near-perfect-gas” models that are defined and unified here, for the first time, through a shared Helmholtz free energy. This generalization from perfect to near-perfect-gas models preserves some perfect-gas properties (e.g., constant adiabatic exponents) but not others (e.g., constant specific heats) and is constrained by thermodynamic consistency. It is important to be aware of the properties of this class of near-perfect-gas models when using them as surrogates for real nonperfect materials that may have properties that cannot be captured by near-perfect-gas models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the accuracy of compressibility transformations

This study highlights the importance of satisfying the eddy viscosity equivalence below the logarithmic layer, to deriving accurate compressibility transformations. First, we analyze the ability of known transformations to satisfy the eddy viscosity equivalence and show that the accuracy of these transformations is strongly dependent on this ability. Second, in a step-by-step manner, we devise new transformations that satisfy this hypothesis. An approach based on curve fitting of the incompressible Direct Numerical Simulation data for eddy viscosity profiles below the logarithmic layer provides an extremely accurate transformation, which motivates self-contained methods, making use of mixing length formulas in the inner region. It is shown that the accuracy of existing transformations can be significantly improved by applying these ideas, below the logarithmic layer. Motivated by the effectiveness of the formulations derived from eddy viscosity equivalence, we introduce a new integral transformation based on Reynolds number equivalence between compressible and incompressible flows. This approach is based on defining a new compressible velocity scale, which affects the accuracy of transformations. Several choices for the velocity scale are tested, and in each attempt, it is shown that the eddy viscosity equivalence plays a very important role for the accuracy of compressibility transformations.

42 ENGINEERING↗

Analysis of a Heat Balance Experiment in the TREAT Reactor with SPARTA and OpenMC

This article presents the interpretation of a heat balance experiment in TREAT using the SPARTA methodology and an OpenMC model. The study examines the effects of control rods and temperature variations on the effective multiplication factor ($k_{\rm eff}$) and correction factors. Results indicate that both core temperature and rod positions have a linear effect on $k_{\rm eff}$ and correction factors in the upper-west channel. A critical hypothesis of SPARTA was validated: reactivities and correction factors can be obtained by linearly combining individual effects. The previously considered negligible effect of temperature on correction factors was found to be significant, necessitating updates to the SPARTA methodology. By fitting an exponential curve to mean temperature data, SPARTA's accuracy and convergence have been improved. Comparisons between SPARTA and OpenMC showed consistent reactivity removal due to temperature increase. While SPARTA effectively interprets complex transients, simpler methods may suffice when correction factors have minimal influence on reactivity calculations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Application of machine learning and artificial intelligence to extend EFIT equilibrium reconstruction

Recent progress in the application of machine learning (ML)/artificial intelligence (AI) algorithms to improve the Equilibrium Fitting (EFIT) code equilibrium reconstruction for fusion data analysis applications is presented. A device-independent portable core equilibrium solver capable of computing or reconstructing equilibrium for different tokamaks has been created to facilitate adaptation of ML/AI algorithms. A large EFIT database comprising of DIII-D magnetic, motional Stark effect, and kinetic reconstruction data has been generated for developments of EFIT model-order-reduction (MOR) surrogate models to reconstruct approximate equilibrium solutions. Furthermore, a neural-network MOR surrogate model has been successfully trained and tested using the magnetically reconstructed datasets with encouraging results. Other progress includes developments of a Gaussian process Bayesian framework that can adapt its many hyperparameters to improve processing of experimental input data and a 3D perturbed equilibrium database from toroidal full magnetohydrodynamic linear response modeling using the Magnetohydrodynamic Resistive Spectrum - Feedback (MARS-F) code for developments of 3D-MOR surrogate models.

Gaussian process↗

HOD-dependent systematics in Emission Line Galaxies for the DESI 2024 BAO analysis

The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we estimate the level of systematics induced in the DESI BAO analysis due the assumed Halo Occupation Distribution (HOD) model for the Emission Line Galaxy (ELG) tracer. We make use of mock galaxy catalogs constructed by fitting various HOD models to early DESI data, namely the One-Percent survey data. Our analysis includes typical HOD models for the ELG tracer used in the literature as well as extensions to the baseline models. Among the extensions, we consider various recipes for galactic conformity and assembly bias. We use 25 AbacusSummit simulations under the ΛCDM cosmology for each HOD model and perform independent analyses in Fourier space and in configuration space. To recover the BAO signal from our mocks we perform BAO reconstruction and apply the control variates technique to reduce sample variance noise. Our BAO analyses can recover the isotropic BAO parameter α iso within 0.1% and the Alcock Paczynski parameter α AP within 0.3%. Overall, we find that the systematic error due to the HOD dependence is below 0.17%, with the Fourier space analysis being more robust against the HOD systematics. We conclude that our analysis pipeline is robust enough against the HOD systematics for the ELG tracer in the DESI 2024 BAO analysis, for the assumptions made.

79 ASTRONOMY AND ASTROPHYSICS↗