Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transfer curve”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Observations and Parameterizations of Particle Size Distributions in Deep Tropical Cirrus and Stratiform Precipitation Clouds: Results from In-Situ Observations in TRMM Field Campaigns

In this study, we report on the evolution of particle size distributions (PSDs) and habits as measured during slow, Lagrangian-type spiral descents through deep subtropical and tropical cloud layers in Florida, Brazil, and Kwajalein, Marshall Islands, most of which were precipitating. The objective of the flight patterns was to learn more about how the PSDs evolved in the vertical and to obtain information of the vertical structure of microphysical properties. New instrumentation yielding better information on the concentrations of particles in the size (D) range between 0.2 and 2 cm, as well as improved particle imagery, produced more comprehensive observations for tropical stratiform precipitation regions and anvils than have been available previously. Collocated radar observations provided additional information on the vertical structure of the cloud layers sampled. Most of the spirals began at cloud top, with temperatures (T) as low as -50 C, and ended at cloud base or below the melting layer (ML). The PSDs broadened from cloud top towards cloud base, with the largest particles increasing in size from several millimeters at cloud top to one centimeter or larger towards cloud base. Some continued growth was noted in the upper part of the ML. Concentrations of particles less than 1 mm in size decreased with decreasing height. The result was a consistent change in the PSDs in the vertical. Similarly, systematic changes in the size dependence of the particle cross-sectional area was noted with decreasing height. Aggregation-as ascertained from both the changes in the PSDs and evolution of particle habits as observed in high detail with the cloud particle imager (CPI) probe-was responsible for these trends. The PSDs were generally well-represented by gamma distributions of the form N = N0 gamma D microns e- lambda gamma D that were fitted to the PSDs over 1-km horizontal intervals throughout the spirals. The intercept (N0 gamma), slope (lambda gamma), and dispersion (microns) values were derived for each PSD. Exponential curves (N = N0e- lambdaD; micron = 0) were also fitted to the distributions. The lambda gamma values for given spirals varied systematically with temperature as did the values of lambda (exponential), and the data generally conformed to values found in previous studies involving exponential fits to size distributions in mid-latitude frontal and cirrus layers. Considerable variability often noted in the PSD properties during the loops of individual spirals was manifested primarily in large changes in N0 gamma and N0, but micron, lambda gamma and lambda remained fairly stable. Temperature is not found to be the sole factor controlling lambda gamma or lambda but is a primary one. Direct relationships were found between lambda gamma and N0 gamma or lambda gamma and micron for the gamma distributions and lambda and N0 for the exponential. The latter relationship was not found as distinctly in earlier studies; observed PSDs in this study had better fidelity with less scatter. The micron values changed monotonically with T over the range of temperatures and were directly related to N0 gamma or lambda gamma, thereby reducing the number of variables in the PSD functional equation to two. In the upper part of the ML, N0 and lambda continued to decrease, and in the lower part these values began to increase as the largest particles melted. We developed general expressions relating various bulk microphysical, radar, and radiative transfer-related variables to N0 gamma and lambda gamma, useful for both tropical and mid-latitude clouds. These relationships facilitate the specification of a number of bulk properties in cloud and climate models. The results presented in this paper apply best to temperatures between 0 and -40 C, for which the measured radar reflectivities fall in the range of 0 to 25 dBZe.

Heymsfield, Andrew J.↗
DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Discovery of Orbital Decay in SMC X-1

We report on the results of three observations of the binary X-ray pulsar SMC X-1 with the Ginga satellite. Timing analyses of the 0.71 s X-ray pulsations yield Doppler delay curves which, in turn, enable the most accurate determination of the SMC X-1 orbital parameters available to date. Epochs of phase zero for the 3.9 day orbit were determined for 1987 May, 1988 August, and 1989 August with accuracies of 13, 0.6, and 3 s, respectively. These epochs are combined with two previous determinations of the orbital epoch to yield the rate of change in the orbital period dot-P(orb)/P(orb) = ( 3.36 +/- 0.02) x 10(exp -6) yr(exp -1). An interpretation of the orbital decay is made in the context of tidal evolution, with consideration of the influence of the increasing moment of inertia of the companion star due to its nuclear evolution. We find that, while the orbital decay is probably driven by tidal interactions, the asynchronism between the orbit and the rotation of the companion star is most likely maintained by the evolutionary expansion of the companion star (Sk 160) rather than via the Darwin instability. In this case Sk 160 is likely to be in the hydrogen shell burning phase of its evolution. Finally, a discussion is presented of the relation among the time scales for stellar evolution (less than 10(exp 7) yr), orbital decay (3 x 10(exp 5) yr), and neutron-star spin-up in the SMC X-1 system (2000 yr). In particular, we present the result of a self-consistent calculation for the histories of the spin of the neutron star and the mass transfer in this system. A plausible case can be made for the spin-up time scale being directly related to the lifetime of the luminous X-ray phase which will end in a common-envelope phase within a time of less than approx. 10(exp 4) yr.

