Search NASA⌕ Search

SEARCH · Search NASA

Results for “model simplification”

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

Thermo-hydraulic steam pipe models for district heating simulations: Simplifications to balance accuracy and simulation speed

Steam piping networks are essential for optimizing performance in industrial processes and district heating systems. However, dynamic models that balance thermo-hydraulic accuracy with computational efficiency remain limited. In response, this paper presents a new discretized steam pipe model based on the plug flow approach, capturing key thermo-hydraulic behaviors while simplifying steam phase change processes. Implemented in Modelica, the model accurately calculates temperature and pressure distributions along steam pipelines. To improve computational efficiency for district-scale simulations, five model simplifications are introduced: lumped thermo-hydraulic functions, empirical correlations, fluid state approximations, steady-state dynamics and inclusion of flow derivatives. These simplified models achieve 85%-98% accuracy in predicting pressure drop and condensation losses, including dynamic condensate behavior during pipe warm-up—a factor often overlooked in existing models. The models support diverse network configurations, scaling effectively to systems with multiple distribution pipes and connected building loads. Discrete models provide detailed insights but exhibit a cubic increase in simulation time as the network scales by N connected building O(N 2.42 ). In contrast, lumped models simulate 10–28 times faster than discrete, offering quadratic scaling of simulation time O(N 1.73 ). However, they still require 6 times more computation time than a lossless network, highlighting the inherent computational challenges of modeling compressible fluid flow. In conclusion, the steady-state lumped variant, with its near-linear scalability in computational time O(N 1.01 ), emerges as an efficient solution for preliminary design evaluations and extensive parametric studies.

15 GEOTHERMAL ENERGY↗

All-Atom Modeling and Simulation of Biopolymer Interface: Dual Role of Antifouling Polymer Brushes

Antifouling polymer brushes are well-known for their exceptional resistance to unwanted protein adsorption. While experimental studies have extensively characterized protein–polymer brush interactions, computational investigations remain limited, largely due to the challenges in accurately modeling and integrating these complex, multicomponent systems without resorting to oversimplification. To address this challenge, this study presents one of the most comprehensive and realistic model systems to date, comprising the substrate, grafted polymers, and proteins. Our work interprets the interactions between polymer brush and protein based on realistic modeling without simplification. In particular, this study utilizes molecular modeling and simulation of polycationic and polyzwitterionic brushes─poly(dimethylaminoethyl methacrylate) (PDMAEMA), poly(2-(N-oxide-N,N-dimethylamino)ethyl methacrylate) (PNOMA), and poly(2-(N-3-sulfopropyl-N,N-dimethylammonium)ethyl methacrylate) (PSBMA)─grafted onto α-quartz substrates via polymerization initiator linkers. The brush models were developed to closely replicate experimentally synthesized samples and to provide detailed insights into structural and dynamical changes at the molecular level during protein adsorption. Using steered molecular dynamics simulations, we show that the PSBMA brush, due to its high local density, exhibits the greatest resistance to protein insertion. Cα root-mean-square deviation and interaction pattern analyses further reveal that the PSBMA brush also induces the most significant destabilization of lysozyme, while the PDMAEMA brush enhances protein stability through ion-mediated interactions. The PNOMA brush, while requiring the lowest force for protein adsorption, induces greater protein destabilization than the PDMAEMA brush, primarily due to electrostatic repulsion caused by a short carbon spacer length. Hydration analysis reveals that both the PSBMA brush and the lysozyme interacting with it exhibit the most rapid dehydration, attributed to the brush’s high local chain density, which results in the greatest lysozyme destabilization and the highest adsorption force. These findings highlight the dual role of antifouling polymer brushes: resisting protein adsorption and modulating protein structural dynamics. In conclusion, this study provides valuable insights for the rational design of next-generation antifouling materials and offers a framework for realistic model development in complex multicomponent systems.

Adsorption↗

Simplification of the Grid Model and its Impact on the Analysis of Electrical Power Systems

