Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vessel 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analysis of Postulated Accident Scenarios for the Transformational Challenge Reactor

The Transformational Challenge Reactor (TCR) is a helium-cooled, yttrium-hydride-moderated reactor that was designed for the U.S. Department of Energy Office of Nuclear Energy. A key objective of the TCR was to employ advanced manufacturing techniques in a nuclear system and demonstrate their potential for revolutionizing the nuclear reactor design process. One purpose of the present work is to demonstrate the safety of the TCR under postulated accidents. Based on RELAP5-3D and COMSOL analyses, the TCR remained below all current safety limits and far below the expected failure limits for the core materials. Another purpose of this work is to provide useful insights and recommendations regarding the application of RELAP5-3D to gas-cooled or other advanced reactors. A novel approach was implemented for simultaneously modeling conduction and radiation in RELAP5-3D, which was found to provide reasonable predictions of radial core, vessel, and ex-vessel heat transfer during postulated events. A multicode approach was also applied, in which high-fidelity COMSOL calculations were used to tune the radial heat transfer parameters in RELAP5-3D. The tuned RELAP5-3D model demonstrated comparable peak temperature predictions as COMSOL, despite a coarse treatment of the core in RELAP5-3D consisting of only two lumped heat structures. This high-fidelity tuning approach enabled enhanced accuracy as well as minimal complexity within the RELAP5-3D model, even for complex fuel geometric designs as in the TCR. Finally, investigations were made into the potential for flow reversal during a pressurized loss-of-forced-flow event in the TCR. The TCR is designed with downward helium flow through the core during normal operation. The RELAP5-3D model predicted that this downward flow would persist, without flow reversal, up to several days after the circulator trip. This was attributed to natural circulation hysteresis effects as have been noted in similar thermofluidic systems. Although flow stagnation and eventual reversal did not lead to unsafe TCR conditions, interesting spatial effects were observed which may have safety relevance for other reactor system designs and coolant types that are designed for downward core flow during normal operation, warranting closer investigation of the flow reversal phenomenon.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Exposed Phosphatidylserine as a Biomarker for Clear Identification of Breast Cancer Brain Metastases in Mouse Models

Brain metastasis is the most common intracranial malignancy in adults. The prognosis is extremely poor, partly because most patients have more than one brain lesion, and the currently available therapies are nonspecific or inaccessible to those occult metastases due to an impermeable blood–tumor barrier (BTB). Phosphatidylserine (PS) is externalized on the surface of viable endothelial cells (ECs) in tumor blood vessels. In this study, we have applied a PS-targeting antibody to assess brain metastases in mouse models. Fluorescence microscopic imaging revealed that extensive PS exposure was found exclusively on vascular ECs of brain metastases. The highly sensitive and specific binding of the PS antibody enables individual metastases, even micrometastases containing an intact BTB, to be clearly delineated. Furthermore, the conjugation of the PS antibody with a fluorescence dye, IRDye 800CW, or a radioisotope, 125I, allowed the clear visualization of individual brain metastases by optical imaging and autoradiography, respectively. In conclusion, we demonstrated a novel strategy for targeting brain metastases based on our finding that abundant PS exposure occurs on blood vessels of brain metastases but not on normal brain, which may be useful for the development of imaging and targeted therapeutics for brain metastases.

Oncology↗
DOE OSTI2022

Quantifying the reduction in cavitation-induced erosion damage in the Spallation Neutron Source mercury target by means of small-bubble gas injection

Here, a model developed to represent the progress of erosion damage in liquid-metal spallation target vessels was modified to incorporate the effect of gas injection on the erosion rate. The liquid mercury target system for the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory now operates with helium gas injection to reduce target vessel fatigue stress and cavitation-induced erosion damage. Erosion damage is a primary degradation phenomenon affecting the service life of SNS target vessels, and cavitation mitigation techniques, such as small-bubble gas injection, have been implemented to reduce damage and extend target lifetimes. Erosion depths in samples removed from SNS targets after operation were measured using laser line scanning. These measurements confirmed that gas injection reduced erosion damage. However, quantifying the damage reduction due to gas injection was complicated by variations in lifetime, power, and gas injection rates between different targets. In this study, the operating power and gas injection rate of targets were incorporated into an erosion damage prediction model to quantify their effects on erosion damage reduction. Values of a power scaling factor, β, were calculated by comparing modeled with measured erosion damage. These values indicate that the use of gas injection at the SNS reduced damage to a level equivalent to operating targets without gas injection at 35–47% of the actual beam power. To account for the gas injection effect on the cavitation damage, a simple exponential form based on analysis of the scaling factor β was developed to incorporate the gas rate history with a scaling factor γ in the erosion damage modeling.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Non-invasive characterization of complex coronary lesions

