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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Toward Standardized Microscale Tensile Testing for Two‐Photon Polymerization‐Fabricated Materials in Liquid

Two-photon polymerization (TPP) enables the fabrication of intricate 3D microstructures with submicron precision, offering significant potential in biomedical applications like tissue engineering. In such applications, to print materials and structures with defined mechanics, it is crucial to understand how TPP printing parameters impact the material properties in a physiologically relevant liquid environment. Herein, an experimental approach utilizing microscale tensile testing (μTT) for the systematic measurement of TPP-fabricated microfibers submerged in liquid as a function of printing parameters is introduced. Using a diurethane dimethacrylate-based resin, the influence of printing parameters on microfiber geometry is first explored, demonstrating cross-sectional areas ranging from 1 to 36 μm 2 . Tensile testing reveals Young's moduli between 0.5 and 1.5 GPa and yield strengths from 10 to 60 MPa. The experimental data show an excellent fit with the Ogden hyperelastic polymer model, which enables a detailed analysis of how variations in writing speed, laser power, and printing path influence the mechanical properties of TPP microfibers. The μTT method is also showcased for evaluating multiple commercial resins and for performing cyclic loading experiments. Collectively, this study builds a foundation toward a standardized microscale tensile testing framework to characterize the mechanical properties of TPP printed structures.

mechanical characterization↗
DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗
DOE OSTI2022

Neutrino masses and magnetic moments of electron and muon in the Zee Model

We explore parameter space in the Zee Model to resolve the long-standing tension of the electron and muon anomalous magnetic moment (AMM). The model comprises a second Higgs doublet and a charged singlet at electroweak scale and generates Majorana neutrino masses at one-loop level; the neutral partner of the SU(2) L doublet contributes to the AMM of electron and muon via one loop and two-loop corrections. We propose two minimal flavor structures that can explain these anomalies while fitting the neutrino oscillation data. We find that the neutral Higgs resides in the mass range of roughly 10–300 GeV or 1–30 GeV, depending on the flavor structures. The model is consistent with constraints from colliders, electroweak precision data, and lepton flavor violation. To be comprehensive, we examine the constraints from the electric dipole moment (EDM) and find a region of parameter space that gives a sizable contribution to muon EDM while simultaneously giving corrections to muon AMM. In addition to the light scalar, the two charged scalars with masses as low as 100 GeV can induce nonstandard neutrino interactions ε ee as large as 8%, potentially hinting at new physics. We also investigate the projected capability of future lepton colliders to probe the currently allowed parameter space consistent with both electron and muon AMMs via direct searches in the ℓ + ℓ – → ℓ + ℓ – (H→ℓ + ℓ – )channel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Naturally light Dirac and pseudo-Dirac neutrinos from left-right symmetry

We develop a class of left-right symmetric theories based on the gauge group SU(3) c x SU(2) L x SU(2) R x U(1) with a generalized seesaw mechanism for generating the charged fermion masses. Neutrinos are naturally Dirac particles in this setup with their small masses arising from two-loop quantum corrections. We evaluate these two-loop diagrams exactly and analyze the flavor structure of the lepton sector. We find excellent fits to neutrino oscillation data, independent of the right-handed gauge symmetry breaking scale. We also explore the possibility that neutrinos are pseudo-Dirac particles in this framework, with the tiny mass splittings between active and sterile neutrinos arising from Planck-induced corrections and find possible realizations. These models can be tested in the near future with precision cosmological measurements of ΔN eff in CMB which is predicted to be ≃ 0.14. This class of models allows for a solution to the strong CP problem via parity symmetry without the need for an axion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A nonparametric approach to modeling dependence among production parameters for energy cost estimation

