Search NASA⌕ Search

SEARCH · Search NASA

Results for “box 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 307 records · Page 17

Field Work Proposal ERKJ358: Black-box training for scientific machine learning models (Final Report)

The overarching goal of this project is to develop a scalable black-box training capability for scientific machine learning (SciML) problems that are non-trainable with existing automatic differentiation (AD)-based algorithms. AD assumes that a loss function can be decomposed into a sequence of elementary operations whose derivatives are known. This assumption is violated when the loss function includes a black-box physical model (e.g., a legacy simulator). The current strategy, converting a black-box simulator to an AD-enabled code via differential programming, is inflexible and time-, labor-consuming. Thus, black-box optimization is a main workhorse for training SciML models, e.g., in scientific reinforcement learning, hyper-parameter fine tuning, designing SciML models with adversarial robustness, etc.

97 MATHEMATICS AND COMPUTING↗

Continental-scale river flow in climate models

The hydrologic cycle is a major part of the global climate system. There is an atmospheric flux of water from the ocean surface to the continents. The cycle is closed by return flow in rivers. In this paper a river routing model is developed to use with grid box climate models for the whole earth. The routing model needs an algorithm for the river mass flow and a river direction file, which has been compiled for 4 deg x 5 deg and 2 deg x 2.5 deg resolutions. River basins are defined by the direction files. The river flow leaving each grid box depends on river and lake mass, downstream distance, and an effective flow speed that depends on topography. As input the routing model uses monthly land source runoff from a 5-yr simulation of the NASA/GISS atmospheric climate model (Hansen et al.). The land source runoff from the 4 deg x 5 deg resolution model is quartered onto a 2 deg x 2.5 deg grid, and the effect of grid resolution is examined. Monthly flow at the mouth of the world's major rivers is compared with observations, and a global error function for river flow is used to evaluate the routing model and its sensitivity to physical parameters. Three basinwide parameters are introduced: the river length weighted by source runoff, the turnover rate, and the basinwide speed. Although the values of these parameters depend on the resolution at which the rivers are defined, the values should converge as the grid resolution becomes finer. When the routing scheme described here is coupled with a climate model's source runoff, it provides the basis for closing the hydrologic cycle in coupled atmosphere-ocean models by realistically allowing water to return to the ocean at the correct location and with the proper magnitude and timing.

Miller, James R.↗

POD Model Reconstruction for Gray-Box Fault Detection

Proper orthogonal decomposition (POD) is the mathematical basis of a method of constructing low-order mathematical models for the "gray-box" fault-detection algorithm that is a component of a diagnostic system known as beacon-based exception analysis for multi-missions (BEAM). POD has been successfully applied in reducing computational complexity by generating simple models that can be used for control and simulation for complex systems such as fluid flows. In the present application to BEAM, POD brings the same benefits to automated diagnosis. BEAM is a method of real-time or offline, automated diagnosis of a complex dynamic system.The gray-box approach makes it possible to utilize incomplete or approximate knowledge of the dynamics of the system that one seeks to diagnose. In the gray-box approach, a deterministic model of the system is used to filter a time series of system sensor data to remove the deterministic components of the time series from further examination. What is left after the filtering operation is a time series of residual quantities that represent the unknown (or at least unmodeled) aspects of the behavior of the system. Stochastic modeling techniques are then applied to the residual time series. The procedure for detecting abnormal behavior of the system then becomes one of looking for statistical differences between the residual time series and the predictions of the stochastic model.

Park, Han↗

Comparison of equivalent plate and finite element analysis of a realistic aircraft structural configuration

This paper presents a comparison study to assess how well a realistic aircraft structure can be represented with an equivalent plate model. A wing box and full wing structure of a high-speed civil transport aircraft configuration are analyzed. Static and vibration results from the equivalent plate analysis are compared with results from a finite element model. The results from these comparisons show that, for the simpler wing box models, the static displacements and first vibration frequency agree to within 1.5 percent. The full wing models only agree to within 10 percent, for both models, the higher frequencies vary by as much as 26 percent. These differences are due to inadequate refinement of the equivalent plate model.

Kao, Pi-Jen↗

