Search NASA⌕ Search

SEARCH · Search NASA

Results for “extended least squares”

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 37 records · Page 2

Interface PINNs (I-PINNs): A physics-informed neural networks framework for interface problems

Here, we present a novel physics-informed neural networks (PINNs) framework for modeling interface problems, termed Interface PINNs (I-PINNs). I-PINNs uses different neural networks for any two subdomains separated by a sharp interface such that the neural networks differ only through their activation functions while the other parameters remain identical. The performance of I-PINNs, conventional PINNs, and other existing domain-decomposition PINNs methods such as extended PINNs (XPINNs) and multi-domain PINN (M-PINN) is compared through several one-dimensional, two-dimensional, and three-dimensional benchmark elliptic interface problems. The results demonstrate that I-PINNs provides a root-mean-square-error accuracy, at least two orders of magnitude better than conventional PINNs and XPINNs at approximately one-tenth of the computational cost of conventional PINNs and half the cost of XPINNs. Additionally, while I-PINNs and M-PINN provide comparable accuracies, M-PINN is found to be approximately 50% more expensive.

42 ENGINEERING↗

Operando structural and chemical evolutions of TiS 2 in Na-ion batteries

Titanium disulfide (TiS 2 ) with high electric conductivity, fast rate capability, and good cycling performance is a promising candidate for electrode material in sodium (Na)-ion batteries. Despite the well-studied electrochemical behaviors of TiS 2 in Li-ion batteries, the detailed reaction mechanism of TiS 2 in Na-ion batteries is not yet fully understood due to a more complex multi-phase conversion process. In this work, reactions of TiS 2 in Na-ion batteries are investigated via a multi-modal synchrotron approach: operando X-ray Absorption Spectroscopy (XAS) – including X-ray Absorption Near-Edge Structure (XANES) and Extended X-ray Absorption Fine Structure (EXAFS) – and ex situ X-ray Powder Diffraction (XPD), coupled with computational modeling. Operando XANES spectra indicate that the redox reactions occur in both Ti and S during the electrochemically driven phase transformation. Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS) analysis of XAS suggests that different numbers of components are involved in the lithiation and sodiation of TiS 2 , with the sodiation including at least one intermediate phase in addition to the starting material and the final sodiation product. Ex situ XPD and Rietveld refinement further determined and quantified the unknown phases, showing that three phases, TiS 2 , Na 0.55 TiS 2 , and NaTiS 2 , participate in the sodiation of TiS 2 . Operando EXAFS results show the changes in the Ti–Ti coordination number and interatomic distance. This explains the coulombic efficiency decay due to the incomplete recovery of the coordination number of Ti after cycling. Overall, this work reveals the reaction mechanism occurring in Na–TiS 2 batteries with a greater quantitative understanding of the structural evolution. By combining the multi-modal synchrotron approach and computational work, this study provides a framework for studying a broader range of electrochemically driven phase-transformation systems towards advanced energy storage and conversion applications.

25 ENERGY STORAGE↗

The Area Localized Coupled Model for Analytical Mean Flow Prediction in Arbitrary Wind Farm Geometries

This work introduces the area localized coupled (ALC) model, which extends the applicability of approaches that couple classical wake superposition models and atmospheric boundary layer models to wind farms with arbitrary layouts. Coupling wake and top–down boundary layer models is particularly challenging since the latter requires averaging over planform areas associated with turbine-specific regions of the flow that need to be specified. The ALC model uses Voronoi tessellation to define this local area around each turbine. A top–down description of a developing internal boundary layer is then applied over Voronoi cells upstream of each turbine to estimate the local mean velocity profile. Coupling between the velocity at hub-height based on this localized top–down model and a wake model is achieved by enforcing a minimum least-square-error in mean velocity in each cell. The wake model in the present implementation takes into account variations in wind farm inflow velocity and represents the wake profile behind each turbine as a super-Gaussian function that smoothly transitions between a top-hat shape in the region immediately following the turbine to a Gaussian profile downstream. Detailed comparisons to large-eddy simulation (LES) data from two different wind farms demonstrate the efficacy of the model in accurately predicting both wind farm power output and local turbine hub-height velocity for different wind farm geometries. These validations using data generated from two different LES codes demonstrate the model's versatility with respect to capturing results from different simulation setups and wind farm configurations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Remote Americium Detection Using an Optical Sensor: A D-Optimal Strategy for Efficient PLS-Based Modeling

