Progress on PRONGHORN Application to NGNP Related Problems

We are developing a multiphysics simulation tool for Very High-Temperature gascooled Reactors (VHTR). The simulation tool, PRONGHORN, takes advantages of the Multiphysics Object-Oriented Simulation library, and is capable of solving multidimensional thermal-fluid and neutronics problems implicitly in parallel. Expensive Jacobian matrix formation is alleviated by the Jacobian-free Newton-Krylov method, and physics-based preconditioning is applied to improve the convergence. The initial development of PRONGHORN has been focused on the pebble bed corec concept. However, extensions required to simulate prismatic cores are underway. In this progress report we highlight progress on application of PRONGHORN to PBMR400 benchmark problems, extension and application of PRONGHORN to prismatic core reactors, and progress on simulations of 3-D transients

A Numerical Model for Coupling of Neutron Diffusion and Thermomechanics in Fast Burst Reactors

We develop a numerical model for coupling of neutron diffusion adn termomechanics in order to stimulate transient behavior of a fast burst reactor. The problem involves solving a set of non-linear different equations which approximate neutron diffusion, temperature change, and material behavior. With this equation set we will model the transition from a supercritical to subcritical state and possible mechanical vibration

Jacobian Free-Newton Krylov Discontinuous Galerkin Method and Physics-Based Preconditioning for Nuclear Reactor Simulations

Multidimensional, higher-order (2nd and higher) numerical methods have come to the forefront in recent years due to significant advances of computer technology and numerical algorithms, and have shown great potential as viable design tools for realistic applications. To achieve this goal, implicit high-order accurate coupling of the multiphysics simulations is a critical component. One of the issues that arise from multiphysics simulation is the necessity to resolve multiple time scales. For example, the dynamical time scales of neutron kinetics, fluid dynamics and heat conduction significantly differ (typically >1010 magnitude), with the dominant (fastest) physical mode also changing during the course of transient [Pope and Mousseau, 2007]. This leads to the severe time step restriction for stability in traditional multiphysics (i.e. operator split, semi-implicit discretization) simulations. The lower order methods suffer from an undesirable numerical dissipation. Thus implicit, higher order accurate scheme is necessary to perform seamlessly-coupled multiphysics simulations that can be used to analyze the “what-if” regulatory accident scenarios, or to design and optimize engineering systems

Comparative Analysis of Natural Convection Flows Simulated by both the Conservation and Incompressible Forms of the Navier-Stokes Equations in a Differentially-Heated Square Cavity

This report illustrates a comparative study to analyze the physical differences between numerical simulations obtained with both the conservation and incompressible forms of the Navier-Stokes equations for natural convection flows in simple geometries. The purpose of this study is to quantify how the incompressible flow assumption (which is based upon constant density advection, divergence-free flow, and the Boussinesq gravitational body force approximation) differs from the conservation form (which only assumes that the fluid is a continuum) when solving flows driven by gravity acting upon density variations resulting from local temperature gradients. Driving this study is the common use of the incompressible flow assumption in fluid flow simulations for nuclear power applications in natural convection flows subjected to a high heat flux (large temperature differences). A series of simulations were conducted on two-dimensional, differentially-heated rectangular geometries and modeled with both hydrodynamic formulations. From these simulations, the selected characterization parameters of maximum Nusselt number, average Nusselt number, and normalized pressure reduction were calculated. Comparisons of these parameters were made with available benchmark solutions for air with the ideal gas assumption at both low and high heat fluxes. Additionally, we generated body force, velocity, and divergence of velocity distributions to provide a basis for further analysis. The simulations and analysis were then extended to include helium at the Very High Temperature gas-cooled Reactor (VHTR) normal operating conditions. Our results show that the consequences of incorporating the incompressible flow assumption in high heat flux situations may lead to unrepresentative results. The results question the use of the incompressible flow assumption for simulating fluid flow in an operating nuclear reactor, where large temperature variations are present. The results show that the use of the incompressible flow assumption with the Boussinesq gravitational body force approximation should be restricted to flows where the density change of a fluid particle along a pathline is negligible

Comparison of a New Length Measurement Technique for Inflatable Penile Prosthesis Implantation to Standard Techniques: Outcomes and Patient Satisfaction

