Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

In Situ X-ray Radiography and Computational Modeling to Predict Grain Morphology in $\beta$-Titanium during Simulated Additive Manufacturing

The continued development of metal additive manufacturing (AM) has expanded the engineering metallic alloys for which these processes may be applied, including beta-titanium alloys with desirable strength-to-density ratios. To understand the response of beta-titanium alloys to AM processing, solidification and microstructure evolution needs to be investigated. In particular, thermal gradients (Gs) and solidification velocities (Vs) experienced during AM are needed to link processing to microstructure development, including the columnar-to-equiaxed transition (CET). In this work, in situ synchrotron X-ray radiography of the beta-titanium alloy Ti-10V-2Fe-3Al (wt.%) (Ti-1023) during simulated laser-powder bed fusion (L-PBF) was performed at the Advanced Photon Source at Argonne National Laboratory, allowing for direct determination of Vs. Two different computational modeling tools, SYSWELD and FLOW-3D, were utilized to investigate the solidification conditions of spot and raster melt scenarios. The predicted Vs obtained from both pieces of computational software exhibited good agreement with those obtained from in situ synchrotron X-ray radiography measurements. The model that accounted for fluid flow also showed the ability to predict trends unobservable in the in situ synchrotron X-ray radiography, but are known to occur during rapid solidification. A CET model for Ti-1023 was also developed using the Kurz–Giovanola–Trivedi model, which allowed modeled Gs and Vs to be compared in the context of predicted grain morphologies. Both pieces of software were in agreement for morphology predictions of spot-melts, but drastically differed for raster predictions. The discrepancy is attributable to the difference in accounting for fluid flow, resulting in magnitude-different values of Gs for similar Vs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computational model development and validation of fuel dispersal phenomena

As the US nuclear industry is renewing efforts to extend the rod average burnup limits in existing pressurized water reactors, it becomes imperative to investigate the safety implications associated with high burnup fuel fragmentation and the potential dispersal of fragmented fuel into the reactor coolant system. Additionally, as there is a growing interest in transitioning to new nuclear fuel designs to accommodate higher fuel burnup levels, it becomes essential to conduct studies on fuel behavior during design basis accident scenarios, particularly during a postulated loss-of-coolant accident, in order to ensure the safe operation of light-water reactors. The present study develops a computational model to simulate the complex three-phase flow of the fuel dispersal phenomena that occur following a breach in the cladding, when the stream of solid fuel particles and carrier fission gas stream may interact with the surrounding liquid medium. Furthermore, the proposed framework considers the solid phase as discrete Lagrangian particles while the interpenetrating gas–liquid continuum is modeled using a Eulerian framework. The developed simulation framework is validated with experimental results performed in a separate-effect test facility by comparing predicted particle settlement locations with the experimental observations reported in the literature. Key highlights of this study include the interpretation of the high-pressure boundary conditions and the three-phase flow coupling strategy. Also discussed are the gas–liquid hydrodynamics and the behavior of particles as they are transported through the gas and liquid phases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Optimal size of the block in block GMRES on GPUs: computational model and experiments

The block version of GMRES (BGMRES) is most advantageous over the single right hand side (RHS) counterpart when the cost of communication is high while the cost of floating point operations is not. This is the particular case on modern graphics processing units (GPUs), while it is generally not the case on traditional central processing units (CPUs). Here, in this paper, experiments on both GPUs and CPUs are shown that compare the performance of BGMRES against GMRES as the number of RHS increases, with a particular focus on GPU performance. The experiments indicate that there are many cases in which BGMRES is slower than GMRES on CPUs, but faster on GPUs. Furthermore, when varying the number of RHS on the GPU, there is an optimal number of RHS where BGMRES is clearly most advantageous over GMRES. A computational model for the GPU is developed using hardware specific parameters, providing insight towards how the qualitative behavior of BGMRES changes as the number of RHS increase, and this model also helps explain the phenomena observed in the experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

NANO.PTML model for read-across prediction of nanosystems in neurosciences. computational model and experimental case of study