Dependence existing among production parameters can significantly affect energy costs. For material removal processes, processing times and material amounts to be cut have positive dependence, and the dependence can affect energy costs by varying a peak electricity load. Existing methods to model the dependence of production parameters are based on parametric approaches that cannot accurately represent the dependence due to the rigid nature of parametric models. To address this issue, this study proposes a method to quantify dependence among manufacturing parameters through the application of the empirical copula to the dependence between the milling processing time and the amount of volume to be cut. The proposed method is illustrated by a case study of a manufacturing facility consisting of milling machines; a total of 27 scenarios are simulated for energy cost estimations. The case study clearly shows that the proposed method can capture the dependence between the milling parameters more accurately than a conventional dependence measure from parametric correlation models. Furthermore, the findings from this study would be useful to estimate energy costs in that the proposed method provides a better fit to real production data than conventional parametric approaches can.

42 ENGINEERING↗
DOE OSTI2024

Micro–macro finite element modeling method for rub response in abradable coating materials

Gas turbine engines experience “rub” when the rotating blades come in contact with a static abradable coating. This results in extreme strain rates and dynamics inside a high-temperature/high-pressure environment. Current rub models are phenomenological and do not reflect the underlying microstructures, thus limiting their prediction accuracy. In this work, a microstructure-informed, reduced order modeling framework is introduced for simulating abradable coating “rub" behavior. This framework comprises a microscale model constructed based on digitized abradable microstructure and explicitly simulates the mechanical behavior of each constituent phases and their interactions. After calibration and validation with experiment data, the calibrated microscale model is used to generate data across a vast range of applied strain rates and temperature with various load paths. Then, the virtually generated data are used to fit the macroscopic-reduced order model, which enables fast component scale rub simulation without compromising the integrity of the complex material behavior. In conclusion, the proposed effort will address the technical challenge of predicting abradable material behavior during rub through the application of multiscale modeling from microstructure to engines behavior, effectively reducing the development costs and time of new abradable material for better “rub” properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Residual Stress in Cold Spray SS304L Measured Via Neutron Diffraction and Comparison of Analytical Models to Predict the Residual Stress

Here, this study employs neutron diffraction to investigate the relationship between residual stress and coating thickness in cold sprayed 304L austenitic stainless steel. Results show that shot peening predominantly impacts the residual stress profile, leading to substantial in-plane compressive force. The impact of laser heating, a widely used method to alter cold spray's microstructural properties, on the coating's residual stress is also analyzed. The findings indicate that the maximum compressive residual stress in the in-plane component is mainly independent of coating thickness, which suggests that the material properties determine the maximum residual stress. The cold sprayed deposits possessed compressive, nearly biaxial strain and stresses. After laser heating, these stresses were replaced by tensile residual stresses. Two analytical models, the Tsui and Clyne and the Boruah models, for predicting residual stresses are also evaluated, and both models provide reasonable fits to the experimental data. At this point, the deviations between the experimental results and the models are principally caused by the inability of the current models to address plastic deformation and relaxation, and the residual stresses generated by thermal gradients.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrometallurgical Recycling of Black Mass of Spent Lithium-Ion Batteries Using Methanesulfonic Acid: Leaching, Kinetic Studies, and Potential for Total Recovery of Valuable Components

Methanesulfonic acid (MSA) exhibits several advantageous properties rendering it a promising candidate for circular hydrometallurgical processes. These properties include a high acidity (pKa = - 1.9) comparable to that of classical mineral acids as well as biodegradability, high stability, and high solubility of metal-MSA complexes in aqueous solutions. In this study, MSA was employed as a lixiviant for the leaching of metals (lithium, nickel, cobalt, and manganese) from the black mass of spent lithium-ion batteries (LIBs). The effect of various parameters, including MSA concentration, H 2 O 2 concentration, temperature, and pulp density, was systematically investigated. Under the optimized conditions (1.5 M MSA, 0.2 M H 2 O 2 , 60 °C, and 50 g/L pulp density), quantitative leaching of lithium was achieved within 30 min, while for nickel and cobalt it was after 2 h, and 4 h for manganese leaching. The leaching kinetics of Li, Ni, Co, and Mn were studies using the shrinking particle models (SPM) and the Avrami model. The results indicated that the Avrami model provided the best fit to the kinetic data, with apparent activation energies of 46.81 kJ/mol for Li, 58.61 kJ/mol for Ni, 59.69 kJ/mol for Co, and 58.86 kJ/mol for Mn, within the temperature range of 25-70 °C (except for Li, which was analyzed in the range of 25-60 °C), consistent with chemical reaction control. Subsequently, residual contaminants in the leaching residue were eliminated through pyrolysis. The quantitative leaching of metals in MSA solution (a green lixiviant), combined with the pyrolytic treatment of leaching residues, represents a circular strategy for the total recovery of valuable components from the black mass of spent LIBs.

