Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Equation of State and Thermometry of the 2D SU ( N ) Fermi-Hubbard Model

We characterize the equation of state (EoS) of the SU ( N > 2 ) Fermi-Hubbard Model (FHM) in a two-dimensional single-layer square optical lattice. We probe the density and the site occupation probabilities as functions of interaction strength and temperature for N = 3 , 4, and 6. Our measurements are used as a benchmark for state-of-the-art numerical methods including determinantal quantum Monte Carlo and numerical linked cluster expansion. By probing the density fluctuations, we compare temperatures determined in a model-independent way by fitting measurements to numerically calculated EoS results, making this a particularly interesting new step in the exploration and characterization of the SU ( N ) FHM. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Hybrid Modeling of Three-Phase Grid-Supporting Inverters for Dynamic Studies

Grid technologies connected by power electronic converter (PEC) interfaces continually implement grid support functions mandated by grid codes and standards. The transition to converter-based generation demands precise PEC models to assess system dynamics, which have been previously overlooked in conventional power systems. This study proposes a hybrid method for analyzing grid-connected three-phase PEC dynamics with the IEEE standard 1547-2018 Volt-VAr mode that combines physics and data-driven techniques. The physics model reflects the PEC’s internal behavior, whereas the data-driven modeling technique evaluates the grid-supporting capabilities of the smart PEC. The system identification approach is used to generate dynamic PEC models based on changing grid voltage and measured current injected into the grid by the PEC. In the Volt-VAr support mode, a detailed topological model including switches is utilized to compare the goodness-of-fit of the extracted hybrid dynamic model. The results demonstrate that the hybrid PEC model in the Volt-VAr mode accurately matches the dynamics with the topological model.

Subedi, Sunil↗
DOE OSTI2023

Spatiotemporal measurements of striations in a glow discharge’s positive column using laser-collisional induced fluorescence

Here we have observed the behavior of striations caused by ionization waves propagating in low-pressure helium DC discharges using the non-invasive laser-collision induced fluorescence (LCIF) diagnostic. To achieve this, we developed an analytic fit of collisional radiative model (CRM) predictions to interpret the LCIF data and recover quantitative two-dimensional spatial maps of the electron density, n e , and the ratios of LCIF emission states that can be correlated with T e with the use of accurate distribution functions at localized positions within striated helium discharges at 500 mTorr, 750 mTorr, and 1 Torr. To our knowledge, these are the first spatiotemporal, laser-based, experimental measurements of n e in DC striations. The n e and 447:588 ratio distributions align closely with striation theory. Constriction of the positive column appears to occur with decreased gas pressure, as shown by the radial n e distribution. We identify a transition from a slow ionization wave to a fast ionization wave between 750 mTorr and 1 Torr. These experiments validate our analytic fit of n e , allowing the implementation of an LCIF diagnostic in helium without the need to develop a CRM.

42 ENGINEERING↗
DOE OSTI2021

Identification of crystal plasticity model parameters by multi-objective optimization integrating microstructural evolution and mechanical data

Crystal plasticity models evolve a polycrystalline yield surface using meso-scale descriptions of deformation mechanisms. The activation of deformation mechanisms is governed by crystallography and a set of model parameters, which are typically calibrated through the fitting of mechanical data such as stress–strain curves and elastic lattice strains. Microstructural data such as phase fractions and texture evolution are used for verifying crystal plasticity parameters. In this study, we use a multi-objective genetic algorithm to identify hardening parameters from flow stress curves with an option to incorporate texture into the optimization approach. Robust, generalized objective functions are developed and used to identify sets of parameters pertaining to dislocation density-based hardening laws in visco-plastic and elasto-plastic self-consistent (VPSC and EPSC) homogenization models. First, the parameters are identified for pure Nb directly from texture using an objective function based on generalized spherical harmonics. Since texture evolution is driven by the relative contribution of active slip systems, the parameters governing the evolution of slip resistance ratios can be recovered from fitting discrete textures at a series of strains. Next, a comprehensive set of load reversal data for dual phase (DP) 780 steel is used to fit a hardening law and a back-stress law in EPSC. Finally, parameters pertaining to a complex hardening law for the evolution of slip and twinning in pure α-Ti are identified. Remarkably, using texture as an objective in combination with stress–strain objectives constrains the model of Ti to fully reproduce not only stress–strain and texture evolution but also hierarchical twinning measurements as a function of initial grain size and texture. Furthermore, given an appropriate model fit to representative experimental texture evolution, underlying twin volume fractions contributing to texture evolution can be predicted.