Abstract Neurodegenerative diseases involve progressive neuronal death. Traditional treatments often struggle due to solubility, bioavailability, and crossing the Blood-Brain Barrier (BBB). Nanoparticles (NPs) in biomedical field are garnering growing attention as neurodegenerative disease drugs (NDDs) carrier to the central nervous system. Here, we introduced computational and experimental analysis. In the computational study, a specific IFPTML technique was used, which combined Information Fusion (IF) + Perturbation Theory (PT) + Machine Learning (ML) to select the most promising Nanoparticle Neuronal Disease Drug Delivery (N2D3) systems. For the application of IFPTML model in the nanoscience, NANO.PTML is used. IF-process was carried out between 4403 NDDs assays and 260 cytotoxicity NP assays conducting a dataset of 500,000 cases. The optimal IFPTML was the Decision Tree (DT) algorithm which shown satisfactory performance with specificity values of 96.4% and 96.2%, and sensitivity values of 79.3% and 75.7% in the training (375k/75%) and validation (125k/25%) set. Moreover, the DT model obtained Area Under Receiver Operating Characteristic (AUROC) scores of 0.97 and 0.96 in the training and validation series, highlighting its effectiveness in classification tasks. In the experimental part, two samples of NPs (Fe 3 O 4 _A and Fe 3 O 4 _B) were synthesized by thermal decomposition of an iron(III) oleate (FeOl) precursor and structurally characterized by different methods. Additionally, in order to make the as-synthesized hydrophobic NPs (Fe 3 O 4 _A and Fe 3 O 4 _B) soluble in water the amphiphilic CTAB (Cetyl Trimethyl Ammonium Bromide) molecule was employed. Therefore, to conduct a study with a wider range of NP system variants, an experimental illustrative simulation experiment was performed using the IFPTML-DT model. For this, a set of 500,000 prediction dataset was created. The outcome of this experiment highlighted certain NANO.PTML systems as promising candidates for further investigation. The NANO.PTML approach holds potential to accelerate experimental investigations and offer initial insights into various NP and NDDs compounds, serving as an efficient alternative to time-consuming trial-and-error procedures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Unraveling the Mechanics: Computational Modeling of Residual Stresses in U10Mo Fuel with Aluminum Cladding

LEU (Low Enriched Uranium) plate-type fuel elements consist of a high-density, low-enrichment U–Mo alloy-based fuel foil encapsulated in an aluminum cladding and is fabricated through the Hot Isostatic Pressing (HIP) technique. Understanding the mechanism of possible failure modes of nuclear fuel is critical to mitigate potential consequences. One of the major contributing factors in various failure modes is stress and it greatly affects the integrity under temperature, pressure, and irradiation. These stresses could originate from various sources, including manufacturing, and operation, and can lead to deformation, cracking, or even complete failure of fuel elements. This study focuses on the development of a computational model that accurately predicts the residual stresses generated in the U10Mo fuel during fabrication. It has been observed that cladding creep has a substantial impact on the residual stresses in the U10Mo fuel post the HIP fabrication process. Furthermore, during the HIP bonding process the fuel plate system is heated to a temperature of 560 oC and as a result, the aluminum cladding transitions from (Al 6061-T6 to -O). This presents a challenge in capturing the change in material properties accurately. Therefore, a temperature dependent creep model such as hyperbolic sine creep model is considered to estimate the creep properties of the Aluminum cladding. The proposed calibrated creep model accurately predicts the residual stresses in the U10Mo fuel foil and agrees well with the experiments.

42 ENGINEERING↗
DOE OSTI2024

Establishing a clostridia foundry for biosystems design by integrating computational modeling, systems-level analyses, and cell-free engineering technologies (Final Report)

Rapid population growth, a rise in global living standards, and economic competitiveness have intensified the need for sustainable, low-cost biofuels and bioproducts production. Industrial biotechnology using microbial cell factories – and waste and/or renewable feedstocks – is one of the most attractive approaches for addressing this need, particularly when large-scale chemical synthesis is untenable. Unfortunately, designing, building, and optimizing biosynthetic pathways in cells remains a complex challenge. With support from the Department of Energy, we worked to address this challenge in a new interdisciplinary venture that established the world’s first clostridial Foundry for Biosystems Design (cBioFAB). Working both in vitro (cell-free) and in vivo, the goal of this project was to interweave and advance state-of-the-art computational modeling, genome editing, omics measurements, systems-biology analyses, and cell-free technologies to expand the set of platform organisms that meet DOE bioenergy goals. Specifically, we manufactured fuel and chemical intermediates via existing and de novo pathways. This report covers the outcomes of our research project.

09 BIOMASS FUELS↗
DOE OSTI2024

Computational Modeling, Simulation, and Potential Applications of Optical Stochastic Cooling