25 ENERGY STORAGE↗
DOE OSTI2020

Grain boundary segregation beyond the dilute limit: Separating the two contributions of site spectrality and solute interactions

Solute segregation at grain boundaries (GBs) is known to have a profound impact on material properties, and as such is becoming routinely used as an element in alloy design. Beyond the dilute limit, the extent of solute GB segregation is known to be concentration dependent. Using hybrid Monte Carlo/Molecular Statics simulations of Mg segregation in Al, in this work we decouple the two contributions to this composition dependence: (i) spectrality of atomic environments at the boundary and (ii) solute-solute interactions. Although only contribution (ii) is typically considered in the literature, we argue that both contributions are equally important to understand concentration dependence and correctly quantify GB solute segregation in a binary alloy. Finally, a thermodynamic segregation isotherm is outlined that accounts for both the spectrality of grain boundary sites and solute-solute interactions. Unlike classical isotherms like those of McLean or Fowler-Guggenheim, which can be successfully fitted to GB segregation data only over a limited range of composition and temperature, our proposed model is shown to be accurate across the composition-temperature space.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Results of an interlaboratory study on the working curve in vat photopolymerization

The working curve informs resin properties and print parameters for stereolithography, digital light processing, and other photopolymer additive manufacturing (PAM) technologies. First demonstrated in 1992, the working curve measurement of cure depth vs radiant exposure of light is now a foundational measurement in the field of PAM. Despite its widespread use in industry and academia, there is no formal method or procedure for performing the working curve measurement, raising questions about the utility of reported working curve parameters. Here, an interlaboratory study (ILS) is described in which 24 individual laboratories performed a working curve measurement on an aliquot from a single batch of PAM resin. The ILS reveals that there is enormous scatter in the working curve data and the key fit parameters derived from it. Further, the measured depth of light penetration D p varied by as much as 7x between participants, while the critical radiant exposure for gelation E c varied by as much as 70x. This significant scatter is attributed to a lack of common procedure, variation in light engines, epistemic uncertainties from the Jacobs equation, and the use of measurement tools with insufficient precision. The ILS findings highlight an urgent need for procedural standardization and better hardware characterization in this rapidly growing field.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Uranium hydride corrosion. I New insights into the Condon-Kirkpatrick model of uranium hydride formation rate

An analytic solution to the uranium hydride formation rate equation of Condon and Kirkpatrick (CK) has been developed, which is directly solved via algebraic operations and requires only gaseous hydrogen (H 2 ) pressure and temperature as inputs. The solution does not depend on the specific functional forms of the hydriding and dehydriding rate constants, nor of key physical parameters. Furthermore, a modified form of the CK model has also been developed, where replacement of the hydrogen solubility function results in an improved fit to published experimental data as H 2 pressure decreases toward the corrosion threshold value.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Colocalized Raman spectroscopy – scanning electrochemical microscopy investigation of redox flow battery dialkoxybenzene redoxmer degradation pathways

