Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 739 records · Page 41

Physics-informed machine learning modeling for predictive control using noisy data

Due to the occurrence of over-fitting at the learning phase, the modeling of chemical processes via artificial neural networks (ANN) by using corrupted data (i.e., noisy data) is an ongoing challenge. Therefore, this work investigates the effect of both Gaussian and non-Gaussian noise on the performance of process-structure based recurrent neural networks (RNN) models, which take the form of partially-connected RNN models in this work, that are used to approximate a class of multi-input-multi-outputs nonlinear systems. Furthermore, two different techniques, specifically Monte Carlo dropout and co-teaching, are utilized in the development of partially-connected RNN models. Here, these two techniques are employed to reduce the over-fitting in ANNs when noisy data is used in the training process and, hence, to improve the open-loop accuracy as well as the closed-loop performance under a Lyapunov-based model predictive controller (MPC). Aspen Plus Dynamics, a well-known high-fidelity process simulator, is used to simulate a large-scale chemical process application in order to demonstrate the anticipated improvements in both open-loop approximation and closed-loop controller performance in the presence of Gaussian and non-Gaussian noise in the data set using physics-informed RNNs.

97 MATHEMATICS AND COMPUTING↗

Adsorption of terbium (III) on DGA and LN resins: Thermodynamics, isotherms, and kinetics

Two commercially available extraction chromatography (EXC) resins containing N,N,N’,N’-tetra-n-octyldiglycolamide (DGA Resin, Normal, 50 – 100 μm) and Bis(2-ethylhexyl) phosphate (LN Resin, 100 – 150 μm) were used as adsorbents to study fundamental adsorption properties such as thermodynamic values, equilibrium isotherms, and kinetic uptake models for terbium(III) adsorption. Weight distribution ratios (D w ) for terbium on DGA and LN resins were measured using a [ 160 Tb]Tb 3+ radiometric tracer in nitric acid as a function of acidity, temperature, initial analyte concentration, and equilibrium time. The D w values showed increasing binding affinity for DGA resin at high nitric acid concentrations and decreasing binding affinity for LN resins. Thermodynamic studies for DGA and LN resins revealed that the Gibbs free energy (ΔG) increased consistently with temperature. To model equilibrium data, increasingly higher parameter equilibrium isotherm models (Henry (1) < Langmuir, Freundlich (2) < Redlich-Peterson (3) < Fritz-Schluender (4)) were compared on their root mean squared errors (RMSE) and adjusted determination coefficients to determine the most applicable model. In all cases, the empirical four-parameter Fritz-Schluender isotherm demonstrated a superior fit. Similar comparisons for reaction-based kinetic models (Pseudo-first-order < Pseudo-second-order < Pseudo-n-order) revealed that the higher-order PNO model yielded a superior fit of kinetic data for both resins. Furthermore, in some cases, adsorption isotherms and kinetic models could also be modeled by a lower-order model with minimal change in error parameters. Weber-Morris plots revealed that two linear sections are observed for each resin, where the first linear segment is attributed to fast (film diffusion) adsorption of terbium, followed by slower intraparticle diffusion of terbium through the pores as the rate-limiting step. Based on the Weber-Morris plot, both film and intraparticle diffusion are involved in controlling the kinetic rate of adsorption for DGA and LN resins.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Group additivity values for entropy and heat capacities of C 2 –C 8 alkanes, alkyl hydroperoxides, and their radicals

