Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Robustness”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Resource-Adaptive Federated Text Generation with Differential Privacy

In cross-silo federated learning (FL), sensitive text datasets remain confined to local organizations due to privacy regulations, making repeated training for each downstream task both communication-intensive and privacy-demanding. A promising alternative is to generate differentially private (DP) synthetic datasets that approximate the global distribution and can be reused across tasks. However, pretrained large language models (LLMs) often fail under domain shift, and federated finetuning is hindered by computational heterogeneity: only resource-rich clients can update the model, while weaker clients are excluded, amplifying data skew and the adverse effects of DP noise. We propose a flexible participation framework that adapts to client capacities. Strong clients perform DP federated finetuning, while weak clients contribute through a lightweight DP voting mechanism that refines synthetic text. To ensure the synthetic data mirrors the global dataset, we apply control codes (e.g., labels, topics, metadata) that represent each client’s data proportions and constrain voting to semantically coherent subsets. This two-phase approach requires only a single round of communication for weak clients and integrates contributions from all participants. Experiments show that our framework improves distribution alignment and downstream robustness under DP and heterogeneity.

Wang, Jiayi [ORNL]↗
DOE OSTI2017

Thermo-mechanical assessment of full SiC/SiC composite cladding for LWR applications with sensitivity analysis

SiC/SiC composites are considered among leading candidates for accident tolerant fuel cladding in light water reactors. However, when SiC-based materials are exposed to neutron irradiation, they experience significant changes in dimensions and physical properties. Under a large heat flux application (i.e. fuel cladding), the non-uniform changes in the dimensions and physical properties will lead to build-up of stresses in the structure over the course of time. To ensure reliable and safe operation of such a structure it is important to assess its thermo-mechanical performance under in-reactor conditions of irradiation and elevated temperature. In this work, the foundation for 3D thermo-mechanical analysis of SiC/SiC cladding is put in place and a set of analyses with simplified boundary conditions has been performed. The analyses were carried out with two different codes that were benchmarked against one another and prior results in the literature. A constitutive model is constructed and solved numerically to predict the stress distribution and variation in the cladding under normal operating conditions. The dependence of dimensions and physical properties variation with irradiation and temperature has been incorporated. These robust models may now be modified to take into account the axial and circumferential variation in neutron and heat flux to fully account for 3D effects. Here, the results from the simple analyses show the development of high tensile stresses especially in the circumferential and axial directions at the inner region of the cladding. Based on the results obtained, design guidelines are recommended. For lack of certainty in or tailor-ability for the physical and mechanical properties of SiC/SiC composite material a sensitivity analysis is conducted. The analysis results establish a precedence order of the properties based on the extent to which these properties influence the temperature and the stresses.

Accident tolerant fuel↗
DOE OSTI2021

Measuring magnetic fields in laser-driven coils with dual-axis proton deflectometry

Abstract By driving hot electrons between two metal plates connected by a wire loop, high power lasers can generate multi-tesla, quasi-static magnetic fields in miniature coil targets. Many experiments involving laser-coil targets rely on proton deflectometry directed perpendicular to the coil axis to extract a measurement of the magnetic field. In this paper, we show that quantitative measurements using perpendicular probing are complicated by the presence of GV m −1 electric fields in the target that develop on sub-ns timescales. Probing parallel to the coil axis with fiducial grids is shown to reliably separate the electric and magnetic field measurements, giving current estimates of I ≈ 5 kA in 1 mm- and 2 mm-diameter wire loops. An analytic model of proton deflection in electric and magnetic fields is used to benchmark results from the particle-in-cell code and help deconvolve the magnetic and electric field deflections. Results are used to motivate a new experimental scheme that combines a single-plate target with axial proton probing and direct current measurements. This scheme has several important advantages over the traditional target and diagnostic set-up, enabling the robust measurement of coil magnetic fields and plasma properties, as well as making it easier to validate different theoretical models at a range of laser intensities.

Bradford, P. (ORCID:0000000296186342)↗
DOE OSTI2023

Fully Implicit Conjugate Heat Transfer Analysis of the ARC-Class Vacuum Vessel