Introduction.  Within a study evaluating the redesigned AMS 700MS inflatable penile prosthesis (IPP) (American Medical Systems, Minnetonka, MN, USA), one site used new length measurement technique (NLMT), a more aggressive dilation and measurement of the corpora cavernosa on a stretched penis, to address penile shortening. Aim.  To compare cylinder size and patient satisfaction, between a NLMT and traditional sizing for IPP implantation. Methods.  Fourteen men received IPPs using NLMT, and 55 with traditional sizing. Nationwide sales data from 2005 to 2008 for AMS 700 IPPs was obtained from AMS for comparison; additional surveys captured patient satisfaction. Main Outcome Measure.  Demographic data, cylinder sizes, and patient satisfaction were compared between the NLMT and standard techniques. Results.  The Fisher's exact test ( P  21 cm long and 28.6% (4) received cylinders <21 cm long, as compared with 12.7% (7) and 87.3% (48), respectively, for patients implanted by traditional techniques. There were ethnic differences between the samples: 42.9% (6) NLMT patients were of African‐American descent, as compared with 10.9% (6) in the standard technique group. However, longer cylinders were utilized more often, with 83.3% (5) of African‐Americans treated using the NLMT; as compared with 33.3% (2) of the standard technique group. Nationwide data reveal 12.3% of patients routinely receive 21 cm cylinders. At 6 months postimplantation, patient satisfaction with NLMT was no different than standard techniques. There were no distal erosions, complications, infections, or pain concerns reported through 24 months among the NLMT patients. Conclusions.  The NLMT resulted in a larger number of subjects implanted with larger cylinders. Satisfaction with performance and complication rates for NLMT patients was comparable to those implanted using standard techniques. Henry G, Houghton L, Culkin D, Otheguy J, Shabsigh R, and Ohl DA. Comparison of a new length measurement technique for inflatable penile prosthesis implantation to standard techniques: Outcomes and patient satisfaction. J Sex Med 2011;8:2640–2646.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/86820/1/j.1743-6109.2011.02340.x.pd

Hyperoxemia and excess oxygen use in early acute respiratory distress syndrome : Insights from the LUNG SAFE study

Publisher Copyright: © 2020 The Author(s). Copyright: Copyright 2020 Elsevier B.V., All rights reserved.Background: Concerns exist regarding the prevalence and impact of unnecessary oxygen use in patients with acute respiratory distress syndrome (ARDS). We examined this issue in patients with ARDS enrolled in the Large observational study to UNderstand the Global impact of Severe Acute respiratory FailurE (LUNG SAFE) study. Methods: In this secondary analysis of the LUNG SAFE study, we wished to determine the prevalence and the outcomes associated with hyperoxemia on day 1, sustained hyperoxemia, and excessive oxygen use in patients with early ARDS. Patients who fulfilled criteria of ARDS on day 1 and day 2 of acute hypoxemic respiratory failure were categorized based on the presence of hyperoxemia (PaO2 > 100 mmHg) on day 1, sustained (i.e., present on day 1 and day 2) hyperoxemia, or excessive oxygen use (FIO2 ≥ 0.60 during hyperoxemia). Results: Of 2005 patients that met the inclusion criteria, 131 (6.5%) were hypoxemic (PaO2 < 55 mmHg), 607 (30%) had hyperoxemia on day 1, and 250 (12%) had sustained hyperoxemia. Excess FIO2 use occurred in 400 (66%) out of 607 patients with hyperoxemia. Excess FIO2 use decreased from day 1 to day 2 of ARDS, with most hyperoxemic patients on day 2 receiving relatively low FIO2. Multivariate analyses found no independent relationship between day 1 hyperoxemia, sustained hyperoxemia, or excess FIO2 use and adverse clinical outcomes. Mortality was 42% in patients with excess FIO2 use, compared to 39% in a propensity-matched sample of normoxemic (PaO2 55-100 mmHg) patients (P = 0.47). Conclusions: Hyperoxemia and excess oxygen use are both prevalent in early ARDS but are most often non-sustained. No relationship was found between hyperoxemia or excessive oxygen use and patient outcome in this cohort. Trial registration: LUNG-SAFE is registered with ClinicalTrials.gov, NCT02010073publishersversionPeer reviewe

Changes in Invasive Pneumococcal Disease Caused by Streptococcus pneumoniae Serotype 1 following Introduction of PCV10 and PCV13: Findings from the PSERENADE Project