Here, this paper analyzes the impact of the use of Kron reduction on the state variables of a three-phase electrical system, even when it does not meet the necessary conditions for its application. Reduction is applied to a power line model to eliminate the equation corresponding to the neutral conductor of the line. The ATP program is used to model and simulate the behavior of an electrical system considering different degrees of disequilibrium as a reference for the comparison of results. The results show that under certain conditions of disequilibrium the Kron reduction can lead to significant errors in the state variables of the system.

Electric Power Systems↗

The impact of simplified window and exhaust fan assumptions on indoor air quality in multifamily buildings

In residential buildings, the indoor air quality can be significantly affected by ventilation measures initiated by occupants, including the operation of windows and in-unit exhaust fans in kitchens and bathrooms. Many simulations simplify these factors by disregarding window opening behaviors and using fixed schedules for exhaust fan operation across all residential units. To estimate the impact of these simplifications in the U.S. context, this study used coupled CONTAM and EnergyPlus models to simulate airflow and contaminant transport in multifamily buildings. The coupled models parametrically varied climate zone, building airtightness, and mechanical ventilation system types. The study conducted a sensitivity analysis on two key occupant behaviors: (1) operating kitchen and bathroom exhausts on different schedules in individual dwelling units, and (2) scheduling open windows on ground and top floors. The simplified assumptions (i.e. uniform in-unit exhaust fan operation and window operation) had a minimal impact on inter-unit air flow and contaminant transport simulations across a broad range of building air leakage and mechanical ventilation system types. These findings suggest that for buildings with tight construction it is reasonable for most modelling and simulation efforts to ignore the effects of non-uniform exhaust fan operation and window opening.

Occupant behavior↗

Multi-Scale Thermal Analysis for Design of SiC-Based Medium Voltage Motor Drive: Preprint

This paper describes a multi-scale thermal analysis approach for the design of an air-cooled 1.7-kV SiC MOSFET-based medium-voltage variable-speed motor drive. The National Renewable Energy Laboratory supported the project objectives by providing thermal management design support to the project partner The Ohio State University (OSU). The scope of the models and limited resources required efficient and flexible thermal models to be developed. Two modeling techniques are described that significantly reduced model run time and enabled more complex models to be run faster while retaining needed accuracy. The first technique uses the effectiveness-NTU method to extract convection boundary conditions from a computational fluid dynamics (CFD) model that can be applied to a fast-running finite element analysis model. The second is a porous media technique that enables system-level CFD simulations that incorporate effects from heat exchangers (e.g., pin fin heat sinks) that run in a fraction of the time required for fully resolved CFD simulations. The multi-scale approach to the thermal analysis enabled fast and accurate simulation for the converter design ranging from the die level up to the full system with 36 submodules. The modeling results were validated against experimental data from system tests performed by OSU.

DIRECT ENERGY CONVERSION↗

Assessing Uncertainties and Approximations in Solar Heating of the Climate System

In calculating solar radiation, climate models make many simplifications, in part to reduce computational cost and enable climate modeling, and in part from lack of understanding of critical atmospheric information. Whether known errors or unknown errors, the community's concern is how these could impact the modeled climate. The simplifications are well known and most have published studies evaluating them, but with individual studies it is difficult to compare. Here, we collect a wide range of such simplifications in either radiative transfer modeling or atmospheric conditions and assess potential errors within a consistent framework on climate-relevant scales. We build benchmarking capability around a solar heating code (Solar-J) that doubles as a photolysis code for chemistry and can be readily adapted to consider other errors and uncertainties. The broad classes here include: use of broad wavelength bands to integrate over spectral features; scattering approximations that alter phase function and optical depths for clouds and gases; uncertainty in ice-cloud optics; treatment of fractional cloud cover including overlap; and variability of ocean surface albedo. We geographically map the errors in W m -2 using a full climate re-creation for January 2015 from a weather forecasting model. For many approximations assessed here, mean errors are ~2 W m -2 with greater latitudinal biases and are likely to affect a model’s ability to match the current climate state. Combining this work with previous studies, we make priority recommendations for fixing these simplifications based on both the magnitude of error and the ease or computational cost of the fix.

58 GEOSCIENCES↗

Analysis of AP1000 Small-Break Loss-of-Coolant Accident Using Reactor Transient Simulator