The coupled simulation of fusion reactor blankets including neutronics, thermal-hydraulics and thermo-mechanics is expected to speed up the design cycle of fusion reactor design concepts. In this work we demonstrate tight implicit coupling of conjugate heat transfer using the open-source Computational Fluid Dynamics software OpenFOAM for thermo-fluid mechanics and Diablo for thermo-solid mechanics. The heat transfer analysis is augmented by volumetric energy deposition from neutronic calculations using the Monte Carlo N-particle code on both solid and fluid parts of the vacuum vessel. An additional heat flux is imposed on the first wall estimated from the design power of the reactor. The tight coupling is realized through the open-source coupling library, preCICE, and tested on the vacuum vessel of the affordable, robust, compact reactor design by Commonwealth Fusion Systems. The features of the coupling and the influence of different coupling parameters such as coupling schemes, acceleration techniques and convergence criterion are discussed. The coupled simulation results are compared to a thermal-hydraulics simulation which includes only the fluid domains (the liquid immersion molten salt blanket and cooling channel) to demonstrate usefulness of a coupled simulation. Further analysis is performed to identify regions of hot spots for subsequent design improvement. This introduces the outline for integrating conjugate electromagnetics and fluid/solid mechanics (e.g., allow for deformation of the cooling channel walls) with our present approach for future analysis.

Sircar, Arpan↗
DOE OSTI2025

Probability Density Function for the spatial and intensity distribution of neutron-induced defects in Silicon

The ability to model semiconductor device degradation under neutron irradiation depends upon having a robust modeling capability for the neutron-induced collision cascades as well as a means to analytically fit the resulting probability distributions of defect production and ionizing energy deposition for purposes of extrapolation to low-probability, high-consequence scenarios. In this paper, the widely-utilized binary collision approximation codes MARLOWE and SRIM are deployed in conjunction with a critical examination of their parameterizations as benchmarked against higher-fidelity molecular dynamics simulations. A simple 3-parameter form described by the Generalized Logistic Distribution is shown to be a good fit to Frenkel pair and ionization intensity distributions in bulk silicon. The BCA codes are then applied to simulate cascades in 5 nm layers of a representative gate-all-around nanosheet transistor, where joint probability distributions of threshold levels of damage to multiple layers are evaluated.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The Variable Infiltration Capacity model version 5 (VIC-5): infrastructure improvements for new applications and reproducibility

The Variable Infiltration Capacity (VIC) model is a macroscale semi-distributed hydrologic model. VIC development began in the early 1990s and the model has since been used extensively for basin- to global-scale applications that include hydrologic dataset construction, trend analysis of hydrologic fluxes and states, data evaluation and assimilation, forecasting, coupled climate modeling, and climate change impact assessment. Ongoing operational applications of the VIC model include the University of Washington's drought monitoring and forecasting systems and NASA's Land Data Assimilation System. This paper documents the development of VIC version 5 (VIC-5), which includes a major reconfiguration of the legacy VIC source code to support a wider range of modern hydrologic modeling applications. The VIC source code has been moved to a public GitHub repository to encourage participation by the broader user and developer communities. The reconfiguration has separated the core physics of the model from the driver source code, whereby the latter is responsible for memory allocation, preprocessing and post-processing, and input–output (I–O). VIC-5 includes four drivers that use the same core physics modules, but which allow for different methods for accessing this core to enable different model applications. Finally, VIC-5 is distributed with robust test infrastructure, components of which routinely run during development using cloud-hosted continuous integration. The work described here provides an example to the model development community for extending the life of a legacy model that is being used extensively. The development and release of VIC-5 represents a significant step forward for the VIC user community in terms of support for existing and new model applications, reproducibility, and scientific robustness.

42 ENGINEERING↗
DOE OSTI2017

Orbits for 18 Visual Binaries and Two Double-line Spectroscopic Binaries Observed with HRCAM on the CTIO SOAR 4 m Telescope, Using a New Bayesian Orbit Code Based on Markov Chain Monte Carlo

We present orbital elements and mass sums for 18 visual binary stars of spectral types B to K (five of which are new orbits) with periods ranging from 20 to more than 500 yr. For two double-line spectroscopic binaries with no previous orbits, the individual component masses, using combined astrometric and radial velocity data, have a formal uncertainty of ∼0.1 M{sub ⊙}. Adopting published photometry and trigonometric parallaxes, plus our own measurements, we place these objects on an H-R diagram and discuss their evolutionary status. These objects are part of a survey to characterize the binary population of stars in the Southern Hemisphere using the SOAR 4 m telescope+HRCAM at CTIO. Orbital elements are computed using a newly developed Markov chain Monte Carlo (MCMC) algorithm that delivers maximum-likelihood estimates of the parameters, as well as posterior probability density functions that allow us to evaluate the uncertainty of our derived parameters in a robust way. For spectroscopic binaries, using our approach, it is possible to derive a self-consistent parallax for the system from the combined astrometric and radial velocity data (“orbital parallax”), which compares well with the trigonometric parallaxes. We also present a mathematical formalism that allows a dimensionality reduction of the feature space from seven to three search parameters (or from 10 to seven dimensions—including parallax—in the case of spectroscopic binaries with astrometric data), which makes it possible to explore a smaller number of parameters in each case, improving the computational efficiency of our MCMC code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Creating Continuous Integration Infrastructure for Software Development on U.S. Department of Energy High-Performance Computing Systems

