Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME FUNCTION”

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,387 records · Page 77

Spectral Filtering Induced by Non-Hermitian Evolution with Balanced Gain and Loss: Enhancing Quantum Chaos

Here, the dynamical signatures of quantum chaos in an isolated system are captured by the spectral form factor, which exhibits as a function of time a dip, a ramp, and a plateau, with the ramp being governed by the correlations in the level spacing distribution. While decoherence generally suppresses these dynamical signatures, the nonlinear non-Hermitian evolution with balanced gain and loss (BGL) in an energy-dephasing scenario can enhance manifestations of quantum chaos. In the Sachdev-Ye-Kitaev model and random matrix Hamiltonians, BGL increases the span of the ramp, lowering the dip as well as the value of the plateau, providing an experimentally realizable physical mechanism for spectral filtering. The chaos enhancement due to BGL is optimal over a family of filter functions that can be engineered with fluctuating Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Advanced LIGO constraints on neutron star mergers and r-process sites

The role of compact binary mergers as the main production site of r-process elements is investigated by combining stellar abundances of Eu observed in the Milky Way, galactic chemical evolution (GCE) simulations, and binary population synthesis models, and gravitational wave measurements from Advanced LIGO. We compiled and reviewed seven recent GCE studies to extract the frequency of neutron star–neutron star (NS–NS) mergers that is needed in order to reproduce the observed [Eu/Fe] versus [Fe/H] relationship. We used our simple chemical evolution code to explore the impact of different analytical delay-time distribution functions for NS–NS mergers. We then combined our metallicity-dependent population synthesis models with our chemical evolution code to bring their predictions, for both NS–NS mergers and black hole–neutron star mergers, into a GCE context. Finally, we convolved our results with the cosmic star formation history to provide a direct comparison with current and upcoming Advanced LIGO measurements. When assuming that NS–NS mergers are the exclusive r-process sites, and that the ejected r-process mass per merger event is 0.01 M${}_{\odot }$, the number of NS–NS mergers needed in GCE studies is about 10 times larger than what is predicted by standard population synthesis models. Here, these two distinct fields can only be consistent with each other when assuming optimistic rates, massive NS–NS merger ejecta, and low Fe yields for massive stars. For now, population synthesis models and GCE simulations are in agreement with the current upper limit (O1) established by Advanced LIGO during their first run of observations. Upcoming measurements will provide an important constraint on the actual local NS–NS merger rate, will provide valuable insights on the plausibility of the GCE requirement, and will help to define whether or not compact binary mergers can be the dominant source of r-process elements in the universe.

79 ASTRONOMY AND ASTROPHYSICS↗

A generalized computer code for developing dynamic gas turbine engine models (DIGTEM)

This paper describes DIGTEM (digital turbofan engine model), a computer program that simulates two spool, two stream (turbofan) engines. DIGTEM was developed to support the development of a real time multiprocessor based engine simulator being designed at the Lewis Research Center. The turbofan engine model in DIGTEM contains steady state performance maps for all the components and has control volumes where continuity and energy balances are maintained. Rotor dynamics and duct momentum dynamics are also included. DIGTEM features an implicit integration scheme for integrating stiff systems and trims the model equations to match a prescribed design point by calculating correction coefficients that balance out the dynamic equations. It uses the same coefficients at off design points and iterates to a balanced engine condition. Transients are generated by defining the engine inputs as functions of time in a user written subroutine (TMRSP). Closed loop controls can also be simulated. DIGTEM is generalized in the aerothermodynamic treatment of components. This feature, along with DIGTEM's trimming at a design point, make it a very useful tool for developing a model of a specific turbofan engine.

Daniele, C. J.↗

A generalized computer code for developing dynamic gas turbine engine models (DIGTEM)

This paper describes DIGTEM (digital turbofan engine model), a computer program that simulates two spool, two stream (turbofan) engines. DIGTEM was developed to support the development of a real time multiprocessor based engine simulator being designed at the Lewis Research Center. The turbofan engine model in DIGTEM contains steady state performance maps for all the components and has control volumes where continuity and energy balances are maintained. Rotor dynamics and duct momentum dynamics are also included. DIGTEM features an implicit integration scheme for integrating stiff systems and trims the model equations to match a prescribed design point by calculating correction coefficients that balance out the dynamic equations. It uses the same coefficients at off design points and iterates to a balanced engine condition. Transients are generated by defining the engine inputs as functions of time in a user written subroutine (TMRSP). Closed loop controls can also be simulated. DIGTEM is generalized in the aerothermodynamic treatment of components. This feature, along with DIGTEM's trimming at a design point, make it a very useful tool for developing a model of a specific turbofan engine.