Non-aqueous redox flow batteries offer high voltages for grid-level energy storage technologies. However, decomposition of the redoxmers - redox-active molecules that make up the anolyte and catholyte in the negative and positive cell compartments - is an important challenge to overcome for long-term storage. Here, we present a spectroelectrochemical study of the catholyte candidate 2,3-dimethyl-1,4-dialkoxybenzene (C7) and its decomposition mechanisms in the presence of a model nucleophilic base, pyridine. We utilize colocalized Raman microscopy and scanning electrochemical microscopy (Raman-SECM) to quantify the chemical rates of decom-position and qualitatively identify the reaction intermediates and products. A detailed study on how the Raman-SECM parameters (electrode distance, substrate electrode, and laser focal height) influence the Raman signal of the charged catholyte C7 ∙+ is presented. Using optimized conditions, we monitored the deprotonation of C7 ∙+ via Raman spectroscopy and the subsequent hydrogen abstraction from solvent molecules to regenerate C7 via electrochemistry. Finite element modeling was used to fit electrochemical and spectroscopic data, quantifying the deprotonation rates as k dep = 2000 and 700 L mol -1 s -1 and abstraction rates as k abs = 0.5 and 0.2 s - 1 in propylene carbonate and acetonitrile solvents, respectively. Our results show the value of spatiotemporal reso-lution in evaluating the chemical and electrochemical behavior of materials for redox flow batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Isotherm modeling and techno-economic analysis of a TSA moving bed process using a tetraamine-appended MOF for NGCC applications

Tetraamine-appended metal-organic frameworks (MOF) are a new family of amine-functionalized MOF materials that show potential for CO 2 capture from flue gas conditions relevant to natural gas combined cycle (NGCC) power plants. This work presents isotherm modeling of the tetraamine-appended MOF Mg 2 (dobpdc)(3-4-3)(dobpdc 4- =4,4'-dioxidobiphenyl-3,3'dicarboxylate;3-4-3=N,N'bis(3-aminopropyl)-1,4-diaminobutane), process modeling, scale up, and a techno-economic optimization of a moving bed Temperature Swing Adsorption (TSA) process for carbon capture using this sorbent. This MOF exhibits a unique two-step CO 2 adsorption profile in three different pressure ranges. Thus, arctangent-based logistic functions and, quadratic and Langmuir models were employed to represent such isotherm behavior. The results of the isotherm model show good fitting vs the experimental data with an RMSE of 0.41. Here, to model the carbon capture process, the isotherm model was embedded into a moving bed contactor model, and this was used to simulate a TSA CO 2 capture cycle and evaluate cost-optimal designs considering flue gas from a -650 MWe NGCC power plant. The capital cost model consists of CAPEX correlations for reactors, compressors, ducting, etc., while the operating costs include steam, water, chemicals, and electricity (following NETL's quality guidelines for energy systems studies). A techno-economic optimization of the capture system was performed by using NETL's Framework for the Optimization and Quantification of Uncertainty and Surrogates tool (FOQUS). Results suggest that a moving bed carbon capture system with tetraamine-appended MOF can be competitive compared to conventional MEA solvent-based capture processes for NGCC plants, when a heat recovery efficiency of at least 40% is achieved, and MOF materials can be produced at a cost below -$9/kg.

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

Using Thermodynamic Availability to Predict the Transitional Film Reynolds Number between the Droplet and Jet Modes in Falling Liquid between Horizontal Tubes

The objective of this paper is to propose and validate a model for predicting the transitional film Reynolds number between the droplet and jet modes of falling film flows. The proposed model relies on thermodynamic availability with no fitting to mode transition data. The model directly predicts the numerical values of the transitional film Reynolds numbers at different operating conditions. The results of the model are validated against 43 experimental data points from the literature, which cover 10 working fluids and a wide range of tube geometries. The accuracy of the proposed model in predicting the validation data is -on average- more than 20 % better than that of the widely used Re = a (Ga) b and Re = a (Ga) b (S/ξ) c models. The model also achieves low computation time (less than 30 seconds for the 43 validation data points). The model presented in this paper has the potential of integration with heat transfer, mass transfer, and pressure drop models for optimizing the design and operation of falling film heat and mass exchangers.

42 ENGINEERING↗
DOE OSTI2023

