Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized Linear Model”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Role of Sink Density in Nonequilibrium Chemical Redistribution in Alloys

Nonequilibrium chemical redistribution in open systems submitted to external forces, such as particle irradiation, leads to changes in the structural properties of the material, potentially driving the system to failure. Such redistribution is controlled by the complex interplay between the production of point defects, atomic transport rates, and the sink character of the microstructure. In this work, we analyze this interplay by means of a kinetic Monte Carlo (KMC) framework with an underlying atomistic model for the Fe-Cr model alloy to study the effect of ideal defect sinks on Cr concentration profiles, with a particular focus on the role of interface density. We observe that the amount of segregation decreases linearly with decreasing interface spacing. Within the framework of the thermodynamics of irreversible processes, a general analytical model is derived and assessed against the KMC simulations to elucidate the structure-property relationship of this system. Interestingly, in the kinetic regime where elimination of point defects at sinks is dominant over bulk recombination, the solute segregation does not directly depend on the dose rate but only on the density of sinks. Furthermore, this model provides new insight into the design of microstructures that mitigate chemical redistribution and improve radiation tolerance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2018

Influence of the starting temperature of calorimetric measurements on the accuracy of determined magnetocaloric effect

Availability of a restricted heat capacity data range has a clear influence on the accuracy of calculated magnetocaloric effect, as confirmed by both numerical simulations and experimental measurements. Simulations using the Bean-Rodbell model show that, in general, the approximated magnetocaloric effect curves calculated using a linear extrapolation of the data starting from a selected temperature point down to zero kelvin deviate in a non-monotonic way from those correctly calculated by fully integrating the data from near zero temperatures. However, we discovered that a particular temperature range exists where the approximated magnetocaloric calculation provides the same result as the fully integrated one. These specific truncated intervals exist for both first and second order phase transitions and are the same for the adiabatic temperature change and magnetic entropy change curves. Here, the effect of this truncated integration in real samples was confirmed using heat capacity data of Gd metal and Gd 5 Si 2 Ge 2 compound measured from near zero temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

GaAs Thermophotovoltaic Patterned Dielectric Back Contact Devices with Improved Sub-Bandgap Reflectance

We demonstrate GaAs thermophotovoltaic (TPV) devices with a patterned dielectric back contact (PDBC) architecture, featuring a dielectric spacer between the semiconductor and back metal contact over most of the back surface for high reflectance, and metal point contacts over a smaller area for electrical conduction. In the TPV application, high sub-bandgap reflectance is needed to reflect unused sub-bandgap photons to the thermal emitter to minimize energy losses in this portion of the thermal spectrum. We explore different PDBC fabrication processes with SU-8 and SiO2 dielectric spacer layers to maximize sub-bandgap reflectance while minimizing series resistance to increase TPV conversion efficiency. We successfully demonstrate GaAs SU-8 PDBC TPV devices with 2200 degrees C blackbody-weighted sub-bandgap reflectance of 94.9% and 96.5% with and without a front metal grid, respectively. This is 0.7% and 2.3% (absolute) higher than the mean sub-bandgap reflectance of 94.2% for GaAs baseline TPV devices with 100% Au back contact with front metal grid. Lower sub-bandgap reflectance in TPV devices with front grids indicates the front grid induces light scattering leading to additional parasitic absorption in the TPV device. We also show that for higher contact coverage fractions, the PDBC reflectance cannot in general be treated by a linear interpolation using simple 1D transfer matrix method modeling and should be treated instead as a diffraction grating by solving Maxwell's equations in 3D.

energy storage↗
DOE OSTI2014

Quantum, classical, and hybrid QM/MM calculations in solution: General implementation of the ddCOSMO linear scaling strategy

We present the general theory and implementation of the Conductor-like Screening Model according to the recently developed ddCOSMO paradigm. The various quantities needed to apply ddCOSMO at different levels of theory, including quantum mechanical descriptions, are discussed in detail, with a particular focus on how to compute the integrals needed to evaluate the ddCOSMO solvation energy and its derivatives. The overall computational cost of a ddCOSMO computation is then analyzed and decomposed in the various steps: the different relative weights of such contributions are then discussed for both ddCOSMO and the fastest available alternative discretization to the COSMO equations. Finally, the scaling of the cost of the various steps with respect to the size of the solute is analyzed and discussed, showing how ddCOSMO opens significantly new possibilities when cheap or hybrid molecular mechanics/quantum mechanics methods are used to describe the solute.

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