The Exascale Computing Project (ECP) software deployment effort developed and advanced DevOps capabilities. One goal was to enable robust continuous integration (CI) workflows that span the protected high performance computing (HPC) environments found within many of the Department of Energy’s (DOE) national laboratories. This article highlights several challenges encountered with enabling automation, such as charging models for CI jobs, and meeting individualized security requirements that revolve around strongly associating running code with a human identity. Here, it also describes how the Jacamar CI tool evolved to meet latter requirements and became a key aspect of the solutions currently offered. Derived from this experience, we offer a conceptual framework for understanding current and future CI challenges at DOE facilities and offer suggestions for long-term solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Developing an Automated Uncertainty Quantification Tool to Improve Watershed-Scale Predictions of Water and Nutrient Cycling

Managing the flow of water, nutrients, and contaminants in watersheds is vital to addressing pressing issues related to water scarcity, access to clean drinking water, energy production, resilience to natural and anthropogenic perturbations, and ecological restoration. Decisions about the management of watersheds critically depend on the accuracy with which the flow of water and chemicals through the watershed can be predicted by computer models. Prediction uncertainty can be reduced by matching the model to data, which are collected in the field at great expense. The contribution of watershed characterization data to reducing uncertainty of relevant model predictions can be evaluated in a so-called data-worth analysis, which provides transparent, quantitative metrics about a data set’s value for the support of relevant watershed management objectives. To achieve this goal, we developed a software package that implements the data-worth analysis approach for use with state-of-the-art watershed models. The purpose of the proposed data-worth analysis is to help decision-makers allocate resources for watershed characterization such that the uncertainty in model predictions can be significantly reduced, which leads to better, more effective management decisions. At the same time, watershed characterization costs can be reduced. The specific technical objectives of this SBIR/STTR Phase II project were to develop a framework and associated software toolsets that implement the uncertainty quantification and data-worth analysis approach for use with state-of-the-art watershed models. This goal was achieved by (A) developing a user-friendly, robust software package that is accessible to a wide audience, including watershed managers, policy-makers, and public stakeholders; (B) by demonstrating application of the prototype on several use cases that are representative of complex watershed management challenges spanning a range of scales and that consider different open-source, DOE-based codes and other modeling platforms; and (C) by gathering information about the needs and requirements from potential users to help guide future developments, ensuring that the final product will be commercially viable. The developed software consists of a graphical user interface that guides the user through a sequence of analysis steps, supported by toolsets that leverage state-of-the-art computational simulation-optimization capabilities. A prototype of the software runs on multiple platforms (PC, Mac, multi-processor Linux environment), is linked to diverse watershed simulators (e.g., ECOSYS, TOUGH2, TOUGHREACT, Amanzi-ATS), performs multiple analysis tasks (predictive simulations, sensitivity analysis, uncertainty analysis, automatic parameter estimation, and data-worth analysis, multicomponent geothermometry), and is readily extensible to include external simulators and analysis tools. The software is being commercialized and will be continually updated to address user needs.

58 GEOSCIENCES↗
DOE OSTI2026

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2022

Guidelines for Publicly Archiving Terrestrial Model Data to Enhance Usability, Intercomparison, and Synthesis

Scientific communities are increasingly publishing data to evaluate, accredit, and build on published research. However, guidelines for curating data for publication are sparse for model-related research, limiting the usability of archived simulation data. In particular, there are no established guidelines for archiving data related to terrestrial models that simulate land processes and their coupled interactions with climate. Terrestrial modelers have a unique set of challenges when publishing data due to the diversity of scientific domains, research questions, and the types and scales of simulations. Researchers in the U.S. Department of Energy’s (DOE) projects use a variety of multiscale models to advance robust predictions of terrestrial and subsurface ecosystem processes. Here, we synthesize archiving needs for data associated with different DOE models, and provide guidelines for publishing terrestrial model data components following FAIR (Findable, Accessible, Interoperable, Reusable) principles. The guidelines recommend archiving model inputs and testing data used in final simulation runs along with associated codes, workflow scripts, and metadata in public repositories. Researchers should consider archiving model outputs if they are within the storage limits of the repository. We also provide considerations for how to bundle files into different data publications with citable digital object identifiers. Finally, we identify repository features and tools that would enable storage and reuse of model data. Given the diversity of DOE terrestrial models, these guidelines are transferable to other model types and will enable efficient reuse of simulation data for purposes such as model intercomparisons, initialization, benchmarking, synthesis, and comparisons with field observations.