Daniele, C. J.↗

10 000-hr Cyclic Oxidation Behavior of 68 High-Temperature Co-, Fe-, and Ni- Base Alloys Evaluated at 982 deg. C (1800 deg. F)

Power systems with operating temperatures in the range of 815 to 982 C (1500 to 1800 F) frequently require alloys that can operate for long times at these temperatures. A critical requirement is that these alloys have adequate oxidation resistance. The alloys used in these power systems require thousands of hours of operating life with intermittent shutdown to room temperature. Intermittent power plant shutdowns, however, offer the possibility that the protective scale will tend to spall (i.e., crack and flake off) upon cooling, increasing the rate of oxidative attack in subsequent heating cycles. Thus, it is critical that candidate alloys be evaluated for cyclic oxidation behavior. It was determined that exposing test alloys to ten 1000-hr cycles in static air at 982 10 000-hr Cyclic Oxidation Behavior of 68 High-Temperature Co-, Fe-, and Ni-Base Alloys Evaluated at 982 C (1800 F) could give a reasonable simulation of long-time power plant operation. Iron- (Fe-), nickel- (Ni-), and cobalt- (Co-) based high-temperature alloys with sufficient chromium (Cr) and/or aluminum (Al) content can exhibit excellent oxidation resistance. The protective oxides formed by these classes of alloys are typically Cr2O3 and/or Al2O3, and are usually influenced by their Cr, or Cr and Al, content. Sixty-eight Co-, Fe-, and Ni-base high-temperature alloys, typical of those used at this temperature or higher, were used in this study. At the NASA Lewis Research Center, the alloys were tested and compared on the basis of their weight change as a function of time, x-ray diffraction of the protective scale composition, and the physical appearance of the exposed samples. Although final appearance and x-ray diffraction of the final scale products were two factors used to evaluate the oxidation resistance of each alloy, the main criterion was the oxidation kinetics inferred from the specific weight change versus time data. These data indicated a range of oxidation behavior including parabolic (typical of isothermal oxidation), paralinear, linear, and mixed-linear kinetics.

Barrett, Charles A.↗

Effect of Delayed Energy Release on Power Normalization in Reactor Depletion Calculations

Reactor depletion calculations typically assume that the thermal power of the reactor is directly proportional to the instantaneous fission rate within the core. However, approximately 7% of the sensible energy released per fission comes from the radioactive decay of fission products, which occurs alter the fission event. Thus, the assumed direct relationship between reactor power and fission rate implies that all fission products have reached secular equilibrium within the core. In situations where fission product concentrations have not reached equilibrium (e.g., shortly after reactor startup or a change in reactor power during operation), the delayed energy release due to fission product decay creates a nonlinear relationship between reactor power and fission rate. This non-linear relationship, in turn, affects the neutron flux levels and fuel consumption rate required to maintain the desired thermal power output for the reactor. For example, during initial startup of a commercial light water reactor the apparent energy release per fission is ∼7{sup %} lower than the typical value of ∼200 MeV/fission due to the potential energy stored in radioactive fission products. As a consequence, the flux level (and, therefore, uranium consumption rate) in the reactor must be correspondingly higher to preserve the desired reactor power. As the reactor operates and fission product inventories increase, the decay heat produced by the reactor will rise, producing a corresponding increase in the apparent energy release per fission. The characteristic time-dependent behavior of energy release due to the decay of fission products (i.e., decay heat) following reactor startup, shutdown, or changes in reactor power is well-known and widely-used in reactor safety analyses. However, to our knowledge, the effects of time-dependent decay heating rates on reactor depletion calculations have not been thoroughly investigated. In this paper, we present a simple methodology for explicitly accounting for delayed energy release due to the decay of fission products in reactor depletion calculations. In the proposed method, the balance equations for fission product decay energy are incorporated into the depletion coefficient matrix, using a strategy similar to that described by Isotalo in 2016. The modified system of differential equations can be solved via a matrix exponential approximation (or other suitable numerical method) to yield solutions for both nuclide densities and decay energy emission densities as a function of time. Numerical results from the method are presented for simplified (one-group point (0-D) and 1-D) representations of a commercial pressurized water reactor (PWR), in order to estimate the effect of delayed energy release on fuel consumption and flux level immediately following reactor startup.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Logical Majorana zero modes in a nanowire network