Group additivity values for the thermodynamic properties of oxygenated radicals are poorly determined due to the absence of high quality reference data. Here, a set of 58 group additive values (GAV) for the standard entropy and heat capacity of relevance to alkane oxidation is derived from fits to an extensive and accurate database of standard entropies (298.15 K) and heat capacities (300–3000 K) recently calculated with the “STAR-1D” formalism. The 192 species in this database represent the alkanes (RH), alkyl radicals ($\dot{R}$), alkyl hydroperoxides (RO 2 H), alkyl-peroxy (R$\dot{O}$ 2 ) and hydroperoxy-alkyl ($\dot{Q}$OOH) radicals for all of the isomers of C 2 –C 5 alkane fuels and a select number of isomers of C 6 –C 9 species. The STAR-1D thermochemical data for this set of species was previously obtained from a coupling of scaled B2PLYPD3/cc-pVTZ vibrational analyses with scaled ωB97X-D/cc-pVTZ one-dimensional hindered rotor corrections. The 2σ uncertainties in the GAV results relative to the STAR-1D data set are 2.4 cal K –1 mol –1 for the entropies and at most 2.0 cal K –1 mol –1 for the heat capacities in the temperature range 500–800 K. The 2σ fitting uncertainties in the heat capacity gradually reduce at higher temperatures reaching a value of 0.8 cal K –1 mol –1 at 2000 K decreasing to only 0.5 cal K –1 mol –1 at 3000 K. The high degree of accuracy for the GAV representations is obtained through the introduction of various new group terms, together with the re-optimization of existing group terms. Among the full set of 58 GAV terms, 25 include non-next-nearest neighbor interactions (NNI) and β-corrections. The updated GAVs can be applied in the prediction of entropies and heat capacities for a wide range of hydrocarbons and hydroperoxide species and their radicals, which is important to the accurate prediction of fuel reactivity at low-temperatures in the range 600–1000 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗

Qualitative trend analysis based on a mixed-integer representation

Shape constrained spline fitting is a useful method to impose prior knowledge onto flexible semi-parametric models during parameter estimation. Most typically, the function shape is imposed through order restrictions on the regression coefficients. The intended shape is considered known or selected based on heuristic rules. In this study, we present a method to estimate the optimal set of order restrictions to segment a univariate data series into episodes with distinct shapes. This is also known as the qualitative trend analysis (QTA) problem. The obtained solution uses a trade-off between lack-of-fit and model complexity. Further, our practical implementation takes inspiration from the generalized order restricted information criterion (GORIC) for inequality-constrained model selection. From this, one learns (a) that QTA can be formulated as a mixed-integer quadratic program (MIQP) and (b) that the newly proposed mixed order restricted information criterion (MORIC) enables optimal segmentation. This is illustrated through didactic case studies.

42 ENGINEERING↗

Extracting single fiber transverse and shear moduli from off-axis misalignment fiber tensile testing

Small diameter (<100 μm) fibers (e.g. carbon fibers, Kevlar, and fiberglass) and wires (e.g. ultrafine copper and aluminum wires) are frequently used in many different engineering applications, such as for light weighting structures, electromagnetic shielding for aircraft/infrastructure/EVs, vibration damping, biological sensors, aerospace electrical devices, and electric windings just to name a few. Due to the manufacturing process, the fibers and wires are pulled and stretched to produce a preferential alignment. Therefore, thin fibers and wires typically display different properties along the length of the fiber as opposed to their cross section and many fibers/wires are considered transversely isotropic. The axial properties of fibers/wires can be ascertained via tensile testing of single-filaments or fiber tows, but the radial properties require much more effort to measure. Knowing these properties is important for the accurate prediction of micromechanical models and manipulation of fibers during micromanufacturing. In this paper, a new technique was developed to determine the transverse/shear moduli and strength of a material by conducting tensile tests of the material at increasing misalignment angles from the tensile axis. Due to the transversely isotropic nature of the material, the transverse/shear moduli and strength influence the experimental results recorded by the test machine to different degrees based on the amount of misalignment in the test setup. An equation was derived to determine the influence of each of the material properties based on the misalignment angle by manipulating the stiffness matrix for transversely isotropic materials using the transformation matrices. Then, curve fitted coefficients were used to identify the material properties. Here, the strengths were similarly determined by curve fitting an off-axis Tsai-Hill failure criteria to determine the influence of transverse, shear, and tensile strengths based on the complex loading condition provided by the off-axis tensile test. Zoltek Panex 35 carbon fibers were used to demonstrate this new technique and the determined properties were then compared to those obtained from nanoindentation and from literature. Fracture surfaces provide insight into the different failure mechanisms at various misalignment angles.

36 MATERIALS SCIENCE↗