58 GEOSCIENCES↗
DOE OSTI2021

Resilient Information Architecture Platform for Smart Grid (RIAPS)

A number of emerging trends will substantially alter the operation and control of the electric grid over the next several decades. These trends include ensuring resiliency under severe weather events, increasing integration of renewable electricity generation, supporting changing electricity demand patterns, and the improving cost effectiveness of distributed energy resources. To address these challenges, the future “Smart Grid” management will need to transition from centralized to coordinated distributed control paradigm. Reliable operation of the Smart Grid depends on distributed intelligence realized through software applications that run on distributed computing devices attached to the power system to collect data and collaboratively manage resources. However, much of the existing software for Smart Grid-enabled devices is either proprietary or developed with custom solutions, which limits interoperability among the heterogeneous devices and hinders the ability to manage system-level reliability, security, and resiliency requirements. Additionally, this approach makes Smart Grid applications hard to maintain, evolve, verify, and replace; resulting in high development and deployment costs. Further development of the Smart Grid requires a reusable software base-layer to move from hard-coded functionality to a plug-and-play architecture capable of managing system-level objectives and constraints in addition to providing consistent common services across heterogeneous devices and applications. Vanderbilt University, in collaboration with North Carolina State University and Washington State University has developed a foundation ‘software platform’ for developing and deploying robust, reliable, effective and secure software applications for the Smart Grid. The Resilient Information Architecture Platform for the Smart Grid (RIAPS) provides core services for building effective and powerful smart grid applications. It offers unique services for real-time data dissemination, fault tolerance, and coordination across apps distributed over the network.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Plasma edge and plasma-wall interaction modelling: Lessons learned from metallic devices

Robust power exhaust schemes employing impurity seeding are needed for target operational scenarios in present day tokamak devices with metallic plasma-facing components (PFCs). For an electricity-producing fusion power plant at power density P sep /R > 15 MW/m divertor detachment is a requirement for heat load mitigation. 2D plasma edge transport codes like the SOLPS code as well as plasma-wall interaction (PWI) codes are key to disentangle relevant physical processes in power and particle exhaust. With increased quantitative credibility in such codes more realistic and physically sound estimates of the life-time expectations and performance of metallic PFCs can be accomplished for divertor conditions relevant for ITER and DEMO. An overview is given on the recent progress of plasma edge and PWI modelling activities for (carbon-free) metallic devices, that include results from JET with the ITER-like wall, ASDEX Upgrade and Alcator C-mod. It is observed that metallic devices offer an opportunity to progress the understanding of underlying plasma physics processes in the edge. The validation of models can be substantially improved by eliminating carbon from the experiment as well as from the numerical system with reduced degrees of freedom as no chemical sputtering from amorphous carbon layers and no carbon or hydro-carbon transport are present. With the absence of carbon as the primary plasma impurity and given the fact that the physics of the PWI at metallic walls is less complex it is possible to isolate the crucial plasma physics processes relevant for particle and power exhaust. Finally, for a reliable 2D dissipative plasma exhaust model these are: cross-field drifts, complete kinetic neutral physics, geometry effects (including main-chamber, divertor and sub-divertor structures), SOL transport reflecting also the non-diffusive nature of anomalous transport, as well as transport within the pedestal region in case of significant edge impurity radiation affecting pedestal pressure and hence P sep .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Major Revisions in Arthropod Phylogeny Through Improved Supermatrix, With Support for Two Possible Waves of Land Invasion by Chelicerates