The Westinghouse Electric Company’s Advanced Passive Reactor (AP1000) is characterized by the incorporation of passive safety systems (PSSs) designed to ensure core cooling during transient events. The assessment of PSSs requires evaluation of their performance through a combination of experiments and simulations employing various thermal-hydraulic codes. In addition, detailed evaluation of PSSs for a specific reactor system transient analysis such as loss-of-coolant-accident analysis supports understanding representative integral effects test facility development and the further evolution model development and assessment process. Developing a reactor system code is a complex and time-consuming process that requires significant engineering expertise and effort. It can take several months to even years to complete in the early stages of reactor system design and analysis. However, this process can be expedited through the use of transient simulator models for similar reactor systems, which can be used for lesson learning and training purposes. This study uses the Personal Computer Transient Analyzer (PCTRAN) code. The main advantage of PCTRAN is its ease of use and ability to run faster than real time. This study presents the results obtained for a small-break loss-of-coolant accident (SBLOCA) for two breaks using the full version (licensed) of PCTRAN. The purpose of this investigation is to evaluate the overall system behavior during the postulated SBLOCA event as well as assess the capability of the PCTRAN code to reproduce the system response during transient events. The obtained results were compared with the Westinghouse NOTRUMP system code. The PCTRAN code proved to be reliable in predicting the qualitative behavior of the system in both transient cases. As for the system response, it was found that it is contingent on the activation time of the PSSs. The differences in reactor coolant system pressure between the two codes were attributed to the critical flow model and simplification of mass and energy balance. Despite PCTRAN’s limitations, it can still provide a reasonable prediction of various reactor parameters such as pressure, mass flow rate, and void fraction during a SBLOCA scenario. It is worth noting that PCTRAN currently employs a bulk approach similar to that of the Modular Accident Analysis Program (MAAP) and MELCOR codes. However, the upcoming version of PCTRAN will include an artificial intelligence–based detection and accident prevention system, as well as different models for different reactor components. Consequently, PCTRAN has the potential to be upgraded to match the system thermal-hydraulic codes of the U.S. Nuclear Regulatory Commission and become more widely used in cybersecurity to safeguard nuclear power plants from cyberattacks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Ice storage model-predictive control in an office building with PV: scenario, error and sensitivity analysis

Thermal energy storage (TES) can enable more building-sited renewable electricity generation and lower utility bill costs for buildings owners and occupants, especially when there are high demand and variable time-of-use (TOU) charges. A model predictive control (MPC) strategy can offer additional savings over a schedule-based control with added complexity and reliance on forecasts. Here, this study examines savings for medium office buildings with chiller plants in three locations with building-installed solar photovoltaics (PV) to understand the impact of MPC. Control setpoints are fixed by a schedule-based control or optimized by nonlinear MPC. These control setpoints are actuated within EnergyPlus building models to simulate the utility cost of the chiller plant. NLP solutions can be unstable or unrealistic, but our results show that by regularizing the NLP, the solutions can be reasonably followed by the building model. MPC models make simplifications that lead to errors once the controller is participating in and changing the operation of the building. These errors average 9 % across the cases, showing that the most important parts of the system are represented. The no-thermal load costs are computed to show that the optimization can in some cases achieve both the minimum TOU and minimum monthly demand costs by demand management while reducing TOU energy costs by energy arbitrage. The MPC saves 35–66 % in the annual chiller plant operating costs, which is an additional savings above the schedule by 1–33 %. PV and TES are complementary and mostly independent, but a load with PV often results in better performance for the schedule. Our case study and sensitivity analysis show the importance of modeling and optimization for complex rates, but also the circumstances wherein a simpler strategy achieves the same performance with less potential for error.

14 SOLAR ENERGY↗

Prediction of the thermal-hydraulic-mechanical response of a geological repository at large scale and sensitivity analyses