A fiber-optic visible–near-infrared absorption spectroscopy system in a glove box was demonstrated for remote quantification of Am(III) (0–500 µM) and HNO 3 (0.1–9 M) using partial least squares regression (PLSR) models. The sensor platform, featuring a simple plug-and-play spectrophotometer, can enable noninvasive, real-time monitoring of actinide process solutions. To establish a flexible PLSR model calibration strategy, a D-optimal design developed using Nd(III) in previous studies was successfully extended to an actinide system with Am(III) to effectively minimize sample set size while maintaining robust prediction performance. The results suggest strong spectral similarities between Nd(III) and Am(III) and validate Nd(III) as an effective optical surrogate for trivalent actinide species. This work also supports the generalizability of a D-optimal training set selection approach for two-factor systems. The PLS1 models for Am(III) and HNO 3 outperformed a PLS2 model and maintained reasonable performance in the presence of interfering U(VI). The resulting sensor system and multivariate approach provides a flexible and scalable solution for process monitoring, control, and safety in diverse nuclear applications.

actinide↗

Welding of Haynes 282 to Steels to Enable Modular Rotors for Advanced Ultra Super-Critical Steam Turbines

Steam turbines for an Advanced Ultra Super Critical (AUSC) fossil fired power plant will operate at temperatures well above those of current commercial steam cycles with inlet temperatures more than 760 C , which is beyond the capabilities of alloy steels presently used for steam turbine applications and requires advanced materials. Large components such as steam turbine rotors may be made of nickel based super alloys such as Haynes 282 (H282). However, monolithic forgings of superalloys in the sizes required for large steam turbines rotors can be prohibitively expensive besides many technical challenges. To minimize cost and alleviate the related technical challenges, superalloy use needs to be limited to locations on the steam turbine rotor where strength and temperature requirements cannot be met by conventional steels. This is possible if nickel-based superalloys can be successfully welded to steels and the related technical challenges - machining parts made of dissimilar welded materials, non-destructive examination of such welds for flaw detection; and, material properties of such hybrid components – are sufficiently addressed. In this technology development project, we successfully welded H282 to plates up to ~ 75 mm (~ 3 inches) to a 3.5NiCrMoV steel of similar thickness. Advanced ultrasonic inspection technique called Phased Array Ultrasonic Testing (PAUT) was employed to examine the dissimilar H282-Steel welds, into which flat bottomed side drilled holes (SDH) of various diameters were introduced, to determine the minimum detectable feature sizes; it was shown that with PAUT SDH of dia. down to 0.5 mm could be detected in the base alloys and SDH with dia. down to 2.4 mm could be detected in the weld metal under multiple orientations successfully. An autonomous machining process monitoring system was developed and demonstrated whereby the forces acting on the cutting tool could be actively monitored as the cutting tool transitioned from H282 to Steel across the weld using which the machining parameters can be potentially altered without interruption to extend tool life. This project successfully achieved its objectives of - i. Developing a welding methodology and viable welding geometries, to successfully join H282 to steel 3.5CrMoNiV steel to enable manufacture of modular steam turbine rotors for AUSC applications (conditions of at least 760 °Celsius and 3,100 psia (pounds per square inch absolute pressure) and evaluate the material properties of the welded specimen. ii. Employ the advanced ultrasonic inspection technique to the dissimilar weld metal joint and determine the minimum detectable feature sizes. iii. Develop effective machining techniques to machine such hybrid structures with online tool force monitoring and effect machine state metrics for optimal results.

20 FOSSIL-FUELED POWER PLANTS↗

Optimization of screw design for continuous wet granulation: A case study of metoprolol succinate ER tablets

This study aimed at understanding the effect of screw design on the critical characteristics of granules and tablets of an extended-release (ER) formulation for twin screw granulation process. The screw design parameters assessed included number of kneading elements (KEs) per kneading zone, distance separating kneading zones, staggering angle (SA) of kneading elements and number of sizing elements (SEs). These input variables were varied using a design of experiment (DoE) approach to manufacture granules. Particle size distribution (PSD), flow and bulk properties of the granules, breaking strength and dissolution of tablets manufactured using these granules were characterized. The results of least square fitting showed that KEs, SA, and SEs of the screws significantly (p -values < 0.05) affected the PSD, cohesion, compressibility (CPS), conditioned bulk density (CBD) and permeability of the granules. The KEs and SEs significantly (p -value < 0.05) affected the dissolution, which was attributed to their effects on CPS and CBD of the granules. The distance between kneading zones had no significant effect on granules and tablet characteristics. Finally, these results may be used to further study the interaction of the identified critical screw design parameters with other processing parameters for continuous manufacturing of this ER matrix-based tablet formulation.

60 APPLIED LIFE SCIENCES↗

Macrovoid resolved simulations of transport through HPRO relevant membrane geometries

Modeling the transport properties such as diffusivity and permeability of high pressure reverse osmosis (HPRO) membranes is critical for the selection and manufacture of membranes suitable for operation under high pressure. These properties can be significantly affected by the changes in heterogeneous pore structures due to compaction. Here, the modeling platform presented in this work resolves the two scale porosity in HPRO relevant membranes. A synthetic membrane geometry is constructed based on available experimental visualizations of pore structures. The simulations directly capture the flow channeling that results from a combination of material properties and geometric features of the macrovoids. A parametric study is presented to account for the transition of flow characteristics from a material governed regime to a macrovoid governed regime. Permeability of the membrane is evaluated using the simulation data and compared with an existing model that scales with the square of tortuosity over a range of material properties. The model is found to perform well under a narrow range of tortuosity while deviating from the calculated permeabilities at extreme conditions. The effective permeability of the membrane is found to vary by at least two orders of magnitude between the two flow regimes. It is also observed that tortuosity is a bounded property with its upper limit determined by the macrovoid geometry. Consequently, the tortuosity based correlations fail near a flow regime that is mainly governed by the macrovoids. The modeled permeability can be more than an order of magnitude smaller than the simulation result. A new model based on flux-weighted porosity of a membrane is introduced and its correlation with tortuosity is studied. The model agrees with the simulated data as, in addition to tortuosity, it also accounts for the flux partition within and outside the flow channels. Such correlations enable extending the existing understanding of flow characteristics to enhance predictability of porous media models.

42 ENGINEERING↗

Prediction of Dielectric Constant in Series of Polymers by Quantitative Structure-Property Relationship (QSPR)

This work is devoted to the investigation of dielectric permittivity which is influenced by electronic, ionic, and dipolar polarization mechanisms, contributing to the material’s capacity to store electrical energy. In this study, an extended dataset of 86 polymers was analyzed, and two quantitative structure–property relationship (QSPR) models were developed to predict dielectric permittivity. From an initial set of 1273 descriptors, the most relevant ones were selected using a genetic algorithm, and machine learning models were built using the Gradient Boosting Regressor (GBR). In contrast to Multiple Linear Regression (MLR)- and Partial Least Squares (PLS)-based models, the gradient boosting models excel in handling nonlinear relationships and multicollinearity, iteratively optimizing decision trees to improve accuracy without overfitting. The developed GBR models showed high R2 coefficients of 0.938 and 0.822, for the training and test sets, respectively. An Accumulated Local Effect (ALE) technique was applied to assess the relationship between the selected descriptors—eight for the GB_A model and six for the GB_B model, and their impact on target property. ALE analysis revealed that descriptors such as TDB09m had a strong positive effect on permittivity, while MLOGP2 showed a negative effect. These results highlight the effectiveness of the GBR approach in predicting the dielectric properties of polymers, offering improved accuracy and interpretability.

Ascencio-Medina, Estefania↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

Temporal Ni K-Edge X-ray Absorption Spectroscopy Study Reveals the Kinetics of the Ni Redox Behavior of the Iron-Nickel Oxide Bimetallic OER Catalyst

Operando X-ray absorption spectroscopy (XAS) can be utilized to probe the phase and structural changes of FeNiOx and similar transition metal oxide electrocatalysts during electrocatalytic reactions. However, capturing the temporal changes occurring in the operando chemistry of electrocatalysts has been little studied. Here, in this work, we successfully capture the time-resolved changes at the Ni K-edge for both the X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) regions of an aqueous phase synthesized FeNiO x bimetallic nanoparticle electrocatalyst. During a stepped voltage experiment, the temporal change in the Ni K-edge was observed as the applied voltage was stepped from 0.7 to 0.8 V versus the silver/silver chloride reference electrode, and the change observed is associated with the Ni redox transition from the 2+ to 3+/4+ oxidation state. Individual XAS spectra were obtained in 90 s, and a total of 9 scans post voltage step showed a unique transition from a hydroxide to an oxyhydroxide phase. The shift in absorption edge energy position and the changes in spectral shape for individual scans explain the typically observed broadening of the Ni K-edge XANES spectrum during time-averaged operando XAS, and a kinetic analysis revealed a first-order observed rate constant, |k obs |, of 0.00426 s –1 and a half-life (t 1/2 ) of 163 s. Multivariate curve resolution-alternating least squares analysis followed by linear combination fitting analysis for time-averaged versus time-resolved Ni K-edge changes upon voltage step show distinct differences in estimated contributions from Ni 2+ oxide/hydroxide phases. Detailed EXAFS modeling shows the phase transition from hydroxide to oxyhydroxide on top of the unchanged metallic core.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal neutron imaging detector with wavelength shifting fiber and digital readout