We present a scheme to use physical Majorana quasizero modes at each junction of a two-dimensional nanowire network to build a logical Majorana zero mode, the location of which is controllable through gate voltages. The wire network is a way to realize a recent proposal to imprint a Kekulé vortex pattern on a honeycomb lattice via gate voltages. We show that a specific type of junction, other than a naive Y or T junction, is needed to realize, without breaking time-reversal symmetry, an artificial “graphene” system with Majorana fermions instead of complex ones at each site. The junction we propose (i) traps exactly one physical Majorana (quasi)zero mode at each site of either a brick wall or honeycomb lattice and (ii) allows this mode to hybridize with all three neighboring sites. Using a lattice of these junctions and starting from an electronic, tight-binding model for the wires, we imprint the voltage patterns corresponding to the Kekulé vortex and observe the emergence of the logical Majorana zero mode at the vortex core. By suitable adiabatic modulations of gate voltages as a function of time, the positions of individual Kekulé vortices can be moved, and in this manner, non-Abelian braiding of the logical Majoranas bound to the vortex cores can be performed. We provide the range of parameters where these excitations could be realized experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Propagating tides in the mesosphere

A preliminary search has begun for evidence of tides in the 1-hr average line-of-sight mesospheric velocity data from the Urbana radar in the period 1978 to 1982, inclusive. Observations are restricted to the southeasterly component of those velocities. Since observations are only available for a fraction of a day due to the absence of night-time ionization in the Urbana mesosphere, it was decided to adopt an unusual procedure in the search; namely, to perform a Fourier analysis in the vertical direction and look for rotation in phase of vectors representing spatial frequency components. Propagating tidal modes would then show as vectors with a net rotation corresponding to their downward phase velocity. Five year monthly averages of hourly mean horizontal velocities inferred from the Urbana data are given. Consistent diurnal variation is seen for a number of months. These data were analyzed for vertical spatial periods of 3, 4.5, 6, 9, 12, and 24 km. When plotted as a function of time of day, many of the phasors tended to show a net rotation. The sense of rotation of a spinning vector in the complex plane can be determined objectively by computing the signed area swept out by the vector in saving from point to point. This calculation, indicated downward motion in 8 out of 12 months for the 9-km component and in 11 out of 12 months for the 24-km component. A comparison of the magnitudes of the 6 modes showed that the primary component was the 24-km component, a result confirmed by the tidal models of Forbes.

Bowhill, S. A.↗

Evolution of inclined galactic gas disks. I - A cloud-fluid approach

In this paper, numerical models are developed to follow the time-dependent behavior of inclined gas disks in nonspherical galaxies. Differential equations are developed which describe the time rate of change of the orbits of mass elements composing a thin, smoothly warping and twisting gas disk. These equations include multiple moments due to the nonsphericity of the galactic mass distribution, Coriolis effects arising from determination of orbits in a reference frame which tumbles with the surface figure of the galaxy, and viscous forces arising from cloud-cloud collisions. A cloud-fluid approach is used in determining the viscous forces. Physical arguments are used to determine the maximum permissible local value of the coefficient of kinematic viscosity. For small initial inclinations, analytic expressions are obtained for the inclination as a function of time and for the settling time, and the features of the solution are discussed.

Steiman-Cameron, Thomas Y.↗

Methanol Droplet Combustion in Oxygen-Inert Environments in Microgravity

