Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional rate functions”

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 1,477 records · Page 82

Entrapped Single Tungstate Site in Zeolite for Cooperative Catalysis of Olefin Metathesis with Brønsted Acid Site

Industrial olefin metathesis catalysts generally suffer from low reaction rates and require harsh reaction conditions for moderate activities. This is due to their inability to prevent metathesis active sites (MAS) from aggregation and their intrinsic poor adsorption and activation of olefin molecules. Here, isolated tungstate species as single molecular MAS is immobilized inside zeolite pores by Bronsted acid sites (BAS) on the inner surface. It is demonstrated for the first time that unoccupied BAS in atomic proximity to MAS enhance olefin adsorption and greatly facilitate the formation of metallocycle intermediates in a stereospecific manner. Thus, effective cooperative catalysis takes place over the BAS-MAS pair. In consequence, for the cross-metathesis of ethene and trans-2-butene to propene, under the same mild reaction conditions, the propene production rate over WOx/USY is ca. 7,300 times that over the industrial WO3/SiO2 based catalyst. A propene yield up to 79% (80% selectivity) without observable deactivation was obtained over WOx/USY for a wide range of reaction conditions.

Brønsted acid sites↗

Comparison of the Arrhenius parameters between conventional hydrothermal and microwave-assisted synthesis methods for tin oxide nanoparticles

Microwave (MW) irradiation has emerged as a powerful tool for accelerating materials synthesis, yet the origins of its specific influence on reaction kinetics remain elusive. While multiple studies have attributed the observed enhancements in reaction rates under MW heating to reduced activation energies, other accounts have suggested modifications to the Arrhenius pre-exponential factor as the predominant cause. Distinguishing between these parameters in modern applications of MW processing in nanomaterials requires experimental approaches capable of resolving the dynamic and nuanced structural kinetics that govern MW-assisted chemistry. Here, we combine in-situ synchrotron X-ray total scattering with pair distribution function (PDF) analysis to track the structural evolution of SnO 2 nanoparticles synthesized via MW-assisted and conventional hydrothermal conditions. Avrami modeling and Arrhenius analysis suggest that although MW irradiation yields a higher apparent activation energy, the enhanced crystallization is better explained by a pre-exponential factor several orders of magnitude larger than that of conventional heating. These findings suggest that the MW field induces a higher frequency of successful molecular rearrangements rather than lowering the intrinsic activation barrier. The results contribute further insights clarifying the role of the applied MW field for materials design where MW-specific effects can be deliberately harnessed. Furthermore, this work presents a framework promoting the utility of in-situ PDF characterization coupled with kinetic analysis for developing more sophisticated descriptions of nanoscale transformations for MW-driven reaction kinetics.

36 MATERIALS SCIENCE↗

Black hole superradiance with dark matter accretion

Studies of black hole superradiance often focus on the growth of a cloud in isolation, accompanied by the spin-down of the black hole. Here, in this paper, we consider the additional effect of the accretion of matter and angular momentum from the environment. We show that, in many cases, the black hole evolves by drifting along the superradiance threshold, in which case the evolution of its parameters can be described analytically or semianalytically. We quantify the conditions under which accretion can serve as a mechanism to increase the cloud-to-black hole mass ratio, beyond the standard maximum of about 10%. This occurs by a process we call oversuperradiance, whereby accretion effectively feeds the superradiance cloud, by way of the black hole. We give two explicit examples: accretion from a vortex expected in wave dark matter and accretion from a baryonic disk. In the former case, we estimate the accretion rate by using an analytical fit to the asymptotic behavior of the confluent Heun function. Level transition, whereby one cloud level grows while the other shrinks, can be understood in a similar way.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulation of PV Variability as a Function of PV Generation and Plant Size

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear-and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute time-scale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

irradiance↗

Relative humidity threshold for oxygen generation by high-purity plutonium dioxides in 3013 containers