Conventional invasive diagnostic imaging techniques do not adequately resolve complex Type B and C coronary lesions, which present unique challenges, require personalized treatment and result in worsened patient outcomes. These lesions are often excluded from large-scale non-invasive clinical trials and there does not exist a validated approach to characterize hemodynamic quantities and guide percutaneous intervention for such lesions. This work identifies key biomarkers that differentiate complex Type B and C lesions from simple Type A lesions by introducing and validating a coronary angiography-based computational fluid dynamic (CFD-CA) framework for intracoronary assessment in complex lesions at ultrahigh resolution. Among 14 patients selected in this study, 7 patients with Type B and C lesions were included in the complex lesion group including ostial, bifurcation, serial lesions and lesion where flow was supplied by collateral bed. Simple lesion group included 7 patients with lesions that were discrete, <10mm long and readily accessible. Intracoronary assessment was performed using CFD-CA framework and validated by comparing to clinically measured pressure-based index, such as FFR. Local pressure, endothelial shear stress (ESS) and velocity profiles were derived for all patients. We validates the accuracy of our CFD-CA framework and report excellent agreement with invasive measurements (n=14,R 2 =0.6,p=0.0013). Ultra-high resolution achieved by the model enable physiological assessment in complex lesions and quantify hemodynamic metrics in all vessels up to 1mm in diameter. Importantly, we demonstrate that in contrast to traditional pressure-based metrics, there is a significant difference in the intracoronary hemodynamic forces, such as ESS, in complex lesions compared to simple lesions at both resting and hyperemic physiological states [n = 14, p=0.03]. Higher ESS was observed in the complex lesion group (7.0±4.7 Pa) than in simple lesion group (4.8±3.6 Pa). Complex coronary lesions have higher ESS compared to simple lesions, such differential hemodynamic evaluation can provide much the needed insight into the increase in adverse outcomes for such patients and has incremental prognostic value over traditional pressure-based indices, such as FFR.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Initial development of a high temperature life prediction method directly accounting for variability in material properties

This report presents the initial development of a life prediction method that accounts for the variability of the material properties in Grade 91 steel. We account for material variability by both fitting a variable 3-parameter Weibull distribution to the rupture data of the material and by quantifying the effect that the creep deformation properties have on stresses via a Monte Carlo approach. To ensure a reasonable time frame when applying Monte Carlo to the creep stress, we propose a methodology that models the material as an extremely viscous Stokes flow with a non-Newtonian viscosity, therefore solving the problem in a single step. We use the final probabilistic model for structural failure developed in this document to evaluate a Flat Head Vessel designed with the Section III, Division 5 ASME Boiler Pressure and Vessel Code rules with a 100,000h design life and assess its probability of premature failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

R-CERT Project Analysis Report

The R-CERT project aims to determine if PBX 9501 will experience a deflagration to detonation transition (DDT) while in a high pressure environment. Combined efforts in design and analysis between multiple groups at Los Alamos National Laboratory (LANL) were required to finalize the test design. The proposed design of the case assembly (i.e. the pressure vessel that would contain the high explosive) was analyzed and iteratively modified until modeling efforts predicted the case would maintain pressure up to 5000 psi under quasi-static conditions. To prevent the case from being excessively robust beyond its design requirement, grooves were added to the case geometry with depths and radii of 4 mm, thus weakening the case while maintaining much of the mass of the proposed design.

42 ENGINEERING↗
DOE OSTI2020

Update on Full-Scale Retrieval and Transfer System Testing for Pneumatic Conveyance of High-Level Calcine Waste - 20241