Structural fluctuation of proteins induced by thermodynamic perturbation

A theory to describe structural fluctuations of protein induced by thermodynamic perturbations, pressure, temperature, and denaturant, is proposed. The theory is formulated based on the three methods in the statistical mechanics: the generalized Langevin theory, the linear response theory, and the three dimensional interaction site model (3D-RISM) theory. The theory clarifies how the change in thermodynamic conditions, or a macroscopic perturbation, induces the conformational fluctuation, which is a microscopic property. The theoretical results are applied, on the conceptual basis, to explain the experimental finding by Akasaka et al., concerning the NMR experiment which states that the conformational change induced by pressure corresponds to structural fluctuations occurring in the ambient condition. A method to evaluate the structural fluctuation induced by pressure is also suggested by means of the 3D-RISM and the site-site Kirkwood-Buff theories.

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

Z{sup '} physics with early LHC data

We discuss the prospects for setting limits on or discovering spin-1 Z{sup '} bosons using early LHC data at 7 TeV. Our results are based on the narrow width approximation in which the leptonic Drell-Yan Z{sup '} boson production cross section only depends on the Z{sup '} boson mass together with two parameters c{sub u} and c{sub d}. We carefully discuss the experimental cuts that should be applied and tabulate the theoretical next-to-next-to-leading order corrections which must be included. Using these results the approach then provides a safe, convenient, and unbiased way of comparing experiment to theoretical models which avoids any built-in model-dependent assumptions. We apply the method to three classes of perturbative Z{sup '} boson benchmark models: E{sub 6} models, left-right symmetric models, and sequential standard models. We generalize each class of model in terms of mixing angles which continuously parametrize linear combinations of pairs of generators and lead to distinctive orbits in the c{sub u}-c{sub d} plane. We also apply this method to the strongly coupled four-site benchmark model in which two Z{sup '} bosons are predicted. By comparing the experimental limits or discovery bands to the theoretical predictions on the c{sub u}-c{sub d} plane, we show that the LHC at 7 TeV with integrated luminosity of 500 pb{sup -1} will greatly improve on current Tevatron mass limits for the benchmark models. If a Z{sup '} is discovered our results show that measurement of the mass and cross section will provide a powerful discriminator between the benchmark models using this approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Data-Driver Interpretation of Ultradeep Azimuthal Propagation

We report the interpretation of geological features around the borehole stresses on the search of true earth model parameters. Different from the conventional logging-while-drilling (LWD) technologies, such as wireline induction and propagation resistivity tools, the azimuthal resistivity measurement acquires the ability of directional sensitivity. Within this decade, this type of tools has been greatly enhanced with over hundred feet depth-of-detection (DoD). To resolve the interpretation of geological features within 100 feet from the borehole, the multiple-layer structure must be considered. However, the traditional methods can never handle the highly nonlinear inverse problems especially facing on the complicated earth structure. Meanwhile, no justifiable uncertainty of inverse parameters can be provided. The statistical methods have been adopted to search for the global statistical distribution and uncertainties of earth model parameters. However, the most attempts are still stuck in non-uniqueness of the solutions since the structure of earth model has been predefined, where the number of layers is assumed, which introduces another ambiguity into the model interpretation. This paper develops a general trans-dimensional Bayesian methodology to solve non-linear inverse problem of earth layer resistivity and depth to the layer boundaries, from the interpretation of ultra-deep azimuthal resistivity measurements. The presented work is a data-driven approach to search out the multi-dimensional parameter spaces. Finally, the simulations tell that the algorithm can efficiently solve the earth model interpretation within a hundred feet from borehole, and examine the uncertainty of earth model parameters.

58 GEOSCIENCES↗
DOE OSTI2021

Random circuit block-encoded matrix and a proposal of quantum LINPACK benchmark

The LINPACK benchmark reports the performance of a computer for solving a system of linear equations with dense random matrices. Although this task was not designed with a real application directly in mind, the LINPACK benchmark has been used to define the list of TOP500 supercomputers since the debut of the list in 1993. We propose that a similar benchmark, called the quantum LINPACK benchmark, could be used to measure the whole machine performance of quantum computers. The success of the quantum LINPACK benchmark should be viewed as the minimal requirement for a quantum computer to perform a useful task of solving linear algebra problems, such as linear systems of equations. We propose an input model called the Random Circuit Block-Encoded Matrix (RACBEM), which is a proper generalization of a dense random matrix in the quantum setting. The RACBEM model is efficient to be implemented on a quantum computer and can be designed to optimally adapt to any given quantum architecture, with relying on a black-box quantum compiler. Besides solving linear systems, the RACBEM model can be used to perform a variety of linear algebra tasks relevant to many physical applications, such as computing spectral measures, time series generated by a Hamiltonian simulation, and thermal averages of the energy. We implement these linear algebra operations on IBM Q quantum devices as well as quantum virtual machines, and demonstrate their performance in solving scientific computing problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Application of the OPTEX method for computing reflector parameters

We are investigating the OPTEX reflector model for obtaining few-group reflector parameters consistent with a reference power distribution in the core. In our study, the reference power distribution is obtained using a 142,872-region calculation defined over a 2D eighth-of-core pressurized water reactor and performed with the method of characteristics. The OPTEX method is based on generalized perturbation theory and uses an optimization algorithm known as parametric linear complementarity pivoting. The proposed model leads to few-group diffusion coefficients or P1-weighted macroscopic total cross sections that can be used to represent the reflector in full-core calculations. These few-group parameters can be spatially heterogeneous in order to correctly represent steel baffles present in modern pressurized water reactors. The optimal reflector parameters are compared to those obtained with a flux-volume weighting of the reflector cross sections recovered from the reference calculation. Important improvements in full-core power distribution are observed when the optimal parameters are used. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Nuclear Data Adjustment for Nonlinear Applications in the OECD/NEA WPNCS SG14 Benchmark—A Bayesian Inverse UQ-Based Approach for Data Assimilation

The Organisation for Economic Co-operation and Development Working Party on Nuclear Criticality Safety has proposed a benchmark exercise to assess the performance of current nuclear data adjustment techniques applied to nonlinear applications and experiments with low correlation to applications. This work introduces Bayesian inverse uncertainty quantification (IUQ) employing scientific machine learning surrogate models as a method for nuclear data adjustments in this benchmark, and compares IUQ to the more traditional methods of generalized linear least squares (GLLS) and Monte Carlo Bayes (MOCABA). Posterior predictions from IUQ showed agreement with GLLS and MOCABA for linear applications. Here, when comparing GLLS, MOCABA, and IUQ posterior predictions to computed model responses using adjusted parameters, we observe that the GLLS predictions failed to replicate the computed response distributions for nonlinear applications, while MOCABA showed near agreement, and IUQ used the computed model responses directly. We also discuss observations on why experiments with low correlation to applications can be informative to nuclear data adjustments and identify some properties useful in selecting experiments for inclusion in nuclear data adjustment. Performance in this benchmark indicates potential for Bayesian IUQ in nuclear data adjustments.

Bayesian calibration↗
DOE OSTI2025

Using system identification in modeling the yaw response of tail fins for small wind turbines with bearing friction

Here, we describe three main developments of our previous study of the nonlinear yaw dynamics of the tail fins for small wind turbines [Khedr et al., J. Renewable Sustainable Energy 16, 053305 (2024)]. First, the model constants derived from archived computational and experimental studies are adjusted by employing system identification (SI) to maximize the model's agreement with wind tunnel tests. This adjustment was done for high wind speeds, where yaw bearing friction can be ignored. When starting a turbine at low wind speed, however, friction can become important. Our second development is to implement a model for the frictional resistive torque and use SI to maximize its accuracy. These developments used wind tunnel experiments on generic delta, elliptical, and rectangular planforms that were described by Khedr et al. [J. Renewable Sustainable Energy 16, 053305 (2024)]. Since the aerodynamic and friction models employ a large number of constants, we describe ways to constrain the values using linearized solutions of the response equations for small and large yaw angles. Third, we test the generality of the aerodynamic and friction modeling using the complex planform from a commercial small turbine, for which limited theoretical and computational guidance is available in selecting the model constants. The model implemented with SI is shown to provide an accurate description of the yaw response of the complex planform. Guidelines are given for the use of wind tunnel tests to determine the model constants for tail fins of any planform.

17 WIND ENERGY↗
DOE OSTI2024

FARM User Guidance and Instructions

The FARM (Feasible Actuator Range Modifier) module is a component of the RAVEN-based FORCE framework for the analysis of Integrated Energy Systems (IES). FARM aids HERON in the solution of the power dispatch problem by evaluating feasible set-point signals to be issued to the control systems of the different IES unit components. Set-points need to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To enforce all these limits, a two-stage approach is adopted. First, the power dispatcher algorithm in HERON module estimates set-points that meet the constraints on the production variables, e.g., power levels and power ramp rate limits. These constraints are called explicit constraints. Then, if necessary, FARM adjusts these set-points to ensure the respect of the limits on the process variables of interest, given the knowledge of the system dynamics acquired through machine learning algorithms. These constraints are called implicit constraints. From this standpoint, FARM constitutes a bridge between the HERON power dispatcher that adopts a simplified description of the IES unit (low-resolution physics) and the HYBRID high-fidelity models (high-resolution physics). In this report, the general workflow and the software requirements of FARM module are summarized, and the detailed instructions for installing FARM software, running built-in example cases, deriving Linear Parameter-Varying (LPV) state-space models, and using FARM for user-defined power dispatch problems are provided.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Linear Time-Invariant Models of a Large Cumulus Ensemble

Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Next-Level Energy Management in Manufacturing: Facility-Level Energy Digital Twin Framework Based on Machine Learning and Automated Data Collection

This research introduces an energy prediction framework at the facility level supported by automated data collection and machine learning models. It investigates whether reducing the prediction time scale allows for applying more complex machine learning techniques and if those techniques improve the prediction accuracy. The primary advantages of this framework lie in its automation of the energy prediction process and its provision of real-time energy data suitable for use in energy dashboards or digital twins. A sitewide dataset was created by combining 15 min energy and daily production data of five shops—assembly, battery, body (electric), body (gas), and paint—from a globally recognized electric vehicle manufacturer. Various machine learning models were evaluated on daily, weekly, and monthly datasets, including, in increasingly complex order: naïve, simple linear regression, net regularized generalized linear regression, principal component regression, k-nearest neighbor, random forest, and Bayesian regularized neural network. Compared to the current state-of-the-art energy consumption prediction for the industrial facility level, this research investigates more complex models and smaller time intervals for higher accuracy. The findings revealed that the more complex monthly models require a minimum of a year and a half of data to operate, while weekly models demand a year of data to achieve improved accuracy. Daily models can operate with only six months of data but exhibit poor performance due to reduced prediction accuracy of production. Key challenges identified include access to reliable, high-quality energy and production data and the initial demand for human labor.

digital twin↗
DOE OSTI2021

Topographical shifts in fine root lifespan in a mixed, mesic temperate forest

Root lifespan, often is estimated in landscape- and ecosystem-level carbon models using linear approximations. In water manipulation experiments, fine root lifespan can vary with soil water content. Soil water content is generally structured by complex topography, which is largely unaccounted for in landscape- and ecosystem-scale carbon models. Topography governs the range of soil water content experienced by roots which may impact their lifespan. We hypothesized that root lifespan varied nonlinearly across a temperate, mesic, forested catchment due to differences in soil water content associated with topographic position. We expected regions of the landscape that were too wet or too dry would have soils that were not optimal for roots and thus result in shorter root lifespans. Specifically, we hypothesized that root lifespan would be longest in areas that consistently had soil water content in the middle of the soil water content spectrum, while in soils at either very low or very high soil water content, root lifespan would be relatively short. We tested this hypothesis by collecting and analyzing two years of minirhizotron and soil moisture data in plots widely distributed in the Shale Hills catchment of the Susquehanna-Shale Hills Critical Zone Observatory in Pennsylvania. We found that fine root lifespans were longer in traditionally wetter topographic regions, but detected no short term (biweekly) effect of soil moisture on root lifespan. Additionally, depth in soil, soil series, slope face orientation, and season of birth strongly affected root lifespans across the catchment. In contrast, lifespan was unaffected by root diameter or mycorrhizal association. Failure to account for these variables could result in erroneous estimates of fine root lifespan and, consequentially, carbon flux in temperate forested regions.

54 ENVIRONMENTAL SCIENCES↗