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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Bemnifosbuvir: An HCV NS5B Inhibitor With Multiple Modes of Action

Bemnifosbuvir (BEM) is a potent, pan-genotypic inhibitor targeting the hepatitis C virus (HCV) NS5B polymerase. Its antiviral activity was evaluated in an ascending dose phase I clinical trial involving 30 patients treated once a day for 7 days. After treatment initiation, plasma HCV RNA declined in a biphasic manner with a mean reduction of 2.3 log 10 IU/mL after 24 hours and 4.4 log10 IU/mL by day 7 for the highest dose. Alanine aminotransferase (ALT) also normalized in most patients. We employed a multiscale mathematical model fitted to the HCV RNA and ALT dynamics to quantify the antiviral activity and evaluate the modes of action of BEM. We found that models in which BEM only acted as a typical HCV RNA polymerase inhibitor and reduced the intracellular production of HCV RNA did not fit the data as well as models in which BEM had multiple modes of action, including suppressing viral assembly and secretion and enhancing intracellular HCV RNA degradation. BEM's effectiveness in inhibiting intracellular HCV RNA production increased with dose (150 mg/day: 88.2%, 300 mg/day: 98.8%, 600 mg/day: 99.5%), while inhibition of viral assembly and release was ~95% effective regardless of dose. We observed a dose-dependent enhancement in the degradation of intracellular HCV RNA, with degradation rates 1.5-fold higher in patients receiving 300 mg/day and 2.7-fold higher in those receiving 600 mg/day than in patients receiving 150 mg/day. No significant differences in antiviral activity were detected between HCV genotypes 1b and 3 or between patients with and without compensated cirrhosis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Compressional wave velocity and effective stress in unsaturated soil: Potential application for monitoring moisture conditions in vadose zone sediments

Compressional seismic (P-wave) velocities were measured experimentally at different water potentials and confining pressure states using a variety of sediments collected from the US Department of Energy Hanford Site in southeastern Washington State. P-wave velocity was measured in variably saturated sediments using ultrasonic piezoelectric transducers and water potential was measured with a heat dissipation sensor. We propose a model to relate changes in P-wave velocity to soil stress conditions using an effective stress formulation that includes capillary stress as well as adsorptive stress. We show that compressional wave velocities can be related to water potential and confining stress using a model with five fitted parameters. A single set of parameters provides an excellent fit to data acquired from all of the sediment samples. Under very dry conditions, corresponding to low water potentials, water flow occurs in adsorptive thin films associated with very low unsaturated hydraulic conductivities. The results presented here suggest the utility of seismic methods for establishing low water potential and flow conditions within variably saturated Hanford formation sediments and providing feedback on the performance of mitigation efforts such as soil desiccation and surface infiltration barriers. Seismic methods can provide minimally invasive, three-dimensional information on relevant flow conditions in vadose zone sediments as well as changes that occur over time.

Linneman, Dorothy C.↗
DOE OSTI2021

High-pressure crystal structure and equation of state of ferromagnesian jeffbenite: implications for stability in the transition zone and uppermost lower mantle

Jeffbenite, ideally Mg 3 Al 2 Si 3 O 12 , has been identified as inclusions in super-deep diamonds originating from depths that exceed 300 km. Although Mg-end member jeffbenite has limited stability at upper-mantle conditions, iron-bearing jeffbenite may have broader P–T stability that extends to the transition zone or uppermost lower mantle, incorporating significant amounts of ferric iron. Using synchrotron-based, single-crystal X-ray diffraction (XRD) and synchrotron Mössbauer spectroscopy (SMS) at pressures up to 29 GPa, we report the crystal structure, compressibility, and likely spin transition of iron in ferromagnesian jeffbenite (Mg 2.32 Al 0.03 Fe 2+ 1.28 Fe 3+ 1.77 Si 2.85 O 12 ). High-pressure structure refinements reveal that Fe 3+ substitution for Si in the T2 site, which shares edges with the M2 octahedron, likely stabilizes jeffbenite at high pressure, because it increases the cation-to-cation distance between these sites. Although ferromagnesian jeffbenite does not undergo a structural phase transition below 30 GPa, SMS hyperfine parameters suggest the onset of an electronic spin transition of iron from high-spin (HS) to low-spin (LS) at around 22 GPa, which may increase its stability at high pressures. Pressure–volume data were fit to a third order Birch–Murnaghan equation of state, resulting in V o = 816.54(9), K To = 181.54(1.39), and K' T0 = 2.76(14). These equation of state parameters are applicable to evaluating the encapsulation pressures of super-deep diamonds. The density and bulk modulus of ferromagnesian jeffbenite are similar to or higher than pyrope–almandine, pyrope–majorite, and skiagite–majorite solid solution garnets, further suggesting that jeffbenite may be an important ferric–iron silicate in the deeper parts of the mantle transition zone and uppermost lower mantle. Furthermore, future studies on the influence of temperature and oxidation state on the stability and equations of state of iron-bearing jeffbennite are still needed to determine what role, if any, jeffbenite plays in transition-zone mineralogy.