42 ENGINEERING↗
DOE OSTI2024

Empirical Characterization and Modeling of Cohesive – to – Adhesive Shear Fracture Mode Transition due to Increased Adhesive Layer Thicknesses of Fiber Reinforced Composite Single – Lap Joints

Here, to ensure a strong adhesive bond, most standards and adhesive manufacturers specify a maximum adhesive gap of 1 mm when bonding fiber reinforced composite structures. In manufacturing large components, such as joining two halves of wind turbine blades, meeting this gap tolerance specification is impractical; gaps larger than 10 mm are common in large adhesively bonded composite structures using state-of-the-art manufacturing techniques. Currently, there is a lack of fundamental understanding of the failure mechanics of adhesive gaps larger than 3 mm. To create such understanding, glass fiber - acrylic thermoplastic composite panels bonded using different epoxy adhesives within single-lap joint samples with adhesive thicknesses of 0.1 mm, 0.3 mm, 1 mm, 3 mm, 5 mm, and 10 mm were sheared to failure. A transition from cohesive to adhesive failure was observed to occur about 1 mm to 3 mm joint thicknesses. Plotting the shear stress normalized by the ratio of the joint width to thickness as a function of the joint thickness normalized by the joint length is shown to result in the ability to fit simple empirically derived models of the cohesive-to-adhesive failure transition, regardless of the adhesive. Furthermore, using these normalized variables, all the observed cohesively failed specimens collapse to a single master curve, as do the adhesively failed specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scaling up carboxylic acid production from cheese whey and brewery wastewater via methane-arrested anaerobic digestion

In a circular economy, organic waste streams are valuable resources for sustainably producing chemicals and fuels. This study investigates a new methane-arrested anaerobic digestion (MAAD) process that converts high-strength cheese whey and brewery wastewater into carboxylic acids. The process was developed and optimized under various bench-scale semi-continuous (fed-batch) operating conditions (e.g., retention time, organic loading rate, pH, and feed/harvest frequency). The MAAD responses to various control-failure scenarios were also systematically investigated. The highest total acid productivity was 26g/(L liq ·d) with a substrate conversion of 0.79g COD digested /g COD fed at a hydraulic retention time (HRT) of approximately 2 d in a 14-L digester. The most stable conditions for digester operation (HRT 3 d at pH 6.0 and 40°C) were selected for process scale-up to 100gal (~380 L). Semi-continuous, pilot-scale MAAD successfully produced a total acid concentration of 40.6+/-1.1 g/L with 8.1% acetic acid, 45.1% butyric acid, and 44.5% lactic acid. The links between wastewater characteristics, operation mode, digester scale, and microbial community structure were statistically analyzed. The results show genera Sporolactobacillus and Clostridium positively correlate with butyric acid production (Pearson correlation coefficient>0.5). Moreover, four kinetic models were developed and fit to batch MAAD experimental datasets (R 2 >95%) and were successfully applied to predict the total acid production in both bench-and pilot-scale semi-continuous MAADs. In conclusion, this study shows MAAD has the potential for industrial-scale applications and is a robust platform to valorize low- or negative-value waste streams into high-value bioproducts.

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

Parasitic Effects in Impedance Spectrum of PEM Water Electrolysis Cells: Case Study of High-Frequency Inductive Effects

Electrochemical impedance spectroscopy (EIS) is a powerful tool to characterize and distinguish electrochemical, electrical, and diffusive processes in an electrolysis cell. The EIS response with small impedance (mΩ) depends on components and materials of the cell, but often also on the electrical setup, cables, and connectors, leading to incorrect conclusions about the performance of the electrolysis cells. These parasitic effects are assessed via a short-circuit measurement and confirmed to be external from the membrane electrode assembly (MEA). In the setup, the inductive characteristic is described by a modified inductive reactance and must be included in the equivalent circuit model (ECM) which is fitting the cell EIS spectra in operation (in situ). If points of the spectra showing inductive characteristic (below x-axis in Nyquist plot) are excluded from the ECM evaluation, an artificially increased ohmic resistance at high frequencies is obtained. In addition, considerations on other components of the ECM can be misleading and significantly incorrect. The only way to avoid that inductive behavior caused by the external setup is incorrectly assigned to the ECM describing properties of the MEA is to include the inductive effects in the ECM analysis of in situ EIS measurements.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Spectral line-shape in collinear laser spectroscopy after atomic charge exchange