Levine, A.↗
NASA NTRS1994

Charged-particle mutagenesis II. Mutagenic effects of high energy charged particles in normal human fibroblasts

The biological effects of high LET charged particles are a subject of great concern with regard to the prediction of radiation risk in space. In this report, mutagenic effects of high LET charged particles are quantitatively measured using primary cultures of human skin fibroblasts, and the spectrum of induced mutations are analyzed. The LET of the charged particles ranged from 25 KeV/micrometer to 975 KeV/micrometer with particle energy (on the cells) between 94-603 MeV/u. The X-chromosome linked hypoxanthine guanine phosphoribosyl transferase (hprt) locus was used as the target gene. Exposure to these high LET charged particles resulted in exponential survival curves; whereas, mutation induction was fitted by a linear model. The Relative Biological Effect (RBE) for cell-killing ranged from 3.73 to 1.25, while that for mutant induction ranged from 5.74 to 0.48. Maximum RBE values were obtained at the LET of 150 keV/micrometer. The inactivation cross-section (alpha i) and the action cross-section for mutant induction (alpha m) ranged from 2.2 to 92.0 micrometer2 and 0.09 to 5.56 x 10(-3) micrometer2, respectively. The maximum values were obtained by 56Fe with an LET of 200 keV/micrometer. The mutagenicity (alpha m/alpha i) ranged from 2.05 to 7.99 x 10(-5) with the maximum value at 150 keV/micrometer. Furthermore, molecular analysis of mutants induced by charged particles indicates that higher LET beams are more likely to cause larger deletions in the hprt locus.

NASA Center JSC↗
DOE OSTI2025

Nonthermal Ionization of Kilonova Ejecta: Observable Impacts

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass (≲0.05 M ⊙ ) and high ejecta velocity (≳0.1c). We show that, even at very early times (≲2 days), regions of ejecta fall below critical density and temperature thresholds at which nonlocal thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized r-process elements. We find that incorporating ionization from high-energy electrons produces an “inverted” and “blended” ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta, and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light-curve properties of AT 2017gfo, such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity (∼0.3 c) ejecta components, our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

High energy astrophysics↗
DOE OSTI2023

First Impressions: Early-time Classification of Supernovae Using Host-galaxy Information and Shallow Learning

Substantial effort has been devoted to the characterization of transient phenomena from photometric information. Automated approaches to this problem have taken advantage of complete phase coverage of an event, limiting their use for triggering rapid follow-up of ongoing phenomena. In this work, we introduce a neural network with a single recurrent layer designed explicitly for early photometric classification of supernovae (SNe). Our algorithm leverages transfer learning to account for model misspecification, host-galaxy photometry to solve the data-scarcity problem soon after discovery, and a custom weighted loss to prioritize accurate early classification. We first train our algorithm using state-of-the-art transient and host-galaxy simulations, then adapt its weights and validate it on the spectroscopically confirmed SNe Ia, SNe II, and SNe Ib/c from the Zwicky Transient Facility Bright Transient Survey. On observed data, our method achieves an overall accuracy of 82% ± 2% within 3 days of an event’s discovery, and an accuracy of 87% ± 5% within 30 days of discovery. At both early and late phases, our method achieves comparable or superior results to the leading classification algorithms with a simpler network architecture. These results help pave the way for rapid photometric and spectroscopic follow-up of scientifically valuable transients discovered in massive synoptic surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Radiation hydrodynamics modelling of kilonovae with SNEC

ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Control of low-speed turbulent separated flow over a backward-facing ramp