58 GEOSCIENCES↗
DOE OSTI2024

Stochastic average model methods

We consider the solution of finite-sum minimization problems, such as those appearing in nonlinear least-squares or general empirical risk minimization problems. We are motivated by problems in which the summand functions are computationally expensive and evaluating all summands on every iteration of an optimization method may be undesirable. Here we present the idea of stochastic average model (SAM) methods, inspired by stochastic average gradient methods. SAM methods sample component functions on each iteration of a trust-region method according to a discrete probability distribution on component functions; the distribution is designed to minimize an upper bound on the variance of the resulting stochastic model. We present promising numerical results concerning an implemented variant extending the derivative-free model-based trust-region solver POUNDERS, which we name SAM-POUNDERS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

In-situ irradiation-induced studies of grain growth kinetics of nanocrystalline UO 2

The thermal conductivity of UO 2 fuel needs to be high enough to dissipate the heat generated from the fission reaction. Since grain size affects thermal conductivity and grain size can evolve with irradiation, it is critical to understand in-reactor UO 2 grain growth. Most studies of grain growth in UO 2 are based on thermally driven processes at elevated temperatures. However, studies have shown that grain growth can occur even at cryogenic temperatures by ballistic processes. Such irradiation-induced grain growth in UO 2 is yet to be studied. Advanced in-situ Kr ion irradiation and transmission electron microscopy were systematically performed on nanocrystalline UO 2 thin films at temperatures ranging from 50 K to 1073 K; grain growth was observed at all temperatures. A combination of manual and machine learning techniques was used to measure and plot grain size evolution against irradiation fluence at various irradiation temperatures. The machine learning method has significantly improved the analysis efficiency and reduced human labors. The grain diameter data were fitted using classical grain growth and thermal spike models to describe grain growth kinetics with and without irradiation effect. Grain growth during low temperature irradiation (≤ 475 K) can be well described by the thermal spike model. Above 475 K, there were additional thermally assisted processes that further accelerate the grain growth. At the highest irradiation temperature about 1075 K, both irradiation-induced dislocation loops and cavities/bubbles were observed to form in the UO 2 . In this report, the effects of irradiation-induced defects on grain growth kinetics are discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Improved natural convection heat transfer correlations for reactor cavity cooling systems of high-temperature gas-cooled reactors: From computational fluid dynamics to Pronghorn

The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.

42 ENGINEERING↗
DOE OSTI2023

Piecewise linear approximation with minimum number of linear segments and minimum error: A fast approach to tighten and warm start the hierarchical mixed integer formulation

In several areas of economics and engineering, it is often necessary to fit discrete data points or approximate nonlinear functions with continuous functions. Piecewise linear (PWL) functions are a convenient way to achieve this. PWL functions can be modeled in mathematical problems using only linear and integer variables. Moreover, there is a computational benefit in using PWL functions that have the least possible number of segments. This work proposes a novel hierarchical mixed integer linear programming (MILP) formulation that identifies a continuous PWL approximation with minimum number of linear segments for a given target maximum error. The proposed MILP formulation also identifies the solution with the least maximum error among the solutions with minimum number of segments. Then, this work proposes a fast iterative algorithm that identifies non necessarily continuous PWL approximations by solving O(S log N) linear programming (LP) problems, where N is the number of data points and S is the minimum number of segments in the non necessarily continuous case. This work demonstrates that tight bounds for the MILP problem can be derived from these approximations. Next, a fast algorithm is introduced to transform a non necessarily continuous PWL approximation into a continuous one. Finally, the tight bounds and the continuous PWL approximations are used to tighten and warm start the MILP problem. The tightened formulation is shown in experimental results to be more efficient, especially for large data sets, with a solution time that is up to two orders of magnitude less than the existing literature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗
DOE OSTI2024

A model to separate conduction and radiation in high temperature laser flash measurements for semi-transparent materials

Laser flash analysis (LFA) is a widely used technique in the scientific and industrial communities for measuring the thermal diffusivity and conductivity of materials. However, traditional fitting methods based solely on heat conduction do not account for radiation heat transfer in transparent and semi-transparent materials, which becomes significant at high temperatures. In this study, we have developed a coupled radiation-conduction model that considers the optical properties of materials to accurately extract thermal diffusivity from LFA measurements. We conducted LFA measurements on fused silica glass samples of varying thicknesses and within the temperature range of 25 to 1000 °C and fitted the data with the model. Since the coupled model requires optical properties as input parameters, which are often unavailable for new materials, we have also proposed a modified transparent model. Compared to the coupled model, results from the modified transparent model generally overestimate thermal conductivity, but the errors are acceptable in silica glass with small thickness. Here, we also investigated the applicability of the modified transparent model with materials of different optical properties.