Collinear laser spectroscopy experiments on fast, neutral beams have been extensively used for studies on short-lived radioactive nuclei, taking advantage of its high sensitivity. The resulting resonance line-shape is known to show significant distortion, due to the energy exchange during the charge-exchange neutralization process, which can cause large systematic uncertainty in the determined centroid. A model for the line shape was constructed and simulated to be compared to measured Al, Si, and Ni hyperfine spectra. It is shown that the distortion is caused mainly by the transfer of electron into many different energy levels in the projectile atom and subsequent decays, rather than secondary inelastic collisions, which were often assumed in the line shape analysis before. Furthermore, the model can also be applied to other projectile–alkali pairs, providing a reliable line-shape with less fitting parameters than conventional phenomenological models.

Atom↗
DOE OSTI2025

Harnessing ionic complexity: A modeling approach for hierarchical ionic circuit design

Since the 1950s, soft ionic devices have evolved from individual components to an expanding library of sensors, actuators, signal transmitters, and processors. However, integrating these components into complex, multifunctional systems remains challenging due to the nonintuitive and nonlinear interactions between ionic elements. In this work, we address these fundamental challenges by developing a lumped element model that enables interrogation of the physics that governs ionic circuits, as well as rapid design and optimization. Our model captures features specific to ionic charge carriers, while preserving the hierarchical design flexibility and computational efficiency of traditional circuit modeling. We demonstrate that our model can not only fit individual device behavior but also accurately predict the behavior of larger circuits formed by combining those devices. Additionally, we show how our tool utilizes the intrinsic nonlinearities of ionic systems to enable extended functionality, revealing how factors such as ion enrichment, ion leakage, and polymer charge density influence performance. Lastly, we present a fully ionic power supply, sensor, control system, and actuator for a soft robot that adapts its motion in response to environmental salt, illustrating the tool’s potential to accelerate advancements in chemical sensing, biointerfacing, biomimetic systems, and adaptive materials.

42 ENGINEERING↗
DOE OSTI2024

Application of the immobilized low-activity waste glass corrosion model to the static dissolution of 24 statistically-designed alkali-borosilicate waste glasses

Glass corrosion models that capture the complex mechanisms of the glass-water reaction enable the prediction of nuclear waste glass durability in disposal scenarios. Parameterization of such models is challenging because of the need to capture changes in corrosion behavior with time, reaction conditions, and glass composition. In this report we describe and employ the immobilized low activity waste (ILAW) glass corrosion model (IGCM) in geochemical simulations of static dissolution tests, at two temperatures (40 °C and 90 °C), for a matrix of 24 enhanced low activity waste (eLAW) glasses statistically designed to cover a processable composition space defined by 8 major glass components (Al 2 O 3 , B 2 O 3 , CaO, Na 2 O, SiO 2 , SnO 2 , ZrO 2 , and Others as defined in the text). The IGCM includes a first-order chemical affinity term and an ion-exchange term that represents the net exchange of Na+ ions in the pristine glass with protons in aqueous solution. Constant, time-dependent, and time- and pH-dependent functional forms of the ion-exchange term are evaluated to reproduce the change in corrosion behavior with time in saturated, static dissolution tests. The agreement with measured aqueous concentrations of the main glass components (B, Na, and Si) improved significantly upon addition of a time-dependent term, which therefore constitutes a simple representation of the glass-water reaction progress. Due to the limited changes in pH in the static dissolution tests, past a short initial period of rapid increase, addition of a pH dependent term did not appreciably improve the fits, indicating that comprehensive model parameterization requires more than one type of glass corrosion test to capture a wide range of solution chemistries. IGCM parameters were found to be dependent on glass composition, and the parameter sets generated in this work will enable the development of composition–parameter correlation models that offer the promise of predicting IGCM parameters, and thus glass corrosion behavior, solely based on glass composition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Controlling the pandemic during the SARS-CoV-2 vaccination rollout