High-purity plutonium dioxide with adsorbed water in a sealed container can generate oxygen, but oxygen will not be generated in low water loading conditions. A threshold for oxygen generation has been observed in the hydrogen generation rate normalized by the specific power for high-purity plutonium dioxide. Normalized hydrogen generation rates above this threshold generate oxygen and below the threshold consume oxygen. The threshold rate is determined by the amount of water adsorbed onto the surface of the oxide. The relative humidity of the atmosphere that the oxide is in equilibrium with determines the amount of physisorbed water adsorbed on the surface. We calculate the relative humidity corresponding to the threshold rate as a function of mass and specific surface area of material stored in 3013 containers. Exceeding the calculated relative humidity may result in flammable gas mixtures of hydrogen and oxygen at some time during storage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Surface-slip equations for multicomponent, nonequilibrium air flow

Equations are presented for the surface slip (or jump) values of species concentration, pressure, velocity, and temperature in the low-Reynolds-number, high-altitude flight regime of a space vehicle. These are obtained from closed-form solutions of the mass, momentum, and energy flux equations using the Chapman-Enskog velocity distribution function. This function represents a solution of the Boltzmann equation in the Navier-Stokes approximation. The analysis, obtained for nonequilibrium multicomponent air flow, includes the finite-rate surface catalytic recombination and changes in the internal energy during reflection from the surface. Expressions for the various slip quantities have been obtained in a form which can readily be employed in flow-field computations. A consistent set of equations is provided for multicomponent, binary, and single species mixtures. Expression is also provided for the finite-rate species-concentration boundary condition for a multicomponent mixture in absence of slip.

Gupta, Roop N.↗

Lyapunov function gradient generated robust control in the absence of the nominal stabilizing control

A relevant facet of the application of Lyapunov gradient-generated robust control to unstable linear autonomous plants is explored. It is demonstrated that if the plant, the output, and the nominal stabilizing control satisfy certain conditions, then the robust component alone stabilizes the nominal plant. An example characterized by two zero eigenvalues and two negative real value poles is presented. These results assure that the robust component will fulfill the role of nominal stabilization successfully so long as the possible magnitude of the robust component can overcome the contribution of the instability to positiveness of the Lyapunov rate.

Blackwell, C. C.↗

Application of a Reynolds Stress turbulence model to a supersonic hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows undergoing hydrogen-air chemical reactions. An assumed Beta Probability Density Function is applied to account for the chemical source terms in the conservation equations. An algebraic Reynolds Flux model is used for the fluctuating density-velocity as well as the species mass fraction-velocity correlations. The variances of temperature and species fluctuations are also modelled using an algebraic flux technique. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion processes. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant.

Chandrasekhar, R.↗

Head perturbations during walking while viewing a head-fixed target

Inexpensive, head-fixed computer displays are now available that subjects can wear during locomotion. Our hypothesis is that viewing a head-fixed visual display will change the character- istics of rotational head perturbations during natural walking. Using a 3-axis angular rate sensor, we measured head rotations during natural or treadmill walking, in 10 normal subjects and 2 patients with deficient vestibular function, as they attempted to view (1) a stationary target at optical infinity; and (2) a target at a distance of 20 cm rigidly attached to the head. Normal subjects and patients showed no significant change in the predominant frequency of head rotations in any plane (ranging 0.7-5.7 Hz) during the two different viewing tasks. Mean peak head velocities also showed no difference during the two viewing conditions except in the yaw plane, in which values were greater while viewing the near target. Predominant frequencies of head rotations were similar in the pitch plane during natural or treadmill walking; however, peak velocities of pitch head rotations were substantially greater during natural walking. One vestibular patient showed modest increases of head velocity during natural walking compared with normal subjects. Rotational head perturbations that occur during natural walking are largely unaffected when subjects view a head-fixed target. There is need to study how such perturbations, which induce vestibular eye movements, affect vision of head-fixed displays.

Das, Vallabh E.↗

Pattern Formation in Diffusion Flames Embedded in von Karman Swirling Flows