To better understand the Thermal-Hydraulic-Mechanical (THM) response of the rock caused by the high-level radioactive waste (HLW) released heat in the Callovo-Oxfordian formation (COx) in the near-field and far-field areas, a series of coupled THM modelling is performed. This study presents a case study based on the French geological repository for HLW (Cigéo project) to better assess modelling of deep geological repositories (DGR) within the DECOVALEX-2019 framework. In this study, the proposed coupled THM model, implemented in the platform of COMSOL, is validated by comparing the modelled THM results for a point heat source in an infinite rock mass with an analytical solution. It is also validated by comparing its thermal component for the Base Case with accurate results calculated using another numerical method. To shorten the calculation time, the DGR is modelled using six cuboid blocks representing 162 HLW cells and only six central HLW cells of interest modelled in detail. This simplification is validated. This model can easily be used to model not only the near-field THM response of the rock but also the far-field THM response. The influence of the boundary conditions applied on the gallery wall, HLW cell walls, and external surfaces of the model and the influence of different vertical dimensions of the geometry are studied. The sensitivity analyses of the THM parameters of COx on the THM response at different locations are performed. Finally, the influence of using a 2-dimensional (2D) model to represent a 3-dimensional (3D) repository is also investigated.

58 GEOSCIENCES↗

Model Package Report: Composite Analysis Solid Waste Release Model (CASWR Model)

This document describes the implementation of a solid waste form release model in GoldSim for the Hanford Site Composite Analysis (CA) Update. This Composite Analysis Solid Waste Release model (CASWR model) was designed to generate deterministic radionuclide release rates for Hanford’s Central Plateau solid waste disposal sites using single realizations of release model coefficients. Five generalized waste form types are used for the conceptual model of waste release: surplus reactor block, cement, soil-debris, grouted residual waste, and ancillary equipment. The surplus reactor block waste form consists of radionuclide leaching from surplus production graphite reactor core blocks via unspecified processes (White et al., 1984 as cited in PNNL-15965). The cement waste form represents solidified wastes whose permeability is much lower than that of the surrounding soil. The soil-debris waste form type is defined as unconsolidated waste mixed with soil material. Tanks and canyon complexes comprise the grouted residual waste form such that their solid waste will be grouted and capped with a surface barrier at the completion of their cleanup. Finally, the ancillary equipment waste form constitutes contaminant releases from ancillary and auxiliary waste form residues associated with tank farms at closure. Individual sub-models are implemented to numerically represent a respective waste form within the CASWR Model: Surplus Reactor Block Sub-model, Cement Sub-model, SoilDebris Sub-model, Grouted Residual Waste Sub-model, and Ancillary Equipment Submodel. Advection is assumed to be the primary transport process governing the release of radionuclides in the Ancillary Equipment and Soil-Debris Sub-models. Diffusion is assumed to be the primary release process in the Grouted Residual Waste and Cement Sub-models. An unspecified zero-order release process is considered in the Surplus Reactor Block Sub-model due to the lack of information regarding actual processes involved in irradiated graphite leaching. The Surplus Reactor Block, Cement, and SoilDebris Sub-models were compared against analytical solutions (PNNL-11800, Composite Analysis for Low-Level Waste Disposal in the 200 Area Plateau of the Hanford Site). The agreement between results of these analytical solutions and the corresponding waste form models verified their correct implementation in GoldSim. A 1-D modeling abstraction approach for the Grouted Residual Waste and Ancillary Equipment Sub-models was adopted from existing Performance Assessment (PA) models. Despite the simplifications made in these sub-models, they were found to be appropriate representations of the waste forms, similar to what was used in the Waste Management Area C PA model (RPP-ENV-58782, Performance Assessment of Waste Management Area C, Hanford Site, Washington, Rev. 0). A sensitivity analysis was conducted to identify the most influential parameters in each waste form sub-model. Suggestions for considering pH-dependent and redox-dependent release mechanisms are formulated through the development of a conditional constant approach in GoldSim.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Develop a Fast Analysis Solver for Welding Sequence Optimization