Processing and mechanical characterization of short carbon fiber-reinforced epoxy composites for material extrusion additive manufacturing

Fiber-reinforced polymer composites have been extensively utilized in recent years as feedstock materials for material extrusion additive manufacturing (AM) processes to improve strength, stiffness, and functionality of printed parts over unfilled printed polymers. However, the widespread adoption of AM of fiber-reinforced polymer composites requires a deeper understanding of the process-structure-property relationships in printed components, and such relationships are not well understood yet. Fiber length is critically important to the mechanical performance of short fiber composites, but very few studies to-date have focused on how the fiber length distribution (FLD) evolves during processing of composite feedstocks and how this evolution affects printing behavior and mechanical properties in 3D-printed composites. Here, FLD is measured for carbon fiber reinforced epoxy composites over a wide range of ink compositions and shear mixing times, and the distributions are fit with a Weibull-type distribution function. The effects of FLD on the tradeoff between ink processability, ink rheology, printing behavior and mechanical properties are investigated. Furthermore, the effects of printing parameters (nozzle size and print speed) on mechanical anisotropy and fiber orientation distribution (FOD) in printed composites are explored. Mechanical properties of printed composites are characterized via 3 pt-flexural testing, and microstructure is investigated using optical and scanning electron microscopy (SEM), and x-ray computed tomography. Finally, the fitted Weibull parameters are fed into a composite model that incorporates FLD and FOD, and model predictions are found to be in excellent agreement with experimental observations.

3D printing↗

Data-driven surrogate modeling of hPIC ion energy-angle distributions for high-dimensional sensitivity analysis of plasma parameters' uncertainty

In this work, we present a data-driven strategy for effective construction of a surrogate model in high-dimensional parameter space for the ion energy-angle distribution (IEAD) output of hPIC simulations of plasma-surface interactions. The methodology is based on a bin-by-bin least-squares fitting of the IEAD in the parameter space. The fitting is performed in a transformed coordinate system to normalize the IEAD, and it employs sparse grids for sampling the parameter space to overcome sampling challenges in high dimensions. The surrogate model is significantly cheaper computationally than direct hPIC simulations yet maintains high fidelity to them, providing a fast emulator for hPIC simulations. Sensitivity analysis based on the surrogate model is utilized to characterize the dependence of the ion impact angle and energy moments on the physical parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerating Markov Chain Monte Carlo sampling with diffusion models

Global fits of physics models require efficient methods for exploring high-dimensional and/or multimodal posterior functions. We introduce a novel method for accelerating Markov Chain Monte Carlo (MCMC) sampling by pairing a Metropolis-Hastings algorithm with a diffusion model that can draw global samples with the aim of approximating the posterior. We briefly review diffusion models in the context of image synthesis before providing a streamlined diffusion model tailored towards low-dimensional data arrays. We then present our adapted Metropolis-Hastings algorithm which combines local proposals with global proposals taken from a diffusion model that is regularly trained on the samples produced during the MCMC run. Our approach leads to a significant reduction in the number of likelihood evaluations required to obtain an accurate representation of the Bayesian posterior across several analytic functions, as well as for a physical example based on a global fit of parton distribution functions. Our method is extensible to other MCMC techniques, and we briefly compare our method to similar approaches based on normalising flows. A code implementation can be found at https://github.com/NickHunt-Smith/MCMC-diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A new database website for nuclear level densities