Epitaxial growth and characterization of magnesium gallate (MgGa 2 O 4 ) thin films by pulsed laser deposition

To advance high power device application of MgGa 2 O 4 , it is imperative to achieve crystalline epitaxial thin films and to observe its materials’ characteristics within a wide range of growth conditions. However, there is no published work on the growth of this promising ultrawide bandgap material till this day. In this study, we report the crystalline growth of MgGa 2 O 4 thin films on c-plane (006) sapphire substrates using the pulsed laser deposition (PLD) technique in a broad range of temperature and oxygen pressure. The temperature range for this crystalline growth was within the temperature range of 300°C to 700°C and a pressure range of 1x10 -1 to 1x10 -3 Torr. Within this parameter range, the structural, chemical, optical, and surface property of MgGa 2 O 4 thin films have been investigated. Specifically, the acquired XRD patterns confirmed the film growth along the [111] preferential crystal orientation in the lattice. Here, the rocking curve measurement of the prominent (222) plane showed an increasing trend of the crystallinity with growth temperature and pressure. Furthermore, the XRD phi (φ) scan demonstrated the six-fold rotational symmetry of the MgGa 2 O 4 films and the epitaxial relationship of 30° between the film and the sapphire substrate. The XPS spectra confirmed the presence of +2 and +3 oxidation state for Mg and Ga respectively, in the films. The direct bandgap of MgGa 2 O 4 films was obtained ~5.27±0.03 eV by analyzing the UV-Vis absorbance spectra using the Tauc equation. The SEM images exhibited a granular surface morphology of the MgGa 2 O 4 films with an increasing trend in the grain size from low to high temperature and pressure. The refractive index and thickness of these films were in the range of ~1.90±0.02 and ~70±2.0 nm, respectively determined by fitting the spectroscopic ellipsometry data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The low-temperature heat capacity and thermodynamic properties of greigite (Fe 3 S 4 )

Heat capacity measurements provide important insights into the energetic, thermodynamic, and magnetic properties of materials. Herein we report the heat capacity of greigite (Fe 3 S 4 ) from 1.8 to 300 K. Greigite is a magnetic spinel mineral and through a ferromagnetic magnon term, C fsw = B fsw T 3/2 , ferrimagnetic ordering is observed in the low-temperature heat capacity. Using a set of theoretical fits of the experimental data, we calculate the thermodynamic functions, including the standard entropy ($Δ_{0}^{T}$S m °). Greigite is important in iron sulfide formation and reaction pathways in environmental, ore-forming, and technological settings and previous work has measured enthalpies $ΔH_{r}°$ of formation and decomposition to neighboring phases. In this work, the stability of greigite relative to the elements is demonstrated with a negative Gibbs energy ($ΔG_{r}°$) of formation and the stability relative to decomposition products of pyrrhotite (FeS 1.092 ) and pyrite (FeS 2 ) is demonstrated with a positive Gibbs energy ($ΔG_{r}°$) of decomposition. Values of the standard thermodynamic functions C p,m °, $Δ_{0}^{T}$S m °, $Δ_{0}^{T}$H m °, and Φ m ° are tabulated.

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

Current interrupt method for calculating the electrochemical impedance in a solid oxide electrolysis stack

Here, in this work the time domain response of Solid Oxide Electrolysis Cells (SOEC) to a current interruption was transformed into the frequency domain using a carrier function Laplace transform, which is fit to the experimental data using a MATLAB Complex Nonlinear Least Squares (CNLS) solver. The hardware implementation, consisting principally of a high-speed switch and a fast-logging Analog to Digital Converter (ADC), was assembled and tested using a calibration module to assess the accuracy, repeatability, and speed of acquisition of the prototype device as compared against a calibrated commercial impedance spectrometer. Additionally, the current interrupt device and commercial FRA were used to acquire the impedance spectra of a four cell SOEC stack with a large, 300 cm 2 , active cell area.

SOEC↗