During the shipbuilding manufacturing process, materials are exposed to significant stresses, as induced both thermally and mechanically, that alter the intended design and significantly affect the production schedule, labor hours (fitting, welding, rework, etc.), and material structural performance. The type and magnitude of deformation of a given structure depends on many factors such as the material, thickness and quality of components, the process heat input, preheat and inter-pass temperatures, type and size of welds, welding sequence and direction, location, sequence, and degree of fixturing. Numerical simulations using finite element analysis (FEA) have long been used to analyze welding-induced structural distortion. For large assemblies, transient thermal elastic-plastic analysis (TEPA) can take days or weeks to run, and optimization of welding sequence is not feasible. Simplified analysis methods were developed to reduce computational time. However, it is challenging to use these techniques to fully optimize welding sequencing because of their applied simplifications in modeling weld details. A fast analysis solver that could be used by the shipbuilding industry is being developed for optimizing welding sequences by taking full advantage of modern GPU-based HPC hardware and incorporating patented acceleration schemes. The accelerated processing factors are up to 2200 times greater for large, multi-pass welded structures.

Yang, Yu-Ping↗

Noble gas constraints on spent fuel irradiation histories

Fission gas isotopic compositions are sensitive to a variety of reactor operating parameters that include the neutron flux, neutron fluence at different neutron energies, and operating temperature. Measurements of fission gas isotopic compositions thus have potential to constrain reactor simulations for nuclear forensics and safeguards applications. In this paper, we present Kr and Xe isotope measurements from a suite of samples obtained from locations that span the axial length a fuel pin with a well-characterized irradiation history and compare these data to spatially resolved reactor simulations. Here, we observed positive correlations between fluence sensitive isotopic ratios and burnup and between a flux sensitive ratio and power, although some discrepancies are observed between the measured data and model predictions. These differences may be due to simplifications in the model and/or inaccuracies in the cross sections. A much broader measurement to model comparison is required to better understand the discrepancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Increasing Functional Diversity in a Global Land Surface Model Illustrates Uncertainties Related to Parameter Simplification

Abstract Simulations of the land surface carbon cycle typically compress functional diversity into a small set of plant functional types (PFT), with parameters defined by the average value of measurements of functional traits. In most earth system models, all wild plant life is represented by between five and 14 PFTs and a typical grid cell (≈100 × 100 km) may contain a single PFT. Model logic applied to this coarse representation of ecological functional diversity provides a reasonable proxy for the carbon cycle, but does not capture the non‐linear influence of functional traits on productivity. Here we show through simulations using the Energy Exascale Land Surface Model in 15 diverse terrestrial landscapes, that better accounting for functional diversity markedly alters predicted total carbon uptake. The shift in carbon uptake is as great as 30% and 10% in boreal and tropical regions, respectively, when compared to a single PFT parameterized with the trait means. The traits that best predict gross primary production vary based on vegetation phenology, which broadly determines where traits fall within the global distribution. Carbon uptake is more closely associated with specific leaf area for evergreen PFTs and the leaf carbon to nitrogen ratio in deciduous PFTs.

58 GEOSCIENCES↗

Parametric Study of 239 Pu and 235 U Mass Equivalency in Spherical Geometries

In this technical document, parametric studies were performed to determine if a bounding mass equivalency can be used to simplify multi-material models in MCNP6 for criticality evaluations. The studies explored the simplification of MCNP6 models by substituting materials, particularly 239 Pu for 235 U. The initial hypothesis for the material simplification is based on the subcritical mass limits from the ANSI/ANS- 8.1 standard which states that for water reflected metal spheres the limits are 20.1 kg 235 U and 5.0 kg 239 Pu. Using those masses, it was hypothesized that a 20.1 kg sphere of 235 U and a 5.0 kg sphere of 239 Pu could be considered to have the same multiplication. Given that equivalency, 4kg 235 U could be modeled by 1kg 239 Pu (i.e., a mass ratio of 4:1). The parametric studies compare the materials 239 Pu and 235 U independently and in multi-material spherical models. Critical mass, critical thickness, critical mass ratio and critical thickness ratio values are presented as functions of inner radius. keff values are presented as a function of inner radius at different mass ratios and total masses. Results indicate that the conservative bounding mass equivalency is one. That is, if substituting 239 Pu for 235 U it should be done on a gram for gram basis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unpacking model inadequacy: The quantification of silver release from TRISO fuel by considering empirical and mechanistic approaches