We introduce a new open-access, web-based database (http://nld.ascsn.net), Current Archive of Nuclear Density of Levels (CANDL), that hosts experimental nuclear level density (NLD) datasets from a variety of techniques and energy ranges. Built using the Dash framework in Python, the database is designed to be interactive and user-friendly, allowing researchers to search, visualize, fit, and export NLD data with minimal effort. This resource includes data extracted from evaporation spectra, Oslo method variants, and other experimental techniques that cover excitation energies beyond the neutron resonance region. The database supports on-the-fly fitting with two widely-used phenomenological models—the Constant Temperature (CT) model and the Back-Shifted Fermi Gas (BSFG) model—selected for their simplicity and computational efficiency. Future versions aim to include additional datasets and model types, as well as easy-to-use interfaces to data science techniques. Here, this platform offers a vital tool for the nuclear physics, astrophysics, medicine, and reactor design communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Testing α -attractor quintessential inflation against CMB and low-redshift data

Due to universality and attractor properties, α-attractor quintessential inflation establishes direct relations between inflationary observables such as the scalar tilt n s and the tensor-to-scalar ratio r, and late-time dark energy equation of state parameters w 0 and w a . In this work, we examine three different physically motivated regimes, considering complete freedom in the parameter α, models inspired by supergravity where α takes on values up to α = 7/3, and Starobinsky inflation (α = 1). We investigate the consistency and constraints imposed by Cosmic Microwave Background measurements from the Planck satellite, B-mode polarization data from the BICEP/Keck collaboration, and low-redshift observations. Additionally, we consider small-scale CMB measurements released by the Atacama Cosmology Telescope, which give results approaching the Harrison– Zel’dovich spectrum (n s ≈ 1). Here α-attractors lead to an improved fit over $\Lambda$CDM. For the large-scale CMB measurements, α ≳ 2 models can provide equally good fits as $\Lambda$CDM

79 ASTRONOMY AND ASTROPHYSICS↗

High salinity seawater boiling point elevation: Experimental verification

Boiling point elevation (BPE) for high salinity solutions is a crucial parameter that significantly affects the overall thermal efficiency of desalination, water treatment, and many other industrial processing systems, especially when minimum or zero liquid discharge is a key requirement. This paper addresses the discrepancies among the empirical correlations reported in the literature that were commonly used to determine BPE for seawater solutions. An evaluation of the experimental data with low solution concentration reveals that the data cannot be extrapolated to determine the BPE at high concentration (X > 140 ppt or g/kg). Besides, disparities are established within different equations used to determine the BPE for high salinity seawater. As a result of the minimal availability of experimental data of BPE for X > 140 ppt and the need to eliminate the discrepancies, an experimental setup was built, which can produce reliable and repeatable BPE data for any solution concentration. Boiling points were obtained experimentally for seawater solutions with a wide range of salinity (X between 20 and 280 ppt) and brine solution temperature (t b between 20 and 80°C). The experimental results of boiling point elevations for seawater are further compared with the experimental data and analytical results reported in the literature. The newly obtained experimental results show that the minimum and maximum BPE values for the brine solutions tested are 0.175°C at t b = 80°C and X = 20 ppt, and 11.6°C at t b = 80°C and X = 280 ppt, respectively. Furthermore, the measured BPE data for seawater concentration over 80 ppt are fitted and presented using an appropriate fitting model. The new correlation can accurately predict the BPE for high salinity seawater brines and thus should be very helpful in the design of high-energy-efficient evaporation systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Electron and proton conductivity of Fe-N-C cathodes for PEM fuel cells: A model-based electrochemical impedance spectroscopy measurement

Impedance spectra of a PEM fuel cell with three Fe-N-C cathodes have been measured under the H 2 /N 2 testing regime. The spectra have been fitted using a recently developed physics-based impedance model, which takes into account finite proton ($σ_p$) and electron ($σ_e$) conductivity of the oxygen-free cathode catalyst layer. Fitting allowed to extract numerical data for $σ_p$, $σ_e$, the double layer capacitance, and the inductance of cables used for measuring impedance spectra. The values of σp and σe are close to what previously found for standard Pt/C electrodes, which is found for the first time using PGM-free catalysts. The method enables simultaneous measurement of reference proton and electron conductivity of PEMFC cathode.

25 ENERGY STORAGE↗

Constraint-dependent CTOA determination for stable ductile crack growth

The crack tip opening angle (CTOA) has been used as a reliable fracture toughness to characterize stable crack growth for thin-wall structures in low-constraint conditions. Recently, it is found that CTOA can also be utilized as a robust fracture parameter to describe the arrest fracture toughness for gas transmission pipelines in ductile steels. For better understanding the concept of CTOA, this paper presents a brief review on the CTOA definitions and its test methods. And then, four CTOA models are developed for determining the CTOA toughness using conventional single edge notched bend (SENB) specimens, including the load–displacement linear fit model, the logarithmic load–displacement linear fit model, the stable tearing energy model, and the J-differentiation model. Here fracture test data from SENB specimens in A285 carbon steel are employed to evaluate the proposed CTOA models. The results show that these CTOA models can determine either identical or comparable critical CTOA values for stable ductile crack growth using the SENB specimens. A constraint-corrected CTOA resistance curve and a constraint-corrected critical CTOA equation are also obtained for A285 carbon steel.

42 ENGINEERING↗

The direct and indirect effects of copper on vector-borne disease dynamics

Metal pollution is a growing concern that affects the health of humans and animals globally. Copper is an essential insect micronutrient required for respiration, pigmentation and oxidative stress protection but can also act as a potentially toxic trace element. While several studies have focused on the negative fitness effects of copper on the aquatic larvae of mosquitoes, the effects of larval copper exposure on adult mosquito fitness (i.e., survival and fecundity) and their ability to transmit parasites (i.e., vector competence) remains unclear. Here, using a well-studied model vector-parasite system, the mosquito Aedes aegypti and parasite Dirofilaria immitis, we show that sublethal copper exposure in larval mosquitoes alters adult female fecundity and vector competence. Specifically, mosquitoes exposed to copper had a hormetic fecundity response and mosquitoes exposed to 600 μg/L of copper had significantly fewer infective parasite larvae than control mosquitoes not exposed to copper. Therefore, exposure of mosquito larvae to copper levels far below EPA-mandated safe drinking water limits (1300 μg/L) can impact vector-borne disease dynamics not only by reducing mosquito abundance (through increased larval mortality), but also by reducing parasite transmission risk. Our results also demonstrated that larval copper is retained through metamorphosis to adulthood in mosquitoes, indicating that these insects could transfer copper from aquatic to terrestrial foodwebs, especially in urban areas where they are abundant. To our knowledge this is the first study to directly link metal exposure with vector competence (i.e., ability to transmit parasites) in any vector-parasite system. Additionally, it also demonstrates unequivocally that mosquitoes can transfer contaminants from aquatic to terrestrial ecosystems. These results have broad implications for public health because they directly linking contaminants and vector-borne disease dynamics, as well as linking mosquitoes and contaminant dynamics.

54 ENVIRONMENTAL SCIENCES↗

Neural-based time series forecasting of loss of coolant accidents in nuclear power plants

During the last few years, deep learning in neural networks has demonstrated impressive successes in the areas of computer vision, speech and image recognition, text generation, and many others. However, sensitive engineering areas such as nuclear engineering benefited less from these efficient techniques. In this work, deep learning expert systems are utilized to model and predict time series progression of a design-basis nuclear accident, featuring a loss of coolant accident. Two major findings are accomplished in this work. First, the ability to train expert systems with high accuracy, which could help nuclear power plant operators to figure out plant responses during the accident. Second, building fast, efficient, and accurate deep models to simulate nuclear phenomena, which could be valuable to nuclear computational science. In this work, large amount of time series data is obtained from simulation tools by simulating different conditions of the base-case/nominal accident scenario. Four critical outputs/responses are monitored during the accident (e.g. temperature, pressure, break flow rate, water level). Two approaches are adopted in this work. The first approach is to use feedforward deep neural networks (DNN) to fit all time steps and outputs in a single model. The second approach is to use long short-term memory (LSTM) to fit all time steps together for each reactor response separately. Both DNN and LSTM demonstrate very good performance in predicting the test and base-case scenarios, with accuracy as low as 92% and as high as 99%, where these test scenarios are unknown to the expert systems and are not included in the model training. In addition, both approaches demonstrate a significant reduction in computational costs, as the deep expert system is able to accurately predict the accident 100,000 times faster than the original simulation tool. Given sufficient data, the methodology adopted in this study demonstrates that DNN/LSTM expert systems can be used as a decision support system to model advanced time series phenomena within nuclear power plants with high accuracy and negligible computational costs.

42 ENGINEERING↗