Physics-informed hybrid modeling methodology for building infiltration

Infiltration is responsible for one-third to one-half of the space conditioning load of a typical residential home, but the modeling of infiltration for building energy modeling is either represented by over-simplified equations or dependent on over-generalized rules of thumb. Here, this paper develops a physics-informed data-driven methodology for modeling infiltration using building-specific empirical measurements. The developed hybrid methodology combines machine-learning categorization and grey-box sub-modeling to improve the accuracy and generalization of commonly used grey-box infiltration models. The developed methodology excels at predicting infiltration by improving the ability to predict infiltration under unseen environmental conditions using machine learning algorithms with physical significance. In a case study conducted using the iUnit, a modular studio apartment experimental test facility located at the National Renewable Energy Laboratory, we use empirical airtightness measurements to fit an infiltration model using the developed methodology. We find that the developed methodology can improve the overall model accuracy by 43% and improve extrapolation by 38%, compared with the model based on the common grey-box infiltration equation. We also notice that the selected features can improve the performance of a pure machine-learning model, indicating that our methodology identifies the features with the most physical significance to infiltration modeling.

97 MATHEMATICS AND COMPUTING↗

Linear model decision trees as surrogates in optimization of engineering applications

Machine learning models are promising as surrogates in optimization when replacing difficult to solve equations or black-box type models. This work demonstrates the viability of linear model decision trees as piecewise-linear surrogates in decision-making problems. Linear model decision trees can be represented exactly in mixed-integer linear programming (MILP) and mixed-integer quadratic constrained programming (MIQCP) formulations. Furthermore, they can represent discontinuous functions, bringing advantages over neural networks in some cases. We present several formulations using transformations from Generalized Disjunctive Programming (GDP) formulations and modifications of MILP formulations for gradient boosted decision trees (GBDT). We then compare the computational performance of these different MILP and MIQCP representations in an optimization problem and illustrate their use on engineering applications. Importantly, we observe faster solution times for optimization problems with linear model decision tree surrogates when compared with GBDT surrogates using the Optimization and Machine Learning Toolkit (OMLT).

42 ENGINEERING↗

Power Electronics Based Self-Monitoring and Diagnosing for Photovoltaic Systems

Faults in photovoltaic (PV) systems can seriously affect the efficiency, energy yield, cost, safety, and reliability of PV plants. Condition monitoring of PV plants is, therefore, a very important approach to estimating the health condition of PV modules and power electronics in the system. However, additional hardware for PV system monitoring adds cost to the system's operation; delayed maintenance service also causes additional energy production loss. The Center for Power Electronics Systems (CPES) at the Virginia Polytechnic Institute and State University and Siemens Cooperate Research developed the online impedance measurement for a PV panel self-monitoring and diagnosing technology using the DC-DC converter connected to the panel. Small-signal impedances of a monocrystalline silicon PV panel were modeled and simulated to reflect fault conditions such as the short-circuit, hot-spot, and junction box faults. Modeling and simulation results were validated firstly using a test setup consisting of a solar simulator, a network analyzer, small-signal injectors, and a monocrystalline PV panel rated at 300 W.

14 SOLAR ENERGY↗

A New Higher-Order Composite Theory for Analysis and Design of High Speed Tilt-Rotor Blades

A higher-order theory is developed to model composite box beams with arbitrary wall thicknesses. The theory, based on a refined displacement field, represents a three-dimensional model which approximates the elasticity solution. Therefore, the cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are automatically included in the formulation. The model accurately captures the transverse shear stresses through the thickness of each wall while satisfying all stress-free boundary conditions. Several numerical results are presented to validate the present theory. The developed theory is then used to model the load carrying member of a tilt-rotor blade which has thick-walled sections. The composite structural analysis is coupled with an aerodynamic analysis to compute the aeroelastic stability of the blade. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt-rotor aircraft. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem and a hybrid approximate analysis is used to reduce the computational effort. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt-rotor blade.

McCarthy, Thomas Robert↗

Multidisciplinary Optimization of Tilt Rotor Blades Using Comprehensive Composite Modeling Technique