The Flame Extinguishment (FLEX) experiment that is currently underway in the Combustion Integrated Rack facility onboard the International Space Station is aimed at understanding the effects of inert diluents on the flammability of condensed phase fuels. To this end, droplets of various fuels, including alkanes and alcohols, are burned in a quiescent microgravity environment with varying amounts of oxygen and inert diluents to determine the limiting oxygen index (LOI) for these fuels. In this study we report experimental observations of methanol droplets burning in oxygen-nitrogen-carbon dioxide and oxygen-nitrogen-helium gas mixtures at 0.7 and 1 atmospheric pressures. The initial droplet size varied between approximately 1.5 mm and 4 mm to capture both diffusive extinction brought about by insufficient residence time at the flame and radiative extinction caused by excessive heat loss from the flame zone. The ambient oxygen concentration varied from a high value of 30% by volume to as low as 12%, approaching the limiting oxygen index for the fuel. The inert dilution by carbon dioxide and helium varied over a range of 0% to 70% by volume. In these experiments, both freely floated and tethered droplets were ignited using symmetrically opposed hot-wire igniters and the burning histories were recorded onboard using digital cameras, downlinked later to the ground for analysis. The digital images yielded droplet and flame diameters as functions of time and subsequently droplet burning rate, flame standoff ratio, and initial and extinction droplet diameters. Simplified theoretical models correlate the measured burning rate constant and the flame standoff ratio reasonably well. An activation energy asymptotic theory accounting for time-dependent water dissolution or evaporation from the droplet is shown to predict the measured diffusive extinction conditions well. The experiments also show that the limiting oxygen index for methanol in these diluent gases is around 12% to 13% oxygen by volume.

Nayagam, Vedha↗

Multi-photon UV photolysis of gaseous polycyclic aromatic hydrocarbons: Extinction spectra and dynamics

The extinction spectra of static naphthalene and static biphenylene vapor, each buffered with a noble gas at room temperature, were measured as a function of time in the region between 390 and 850 nm after UV multi-photon laser photolysis at 308 nm. Employing incoherent broadband cavity enhanced absorption spectroscopy (IBBCEAS), the spectra were found to be unstructured with a general lack of isolated features suggesting that the extinction was not solely based on absorption but was in fact dominated by scattering from particles formed in the photolysis of the respective polycyclic aromatic hydrocarbon. Following UV multi-photon photolysis, the extinction dynamics of the static (unstirred) closed gas-phase system exhibits extraordinary quasi-periodic and complex oscillations with periods ranging from seconds to many minutes, persisting for up to several hours. Depending on buffer gas type and pressure, several types of dynamical responses could be generated (classified as types I, II, and III). They were studied as a function of temperature and chamber volume for different experimental conditions and possible explanations for the oscillations are discussed. A conclusive model for the observed phenomena has not been established. However, a number of key hypotheses have made based on the measurements in this publication: (a) Following the multi-photon UV photolysis of naphthalene (or biphenylene), particles are formed on a timescale not observable using IBBCEAS. (b) The observed temporal behavior cannot be described on basis of a chemical reaction scheme alone. (c) The pressure dependence of the system's responses is due to transport phenomena of particles in the chamber. (d) The size distribution and the refractive indices of particles are time dependent and evolve on a timescale of minutes to hours. The rate of particle coagulation, involving coalescent growth and particle agglomeration, affects the observed oscillations. (e) The walls of the chamber act as a sink. The wall conditions (which could not be quantitatively characterized) have a profound influence on the dynamics of the system and on its slow return to an equilibrium state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The eigenmode perspective of NMR spin relaxation in proteins