14 SOLAR ENERGY↗
DOE OSTI2021

Parameterization of the Stoner-Wohlfarth model of magnetic hysteresis

The Stoner-Wohlfarth (SW) is the most used model of magnetic hysteresis, but its computation is time-consuming. We approximate piecewise this model by easy-to-compute analytic functions. Overall, our parametrization is suitable for fast quantitative evaluations and fitting experimental data, which we exemplify. Applicability of this parametrization is as broad as that of the SW model, which can be extended to include more parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Lamellae of α phase in a metastable β -Ti alloy studied by small-angle x-ray scattering

The shapes and orientations of lamellae of the hexagonal close-packed α phase in single crystals of the body-centred cubic β phase in a metastable β -Ti alloy (Timetal LCB) were investigated by synchrotron small-angle x-ray scattering. Several steps of thermal treatment with various annealing times and temperatures were employed. The measured scattering data were fitted to a model of ellipsoidal lamellae of random sizes, and the mean lamella shapes and orientations of their axes with respect to the β -Ti matrix were determined from the fit. The number and relative volume of the lamellae were estimated from the integrated scattered intensity. The elastic energy of the lamellae as well as the energy of defects at the α / β interface (misfit dislocations) were calculated from the parameters of the lamellae. In conclusion, the results provide information on the kinetics of the β $ \rightarrow$ α transformation in metastable β -Ti alloys and allow assessing the degree of thermodynamic equilibrium of the final α / β composite.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Relationship between ductility and the porosity of additively manufactured AlSi10Mg

Additive manufacturing via selective laser melting can result in variable levels of internal porosity both between build plates and within components from the same build. In this work, sample porosity levels were compared to tensile properties for 176 samples spanning eight different build plates. Sample porosity was measured both by Archimedes density, which provided an estimation of overall porosity, and by observation of voids in the fracture surface, which provided an estimation of the porosity at the failure plane. The porosity observed at the fracture surface consistently demonstrated higher porosity than that suggested by Archimedes density. The porosity values obtained from both methods were compared against the mechanical results. Sample porosity appears to have some correlation to the ultimate tensile strength, yield strength, and modulus, but the strongest relationship is observed between porosity and ductility. Three different models were used to relate the fracture surface porosity to the ductility. The first method was a simple linear regression analysis, while the other two models have been used to relate porosity to ductility in cast alloys. It is shown that all three models fit the data well over the observed porosity ranges, suggesting that the models taken from casting theory can extend to additively manufactured metals. Finally, it is proposed that the non-destructive Archimedes method could be used to estimate an approximate sample ductility through the use of correlations realized here. Such a relationship could prove useful for design and for a deeper understanding of the impact of pores on tensile behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Impact of light output on the timing resolution of organic glass scintillator bars in a dual-ended readout configuration

The effectiveness of detector system modeling and validation depends on the quality of characterization performed on the system. For nuclear nonproliferation applications, which require detectors to address complex and variable field conditions, access to accurate system models is essential. This study presents the timing characterization of organic glass scintillator bars used in an imaging system developed at the University of Michigan. 137 Cs and 60 Co gamma-ray sources were used to determine the coincidence time resolution as a function of measured light output for two Compton scatter events between two organic glass scintillator bars. Timing resolution for two events, each with light output ranging from 100 to 1100 keVee were characterized. The resulting data were fit to a parametric function that was validated to agree with experimental results within ±25 ps. Simulations were performed to establish appropriate calibration points for each organic scintillator, thereby ensuring accurate calibration of light output values during analysis. This work successfully characterized the light output dependence of the coincidence timing resolution of a pair of organic glass detectors for use in an imaging system. Depending on the amount of scintillation light output from the events, the full width at half maximum of the measured coincidence time resolution ranged from 653.2 ± 16.1 ps for low light yields, 100 keVee, to 121.0 ± 8.4 ps for higher yields at 1100 keVee. The insights obtained from the timing resolution behavior will allow for accurate simulation capabilities in future security and verification efforts using scatter-based imaging systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Implications of styrenic triblock copolymer gel mechanics on midblock bridging fraction