Deep phylogeny involving arthropod lineages is difficult to recover because the erosion of phylogenetic signals over time leads to unreliable multiple sequence alignment (MSA) and subsequent phylogenetic reconstruction. One way to alleviate the problem is to assemble a large number of gene sequences to compensate for the weakness in each individual gene. Such an approach has led to many robustly supported but contradictory phylogenies. A close examination shows that the supermatrix approach often suffers from two shortcomings. The first is that MSA is rarely checked for reliability and, as will be illustrated, can be poor. The second is that, to alleviate the problem of homoplasy at the third codon position of protein-coding genes due to convergent evolution of nucleotide frequencies, phylogeneticists may remove or degenerate the third codon position but may do it improperly and introduce new biases. We performed extensive reanalysis of one of such “big data” sets to highlight these two problems, and demonstrated the power and benefits of correcting or alleviating these problems. Our results support a new group with Xiphosura and Arachnopulmonata (Tetrapulmonata + Scorpiones) as sister taxa. This favors a new hypothesis in which the ancestor of Xiphosura and the extinct Eurypterida (sea scorpions, of which many later forms lived in brackish or freshwater) returned to the sea after the initial chelicerate invasion of land. Our phylogeny is supported even with the original data but processed with a new “principled” codon degeneration. We also show that removing the 1673 codon sites with both AGN and UCN codons (encoding serine) in our alignment can partially reconcile discrepancies between nucleotide-based and AA-based tree, partly because two sequences, one with AGN and the other with UCN, would be identical at the amino acid level but quite different at the nucleotide level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

The Development of the Generalized Spallation Model [Thesis]

The generalized spallation model (GSM) is a general spallation physics event generator with a physics-rich history dating back to the early 1970s during the development of the CEM. GSM estimates the outcome of high-energy hadron-nucleus and nucleus-nucleus collisions, where collision energies often exceed hundreds of MeV. Such a collision may include a neutron incident on carbon, for example. The thesis presented here outlines the continued development and modernization of GSM, whereby it was migrated from a legacy Fortran code to utilize an advanced and modern Fortran architecture. Verification and validation of GSM also took place throughout the continued development of GSM. The results of this work have transformed GSM into a highly modern and powerful event generator, where numerous benefits have been fostered and precipitated, including but not limited to portability, modularization, flexibility, robustness, and scalability. GSM is a robust and well developed software library that clients and end-users may confidently use, within the scope outlined here, as a result of this work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) Version 4.1 Technical Reference Manual

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 4.1. The two most significant changes for HyRAM+ version 4.1 from HyRAM+ version 4.0 are direct incorporation of unconfined overpressure into the QRA calculations and modification of the models for cryogenic liquid flow through an orifice. In QRA mode, the user no longer needs to input peak overpressure and impulse values that were calculated separately; rather, the unconfined overpressure is estimated for the given system inputs, leak size, and occupant location. The orifice flow model now solves for the maximum mass flux through the orifice at constant entropy while conserving energy, which does not require a direct speed of sound calculation. This does not affect the mass flow for all-gaseous releases; the method results in the same speed of sound for choked flow. However, this method does result in a higher (and more realistic) mass flow rate for a given leak size for liquid releases than was previously calculated.

08 HYDROGEN↗
DOE OSTI2017

SOLVING CONJUGATE HEAT TRANSFER PROBLEM IN TWO-PHASE FLOW USING NEWTON-KRYLOV METHOD

In this paper, we successfully demonstrated an application of the Jacobian-free Newton Krylov (JFNK) method in solving the conjugated two-phase channel flow and heat conduction problem, in a fully coupled and fully implicit manner. Such a problem represents the most fundamental piece of safety analysis codes, for example the heat transfer between fuel rods and reactor coolant during postulated accidental scenarios. The two-phase flow is modeled with the widely accepted one-dimensional two-fluid two-phase flow model, and the heat conduction is modeled with two-dimensional heat conduction model. The coupling between these two physics is realized through the wall heat transfer term in fluid equations and heat flux boundary conditions of heat conduction model. For spatial discretization, both the first-order and second- order upwind schemes are considered for flow equations, while the standard central differencing method is used for the heat conduction model. For time integration, both the first-order backward Euler (BDF1) and second-order BDF2 schemes are considered. Several treatments were also introduced to enhance the robustness of the JFNK method in solving such a problem that is discontinuous in nature because of the discrete closure correlations. These treatments include the re-introduction of a transition region between the slug and annular mist flow regimes, and a simplification in the interfacial drag term. Numerical results have been extensively validated with experimental data in both vertical pipe and rod bundle geometries with a wide range of thermal-hydraulics conditions, including pressure, inlet subcooling, wall heat flux, outlet void fraction, and etc. Further benchmark with the RELAP5-3D code confirms the correct implementation of wall heat transfer models.

42 ENGINEERING↗