With the rising demand for intense particles beams, much research is being conducted in the area of particle beam cooling. One of these methods, called Stochastic Cooling (SC) (developed at CERN in the 1970's), delivered a feedback method to improve the quality and lifetime of circulating proton beams by reducing their 6D phase-space and has been widely implemented in a number of hadron machines. However, traditional stochastic cooling schemes are limited by the bandwidth of microwave frequency systems. Optical Stochastic Cooling (OSC) is a promising extension of the stochastic cooling beam cooling technique. OSC instead uses optical wavelengths which allows for improved control and increased cooling but creates its own technical challenges. This lays out work conducted toward the experimental demonstration of OSC at Fermilab's Integrable Optics Test Accelerator (IOTA) storage ring. This includes the design and characterization of parts of the optical delay system, the development and validation of a high-fidelity computational model of the OSC process, and the investigations into possible applications of the OSC mechanism to advanced beam manipulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Computational modeling and validation of a modified Marple cascade impactor

Oak Ridge National Laboratory (ORNL) is investigating respirable aerosol fractions from failed high-burnup (HBU) commercial spent nuclear fuel (CSNF) rods (Montgomery et al., 2022). The work aims to characterize the source term of radioactive aerosols that could result from mechanical failure and handling of failed fuel rods during storage and/or transportation. The experiment involves capturing and characterizing aerosols from a four-point bending experiment of fueled rod segments at ORNL’s Irradiated Fuels Examination Laboratory (IFEL). For this, a Marple cascade impactor (Series 290) with the ability to collect particles up to aerodynamic equivalent diameters of ~15 μm was tested for feasibility. The commercially available Marple cascade impactor design was modified for the experiment, where the cascade inlet was replaced with a 3D printed inlet nozzle. Collection efficiency curves were then generated with the modified impactor, and the cascade’s performance was tested using ISO dust capture experiments followed by SEM analysis of the collection and computational modeling. The current work highlights the performance testing and experimental validation of this modified Marple Impactor for spent nuclear fuel aerosol collection.

Kumar, Vineet↗
DOE OSTI2022

Computational models of direct and indirect X‐ray breast imaging detectors for in silico trials

Abstract Background To facilitate in silico studies that investigate digital mammography (DM) and breast tomosynthesis (DBT), models replicating the variety in imaging performance of the DM and DBT systems, observed across manufacturers are needed. Purpose The main purpose of this work is to develop generic physics models for direct and indirect detector technology used in commercially available systems, with the goal of making them available open source to manufacturers to further tweak and develop the exact in silico replicas of their systems. Methods We recently reported on an in silico version of the SIEMENS Mammomat Inspiration DM/DBT system using an open‐source GPU‐accelerated Monte Carlo x‐ray imaging simulation code (MC‐GPU). We build on the previous version of the MC‐GPU codes to mimic the imaging performances of two other Food and Drug Administration (FDA)‐approved DM/DBT systems, such as Hologic Selenia Dimensions (HSD) and the General Electric Senographe Pristina (GSP) systems. In this work, we developed a hybrid technique to model the optical spread and signal crosstalk observed in the GSP and HSD systems. MC simulations are used to track each x‐ray photon till its first interaction within the x‐ray detector. On the other hand, the signal spread in the x‐ray detectors is modeled using previously developed analytical equations. This approach allows us to preserve the modeling accuracy offered by MC methods in the patient body, while speeding up secondary carrier transport (either electron–hole pairs or optical photons) using analytical equations in the detector. The analytical optical spread model for the indirect detector includes the depth‐dependent spread and collection of optical photons and relies on a pre‐computed set of point response functions that describe the optical spread as a function of depth. To understand the capabilities of the computational x‐ray detector models, we compared image quality metrics like modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE), simulated with our models against measured data. Please note that the purpose of these comparisons with measured data would be to gauge if the model developed as part of this work could replicate commercially used direct and indirect technology in general and not to achieve perfect fits with measured data. Results We found that the simulated image quality metrics such as MTF, NNPS, and DQE were in reasonable agreement with experimental data. To demonstrate the imaging performance of the three DM/DBT systems, we integrated the detector models with the VICTRE pipeline and simulated DM images of a fatty breast model containing a spiculated mass and a calcium oxalate cluster. In general, we found that the images generated using the indirect model appeared more blurred with a different noise texture and contrast as compared to the systems with direct detectors. Conclusions We have presented computational models of three commercially available FDA‐approved DM/DBT systems, which implement both direct and indirect detector technology. The updated versions of the MC‐GPU codes that can be used to replicate three systems are available in open source format through GitHub.

Sengupta, Aunnasha↗
DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗
DOE OSTI2025

Computational Modeling of a 3D Printed Recuperator and Subsequent Experimental Loop for Supercritical Carbon Dioxide Cycles