Triblock copolymer gels have garnered considerable scientific interest over the past few decades and are used in a variety of applications including consumer cushioning, pressure-sensitive adhesives, and ballistics gels. While many of their applications rely upon understanding mechanical properties, mechanical studies of block copolymer gels are relatively disjointed with each focusing on specific copolymers, which inherently have a fixed molecular weight and composition. Here, the present study examines the quasistatic mechanical response of styrenic triblock copolymer gels composed of seven unique triblock copolymers at various concentrations. Resultant stress-extension data is fitted with the slip-tube network (STN) theory that describes gel mechanics based upon crosslinked network (G c ) and chain entanglement (G e ) modulus contributions. Collectively, modulus contributions imply that midblock bridging is independent of copolymer identity (i.e., molecular weight and block fraction) and dependent on triblock copolymer concentration via F b ∝ φ ABA 0.91±0.39 over the concentration range examined (φ ABA = 0.05–0.39). Gels composed of a complementary diblock-triblock copolymer pair were subsequently analyzed to explore the influence of interlocking loop-loop and loop-bridge pairs within the effective midblock bridging population. These results point to the presence of a physically significant quantity of interlocked loops, but quantitative analysis was inconclusive.

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

OzDES Reverberation Mapping of Active Galactic Nuclei: Final Data Release, Black-Hole Mass Results, & Scaling Relations

Over the last decade, the Australian Dark Energy (OzDES) collaboration has used Reverberation Mapping to measure the masses of high redshift supermassive black holes. Here we present the final review and analysis of this OzDES reverberation mapping campaign. These observations use 6-7 years of photometric and spectroscopic observations of 735 Active Galactic Nuclei (AGN) in the redshift range 0.13-3.85 and bolometric luminosity range 44.3 - 47.5 erg/s. Both photometry and spectra are observed in visible wavelengths, allowing for the physical scale of the AGN broad line region to be estimated from reverberations of the H{̱e̱ṯa̱}̱, MgII and CIV emission lines. We successfully use reverberation mapping to constrain the masses of 62 super-massive black holes, and combine with existing data to fit a power law to the lag-luminosity relation for the H{̱e̱ṯa̱}̱ and MgII lines with a scatter of ~0.25 dex, the tightest yet identified, fit specifically for consistency with high redshift AGN. We fit a similarly constrained relation for CIV, resolving a tension with the low luminosity literature AGN by accounting for selection effects arising from finite survey length. We also examine the impact of emission line width and luminosity (related to accretion rate) in reducing the scatter of these scaling relationships and find no significant improvement over the lag-only approach for any of the three lines. Using these relations, we further estimate the masses and accretion rates of 246 AGN with single epoch methods. We also use these relations to estimate the relative sizes of the H{̱e̱ṯa̱}̱, MgII and CIV emitting regions. In short, we provide a comprehensive benchmark of high redshift AGN reverberation mapping at the close of this most recent generation of surveys, including light curves, time-delays, and a set of significantly improved radius-luminosity relations for use with high-redshift populations.

McDougall, Hugh [Queensland U.]↗
DOE OSTI2023

Benchmarking magnetised three-wave coupling for laser backscattering: analytic solutions and kinetic simulations

Understanding magnetised laser–plasma interactions is important for controlling magneto-inertial fusion experiments and developing magnetically assisted radiation and particle sources. For nanosecond pulses at non-relativistic intensities, interactions are dominated by coherent three-wave interactions, whose nonlinear coupling coefficients became known only recently when waves propagate at oblique angles with the magnetic field. In this paper, backscattering coupling coefficients predicted by warm-fluid theory are benchmarked using particle-in-cell simulations in one spatial dimension, and excellent agreements are found for a wide range of plasma temperatures, magnetic field strengths and laser propagation angles, when the interactions are mediated by electron-dominant hybrid waves. Systematic comparisons between theory and simulations are made possible by a rigorous protocol. On the theory side, the initial boundary value problem of linearised three-wave equations is solved, and the transient-time solutions allow the effects of growth and damping to be distinguished. On the simulation side, parameters are carefully chosen and calibration runs are performed to ensure that comparisons are well controlled. Fitting simulation data to analytical solutions yields numerical growth rates that match theory predictions within error bars. Although warm-fluid theory is found to be valid for a wide parameter range, genuine kinetic effects have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Conductivity Spectroscopy for Investigation and Discovery of Photovoltaic Materials

Conductivity spectroscopy is an extremely powerful set of methods for probing the properties of optoelectronic materials, especially photovoltaics, where photoconductivity is one of the best spectroscopic proxies for performance. Despite this power, they are substantially less commonly used than time-resolved photoluminescence (for instance) because they tend to be more expensive to implement (THz) and/or require specialized knowledge (GHz) to construct instruments, which are not widely available. The goal of this review is to illustrate the utility of these experiments in the discovery and study of photovoltaic absorber materials and simultaneously make them more accessible to the community by providing a central tutorial resource. We provide a comprehensive review of how conductivity spectroscopy has developed over the past decade and been applied in the discovery and development of photovoltaic materials, with a primary focus on emerging solution-processable technologies. Along the way we aim to demystify conductivity spectroscopy with focused tutorial sections that explain the physical models used to fit the data and illustrate how to think about “high-frequency conductivity”.

14 SOLAR ENERGY↗