Increasing adoption of the proposed tristructural isotropic (TRISO) particle fuel for both advanced and existing reactors makes it critical to assess and address any uncertainties and inadequacies of TRISO fission product release models. Model inadequacy stems from simplifications made to the computational model when compared to the experiments. The modeling and simulation efforts conducted using the BISON fuel performance code, along with the experimental campaigns carried out under the Advanced Gas Reactor Fuel Development and Qualification Program, afford a unique opportunity to conduct a rigorous modeling inadequacy assessment within the Bayesian uncertainty quantification (UQ) framework. Here, this study compares the standard Bayesian framework against the Kennedy-O'Hagan (KOH) framework, which explicitly represents modeling inadequacy, in regard to UQ for TRISO silver release models. For this purpose, both the traditional Arrhenius equation fitted to experimental data and the more advanced lower-length-scale (LLS)-informed model, which considers microstructure information, are independently considered. Applying the inverse UQ process on the AGR-2 and -3/4 datasets revealed modeling inadequacy to be the most dominant source of uncertainty. Experimental noise uncertainty is also significant; however, model parameter uncertainty can be considered negligible. Interestingly, both the Arrhenius equation and the LLS-informed model demonstrated similar levels of modeling inadequacy. For the forward predictive UQ, the KOH framework improved both the accuracy and quality of quantified uncertainties in comparison to the standard Bayesian framework. This is true for both the Arrhenius equation and the LLS-informed model. In comparing these modeling approaches, both demonstrated similar performance at the engineering scale, while the LLS-informed model expectedly outperformed the Arrhenius equation at the mesoscale. These conclusions highlight the importance of explicitly accounting for modeling inadequacy in the UQ process, and reinforce the need for continuous refinement of physics-based models in order to address the modeling inadequacy.

Advanced reactors↗

A Simplified Building Modeling Approach for Identifying Whole Building Retrofits

Building Energy Modeling (BEM) is an effective strategy for optimizing new and retrofit building design, evaluating its energy savings potential and rating its energy performance. However, because of the complexity and expense of modeling, only a small percentage of buildings are simulated. Over 80% of buildings are <25,000 ft2, and the cost associated with BEM can be a big deterrent for reaching this subsector, which usually ends up complying with energy codes and qualifying for incentives using prescriptive measures or generalized design guides. Similarly, efficiency considerations for small building retrofits are often identified through deemed measures or Technical Reference Manuals, which though a scalable and convenient approach, can limit innovation in design and optimization of measures. Simplification of the modeling process, where appropriate, can lead to increased energy savings and more informed decision making. This paper will discuss the development of a ruleset for a simplified approach based on the Performance Rating Method (PRM) contained in ASHRAE Standard 90.1. The approach aims to simplify and lower the cost of the energy modeling process in a controlled and documented manner. An anticipated outcome is an increased use of BEM and its effectiveness in the design, retrofit and operation of commercial buildings. In addition, simplification of the modeling process will result in fewer errors from misinterpretations of both program requirements and program intent in simulation software. The paper will discuss the need and advantages of a simplified modeling approach for whole building retrofits and the technical process for achieving the same.

Goel, Supriya↗

Embedded Error Bayesian Calibration of Thermal Decomposition of Organic Materials

Organic materials are an attractive choice for structural components due to their light weight and versatility. However, because they decompose at low temperatures relative to tradiational materials they pose a safety risk due to fire and loss of structural integrity. To quantify this risk, analysts use chemical kinetics models to describe the material pyrolysis and oxidation using thermogravimetric analysis. This process requires the calibration of many model parameters to closely match experimental data. Previous efforts in this field have largely been limited to finding a single best-fit set of parameters even though the experimental data may be very noisy. Furthermore the chemical kinetics models are often simplified representations of the true de- composition process. The simplification induces model-form errors that the fitting process cannot capture. In this work we propose a methodology for calibrating decomposition models to thermogravimetric analysis data that accounts for uncertainty in the model-form and experimental data simultaneously. The methodology is applied to the decomposition of a carbon fiber epoxy composite with a three-stage reaction network and Arrhenius kinetics. The results show a good overlap between the model predictions and thermogravimetric analysis data. Uncertainty bounds capture devia- tions of the model from the data. The calibrated parameter distributions are also presented. In conclusion, the distributions may be used in forward propagation of uncertainty in models that leverage this material.

36 MATERIALS SCIENCE↗