There is a consensus that mass vaccination against SARS-CoV-2 will ultimately end the COVID-19 pandemic. However, it is not clear when and which control measures can be relaxed during the rollout of vaccination programmes. We investigate relaxation scenarios using an age-structured transmission model that has been fitted to age-specific seroprevalence data, hospital admissions, and projected vaccination coverage for Portugal. Our analyses suggest that the pressing need to restart socioeconomic activities could lead to new pandemic waves, and that substantial control efforts prove necessary throughout 2021. Using knowledge on control measures introduced in 2020, we anticipate that relaxing measures completely or to the extent as in autumn 2020 could launch a wave starting in April 2021. Additional waves could be prevented altogether if measures are relaxed as in summer 2020 or in a step-wise manner throughout 2021. We discuss at which point the control of COVID-19 would be achieved for each scenario.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

A machine learning framework for accurate and robust analysis of radiation detector pulses

The microscopic properties of atomic nuclei are used to study various scientific questions. They are essential for understanding the fundamental forces of nature and the chemical evolution of the universe. Detecting decay radiation from radioactive nuclei makes it possible to probe these fundamental nuclear properties. Detector waveform traces may contain additional information about the radiation. Generally, advanced signal processing techniques are needed to extract this additional information, often involving fitting the waveform with model response functions using non-linear least-squares optimization with second-order gradient methods. While this is a powerful technique, it is also computationally expensive, leading to slow processing time, which scales with the volume of data. To address this problem, we have developed a machine learning (ML) approach that infers the characteristics of traces from a model detector response function. In particular, we are interested in classifying whether a single recorded trace consists of one or two pulse constituents and estimating the pulse parameters. Furthermore, our proposed ML method can precisely extract the pulses’ parameters, such as energy and timing information, and accurately classify the pulse multiplicity of a trace. Unlike non-learning-based approaches, our ML approach uses neural networks that are significantly faster at inference, as they do not require any optimization during this stage.

Curve fitting↗
DOE OSTI2021

Predicting bulge to total luminosity ratio of galaxies using deep learning

ABSTRACT We present a deep learning model to predict the r-band bulge-to-total luminosity ratio (B/T) of nearby galaxies using their multiband JPEG images alone. Our Convolutional Neural Network (CNN) based regression model is trained on a large sample of galaxies with reliable decomposition into the bulge and disc components. The existing approaches to estimate the B/T ratio use galaxy light-profile modelling to find the best fit. This method is computationally expensive, prohibitively so for large samples of galaxies, and requires a significant amount of human intervention. Machine learning models have the potential to overcome these shortcomings. In our CNN model, for a test set of 20 000 galaxies, 85.7 per cent of the predicted B/T values have absolute error (AE) less than 0.1. We see further improvement to 87.5 per cent if, while testing, we only consider brighter galaxies (with r-band apparent magnitude <17) with no bright neighbours. Our model estimates the B/T ratio for the 20 000 test galaxies in less than a minute. This is a significant improvement in inference time from the conventional fitting pipelines, which manage around 2–3 estimates per minute. Thus, the proposed machine learning approach could potentially save a tremendous amount of time, effort, and computational resources while predicting B/T reliably, particularly in the era of next-generation sky surveys such as the Legacy Survey of Space and Time (LSST) and the Euclid sky survey which will produce extremely large samples of galaxies.

Grover, Harsh (ORCID:0000000321338142)↗
DOE OSTI2024

Modelling the BAO feature in bispectrum

We investigate how well a simple leading order perturbation theory model of the bispectrum can fit the baryon acoustic oscillation (BAO) feature in the measured bispectrum monopole of galaxies. Previous works showed that perturbative models of the galaxy bispectrum start failing at the wavenumbers of k ∼ 0.1 h Mpc −1 . We show that when the BAO feature in the bispectrum is separated, it can be successfully modelled up to much higher wavenumbers. We validate our modelling on GLAM simulations that were run with and without the BAO feature in the initial conditions. We also quantify the amount of systematic error due to BAO template being offset from the true cosmology. We find that the systematic errors do not exceed 0.3 per cent for reasonable deviations of up to 3 per cent from the true value of the sound horizon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Modeling Molecular Hydrogen in Low-metallicity Galaxies