Here, we developed a large area, digital thermal neutron imaging detector. The detector uses a 6 LiF/ZnS(Ag) neutron-sensitive scintillator combined with wavelength shifting fiber technology. The signals from fiber channels are amplified, integrated, and digitized using individual analog-to-digital converters for each channel. The neutron position is determined from the digitized signal using a least-squares gaussian-fitting algorithm. The detector size was 77 × 38 cm 2 (2926 cm 2 ), it provided 1.3 mm resolution, and its resolution can further be improved. The detector was designed for a powder diffraction neutron scattering beamline and provided substantial improvement of d-space resolution compared with existing detectors. This detector may have broader applications due to its large area and high spatial resolution extended over its large area.

6LiF/ZnS(Ag) scintillator↗

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗

Possible Quantum Paramagnetism in Compressed Sr 2 IrO 4

The result of compression on the magnetic ground state of Sr 2 IrO 4 is studied with x-ray resonant techniques in the diamond anvil cell. The weak interlayer exchange coupling between square-planar 2D IrO 2 layers is easily modified upon compression, with a crossover between magnetic structures around 7 GPa mimicking the effect of an applied magnetic field at ambient pressure. Higher pressures drive an order-disorder magnetic phase transition with no magnetic order detected above 17–20 GPa. The persistence of strong exchange interactions between J eff = 1 / 2 magnetic moments within the insulating IrO 2 layers up to at least 35 GPa points to a highly frustrated magnetic state in compressed Sr 2 IrO 4 , opening the door for realization of novel quantum paramagnetic phases driven by extended 5 d orbitals with entangled spin and orbital degrees of freedom.