Pattern formation is observed in nature in many so-called excitable systems that can support wave propagation. It is well-known in the field of combustion that premixed flames can exhibit patterns through differential diffusion mechanism between heat and mass. However, in the case of diffusion flames where fuel and oxidizer are separated initially there have been only a few observations of pattern formation. It is generally perceived that since diffusion flames do not possess an inherent propagation speed they are static and do not form patterns. But in diffusion flames close to their extinction local quenching can occur and produce flame edges which can propagate along stoichiometric surfaces. Recently, we reported experimental observations of rotating spiral flame edges during near-limit combustion of a downward-facing polymethylmethacrylate disk spinning in quiescent air. These spiral flames, though short-lived, exhibited many similarities to patterns commonly found in quiescent excitable media including compound tip meandering motion. Flame disks that grow or shrink with time depending on the rotational speed and in-depth heat loss history of the fuel disk have also been reported. One of the limitations of studying flame patterns with solid fuels is that steady-state conditions cannot be achieved in air at normal atmospheric pressure for experimentally reasonable fuel thickness. As a means to reproduce the flame patterns observed earlier with solid fuels, but under steady-state conditions, we have designed and built a rotating, porous-disk burner through which gaseous fuels can be injected and burned as diffusion flames. The rotating porous disk generates a flow of air toward the disk by a viscous pumping action, generating what is called the von K rm n boundary layer which is of constant thickness over the entire burner disk. In this note we present a map of the various dynamic flame patterns observed during the combustion of methane in air as a function of fuel flow rate and the burner rotational speed.

Nayagam, Vedha↗

Effects of Transparency on Pilot Trust and Agreement in the Autonomous Constrained Flight Planner

We performed a human-in-the-loop study to explore the role of transparency in engendering trust and reliance within highly automated systems. Specifically, we examined how transparency impacts trust in and reliance upon the Autonomous Constrained Flight Planner (ACFP), a critical automated system being developed as part of NASA's Reduced Crew Operations (RCO) Concept. The ACFP is designed to provide an enhanced ground operator, termed a super dispatcher, with recommended diversions for aircraft when their primary destinations are unavailable. In the current study, 12 commercial transport rated pilots who played the role of super dispatchers were given six time-pressured all land scenarios where they needed to use the ACFP to determine diversions for multiple aircraft. Two factors were manipulated. The primary factor was level of transparency. In low transparency scenarios the pilots were given a recommended airport and runway, plus basic information about the weather conditions, the aircraft types, and the airport and runway characteristics at that and other airports. In moderate transparency scenarios the pilots were also given a risk evaluation for the recommended airport, and for the other airports if they requested it. In the high transparency scenario additional information including the reasoning for the risk evaluations was made available to the pilots. The secondary factor was level of risk, either high or low. For high-risk aircraft, all potential diversions were rated as highly risky, with the ACFP giving the best option for a bad situation. For low-risk aircraft the ACFP found only low-risk options for the pilot. Both subjective and objective measures were collected, including rated trust, whether the pilots checked the validity of the automation recommendation, and whether the pilots eventually flew to the recommended diversion airport. Key results show that: 1) Pilots trust increased with higher levels of transparency, 2) Pilots were more likely to verify ACFPs recommendations with low levels of transparency and when risk was high, 3) Pilots were more likely to explore other options from the ACFP in low transparency conditions and when risk was high, and 4) Pilots decision to accept or reject ACFPs recommendations increased as a function of the transparency in the explanation. The finding that higher levels of transparency was coupled with higher levels of trust, a lower need to verify other options, and higher levels of agreement with ACFP recommendations, confirms the importance of transparency in aiding reliance on automated recommendations. Additional analyses of qualitative data gathered from subjects through surveys and during debriefing interviews also provided the basis for new design recommendations for the ACFP.