We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which includes radiative transfer, a nonequilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H 2 abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sensitive to the metallicity of the gas with a progressively smoother density distribution with decreasing metallicity. In addition to the well-known trend of the H I–H 2 transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the ISM drops drastically at Z ≲ 0.2 Z ⊙ as the formation time of H 2 becomes much longer than a typical lifetime of dense regions of the ISM. We present accurate fitting formulae for both volumetric and projected ƒ H 2 measured on different scales as a function of gas metallicity, UV radiation field, and gas density. We show that when the formulae are applied to the patches in the simulated galaxy, the overall molecular gas mass is reproduced to better than a factor of ≲1.5 across the entire range of metallicities and scales. We also show that the presented fit is considerably more accurate than any of the previous ƒ H 2 models and fitting formulae in the low-metallicity regime. The fit can thus be used for modeling molecular gas in low-resolution simulations and semi-analytic models of galaxy formation in the dwarf and high-redshift regimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Spin Waves and Magnetic Exchange Hamiltonian in CrSBr

Abstract CrSBr is an air‐stable two‐dimensional (2D) van der Waals semiconducting magnet with great technological promise, but its atomic‐scale magnetic interactions—crucial information for high‐frequency switching—are poorly understood. An experimental study is presented to determine the CrSBr magnetic exchange Hamiltonian and bulk magnon spectrum. The A ‐type antiferromagnetic order using single crystal neutron diffraction is confirmed here. The magnon dispersions are also measured using inelastic neutron scattering and rigorously fit the excitation modes to a spin wave model. The magnon spectrum is well described by an intra‐plane ferromagnetic Heisenberg exchange model with seven nearest in‐plane exchanges. This fitted exchange Hamiltonian enables theoretical predictions of CrSBr behavior: as one example, the fitted Hamiltonian is used to predict the presence of chiral magnon edge modes with a spin‐orbit enhanced CrSBr heterostructure.

2D materials↗
DOE OSTI2023

Measurement of the branching fraction of ${\textrm{D}}_{\textrm{s}}^{+}$ → τ + ν τ via τ + → μ+ν μ $\overline{\nu}$τ

Utilizing 7 . 33 fb - 1 of e + e - collision data taken at the center-of-mass energies of 4.128, 4.157, 4.178, 4.189, 4.199, 4.209, 4.219, and 4.226 GeV with the BESIII detector, the branching fraction of the leptonic decay ${D}_s^{+}$ → τ + ν τ via τ + → μ + ν μ $\overline{\nu}$ τ is measured to be ${\mathcal{B}}_{D_s^{+}\to {\tau}^{+}{\nu}_{\tau }}$ = (5 . 37 ± 0 . 17 stat ± 0 . 15 syst )%. Combining this branching fraction with the world averages of the measurements of the masses of τ + and ${D}_s^{+}$ as well as the lifetime of ${D}_s^{+}$, we extract the product of the decay constant of ${D}_s^{+}$ and the c → s Cabibbo-Kobayashi-Maskawa matrix element to be ${f}_{D_s^{+}}$ |V cs | = (246 . 7 ± 3 . 9 stat ± 3 . 6 syst ) MeV. Taking |V cs | from a global fit in the standard model we obtain ${f}_{D_s^{+}}$ = (253 . 4 ± 4 . 0 stat ± 3 . 7 syst ) MeV. Conversely, taking ${f}_{D_s^{+}}$ from lattice quantum chromodynamics calculations, we obtain |V cs | = 0 . 987 ± 0 . 016 stat ± 0 . 014 syst .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Search for $t$-channel scalar and vector leptoquark exchange in the high-mass dimuon and dielectron spectra in proton-proton collisions at $\sqrt{s}=13$ TeV

A search for t-channel exchange of leptoquarks (LQs) is performed in dimuon and dielectron spectra using proton-proton collision data collected at $\sqrt{s}=13$ TeV with the CMS detector at the CERN LHC. The data correspond to an integrated luminosity of 138 fb −1 . Eight scenarios are considered, in which up or down quarks couple to muons or electrons via a scalar or vector LQ exchange, for dilepton invariant masses above 500 GeV. The LQ masses are probed up to 5 TeV, beyond a regime probed by previous pair-production and single-production searches. The differential distributions of dilepton events are fit to templates that model the nonresonant LQ exchange and various standard model background processes. Limits are set on LQ-fermion coupling strengths for scalar and vector LQ masses in the 1–5 TeV range at 95% confidence level, establishing stringent limits on first- and second-generation LQs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