36 MATERIALS SCIENCE↗

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↗

Quantification of Rare Earth Elements in the Parts Per Million Range: A Novel Approach in the Application of Laser-Induced Breakdown Spectroscopy

This work extends a previous percentage level concentration study of the optical emission spectra for six rare earth elements, europium (Eu), gadolinium (Gd), lanthanum (La), praseodymium (Pr), neodymium (Nd), and samarium (Sm), along with the transition metal, yttrium (Y) using laser-induced breakdown spectroscopy (LIBS). The concentration of these six rare earth elements and yttrium has been attempted for the first time systematically down to parts per million (ppm) concentration levels ranging from 30 to 300 ppm. In this study, the authors have developed multivariate models for each element capable of predicting concentration with acceptable to excellent levels of accuracy. Additionally, partial least squares regression coefficients were used to identify key spectral features able to be used in this lower concentration regime. This study has demonstrated that it is conceivable to quantify the six rare earth elements along with yttrium at low concentrations in the parts per million levels.

59 BASIC BIOLOGICAL SCIENCES↗

Extended real-time voltage instability identification method based on synchronized phasor measurements

This paper presents an extended adaptive approach designed to accurately estimate the Thévenin equivalent parameters using phasor measurements at a given bus for measurements lying in any of four quadrants of the PQ-plane. The improvement is achieved by using a new condition to properly update the estimated parameters after an initial guess. Based on an adaptive philosophy, the proposed approach can correctly account for the intrinsic nonlinearities of a power system, can provide a real-time estimation of Thévenin parameters, and does not require network topology knowledge. The method is validated using the Kundur 2-area system, showing estimation improvements compared to the current adaptive approach and the classical recursive least-squares method. The proposed approach is able to estimate both sides of the system with respect to the measurement bus. In addition, a data-driven voltage stability index is developed. To illustrate the performance of the proposed approach in a larger power system, a voltage stability assessment is carried out on the IEEE 39-bus system, considering the action of overexcitation limiters of generators and nonlinear loads. The proposed approach is suitable for applications that require an accurate Thévenin equivalent estimation in real-time, such as for voltage stability assessment. The new approach provides a reliable tool for the system operators to make proper and timely decisions.

42 ENGINEERING↗

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computationally Robust Line Outage Detection and Identification in Three-Phase Networks

Detection and identification of individual phase outages remains a challenging problem due to insufficient metering in three-phase unbalanced power networks. This problem was tackled for the transmission systems in our previous work as documented in [1]. In this paper, this work is extended to detect phase outages in three-phase unbalanced systems using only the sparse estimation method. In addition, further improvements are introduced to increase the estimation accuracy for the virtual power injections at the terminal buses of disconnected lines, once the disconnected line is identified by sparse estimation methods. Simulation results are provided to experimentally validate the increased accuracy in detecting phase outages while decreasing the computational time by using the proposed approach.

Sparse Estimation, Line Outage, Phase Outage, LASS↗