The Idaho National Laboratory processed spent nuclear fuel from 1953 to 1994. The process used solvent extraction to separate uranium, resulting in a liquid waste. A calcination process was used to convert the liquid waste to a high-level radioactive solid powder, referred to as calcine. The calcine is stored in stainless steel bins contained within concrete vaults called the Calcined Solids Storage Facilities. There are a total of seven CSSFs. The Calcine Retrieval Project is tasked with retrieving the calcine from the bins of CSSF 1 and pneumatically transferring it to the bins of CSSF 6. CSSF 1 was constructed to be a final calcine repository with no forethought of retrieval. It represents the greatest technical hurdle to calcine extraction of all the CSSFs. The challenges inherent with accessing the calcine, safely extracting it, and securely transporting it are copious and complex. Once the calcine is removed from CSSF1, that facility may be closed, reducing the footprint of waste storage at the INL. The same process for removing Calcine from CSSF1 may then be utilized to remove the calcine from the remaining CSSFs for final treatment and transport to a long-term storage facility. To overcome these inherent challenges with accessing CSSF 1 and extracting the calcine in a timely manner, a full-scale integrated mockup has been constructed. The mockup includes a full-size model of a concentric group of bins from CSSF 1 and a storage vessel representing CSSF 6. All essential elements of the transport design are represented in the Mockup, including motive and makeup air supply (air compressor), a replica of the CSSF6 cyclone, pre-filter, HEPA filter, vacuum pump, and retrieval control devices. Construction on the Mockup was completed in 2018 and testing began in February 2019. It has already proven very beneficial to the CRP. Experience in running the Mockup have aided in the development of retrieval and transport methodologies and the selection of process parameter values. It has also uncovered some unforeseen technological challenges and will continue to be valuable to the testing team as these hurdles are overcome and the process is proven effective. Operations personnel have also benefitted from its use as a training tool and will continue to use the Mockup as a simulator to reproduce problems encountered during extraction and transfer of Calcine from CSSF 1. The knowledge gained during this process will also be of continued benefit as future work is done to remove various types of calcine from the remaining CSSFs. The methodology of safe and secure transfer of radioactive material may be of benefit to others throughout the U.S. Department of Energy complex. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Characterizing the flux effect on the irradiation embrittlement of reactor pressure vessel steels using machine learning

In-service exposure to high-energy neutrons embrittles reactor pressure vessel (RPV) steels. An increase in the yield stress (Δσ y ) results in a corresponding increase in the brittle to ductile transition temperature (ΔT c ). Most existing models underpredict ΔT c at higher fluence following accelerated irradiations in test reactors. High fluence, up to 10 20 n/cm 2 in some cases, will be reached over extended RPV vessel operation of 80 years, or more, at low service flux. Embrittlement has been extensively studied in accelerated, higher flux test reaction irradiations. However, the use of test reactor data naturally raises the question of flux effects. Here, this study used a machine learning approach trained on a set of hardening data, covering a wide range of flux, fluence, and steel compositions to determine the interactive effects of both irradiation and material variables on Δσ y . The analysis included machine learning-based cross-plots of the variable dependence of Δσ y for six core steels (i.e., CM6, LC, LD, LG, LH, and LI), with controlled differences in their Cu and Ni contents. A primary objective is to evaluate an effective to actual fluence (Φt e /Φt) ratio, as a function of flux, fluence, and steel composition. This is information critical to properly use intermediate flux-high fluence data in calibrating a low flux-high fluence embrittlement model. The predicted Φt e /Φt is reasonably consistent with estimates previously derived from a physics-based solute recombination trap model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A heat and mass transfer model for evaluation of damaged cesium chloride radiological sources

Radiological sources are vital to many applications, and typically contain byproducts from waste streams which decay over time to form complex mixtures of elements. Daughter products resulting from these decay processes can introduce complicating factors to the integrity and safety of the vessel that contains the radiological material. Sources containing cesium 137 ( 137 Cs) are of particular interest, because the decay of this isotope produces barium metal which reacts readily and exothermically with oxygen. This work employs physics-based numerical tools to examine the thermal response of a radiological source containing a mixture of cesium and barium in the event that vessel walls are damaged and atmospheric gases contact source material. A parametric study was conducted to determine the sensitivity of the response to various factors, including vessel geometry, source material age, the degree of vessel damage, and other parameters. Here, it was found that the peak temperatures that occurred within the source material strongly depend on these parameters, particularly vessel geometry and age, which determine whether a breach would be a relatively minor accident or a catastrophic incident. Finally, the model’s sensitivity to uncertain thermophysical properties is discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10 –2 and a vacuumisland overlapwidth of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. Furthermore, these findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

3D coil↗
DOE OSTI2021

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

VERA User's Guide for Ex-core Applications