The relative performance and flow phenomena associated with several devices for controlling turbulent separated flow were investigated at low speeds. Relative performance of the devices was examined for flow over a curved, backward-facing ramp in a wind tunnel, and the flow phenomena were examined in a water tunnel using dye-flow visualization. Surface static pressure measurements and oil-flow visualization results from the wind tunnel tests indicated that transverse grooves, longitudinal grooves, submerged vortex generators, vortex generator jets (VGJ's), Viets' fluidic flappers, elongated arches at positive angle of attack, and large-eddy breakup devices (LEBU's) at positive angle of attack placed near the baseline separation location reduce flow separation and increase pressure recovery. Spanwise cylinders reduce flow separation but decrease pressure recovery downstream. Riblets, passive porous surfaces, swept grooves, Helmholtz resonators, and arches and LEBU's with angle of attack less than or = 0 degrees had no significant effect in reducing the extent of the separation region. Wall-cooling computations indicated that separation delay on a partially-cooled ramp is nearly the same as on a fully-cooled ramp, while minimizing the frictional drag increase associated with the wall cooling process. Dry-flow visualization tests in the water tunnel indicated that wishbone vortex generators in the forward orientation shed horseshoe vortices; wishbone vortex generators oriented in the reverse direction and doublet vortex generators shed streamwise counterrotating vortices; a spanewise cylinder located near the wall and LEBU's at angle of attack = -10 degrees produced eddies or transverse vortices which rotated with the same sign as the mean vorticity in a turbulent boundary layer; and the most effective VGJ's produced streamwise co-rotating vortices. Comparative wind-tunnel test results indicated that transferring momentum from the outer region of a turbulent boundary layer through the action of embedded streamwise vortices is more effective than by transverse vortices for the separation control application studied herein.

Lin, John C.↗
DOE OSTI2022

Magnetic braking saturates: evidence from the orbital period distribution of low-mass detached eclipsing binaries from ZTF

ABSTRACT We constrain the orbital period (Porb) distribution of low-mass detached main-sequence eclipsing binaries (EBs) with light-curves from the Zwicky Transient Facility (ZTF), which provides a well-understood selection function and sensitivity to faint stars. At short periods (Porb ≲ 2 d), binaries are predicted to evolve significantly due to magnetic braking (MB), which shrinks orbits and ultimately brings detached binaries into contact. The period distribution is thus a sensitive probe of MB. We find that the intrinsic period distribution of low-mass (0.1 ≲ M1/M⊙ < 0.9) binaries is basically flat (${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^0$) from Porb = 10 d down to the contact limit. This is strongly inconsistent with predictions of classical MB models based on the Skumanich relation, which are widely used in binary evolution calculations and predict ${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^{7/3}$ at short periods. The observed distributions are best reproduced by models in which the magnetic field saturates at short periods with a MB torque that scales roughly as $\dot{J}\propto P_{\rm orb}^{-1}$, as opposed to $\dot{J} \propto P_{\rm orb}^{-3}$ in the standard Skumanich law. We also find no significant difference between the period distributions of binaries containing fully and partially convective stars. Our results confirm that a saturated MB law, which was previously found to describe the spin-down of rapidly rotating isolated M dwarfs, also operates in tidally locked binaries. We advocate using saturated MB models in binary evolution calculations. Our work supports previous suggestions that MB in cataclysmic variables (CVs) is much weaker than assumed in the standard evolutionary model, unless mass transfer leads to significant additional angular momentum loss in CVs.

El-Badry, Kareem (ORCID:0000000268711752)↗
NASA NTRS2022

High-energy Photon Opacity in the Twisted Magnetospheres of Magnetars

Magnetars are neutron stars characterized by strong surface magnetic fields generally exceeding the quantum critical value of 44.1 TG. High-energy photons propagating in their magnetospheres can be attenuated by QED processes like photon splitting and magnetic pair creation. In this paper, we compute the opacities due to photon splitting and pair creation by photons emitted anywhere in the magnetosphere of a magnetar. Axisymmetric, twisted dipole field configurations embedded in the Schwarzschild metric are treated. The paper computes the maximum energies for photon transparency that permit propagation to infinity in curved spacetime. Special emphasis is given to cases where photons are generated along magnetic field loops and/or in polar regions; these cases directly relate to resonant inverse Compton scattering models for the hard X-ray emission from magnetars and Comptonized soft gamma-ray emission from giant flares. We find that increases in magnetospheric twists raise or lower photon opacities, depending on both the emission locale and the competition between field-line straightening and field strength enhancement. Consequently, given the implicit spectral transparency of hard X-ray bursts and persistent "tail" emission of magnetars, photon splitting considerations constrain their emission region locales and the twist angle of the magnetosphere; these constraints can be probed by future soft gamma-ray telescopes such as COSI and AMEGO. The inclusion of twists generally increases the opaque volume of pair creation by photons above its threshold, except when photons are emitted in polar regions and approximately parallel to the field.

radiation mechanisms: non-thermal↗
DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1998

A Geometry Based Infra-Structure for Computational Analysis and Design

The computational steps traditionally taken for most engineering analysis suites (computational fluid dynamics (CFD), structural analysis, heat transfer and etc.) are: (1) Surface Generation -- usually by employing a Computer Assisted Design (CAD) system; (2) Grid Generation -- preparing the volume for the simulation; (3) Flow Solver -- producing the results at the specified operational point; (4) Post-processing Visualization -- interactively attempting to understand the results. For structural analysis, integrated systems can be obtained from a number of commercial vendors. These vendors couple directly to a number of CAD systems and are executed from within the CAD Graphical User Interface (GUI). It should be noted that the structural analysis problem is more tractable than CFD; there are fewer mesh topologies used and the grids are not as fine (this problem space does not have the length scaling issues of fluids). For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. In most cases, the output from a CAD system could go to Initial Graphics Exchange Specification (IGES) or Standard Exchange Program (STEP) files. The output from Grid Generators and Solvers do not really have standards though there are a couple of file formats that can be used for a subset of the gridding (i.e. PLOT3D data formats). The user would have to patch up the data or translate from one format to another to move to the next step. Sometimes this could take days. Specifically the problems with this procedure are:(1) File based -- Information flows from one step to the next via data files with formats specified for that procedure. File standards, when they exist, are wholly inadequate. For example, geometry from CAD systems (transmitted via IGES files) is defined as disjoint surfaces and curves (as well as masses of other information of no interest for the Grid Generator). This is particularly onerous for modern CAD systems based on solid modeling. The part was a proper solid and in the translation to IGES has lost this important characteristic. STEP is another standard for CAD data that exists and supports the concept of a solid. The problem with STEP is that a solid modeling geometry kernel is required to query and manipulate the data within this type of file. (2) 'Good' Geometry. A bottleneck in getting results from a solver is the construction of proper geometry to be fed to the grid generator. With 'good' geometry a grid can be constructed in tens of minutes (even with a complex configuration) using unstructured techniques. Adroit multi-block methods are not far behind. This means that a million node steady-state solution can be computed on the order of hours (using current high performance computers) starting from this 'good' geometry. Unfortunately, the geometry usually transmitted from the CAD system is not 'good' in the grid generator sense. The grid generator needs smooth closed solid geometry. It can take a week (or more) of interaction with the CAD output (sometimes by hand) before the process can begin. One way Communication. (3) One-way Communication -- All information travels on from one phase to the next. This makes procedures like node adaptation difficult when attempting to add or move nodes that sit on bounding surfaces (when the actual surface data has been lost after the grid generation phase). Until this process can be automated, more complex problems such as multi-disciplinary analysis or using the above procedure for design becomes prohibitive. There is also no way to easily deal with this system in a modular manner. One can only replace the grid generator, for example, if the software reads and writes the same files. Instead of the serial approach to analysis as described above, CAPRI takes a geometry centric approach. This makes the actual geometry (not a discretized version) accessible to all phases of the analysis. The connection to the geometry is made through an Application Programming Interface (API) and NOT a file system. This API isolates the top-level applications (grid generators, solvers and visualization components) from the geometry engine. Also this allows the replacement of one geometry kernel with another, without effecting these top-level applications. For example, if UniGraphics is used as the CAD package then Parasolid (UG's own geometry engine) can be used for all geometric queries so that no solid geometry information is lost in a translation. This is much better than STEP because when the data is queried, the same software is executed as used in the CAD system. Therefore, one analyzes the exact part that is in the CAD system. CAPRI uses the same idea as the commercial structural analysis codes but does not specify control. Software components of the CAD system are used, but the analysis suite, not the CAD operator, specifies the control of the software session. This also means that the license issues (may be) minimized and individuals need not have to know how to operate a CAD system in order to run the suite.

Haimes, Robert↗
NASA NTRS1992

Bypass transition in compressible boundary layers

Transition to turbulence in aerospace applications usually occurs in a strongly disturbed environment. For instance, the effects of free-stream turbulence, roughness and obstacles in the boundary layer strongly influence transition. Proper understanding of the mechanisms leading to transition is crucial in the design of aircraft wings and gas turbine blades, because lift, drag and heat transfer strongly depend on the state of the boundary layer, laminar or turbulent. Unfortunately, most of the transition research, both theoretical and experimental, has focused on natural transition. Many practical flows, however, defy any theoretical analysis and are extremely difficult to measure. Morkovin introduced in his review paper the concept of bypass transition as those forms of transition which bypass the known mechanisms of linear and non-linear transition theories and are currently not understood by experiments. In an effort to better understand the mechanisms leading to transition in a disturbed environment, experiments are conducted studying simpler cases, viz. the effects of free stream turbulence on transition on a flat plate. It turns out that these experiments are very difficult to conduct, because generation of free stream turbulence with sufficiently high fluctuation levels and reasonable homogeneity is non trivial. For a discussion see Morkovin. Serious problems also appear due to the fact that at high Reynolds numbers the boundary layers are very thin, especially in the nose region of the plate where the transition occurs, which makes the use of very small probes necessary. The effects of free-stream turbulence on transition are the subject of this research and are especially important in a gas turbine environment, where turbulence intensities are measured between 5 and 20 percent, Wang et al. Due to the fact that the Reynolds number for turbine blades is considerably lower than for aircraft wings, generally a larger portion of the blade will be in a laminar transitional state. The effects of large free stream turbulence in compressible boundary layers at Mach numbers are examined both in the subsonic and transonic regime using direct numerical simulations. The flow is computed over a flat plate and curved surface. while many applications operate in the transonic regime. Due the nature of their numerical scheme, a non-conservation formulation of the Navier-Stokes equations, it is a non-trivial extension to compute flow fields in the transonic regime. This project aims at better understanding the effects of large free-stream turbulence in compressible boundary layers at mach number both in the subsonic and transonic regime using direct numerical simulations. The present project aims at computing the flow over a flat plate and curved surface. This research will provide data which can be used to clarify mechanisms leading to transition in an environment with high free stream turbulence. This information is useful for the development of turbulence models, which are of great importance for CFD applications, and are currently unreliable for more complex flows, such as transitional flows.

Vandervegt, J. J.↗
DOE OSTI2020

Using Digital and Web-Based Platforms to Enhance Nuclear Power Plant Knowledge Transfer - 20203

As the workforce ages, the loss of experienced personnel can result in a loss of important knowledge. Updating Electric Power Research Institute (EPRI) research and guidance with the latest technology developments and operating experience is an important method to transfer knowledge between the retiring and incoming workforce. Studies suggest that the incoming workforce is accustomed to using digital and/or internet-based platforms to find information and may have a higher comfort with these types of platforms than they have with static, written reports. Additionally, the ease and speed of updating and delivering content through digital formats may offer an effective method for knowledge transfer and combining information from multiple reports into a single access point. EPRI has recently produced an online Decommissioning Hub available to members and is working on replicating the effort in the area of radioactive waste management called the Radwaste Web References. EPRI has extensive guidance in both disciplines. The Decommissioning Hub uses a wiki format to provide a platform that allows for both the updating of that guidance resulting from new research and new operating experiences, and combining related topical information from multiple reports. It also allows users to access the information based on key words and topical searches, and can be a tool for users to input experience to share with others. The Decommissioning Hub is accessed via a web-based portal and includes topical pages; experience summaries; a question and answer capability; and a smart search function of all EPRI decommissioning reports. Use of a wiki format allows users to input new experiences or modify existing experience summaries to add additional detail. EPRI moderates all content. The experience data base has been initially populated with information available in published EPRI reports and unpublished information currently available within EPRI. Moving forward, population of the database will continue with information available within the open literature, new information developed from current decommissioning projects, and user inputs. The web site became active in late 2017, and a major update will be issued in early 2020. Larger scale updates are planned biannually, and smaller-scale updates will occur continuously. Capturing and applying lessons learned from industry experiences is a hallmark of the global nuclear industry. This is even more important in the decommissioning technology area since very few plant staff have relevant experience and thus there is a steep learning curve. Moreover, decommissioning involves a number of complex, non-routine and specialized tasks. The web-accessible and searchable experience database provides an easy-to-use resource for capturing and disseminating practical experiences from completed and ongoing decommissioning projects. The EPRI Radwaste Web References will build upon the philosophy and objectives of the Radwaste Desk References that were published in the 1990's. The purpose of the original Radwaste Desk References was to provide nuclear power plant professionals with a how-to manual for managing radioactive waste. The Radwaste Web References will adopt this original purpose and will be a wiki type web site that provide fundamental practical and theoretical information, operating experience (OE), and information about advanced technologies and methods. As much as practicable, features that allow for interaction between users to share OE, best practices, and benchmarking will be added. Where they exist and are appropriate, the various EPRI Guidelines documents on radwaste topics will be updated and incorporated into the Radwaste Web References. The first Radwaste Web Reference will be on Low Level Waste Characterization and will build upon the EPRI Low Level Waste Characterization Guidelines (EPRI Report TR-107201, 1996) and will be available in 2020. The next Radwaste Web Reference will be on the topic of liquid radioactive waste processing and management, tentatively set for unveiling in 2021. This paper summarizes the EPRI's previous experiences with web-based platforms and future plans in this area. (authors)

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

Involute Working Group – Validation of CFD Turbulence Models for Steady-State Safety Analysis

There are three research reactors in the world having fuel plates curved as circle-involute (a spiral generated around a circle): The Oak Ridge National Laboratory (ORNL) High Flux Isotope Reactor (HFIR) located in Tennessee, U.S.A. [ORNL, 2018]; The Laue-Langevin Institute (ILL) High Flux Reactor (RHF) located in Grenoble, France [ILL, 2018]; and, The Technical University of Munich (TUM) Research Neutron Source Heinz Maier-Leibnitz (FRM II) located in Garching, Germany [TUM, 2018]. All three reactors are currently using Highly Enriched Uranium (HEU, 235U/U ≥ 20 wt. %) as fuel and all three are actively engaged in activities to convert to Low Enriched Uranium (LEU, 235U/U < 20 wt. %) fuel. More information on these reactors can be found in Section two. For various reasons, these reactors have expressed interest in using Computational Fluid Dynamics (CFD) tools to perform their Steady-State Thermal-Hydraulic (SSTH) safety calculations. Using CFD tools for SSTH analysis represents generally a significant departure from traditional methods and requires thorough verification and validation (V&V) to be accepted by regulators. This is why Argonne National Laboratory (ANL) and the involute reactors formed an informal group to help each other in this endeavor. This so-called Involute Working Group (IWG) aims at supporting the qualification of CFD tools for SSTH safety calculations. Activities include benchmarking, code-to-code comparison, Verification and Validation (V&V) as well as technical support. More information on the IWG and the reasons why involute-plate reactors are interested in CFD can be found in Section three. The present study aims at testing three turbulence models commonly implemented in commercial software (namely K-Epsilon, K-Omega and Reynolds-Stress) and assessing their relevance for involute-plate reactor application. An experiment performed by Gambill and Bundy in 1961 was performed to measure heat transfer coefficients in conditions prototypic of the HFIR reactor [Gambill, 1961]. More details on this experiment are provided in Section four. Based on the information available, CFD models of this experiment have been created and the experimental results have been compared to the ones obtained with these CFD models. CFD models and the results of this comparison are provided in Section five.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1998

Theoretical Studies of Important Processes in Planetary and Comet Atmospheres

Using theoretical quantum chemical calculations, I have successfully described the dissociative recombination (DR) of O2(+) leading to the excited S-1 state of atomic oxygen, the upper state of the well known green line emission. The process is described by O2(+) + e(-) yields O(S-1)) + O(D-1) (1) where e(-) is an electron and the product oxygen atoms are both excited. This process is important in the atmospheres of Venus, Mars and Earth. I have shown in prior calculations that only one repulsive potential curve of O2, f(sup l)Sigma(sub u, sup +) can generate O(S-1) from DR of the lowest vibrational levels of O2(+). However, in the prior results, the calculated quantum yield (i.e. the number of O(S-1) atoms produced for every two product atoms) from the O2(+) v = 0 level was smaller than the laboratory and atmospheric measurements by more than an order of magnitude. Including only direct recombination, the calculated quantum yield for O(S-1) is only 0.0016. In a calculation that accounts for both direct and indirect recombination, the quantum yield is 0.0012. The range of experimentally determined quantum yields is between 0.01 and 0.23. Because of this large difference between the theoretical and experimental quantum yields, it was thought for some time that the ionospheric and laboratory O2(+) must be vibrationally excited since for excited levels, theory gave quantum yields that are similar to experimental yields. It was also suggested that some other process was generating O(S-1) but none could be identified. Under current NASA support, I have found that reaction (1) proceeds via an unusual mechanism. The f(sup 1)Sigma(sub u, sup,+) state does not cross the ion between the turning points of the v = 0 level of the O2(+) ground state. The lack of a favorable crossing leads to a very small calculated DR rate coefficient. However, this mechanism assumes that initial electron capture must occur into the repulsive state that leads to )(S-1). Instead, I have found that initial electron capture occurs mostly into the B(sup 3)Sigma(sup u, sub -) state which crosses the ion between the turning points of the v = 0 ion level and has a large DR rate coefficient. The B state dissociates to O(D-1) and )(P-3). After capture, some of the flux is transferred to the f(sup 1)Sigma(sup u, sub +) state via symmetry mixed intermediate Rydberg states. The neutral Rydberg states are a mixture of Sigma(sub u, sup +)-1 and Sigma(sub u, sup -)-3 symmetry.

Guberman, Steven L.↗
NASA NTRS2014

Optimal Stimulus Amplitude for Vestibular Stochastic Stimulation to Improve Sensorimotor Function

Sensorimotor changes such as postural and gait instabilities can affect the functional performance of astronauts when they transition across different gravity environments. We are developing a method, based on stochastic resonance (SR), to enhance information transfer by applying non-zero levels of external noise on the vestibular system (vestibular stochastic resonance, VSR). Our previous work has shown the advantageous effects of VSR in a balance task of standing on an unstable surface. This technique to improve detection of vestibular signals uses a stimulus delivery system that is wearable or portable and provides imperceptibly low levels of white noise-based binaural bipolar electrical stimulation of the vestibular system. The goal of this project is to determine optimal levels of stimulation for SR applications by using a defined vestibular threshold of motion detection. A series of experiments were carried out to determine a robust paradigm to identify a vestibular threshold that can then be used to recommend optimal stimulation levels for SR training applications customized to each crewmember. Customizing stimulus intensity can maximize treatment effects. The amplitude of stimulation to be used in the VSR application has varied across studies in the literature such as 60% of nociceptive stimulus thresholds. We compared subjects' perceptual threshold with that obtained from two measures of body sway. Each test session was 463s long and consisted of several 15s sinusoidal stimuli, at different current amplitudes (0-2 mA), interspersed with 20-20.5s periods of no stimulation. Subjects sat on a chair with their eyes closed and had to report their perception of motion through a joystick. A force plate underneath the chair recorded medio-lateral shear forces and roll moments. First we determined the percent time during stimulation periods for which perception of motion (activity above a pre-defined threshold) was reported using the joystick, and body sway (two standard deviation of the noise level in the baseline measurement) was detected by the sensors. The percentage time at each stimulation level for motion detection was normalized with respect to the largest value and a logistic regression curve fit was applied to these data. The threshold was defined at the 50% probability of motion detection. Comparison of threshold of motion detection obtained from joystick data versus body sway suggests that perceptual thresholds were significantly lower, and were not impacted by system noise. Further, in order to determine optimal stimulation amplitude to improve balance, two sets of experiments were carried out. In the first set of experiments, all subjects received the same level of stimuli and the intensity of optimal performance was projected back on subjects' vestibular threshold curve. In the second set of experiments, on different subjects, stimulation was administered from 20-400% of subjects' vestibular threshold obtained from joystick data. Preliminary results of our study show that, in general, using stimulation amplitudes at 40-60% of perceptual motion threshold improved balance performance significantly compared to control (no stimulation). The amplitude of vestibular stimulation that improved balance function was predominantly in the range of +/- 100 to +/- 400 micro A. We hypothesize that VSR stimulation will act synergistically with sensorimotor adaptability (SA) training to improve adaptability by increasing utilization of vestibular information and therefore will help us to optimize and personalize a SA countermeasure prescription. This combination will help to significantly reduce the number of days required to recover functional performance to preflight levels after long-duration spaceflight.

Goel, R.↗