An optimization procedure is developed for addressing the design of composite tilt rotor blades. A comprehensive technique, based on a higher-order laminate theory, is developed for the analysis of the thick composite load-carrying sections, modeled as box beams, in the blade. The theory, which is based on a refined displacement field, is a three-dimensional model which approximates the elasticity solution so that the beam cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are included automatically in the formulation. The model can accurately capture the transverse shear stresses through the thickness of each wall while satisfying stress free boundary conditions on the inner and outer surfaces of the beam. The aerodynamic loads on the blade are calculated using the classical blade element momentum theory. Analytical expressions for the lift and drag are obtained based on the blade planform with corrections for the high lift capability of rotor blades. The aerodynamic analysis is coupled with the structural model to formulate the complete coupled equations of motion for aeroelastic analyses. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt rotor aircraft. The objective functions include the figure of merit in hover and the high speed cruise propulsive efficiency. Structural, aerodynamic and aeroelastic stability criteria are imposed as constraints on the problem. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem. The search direction is determined by the Broyden-Fletcher-Goldfarb-Shanno algorithm. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt rotor blade.

Chattopadhyay, Aditi↗

Validation and parameterization of a novel physics-constrained neural dynamics model applied to turbulent fluid flow

We report, in fluid physics, data-driven models to enhance or accelerate time to solution are becoming increasingly popular for many application domains, such as alternatives to turbulence closures, system surrogates, or for new physics discovery. In the context of reduced order models of high-dimensional time-dependent fluid systems, machine learning methods grant the benefit of automated learning from data, but the burden of a model lies on its reduced-order representation of both the fluid state and physical dynamics. In this work, we build a physics-constrained, data-driven reduced order model for Navier–Stokes equations to approximate spatiotemporal fluid dynamics in the canonical case of isotropic turbulence in a triply periodic box. The model design choices mimic numerical and physical constraints by, for example, implicitly enforcing the incompressibility constraint and utilizing continuous neural ordinary differential equations for tracking the evolution of the governing differential equation. We demonstrate this technique on a three-dimensional, moderate Reynolds number turbulent fluid flow. In assessing the statistical quality and characteristics of the machine-learned model through rigorous diagnostic tests, we find that our model is capable of reconstructing the dynamics of the flow over large integral timescales, favoring accuracy at the larger length scales. More significantly, comprehensive diagnostics suggest that physically interpretable model parameters, corresponding to the representations of the fluid state and dynamics, have attributable and quantifiable impact on the quality of the model predictions and computational complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

GIST (Grid Impedance Scan Tool) [SWR-22-73]

NREL has developed GIST (Grid Impedance Scan Tool) to scan the impedance response of inverter-based resources such as land-based and offshore wind power plants, solar power plants, etc., from their electromagnetic transient simulation models to evaluate their impact on grid stability. GIST can scan impedance response of inverter-based resources from vendor-supplied black-box simulation models without requiring internal proprietary details, enabling detailed stability analysis of power systems with high levels of inverter-based generation.

Shah, Shahil↗

Model of Inverter-Based Resources

Blackbox modelling for SC analysis is a possible solution. Accuracy can be acceptable even without having vendor control diagrams. Is an NDA required to share a vendor Blackbox model for SC analysis? Differences in VRT detection and injection logics. Angle rotation is not addressed in control logics of most vendors. Current limitation logic during unbalanced faults is not clear. Model IBR as a current source with shunt to improve convergence

IBRs, Model of IBRs, Short Circuit Analysis, Class↗

The impact of global warming on river runoff

A global atmospheric model is used to calculate the annual river runoff for 33 of the world's major rivers for the present climate and for a doubled CO2 climate. The model has a horizontal resolution of 4 x 5 deg, but the runoff from each model grid box is quartered and added to the appropriate river drainage basin on a 2 x 2.5 deg resolution. The computed runoff depends on the model's precipitation, evapotranspiration, and soil moisture storage. For the doubled CO2 climate, the runoff increased for 25 of the 33 rivers, and in most cases the increases coincide with increased rainfall within the drainage basins. There were runoff increases in all rivers in high northern latitudes, with a maximum increase of 47 percent. At low latitudes there were both increases and decreases ranging from a 96 increase to a 43 percent decrease. The effect of the simplified model assumptions of land-atmosphere interactions on the results is discussed.