We developed in recent years the two-body (protein and probe) coupled-rotator slowly relaxing local structure (SRLS) approach for elucidating protein dynamics from NMR spin relaxation. So far we used as descriptors the set of physical parameters that enter the SRLS model. They include the global (protein-related) diffusion tensor, D{sub 1}, the local (probe-related) diffusion tensor, D{sub 2}, and the local coupling/ordering potential, u. As common in analyzes based on mesoscopic dynamic models, these parameters have been determined with data-fitting techniques. In this study, we describe structural dynamics in terms of the eigenmodes comprising the SRLS time correlation functions (TCFs) generated by using the best-fit parameters as input to the Smoluchowski equation. An eigenmode is a weighted exponential with decay constant given by an eigenvalue of the Smoluchowski operator, and weighting factor determined by the corresponding eigenvector. Obviously, both quantities depend on the SRLS parameters as determined by the SRLS model. Unlike the set of best-fit parameters, the eigenmodes represent patterns of motion of the probe-protein system. The following new information is obtained for the typical probe, the {sup 15}N−{sup 1}H bond. Two eigenmodes, associated with the protein and the probe, dominate when the time scale separation is large (i.e., D{sub 2} ≫ D{sub 1}), the tensorial properties are simple, and the local potential is either very strong or very weak. When the potential exceeds these limits while the remaining conditions are preserved, new eigenmodes arise. The multi-exponentiality of the TCFs is associated in this case with the restricted nature of the local motion. When the time scale separation is no longer large, the rotational degrees of freedom of the protein and the probe become statistically dependent (coupled dynamically). The multi-exponentiality of the TCFs is associated in this case with the restricted nature of both the local and the global motion. The effects of local diffusion axiality, potential strength, and extent of mode-coupling on the eigenmode setup are investigated. We detect largely global motional or largely local motional eigenmodes. In addition, we detect mixed eigenmodes associated with correlated/prograde or anti-correlated/retrograde rotations of the global (D{sub 1}) and local (D{sub 2}) motional modes. The eigenmode paradigm is applied to N−H bond dynamics in the β-sheet residue K19, and the α-helix residue A34, of the third immunoglobulin-binding domain of streptococcal protein G. The largest contribution to the SRLS TCFs is made by mixed anti-correlated D{sub 1} and D{sub 2} eigenmodes. The next largest contribution is made by D{sub 1}-dominated eigenmodes. Eigenmodes dominated by the local motion contribute appreciably to A34 and marginally to K19. Correlated D{sub 1} and D{sub 2} eigenmodes contribute exclusively to K19 and do not contribute above 1% to A34. The differences between K19 and A34 are delineated and rationalized in terms of the best-fit SRLS parameters and mode-mixing. It may be concluded that eigenmode analysis is complementary and supplementary to data-fitting-based analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Materials International Space Station Experiment-6 (MISSE-6) Atomic Oxygen Fluence Monitor Experiment

An atomic oxygen fluence monitor was flown as part of the Materials International Space Station Experiment-6 (MISSE-6). The monitor was designed to measure the accumulation of atomic oxygen fluence with time as it impinged upon the ram surface of the MISSE 6B Passive Experiment Container (PEC). This was an active experiment for which data was to be stored on a battery-powered data logger for post-flight retrieval and analysis. The atomic oxygen fluence measurement was accomplished by allowing atomic oxygen to erode two opposing wedges of pyrolytic graphite that partially covered a photodiode. As the wedges of pyrolytic graphite erode, the area of the photodiode that is illuminated by the Sun increases. The short circuit current, which is proportional to the area of illumination, was to be measured and recorded as a function of time. The short circuit current from a different photodiode, which was oriented in the same direction and had an unobstructed view of the Sun, was also to be recorded as a reference current. The ratio of the two separate recorded currents should bear a linear relationship with the accumulated atomic oxygen fluence and be independent of the intensity of solar illumination. Ground hyperthermal atomic oxygen exposure facilities were used to evaluate the linearity of the ratio of short circuit current to the atomic oxygen fluence. In flight, the current measurement circuitry failed to operate properly, thus the overall atomic oxygen mission fluence could only be estimated based on the physical erosion of the pyrolytic graphite wedges. The atomic oxygen fluence was calculated based on the knowledge of the space atomic oxygen erosion yield of pyrolytic graphite measured from samples on the MISSE 2. The atomic oxygen fluence monitor, the expected result and comparison of mission atomic oxygen fluence based on the erosion of the pyrolytic graphite and Kapton H atomic oxygen fluence witness samples are presented in this paper.

Banks, Bruce A.↗

Planetary Geophysics and Tectonics