The Virtual Environment for Reactor Applications, or VERA, allows users to set up models to calculate time-dependent and fully coupled solutions for ex-core quantities of interest such as vessel and coupon fluence, and detector responses for multiple statepoints and cycles. MPACT and COBRA-TF together perform in-core transport calculations with temperature feedback while Shift performs the fluence and detector response calculations in the ex-core region. The in-core region is modeled using VERA’s native input format and the ex-core region is defined using Shift’s general geometry package, also known as Omnibus General Geometry. Fixed source ex-core calculations with Shift can be run in forward mode without advanced variance reduction (VR) methods, or with Consistent Adjoint Driven Importance Sampling (CADIS), which is an automated VR method. This document serves as a guide for setting up inputs, running ex-core calculations and post-processing the results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Spatiotemporal forecasting of the edge localized modes in tokamak plasmas using neural networks

Artificial intelligence techniques have been increasingly adopted by the plasma and fusion science to address problems like plasma reconstruction, surrogate modeling, and tokamak/stellarator optimization. A key focus in sustained fusion research is the prediction and mitigation of edge-localized-modes (ELMs), instabilities that occur in short, periodic bursts and can cause erosion to the tokamak vessel wall. Recent research has demonstrated the power of neural networks in approximating continuous functions. In this work, we build spatiotemporal forecasting models that can predict the onset of ELMs and their evolution at early stages. We leverage recent advances in generative modeling, sequence-to-sequence modeling, and Fourier neural operators to propose architectures and training strategies that can learn to forecast short to long term dynamics of the noisy signals due to ELMs. We benchmark the developed model against a state-of-the-art foundation model using the beam emission spectroscopy (BES) data that captures the plasma fluctuations due to ELMs over a 8 x 8 spatial grid. Our models demonstrate high accuracy, outperforming the baselines, in predicting the evolution of BES signals during ELM events. Furthermore, the developed models exhibit high accuracy in predicting the rapid rise and relaxation of the signals due to ELMs within 30–80 µs.

edge localized modes↗
DOE OSTI2022

Neural net modeling of equilibria in NSTX-U

Neural networks (NNs) offer a path towards synthesizing and interpreting data on faster timescales than traditional physics-informed computational models. In this work we develop two NNs relevant to equilibrium and shape control modeling, which are part of a suite of tools being developed for the National Spherical Torus Experiment-Upgrade for fast prediction, optimization, and visualization of plasma scenarios. The networks include Eqnet, a free-boundary equilibrium solver trained on the EFIT01 (Equilibrium FITtting 01) reconstruction algorithm, and Pertnet, which is trained on the Gspert code and predicts the non-rigid plasma response, a nonlinear term that arises in shape control modeling. The NNs are trained with different combinations of inputs and outputs in order to offer flexibility in use cases. In particular, Eqnet can use magnetic diagnostics as inputs and act as an EFIT-like reconstruction algorithm, or, by using pressure and current profile information the NN can act as a forward Grad–Shafranov equilibrium solver. This forward-mode version is envisioned to be implemented in the suite of tools for simulation of plasma scenarios. The reconstruction-mode version gives some performance improvements compared to the online reconstruction code real-time EFIT, especially when vessel eddy currents are significant. Here, we report strong performance for all NNs indicating that the models could reliably be used within closed-loop simulations or other applications. Some limitations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Assessment of simulated and observed cavitation-induced erosion damage in Spallation Neutron Source target vessels

Cavitation-induced erosion damage in different Spallation Neutron Source (SNS) target designs are simulated using explicit finite element–based techniques and compared with observations of erosion in targets after operation. The efficacy of the previously developed method, called saturation time, was evaluated using erosion-damaged samples from new target designs. A new metric called maximum bubble size was implemented under the rationale that larger cavitation bubbles will collapse more intensely. The maximum cavitation bubble size over 1 ms of simulated time was calculated based on the Rayleigh–Plesset equation for each element integration point and presented as a contour map at the vessel surface for assessing with erosion observations. SNS targets are now operated with helium gas injection to reduce cavitation damage. A simulation method using a material model for the mixture of mercury and gas bubbles was recently developed and used to account for the effect of small gas bubbles on the structural response of the target vessel. Furthermore, this work compares the new method's results with observed cavitation damage. Maps of the calculated maximum bubble size for targets operated with and without gas injection were compared with photographs of erosion damage observed in SNS targets. The patterns in maximum bubble size maps correlated well with observations of erosion patterns in target vessels after service. Advantages and challenges of the maximum bubble size simulation technique are provided, and differences between results from the previous and the newly proposed metric are discussed.

Jiang, Hao↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