Oak Ridge National Laboratory (ORNL), in collaboration with mechanical-thermal energy storage (mTES) provider EarthEn, a US Department of Energy (DOE) Lab-Embedded Entrepreneurship Program (LEEP) recipient at ORNL’s Innovation Crossroads 2023, is utilizing a state-of-the-art patented 3D printing technique to design an additively manufactured (AM) supercritical CO2 (sCO2) recuperator (REC) for EarthEn’s charge/discharge cycle. The AM REC will be printed at ORNL’s Manufacturing Demonstration Facility using Inconel Alloy 718 and tested on a closed-loop, ∼100 kW scale experimental facility that is under construction. The testing will compare the printed design against a commercial-off-the-shelf Printed Circuit Heat Exchanger (PCHE) REC. The design of the sCO2 facility is guided by a Modelica-based system model which is primarily dependent on the open-source TRANSFORM library developed at ORNL and uses the open-source CoolProp library for thermophysical properties of sCO2 via the External Media library. It is envisioned that an iterative process will be followed between the physical loop and the system model wherein the initial experimental data will be used to tune the model, which in turn will be used to guide future loop operation. Simultaneously, the AM REC is being designed using computer-aided design models, and it is also being analyzed for hydraulic and thermomechanical response using commercial computational fluid dynamics software, Simcenter STAR-CCM+, on highperformance computing resources.1

See, Nate [ORNL] (ORCID:0000000178581202)↗
DOE OSTI2025

Combining computational modeling and experimental library screening to affinity-mature VEEV-neutralizing antibody F5

Engineered monoclonal antibodies have proven to be highly effective therapeutics in recent viral outbreaks. However, despite technical advancements, an ability to rapidly adapt or increase antibody affinity and by extension, therapeutic efficacy, has yet to be fully realized. We endeavored to stand-up such a pipeline using molecular modeling combined with experimental library screening to increase the affinity of F5, a monoclonal antibody with potent neutralizing activity against Venezuelan Equine Encephalitis Virus (VEEV), to recombinant VEEV (IAB) E1E2 antigen. We modeled the F5/E1E2 binding interface and generated predictions for mutations to improve binding using a Rosetta-based approach and dTERMen, an informatics approach. The modeling was complicated by the fact that a high-resolution structure of F5 is not available and the H3 loop of F5 exceeds the length for which current modeling approaches can determine a unique structure. A subset of the predicted mutations from both methods were incorporated into a phage display library of scFvs. This library and a library generated by error-prone PCR were screened for binding affinity to the recombinant antigen. Results from the screens identified favorable mutations which were incorporated into 12 human-IgG1 variants. The best variant, containing eight mutations, improved KD from 0.63 nM (parental) to 0.01 nM. While this did not improve neutralization or therapeutic potency of F5 against IAB, it did increase cross-reactivity to other closely related VEEV epizootic and enzootic strains, demonstrating the potential of this method to rapidly adapt existing therapeutics to emerging viral strains.

affinity-maturation↗
DOE OSTI2023

DFO Computational Modeling Project - Catalytic Upgrading of Bio-based Furfural to 1,5-Pentanediol: A New Renewable Monomer for the Coatings Industry (Final Report)

This report summarizes the results of a collaborative effort between Oak Ridge National Laboratory (ORNL) and Pyran™ Inc. to utilize modeling capabilities developed by the CCPC (Consortium for Computational Physics and Chemistry) to assist Pyran in scaling up its proprietary process for thermocatalytic conversion of furfural to 1,5 pentanediol (PDO). Pyran’s bio-based PDO is a direct replacement for petroleum-based PDO and 1,6 hexanediol (HDO) currently used in polymers, additives, coatings, adhesives and sealants. The current market is in excess of $\$$1B/yr. According to Pyran, their production process results in 95+% reduction in fossil CO 2 at lower production costs compared to the existing petroleum-based routes for PDO & HDO. Other chemicals based on intermediates from this process have current markets in excess of $\$$10B/yr.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Discovering Active Subspaces for High-Dimensional Computer Models

Dimension reduction techniques have long been an important topic in statistics, and active subspaces (AS) have received much attention this past decade in the computer experiments literature. The most common approach towards estimating the AS is to use Monte Carlo with numerical gradient evaluation. While sensible in some settings, this approach has obvious drawbacks. Recent research has demonstrated that active subspace calculations can be obtained in closed form, conditional on a Gaussian process (GP) surrogate, which can be limiting in high-dimensional settings for computational reasons. In this paper, we produce the relevant calculations for a more general case when the model of interest is a linear combination of tensor products. These general equations can be applied to the GP, recovering previous results as a special case, or applied to the models constructed by other regression techniques including multivariate adaptive regression splines (MARS). Furthermore, using a MARS surrogate has many advantages including improved scaling, better estimation of active subspaces in high dimensions and the ability to handle a large number of prior distributions in closed form. In one real-world example, we obtain the active subspace of a radiation-transport code with 240 inputs and 9,372 model runs in under half an hour.

97 MATHEMATICS AND COMPUTING↗