Research supported by grant NAGW-1928 has addressed a variety of problems related to planetary evolution. One important focus has been on questions related to the role of chemical buoyancy in planetary evolution with application to both Venus and the Moon. We have developed a model for the evolution of the Moon (Hess and Parmentier, 1995) in which dense, highly radioactive, late stage magma ocean cumulates sink forming a core. This core heats the overlying, chemically layered mantle giving rise to a heated, chemically well-mixed layer that thickens with time. This Mixed layer eventually becomes hot enough and thick enough that its top begins to melt at a pressure low enough that melt is buoyant, thus creating mare basalts from a high pressure source of the correct composition and at an appropriate time in lunar evolution. In work completed during the last year, numerical experiments on convection in a chemically stably stratified fluid layer heated from below have been completed. These results show us how to calculate the evolution of a mixed layer in the Moon, depending on the heat production in the ilmenite- cumulate core and the chemical stratification of the overlying mantle. Chemical stratification of the mantle after its initial differentiation is would trap heat in the deep interior and prevent the rapid rise of plumes with accompanying volcanism. This trapping of heat in the interior can explain the thickness of the lunar lithosphere as a function of time as well as the magmatic evolution. We show that heat transported to the base of the lithosphere at a rate determined by current estimates of radioactivity in the Moon would not satisfy constraints on elastic lithosphere thickness from tectonic feature associated with basin loading. Trapping heat at depth by a chemically stratified mantle may also explain the absence of global compressional features on the surface that previous models predict for an initially hot lunar interior. For Venus, we developed a model in which the chemical buoyancy of crust and a depleted mantle layer stabilizes the lithosphere for long periods of time and provides a mechanism of episodic planetary evolution (Parmentier and Hess, 1992). Continued thickening of a residual depleted mantle layer eventually suppresses pressure release melting and the creation of depleted mantle. Continued cooling then allows the lithosphere to become heavier than the underlying hotter, undepleted mantle. This repeated instability can occur on time scales appropriate for episodic global resurfacing on Venus. We have also examined the role of the gabbro-eclogite phase transformation on crust and lithosphere stability and as a mechanism of crustal recycling in the absence of plate tectonics. Our work thus far concentrates on the scale of instability that would occur due to cooling or crustal thickening associated with horizontal shortening. Whether repeated overturn can explain the evolution of Venus depends in part on whether sufficient heat transfer can occur between overturns and on constraints provided by understanding observed surface features and evolution.

Parmentier, E. M.↗

The Gravitational Wave Signal from Core-collapse Supernovae

We report on gravitational waves (GWs) from a set of 2D multigroup neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Our aim is to systematize the current knowledge about the post-bounce CCSN GW signal and recognize the templatable features that could be used by the ground-based laser interferometers. We demonstrate that, starting from ~400 ms after core bounce, the dominant GW signal represents the fundamental quadrupole (l = 2) oscillation mode (f-mode) of the proto–neutron star (PNS), which can be accurately reproduced by a linear perturbation analysis of the angle-averaged PNS profile. Before that, in the time interval between ~200 and ~400 ms after bounce, the dominant mode has two radial nodes and represents a g-mode. We associate the high-frequency noise in the GW spectrograms above the main signal with p-modes, while below the dominant frequency there is a region with very little power. The collection of models introduced here summarizes the dependence of the CCSN GW signal on the progenitor mass, equation of state, many-body corrections to the neutrino opacity, and rotation. Weak dependence of the dominant GW frequency on the progenitor mass motivates us to provide a simple fit for it as a function of time, which can be used as a prior when looking for CCSN candidates in the LIGO data.

79 ASTRONOMY AND ASTROPHYSICS↗

Real-time crystallization in fluorinated parylene probed by conductivity spectra

Dielectric relaxation spectroscopy experiments were performed at high temperature on fluorinated parylene films during the occurrence of the isothermal crystalline phase transition. For this polymer, since the difference between the glass transition temperature (T{sub g}) and the phase transition temperature (T{sub c}) is very strong (T{sub c} ≥ 4T{sub g}), segmental and dipolar relaxation usually used to probe the crystallization are not shown in the experiment frequency window (10{sup −1} to 10{sup 6} Hz) during the crystallization. The charge diffusion becomes the only electrical marker that allows probing the phase transition. During the transition phase, a continuous decrease of about two orders of magnitude is observed in the conductivity values below an offset frequency (f{sub c}) with a tendency to stabilization after 600 min. Below the offset frequency, the decrease of the normalized conductivity to the initial value as function of time is frequency independent. The same behavior is also observed for the f{sub c} values that decrease from 160 Hz to about 20 Hz. Above the offset frequency, the electronic hopping mechanism is also affected by the phase transition and the power law exponent (n) of the AC conductivity shows a variation from 0.7 to 0.95 during the first 600 min that tend to stabilize thereafter. Accordingly, three parameters (n, f{sub c}, and AC conductivity values for frequencies below f{sub c}) extracted from the AC conductivity spectra in different frequency windows seem suitable to probe the crystalline phase transition.