Miller, James R.↗

SCALE modeling of foil irradiations at WSU’s TRIGA with sensitivity/uncertainty analysis

We report the Washington State University TRIGA reactor was modeled in the SCALE system’s KENO neutron transport code. The model includes 119 fuel rods, a central boral control rod, 4 control blades, 5 irradiation ports, and 20 graphite reflector blocks within a rectangular aluminum grid box. The model accurately simulated the irradiation of depleted uranium and highly enriched uranium foils and predicted reaction rate ratios for 238 U(n,2n) 237 U / 238 U(n,γ) 239 U and 99 Mo / fission, which compared favorably to empirical values previously measured at Los Alamos. Finally, a sensitivity/uncertainty analysis was performed using the continuous-energy TSUNAMI-3D module of SCALE.

237Np↗

A novel implicit hybrid machine learning model and its application for reinforcement learning

A novel methodology to develop implicit hybrid models is presented. PyTorch is used to integrate physics-based equations with machine learning models. Automatic differentiation of the hybrid model is leveraged to solve the implicit equations. Iterative solving enables gradient based updates to the machine learning model. The novel methodology is compared to an explicit hybrid approach on a continuously stirred tank reactor (CSTR). The novel method results in a lower modelling error. Both hybrid models effectively train with noisy data. To test the implicit hybrid model, it is employed as a reinforcement learning (RL) training model. The RL algorithm trained on the hybrid model outperforms real time optimization of the CSTR and performs nearly as well as RL trained directly on the CSTR and a traditional gradient based approach. Training RL directly on the CSTR requires over 60,000 system interactions compared to 6000 historical data points for hybrid model development.

42 ENGINEERING↗

A numerical analysis of shipboard and coastal zone color scanner time series of new production within Gulf Stream cyclonic eddies in the South Atlantic Bight

Eddy-induced upwelling occurs along the western edge of the Gulf Stream between Cape Canaveral, Florida, and Cape Hatteras, North Carolina, in the South Atlantic Bight (SAB). Coastal zone color scanner images of 1-km resolution spanning the period April 13-21, 1979, were processed to examine these eddy features in relation to concurrent shipboard and current/temperature measurements at moored arrays. A quasi-one-dimensional (z), time dependent biological model, using only nitrate as a nutrient source, has been combined with a three-dimensional physical model in an attempt to replicate the observed phytoplankton field at the northward edge of an eddy. The model is applicable only to the SAB south of the Charleston Bump, at approximately 31.5 deg N, since no feature analogous to the bump exists in the model bathymetry. The modeled chlorophyll, nitrate, and primary production fields of the euphotic zone are very similar to those obtained from the satellite and shipboard data at the leading edges of the observed eddies south of the Charleston Bump. The horizontal and vertical simulated fluxes of nitrate and chlorophyll show that only approximately 10% of the upwelled nitrate is utilized by the phytoplankton of the modeled grid box on the northern edge of the cyclone, while approximately 75% is lost horizontally, with the remainder still in the euphotic zone after the 10-day period of the model. Loss of chlorophyll due to sinking is very small in this strong upwelling region of the cyclone. The model is relatively insensitive to variations in the sinking parameterization and the external nitrate and chlorophyll fields but is very sensitive to a reduction of the maximum potential growth rate to half that measured. Given the success of this model in simulating the new production of the selcted upwelling region, other upwelling regions for which measurements or successful models of physical and biological quantities and rates exist could be modeled similarly.

Pribble, J. Raymond↗

On the nature of the cosmic ray positron spectrum

A calculation was made of the flux of secondary positrons above 100 MeV expected for various propagation models. The models investigated were the leaky box or homogeneous model, a disk halo diffusion model, a dynamical halo model, and the closed galaxy model. In each case the parameters of these models were adjusted for agreement with the observed secondary or primary ratios and Be 10 abundance. The positron flux predicted for these models was compared with the available data. The possibility of a primary positron component was considered.

Protheroe, R. J.↗