trust↗

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution, electric propulsion system, and heat rejection system. Specific mass, or α (kg/kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass predictions for each individual component. To inform technology maturation planning activities, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively-cooled heat pipe reactor concepts. Each concept requires specific geometries, fluids, and power conversion interface conditions (temperature, pressure, flow rate) to meet desired performance and mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support concept down-selection along with related technology development planning. The reactor and shield α are a function of several PCS and heat rejection system design choices, and reactor scaling with these parameters must be considered to enable an informed decision on an optimal reactor geometry and working fluid combination.

Nuclear Electric Propulsion↗

A preliminary assessment of olivine phenocrysts from the monogenetic basalt of the McCartys Flow, Zuni-Bandera Volcanic Field, New Mexico

Monogenetic small-volume basaltic volcanoes are the most abundant subaerial volcanic landforms on Earth but are some of the most poorly understood systems. Their short durations, small volumes, and lack of recurrence make monitoring and hazard assessment difficult. The Zuni-Bandera volcanic field in western New Mexico contains small-volume basaltic centers erupting tholeiitic to alkalic basalts. Evidence shows no correlation of magma composition with eruption age, location, or volumetric output, prompting questions about the influence of magma ascent rates, magma storage conditions, and mantle source characteristics on lava compositions. Here, we present olivine major and minor element mineral chemistry from the 3200-year-old McCartys Flow, the youngest tholeiite basalt in the volcanic field. Olivine displays four phenocryst types with unique textures and major and minor element compositions. Multiple olivine types co-exist at the thin section scale. Major and minor element diffusion at frozen melt–phenocryst interfaces was modeled, revealing magma residence times ranging from 3–9 months. Type 3 olivine phenocrysts require step function initial conditions and record diffusion re-equilibration followed by magma mixing. These profiles indicate the magma resided in the reservoir for 10–15 years and accumulated from multiple batches of mixed magmas less than 10 days before the eruption. Our results show that primitive magmas in small-volume monogenetic volcanoes have complex lithospheric magmatic histories and stored in magma bodies influenced by an open system to develop different local chemical environments.

Gary S. Michelfelder↗

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗

Simulation of PV Variability as a Function of PV Generation and Plant Size: Preprint

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear- and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute timescale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

41 EE - Solar Energy Technologies Office (EE-4S)↗

A priori and a posteriori investigations for developing large eddy simulations of multi-species turbulent mixing under high-pressure conditions

A Direct Numerical Simulation (DNS) database was created representing mixing of species under high-pressure conditions. The configuration considered is that of a temporally evolving mixing layer. The database was examined and analyzed for the purpose of modeling some of the unclosed terms that appear in the Large Eddy Simulation (LES) equations. Several metrics are used to understand the LES modeling requirements. First, a statistical analysis of the DNS-database large-scale flow structures was performed to provide a metric for probing the accuracy of the proposed LES models as the flow fields obtained from accurate LESs should contain structures of morphology statistically similar to those observed in the filtered-and-coarsened DNS (FC-DNS) fields. To characterize the morphology of the large-scales structures, the Minkowski functionals of the iso-surfaces were evaluated for two different fields: the second-invariant of the rate of deformation tensor and the irreversible entropy production rate. To remove the presence of the small flow scales, both of these fields were computed using the FC-DNS solutions. It was found that the large-scale structures of the irreversible entropy production rate exhibit higher morphological complexity than those of the second invariant of the rate of deformation tensor, indicating that the burden of modeling will be on recovering the thermodynamic fields. Second, to evaluate the physical effects which must be modeled at the subfilter scale, an a priori analysis was conducted. This a priori analysis, conducted in the coarse-grid LES regime, revealed that standard closures for the filtered pressure, the filtered heat flux, and the filtered species mass fluxes, in which a filtered function of a variable is equal to the function of the filtered variable, may no longer be valid for the high-pressure flows considered in this study. The terms requiring modeling are the filtered pressure, the filtered heat flux, the filtered pressure work, and the filtered species mass fluxes. Improved models were developed based on a scale-similarity approach and were found to perform considerably better than the classical ones. These improved models were also assessed in an a posteriori study. Different combinations of the standard models and the improved ones were tested. At the relatively small Reynolds numbers achievable in DNS and at the relatively small filter widths used here, the standard models for the filtered pressure, the filtered heat flux, and the filtered species fluxes were found to yield accurate results for the morphology of the large-scale structures present in the flow. Analysis of the temporal evolution of several volume-averaged quantities representative of the mixing layer growth, and of the cross-stream variation of homogeneous-plane averages and second-order correlations, as well as of visualizations, indicated that the models performed equivalently for the conditions of the simulations. The expectation is that at the much larger Reynolds numbers and much larger filter widths used in practical applications, the improved models will have much more accurate performance than the standard one.

42 ENGINEERING↗

Estimation of the time for steam generator trip due to cyber intrusions

The time required to trip a pressurized water reactor (PWR) by inserting malicious signals into its steam generator (SG) control system has been studied using the Generic PWR (GPWR) Simulator. A semi-analytical model is developed to approximately reproduce the simulator response and understand the dynamics of the control unit. A series of two proportional-integral controllers determines control action according to preset constants, the readings from the feedwater level sensor, and those from feedwater and steam flowrate transmitters. It is observed that the most important factor that determines whether a trip will occur is how much additional water is added to or withheld from the SG over time compared to normal operating conditions. In order to determine the effects of control action on the SG, changes in mass inventory are considered. This approach models the SG water level as a function of mass inventory and has a backward temporal memory. A Python interface is developed for the GPWR framework to automatically simulate different spoofing scenarios and post-process the related data. We observe that the trip times predominantly depend on flow mismatch and/or level errors. Controller parameters, including the integral time and gain constants, either speed up or slow down the rate of progression to a trip setpoint but do not cause a trip by themselves. The reactor can trip on a high-level signal when the reading crosses above 78%, increased from its reference level of 57%, or a low-level reading when it is below 25%. The present results show roughly how long the operators would have to respond to an attack, given a specific set of spoofing signals within the issue space analyzed. Furthermore, we have generated a simple surface by fitting a combination of exponential functions to the data obtained from the GPWR Simulator. In general, trips on a low level have been observed to occur faster than those on a high level.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Determining the impacts of venoarterial extracorporeal membrane oxygenation on cerebral oxygenation using a one-dimensional blood flow simulator

Venoarterial extracorporeal membrane oxygenation (VA-ECMO) is a mechanical system that provides rapid and short-term support for patients with cardiac failure. In many patients, pulmonary function is also impaired, resulting in poorly-oxygenated cardiac outflow competing against well-oxygenated VA-ECMO outflow, a condition known as North-South syndrome. North-South syndrome is a primary concern because of its potential to cause cerebral hypoxia, which has a critical influence on neurological complications often seen in this patient population. In order to reduce ischemic neurological complications, it is important to understand how clinical decisions regarding VA-ECMO parameters influence blood oxygenation. In our work, we studied the impacts of flow rate and cannulation site on oxygenation using a one-dimensional (1D) model to simulate blood flow. Our model was initially tested by comparing blood flow results to those observed from experimental work in VA-ECMO patients. The 1D model was combined with a two-phase flow model to simulate oxygenation. Additionally, the influence of various other clinician-tunable parameters on oxygenation in the common carotid arteries (CCAs) were tested, including, blood viscosity, cannula position within the insertion artery, heart rate, and systemic vascular resistance (SVR), as well as geometrical changes such as arterial radius and length. Our results indicated that blood oxygenation to the brain strongly depended on the cannula insertion site and the VA-ECMO flow rate with a weaker but potentially significant dependence on arterial radius. During femoral cannulation, VA-ECMO flow rates greater than ~4.9 L/min were needed to perfuse the CCAs. However, axillary and central cannulation began to perfuse the CCAs at significantly lower flow (~1 L/min). These results may help explain the incidence of cerebral hypoxia in this patient population and the common need to change cannulation strategies during treatment to address this clinical problem. While this work describes patient-averaged results, determining these relationships between VA-ECMO parameters and cerebral hypoxia is an important step towards future work to develop patient-specific models that clinicians can use to improve outcomes.

1D blood flow↗