36 MATERIALS SCIENCE↗

IMPLEMENTING PRIMARY CREEP CALCULATIONS DURING STRESS CHANGES AND REVERSALS IN THE FUEL PERFORMANCE CODE FRAPCON-3

Recently, the creep models in FRAPCON-3.4 were used to model the creep tests performed in IFA-699 up through cycle 4. IFA-699 is an instrumented test rig in the Halden Boiling Water Reactor that is being used to perform in-reactor creep experiments under changing stress conditions. The cladding creep tests in IFA-699 were on four fresh fuel segments with Zircaloy-4, E110, M-MDA, and M5™ cladding. These fuel segments were loaded in a test rig and placed in a PWR loop and operated under PWR conditions. Each of the test fuel rods were connected to gas lines in order to apply tensile and compressive stresses to the cladding. Diameters of the cladding were measured by diameter gages at several intervals to measure the creep behavior of the cladding as a function of time. The creep tests were performed at a temperature of around 360°C and fast neutron flux of 1.8 to 2.4×1017 n/m²-s (depends on specimen location). During the four cycles of operation, hoop stresses of -75, -50, +30, and +110 MPa were applied to the specimens. In addition to providing assessment of the FRAPCON-3.4 creep modeling, the creep model comparisons were also used to provide insight into the nature of primary creep following stress changes, whether creep behavior is symmetric under tensile and compressive loading, and the effect of fast fluence on creep rate. It was determined that creep strain should be calculated as a function of effective stress rather than hoop stress. This accounts for the observed difference between creep under compressive and tensile stress for the same hoop stress. It was also determined that the primary creep changes in proportion to the change in stress following stress change or stress reversal. Finally, it was observed that the irradiation creep rate decreases somewhat with increasing fluence. The creep model in FRAPCON-3.4 has already been developed to be a function of effective stress and includes a term that decreases the irradiation creep rate with increasing fluence. A new mathematical formulation has been developed that is suitable for a fuel performance code to calculate the primary creep strain under conditions of changing stress and stress reversal. The creep model in FRAPCON-3.4 and the formulation for calculating primary creep under changing stress conditions will be described in this paper. Comparisons to data will also be provided.

Geelhood, Kenneth J.↗

Characterization of self-magnetic pinch radiographic diode performance on RITS-6 at Sandia National Laboratories. II. Coupling between the inductive voltage adder and the SMP load

The self-magnetic pinch (SMP) diode is a type of radiographic diode used to generate an intense electron beam for radiographic applications. At Sandia National Laboratories, SMP was the diode load for the six-cavity radiographic integrated test stand inductive voltage adder (IVA) driver operated in a magnetically insulated transmission line (MITL). The MITL contributes a flow current in addition to the current generated within the diode itself. Extensive experiments with a MITL of 40Ω load impedance [T. J. Renk et al., Phys. Plasmas 29, 023105 (2022)] indicate that the additional flow current leads to results similar to what might be expected from a conventional high-voltage interface driver, where flow current is not present. However, when the MITL flow impedance was increased to 80 Ω, qualitatively different diode behavior was observed. This includes large retrapping waves suggestive of an initial coupling to low impedance as well as diode current decreasing with time even as the total current does not. A key observation is that the driver generates total current (flow + diode) consistent with the flow impedance of the MITL used. The case is made in this paper that the 80Ω MITL experiments detailed here can only be understood when the IVA-MITL-SMP diode is considered as a total system. Additionally, the constraint of fixed total current plus the relatively high flow impedance limits the ability of the diode (whether SMP or other type) to act as an independent load. An unexpected new result is that in tracking the behavior of the electron strike angle on the converter as a function of time, we observed that the conventional cIV x “Radiographic” radiation scaling (where x~2.2) begins to break down for voltages above 8 MV, and cubic scaling is required to recover accurate angle tracking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