Streptococcus pneumoniae serotype 1 (ST1) was an important cause of invasive pneumococcal disease (IPD) globally before the introduction of pneumococcal conjugate vaccines (PCVs) containing ST1 antigen. The Pneumococcal Serotype Replacement and Distribution Estimation (PSERENADE) project gathered ST1 IPD surveillance data from sites globally and aimed to estimate PCV10/13 impact on ST1 IPD incidence. We estimated ST1 IPD incidence rate ratios (IRRs) comparing the pre-PCV10/13 period to each post-PCV10/13 year by site using a Bayesian multi-level, mixed-effects Poisson regression and all-site IRRs using a linear mixed-effects regression (N = 45 sites). Following PCV10/13 introduction, the incidence rate (IR) of ST1 IPD declined among all ages. After six years of PCV10/13 use, the all-site IRR was 0.05 (95% credibility interval 0.04-0.06) for all ages, 0.05 (0.04-0.05) for <5 years of age, 0.08 (0.06-0.09) for 5-17 years, 0.06 (0.05-0.08) for 18-49 years, 0.06 (0.05-0.07) for 50-64 years, and 0.05 (0.04-0.06) for ≥65 years. PCV10/13 use in infant immunization programs was followed by a 95% reduction in ST1 IPD in all ages after approximately 6 years. Limited data availability from the highest ST1 disease burden countries using a 3+0 schedule constrains generalizability and data from these settings are needed

Global, regional, and national disease burden estimates of acute lower respiratory infections due to respiratory syncytial virus in young children in 2015:a systematic review and modelling study

Background: We have previously estimated that respiratory syncytial virus (RSV) was associated with 22% of all episodes of (severe) acute lower respiratory infection (ALRI) resulting in 55 000 to 199 000 deaths in children younger than 5 years in 2005. In the past 5 years, major research activity on RSV has yielded substantial new data from developing countries. With a considerably expanded dataset from a large international collaboration, we aimed to estimate the global incidence, hospital admission rate, and mortality from RSV-ALRI episodes in young children in 2015. Methods: We estimated the incidence and hospital admission rate of RSV-associated ALRI (RSV-ALRI) in children younger than 5 years stratified by age and World Bank income regions from a systematic review of studies published between Jan 1, 1995, and Dec 31, 2016, and unpublished data from 76 high quality population-based studies. We estimated the RSV-ALRI incidence for 132 developing countries using a risk factor-based model and 2015 population estimates. We estimated the in-hospital RSV-ALRI mortality by combining in-hospital case fatality ratios with hospital admission estimates from hospital-based (published and unpublished) studies. We also estimated overall RSV-ALRI mortality by identifying studies reporting monthly data for ALRI mortality in the community and RSV activity. Findings: We estimated that globally in 2015, 33·1 million (uncertainty range [UR] 21·6–50·3) episodes of RSV-ALRI, resulted in about 3·2 million (2·7–3·8) hospital admissions, and 59 600 (48 000–74 500) in-hospital deaths in children younger than 5 years. In children younger than 6 months, 1·4 million (UR 1·2–1·7) hospital admissions, and 27 300 (UR 20 700–36 200) in-hospital deaths were due to RSV-ALRI. We also estimated that the overall RSV-ALRI mortality could be as high as 118 200 (UR 94 600–149 400). Incidence and mortality varied substantially from year to year in any given population. Interpretation: Globally, RSV is a common cause of childhood ALRI and a major cause of hospital admissions in young children, resulting in a substantial burden on health-care services. About 45% of hospital admissions and in-hospital deaths due to RSV-ALRI occur in children younger than 6 months. An effective maternal RSV vaccine or monoclonal antibody could have a substantial effect on disease burden in this age group

A second-order accurate in time IMplicit–EXplicit (IMEX) integration scheme for sea ice dynamics

Current sea ice models use numerical schemes based on a splitting in time between the momentum and continuity equations. Because the ice strength is explicit when solving the momentum equation, this can create unrealistic ice stress gradients when using a large time step. As a consequence, noise develops in the numerical solution and these models can even become numerically unstable at high resolution. To resolve this issue, we have implemented an iterated IMplicit–EXplicit (IMEX) time integration method. This IMEX method was developed in the framework of an already implemented Jacobian-free Newton–Krylov solver. The basic idea of this IMEX approach is to move the explicit calculation of the sea ice thickness and concentration inside the Newton loop such that these tracers evolve during the implicit integration. To obtain second-order accuracy in time, we have also modified the explicit time integration to a second-order Runge–Kutta approach and by introducing a second-order backward difference method for the implicit integration of the momentum equation. These modifications to the code are minor and straightforward. By comparing results with a reference solution obtained with a very small time step, it is shown that the approximate solution is second-order accurate in time. The new method permits to obtain the same accuracy as the splitting in time but by using a time step that is 10 times larger. Results show that the second-order scheme is more than five times more computationally efficient than the splitting in time approach for an equivalent level of error