Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Decay 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 379 records · Page 21

Ultraviolet and Extreme-Ultraviolet Emissions at the Flare Footpoints Observed by Atmosphere Imaging Assembly

A solar flare is composed of impulsive energy release events by magnetic reconnection, which forms and heats flare loops. Recent studies have revealed a two-phase evolution pattern of UV 1600 A emission at the feet of these loops: a rapid pulse lasting for a few seconds to a few minutes, followed by a gradual decay on timescales of a few tens of minutes. Multiple band EUV observations by the Atmosphere Imaging Assembly further reveal very similar signatures. These two phases represent different but related signatures of an impulsive energy release in the corona. The rapid pulse is an immediate response of the lower atmosphere to an intense thermal conduction flux resulting from the sudden heating of the corona to high temperatures (we rule out energetic particles due to a lack of significant hard X-ray emission). The gradual phase is associated with the cooling of hot plasma that has been evaporated into the corona. The observed footpoint emission is again powered by thermal conduction (and enthalpy), but now during a period when approximate steady-state conditions are established in the loop. UV and EUV light curves of individual pixels may therefore be separated into contributions from two distinct physical mechanisms to shed light on the nature of energy transport in a flare.We demonstrate this technique using coordinated, spatially resolved observations of UV and EUV emissions from the footpoints of a C3.2 thermal flare.

Flares↗

Switching Characteristics of a 4H-SiC Based Bipolar Junction Transistor to 200 C

Static curves and resistive load switching characteristics of a 600 V, 4 A rated, SiC-based NPN bipolar power transistor (BJT) were observed at selected temperatures from room to 200 C. All testing was done in a pulse mode at low duty cycle (approx.0.1 percent). Turn-on was driven by an adjustable base current pulse and turn-off was accelerated by a negative base voltage pulse of 7 V. These base drive signals were implemented by 850 V, gated power pulsers, having rise-times of roughly 10 ns, or less. Base charge sweep-out with a 7 V negative pulse did not produce the large reverse base current pulse seen in a comparably rated Si-based BJT. This may be due to a very low charge storage time. The decay of the collector current was more linear than its exponential-like rise. Switching observations were done at base drive currents (I(sub B)) up to 400 mA and collector currents (I(sub C)) up to 4 A, using a 100 Omega non-inductive load. At I(sub B) = 400 mA and I(sub C) = 4 A, turn-on times typically varied from 80 to 94 ns, over temperatures from 23 to 200 C. As expected, lowering the base drive greatly extended the turn-on time. Similarly, decreasing the load current to I(sub C) = 1 A with I(sub B) = 400 mA produced turn-on times as short as 34 ns. Over the 23 to 200 C range, with I(sub B) = 400 mA and I(sub C) = 4 A, turn-off times were in the range of 72 to 84 ns with the 7 V sweep-out.

Niedra, Janis M.↗

The effects of streamline curvature and swirl on turbulent flows in curved ducts

A technique for improving the numerical predictions of turbulent flows with the effect of streamline curvature is developed. Separated flows, the flow in a curved duct, and swirling flows are examples of flow fields where streamline curvature plays a dominant role. A comprehensive literature review on the effect of streamline curvature was conducted. New algebraic formulations for the eddy viscosity incorporating the kappa-epsilon turbulence model are proposed to account for various effects of streamline curvature. The loci of flow reversal of the separated flows over various backward-facing steps are employed to test the capability of the proposed turbulence model in capturing the effect of local curvature. The inclusion of the effect of longitudinal curvature in the proposed turbulence model is validated by predicting the distributions of the static pressure coefficients in an S-bend duct and in 180 degree turn-around ducts. The proposed turbulence model embedded with transverse curvature modification is substantiated by predicting the decay of the axial velocities in the confined swirling flows. The numerical predictions of different curvature effects by the proposed turbulence models are also reported.

Cheng, Chih-Hsiung↗

The 100-month Swift catalogue of supergiant fast X–ray transients

Supergiant fast X-ray transients (SFXTs) are high mass X-ray binaries (HMXBs) displaying X-ray outbursts that can reach peak luminosities up to 10 38 erg s -1 and spend most of their lives in more quiescent states with luminosities as low as 10 32 -10 33 erg s -1 . During the quiescent states, less luminous flares are also frequently observed with luminosities of 10 34 -10 35 erg s -1 . The main goal of the comprehensive and uniform analysis of the SFXT Swift triggers presented in this paper is to provide tools to predict whether a transient that has no known X-ray counterpart may be an SFXT candidate. These tools can be exploited for the development of future missions exploring the variable X-ray sky through large field-of-view instruments. We examined all available data on outbursts of SFXTs that triggered the Swift/Burst Alert Telescope (BAT) collected between 2005 August 30 and 2014 December 31, in particular those for which broad-band data, including the Swift/X-ray Telescope (XRT) data, are also available. This work complements and extends our previous catalogue of SFXT flares detected by BAT from 2005 February 12 to 2013 May 31, since we now include the additional BAT triggers recorded until the end of 2014 (i.e. beyond the formal first 100 months of the Swift mission). Due to a change in the mission’s observational strategy, virtually no SFXT triggers obtained a broad-band response after 2014. We processed all BAT and XRT data uniformly by using the Swift Burst Analyser to produce spectral evolution dependent flux light curves for each outburst in the sample. The BAT data allowed us to infer useful diagnostics to set SFXT triggers apart from the general γ-ray burst population, showing that SFXTs uniquely give rise to image triggers and are simultaneously very long, faint, and ‘soft’ hard-X-ray transients. We find that the BAT data alone can discriminate very well the SFXTs from other classes of fast transients, such as anomalous X-ray pulsars and soft gamma repeaters. On the contrary, the XRT data collected around the time of the BAT triggers are shown to be decisive for distinguishing SFXTs from, for instance, accreting millisecond X-ray pulsars and jetted tidal disruption events. The XRT observations of 35 (out of 52 in total) SFXT BAT triggers show that in the soft X-ray energy band, SFXTs display a decay in flux from the peak of the outburst of at least three orders of magnitude within a day and rarely undergo large re-brightening episodes, favouring in most cases a rapid decay down to the quiescent level within three to five days (at most).

79 ASTRONOMY AND ASTROPHYSICS↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Solar PV

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) research platform. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence data centers and other variable loads. This dataset entry describes the behavior of a 1.25-MW proton exchange membrane MC250 electrolyzer system, manufactured by Nel Hydrogen , [1] when fed historical data generated by the 430-kW, fixed-axis solar photovoltaic (PV) array located at NLR’s Flatirons Campus. (While the electrolyzer balance of plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack.) Solar PV power output data for the 2020 calendar year were categorized on a daily basis by total energy generation and standard deviation. Each day was then ranked by these metrics, and the 25th, 50th, and 100th percentiles were selected. The 75th percentile day did not exhibit sufficient variability to make for a valuable experiment. A similar process was used for the related historical wind dataset . [2] The historical days in 2020 that represented these percentiles are Dec. 19, March 29, and May 4, respectively. The entire solar day’s power profile was then fed through the MC250 electrolyzer. Due to its length, the 100th percentile day experiment was split into two parts, and the final 3 hours of the solar day were not captured. These final 3 hours contained no spikes or dips of interest and simply represented a slow decay of input solar power. Also, a single timestamp (13:13:47 on Jan. 14, 2026) was lost in the hydrogen system supervisory control and data acquisition. Finally, during the 25th percentile experiment (solar day Dec. 19, 2020) data recording was lost from 11:00:13 to 11:14:45. The roughly 15 minutes of the solar profile were rerun at the end of the experiment and spliced into this time slot during post-processing. The electrolysis system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operation of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The historical solar profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. For more details on the statistical analysis process, see the slide deck “Public Reference Data for Megawatt-Scale Hydrogen Electrolysis: NLR Historical Solar PV Analysis and Profile Generation” accessible with this data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single solar PV electrolysis experiment and is formatted as: {technology}_{percentile}_{scaling factor} For instance, “solarPV-430kW_25_2x.zip” reports the experiment using the 25th percentile solar data from the historical 2020 solar PV dataset, scaled to 200%. Scaling factors were applied to the generated solar PV power output files to more closely match the 1.25-MW capacity of the electrolyzer. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and solar power input. A PDF file detailing the historical solar data statistical analysis used to generate the solar profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all experiments combined into one dataset labeled "combined_solarPV_experiments.csv". [1] nelhydrogen.com/product/mc-series-electrolyser . [2] data.nlr.gov/submissions/316 .

08 HYDROGEN↗

Vortex bursting and associated twist dynamics on helical vortex tubes and vortex rings

The interaction of opposite-signed twist waves on vortex tubes can lead to vortex bursting, a process where the core expands into a double ring-like structure with strong swirling flows. Previous works have studied vortex bursting on rectilinear vortices by axially perturbing the initial core size to generate the twist waves, and observed largely axisymmetric bursting dynamics. In this work, we numerically study bursting on vortical structures with curved centrelines, analysing the interaction between the centreline dynamics, twist wave generation and propagation, and vortex bursting. We focus on axially perturbed helical vortex tubes with small radius-to-pitch ratios up to$0.0625$, as well as vortex rings with a large radius-to-core size ratio$10$, both at a circulation-based Reynolds number$5000$. The results show that though the initial twist wave propagation speeds are relatively unaffected by the curvature and torsion of the centreline, the bursting process is altered significantly compared with rectilinear vortices. The self-induced rotation of the centreline of the helical tube induces a non-axisymmetric distortion of the bursting structure, which rapidly breaks up the vortex core into small-scale helical structures. A similar destabilization of the bursting structure also occurs on vortex rings. The enstrophy increase and accelerated energy decay associated with bursting are predominantly determined by the twist wave strength, rather than the curvature and torsion of the centreline. Combined, our findings imply that bursting could play an important role in transferring and dissipating energy of vortical structures in wakes, and turbulent flows in general.

Mechanics↗

Disk instability and the spectral evolution of the 1992 outburst of the intermediate polar GK Persei

The disk instability model can explain the previous history of dwarf-nova-like outbursts in the intermediate polar GK Per, which occur about once every three years. Disk models that reproduce the recurrence time and outburst light curves suggest that GK Per has a large effective inner disk radius (approx. 30-40 white dwarf radii) truncated by a strong magnetic field (10(exp 7) G). In this context, the effective radius is that of the portion of the disk that participates in the disk thermal instability. The radius derived is larger than the corotation radius, which must be an upper limit on the true dynamical inner radius of the disk. Disk instability models with this large effective inner radius predict that the ultraviolet continuum should be rather flat. Here we compare the predictions of the disk instability model to IUE observations of the 1981 outburst and to IUE and ROSAT observations of the recent 1992 outburst of GK Per. The model disk continuum spectral evolution is consistent with the observed UV and optical spectra, especially at maximum and in the early decay phase of the outburst. The consistency of the model with the observed UV spectra suggests that the effective inner radius of the disk is almost constant, independent of mass accretion rate, and that whatever structure lies between the effective inner radius and the corotation radius neither participates in the disk instability nor radiates substantially in the UV. The related physics of the inner disk region will be briefly discussed.

Kim, S.-W.↗

Exploring unsupervised top tagging using Bayesian inference

Recognizing hadronically decaying top-quark jets in a sample of jets, or even its total fraction in the sample, is an important step in many LHC searches for Standard Model and Beyond Standard Model physics as well. Although there exists outstanding top-tagger algorithms, their construction and their expected performance rely on Montecarlo simulations, which may induce potential biases. For these reasons we develop two simple unsupervised top-tagger algorithms based on performing Bayesian inference on a mixture model. In one of them we use as the observed variable a new geometrically-based observable \tilde{A}_{3} A ̃ 3 , and in the other we consider the more traditional \tau_{3}/\tau_{2} τ 3 / τ 2 N N -subjettiness ratio, which yields a better performance. As expected, we find that the unsupervised tagger performance is below existing supervised taggers, reaching expected Area Under Curve AUC \sim 0.80-0.81 ∼ 0.80 − 0.81 and accuracies of about 69% - − 75% in a full range of sample purity. However, these performances are more robust to possible biases in the Montecarlo that their supervised counterparts. Our findings are a step towards exploring and considering simpler and unbiased taggers.

Alvarez, Ezequiel↗

Observations of gamma-ray line profiles from SN 1987A

Fully resolved gamma-ray line observations from the decay of Co-56 in SN 1987A are presented. Data are from the first two balloon flights of the Gamma-Ray Imaging Spectrometer. On day 433 (after the initial optical sighting) the 847 and 1238 keV lines from Co-56 were observed at 2.3 and 4.3 sigma significance. On day 613, the lines at 847, 1238, and 2599 keV were observed at 4.6, 3.4, and 1.9 sigma, respectively. The combined significance for the three-line complex in both flights is 7.8 sigma. Gaussian profiles yield acceptable least-squares fits to the lines. The line profiles are centered on the red side of the rest energy with typical velocity dispersions of about 3500 km/s FWHM, consistent with an optically thin source, but the line intensities are less than about 30 percent of those produced by the 0.075 solar mass of Co-56 determined from the bolometric light curve.

Tueller, J.↗

RHESSI and Trace Observations of the 21 April 2002 X1.5 Flare

Observations of the X1.5 flare on 21 April 2002 are reviewed using the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and the Transition Region and Coronal Explorer (TRACE). The major findings are as follows: 1. The 3-25 keV X-rays started 54 mins before the EUV (195 A) emission suggesting that the initial energy release heated plasma directly to 220 MK, well above the 1.6 MK needed to produce the Fe XII (195 A) line. 2. Using coaligned 12-25 keV RHESSI and TRACE images, further evidence is found for the existence of hot (15-20 MK) plasma in the 195 A passband. This hot, diffuse emission is attributed to the presence of the Fe XXIV (192 A) line within the TRACE 195 A passband. 3. The 12-25 keV source centroid moves away from the limb with an apparent velocity of approx. 9.9 km/ s, slowing to approx. 1.7 km/ s after 3 hours, its find altitude being approx. 120 Mm after approx. 12 hours. This suggests that the energy release site moves to higher altitudes in agreement with classical flare models. 4. The 50-100 keV emission correlates well with EUV flare ribbons, suggesting thick-target interactions at the footpoints of the magnetic arcade. The 50-100 keV time profile matches the time derivative of the GOES light curve (Neupert effect), which suggests that the same electrons that produced the thick-target hard X-ray emission also heat the plasma seen in soft X-rays. 5. X-ray footpoint emission has an E(sup -3) spectrum down to approx. 10 keV suggesting a lower electron cutoff energy than previously thought. 6. The hard X-ray (25-200 keV) peaks have FWHM durations of approx. 1 min suggesting a more gradual energy release process than expected. 7. The TRACE images reveal a bright symmetric front propagating away from the main flare site at speeds of greater than or = 120 km/ s. This may be associated with fast CME observed several minutes later by LASCO. 8. Dark sinuous lanes are observed in the TRACE images that extend almost radially from the post-flare loop system. This "fan of spines" becomes visible well into the decay phase of the flare and shows evidence for both lateral and downward motions.

Gallagher, Peter T.↗

Cumulative short-lived photofission product yields for nuclear forensic application

Experiments have long been carried out with linear accelerator (linac) produced bremsstrahlung X-rays to study high-energy photon-induced fission, known as “photofission”. One benefit of photofission, when compared with neutron-induced fission, is the ability to investigate fundamental physics of fission at excitation energies lower than what is possible with neutrons [1]. Absorption of an incident neutron results in a compound nucleus with an excitation energy at least equal to the neutron binding energy; thus, the lowest energy region of the fission barrier cannot be studied. Fission product yield distributions are dependent on the incident particle energy and the target nucleus mass. They are generally represented by a double-humped curve, with low and high mass “peaks” and an intermediate mass “valley”. As incident particle energy increases, the excitation energy of the nucleus increases, thereby increasing the overall yield of isotopes in the valley with respect to the two maxima or peaks of the fission product curve [3]. Photofission is most probable for incident photons around 14 MeV due to the giant dipole resonance. Induced photofission around this energy results in a fission product yield distribution that resembles 14 MeV DT neutron fission more than thermal or fast (500 keV) neutron fission. In post-detonation nuclear forensics, the peak-to-valley ratio of the fission product yield curve could be used for characterizing device type by providing an indication of the neutron energy spectrum that had induced fission (e.g., 500 keV vs. 14 MeV). The goal of this work is to improve nuclear data for isotopes of interest in nuclear forensics and special nuclear material (SNM) detection, specifically the fission product yield data via photofission experiments for the isotopes of 238U and 232Th. Photofission has been proposed as a less costly and more flexible production method for the valley isotopes that are representative of DT neutron fission product yield distributions, which would benefit nuclear forensics exercises in testing for analysis techniques. Natural uranium and thorium targets were irradiated with bremsstrahlung photons (nominal endpoint energies of 8, 14, and 20 MeV), generated via impingement of high-energy electrons on a tungsten electron-to-photon radiator. Gamma-ray spectroscopy was used to observe and identify the decay of the resulting fission products. The photofission yields of 102 fission products were measured for the set of the 238U and 232Th targets at the three bremsstrahlung X-ray endpoint energies 8, 14, and 20 MeV. These fragments included masses ranging from A=84 to A=144, with half-lives as short as 1.07 seconds.

07 ISOTOPE AND RADIATION SOURCES↗

Quantification of Residual Water in Spent Fuel Dry Storage Canisters Using Mass Spectrometry

The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal, such as fuel degradation; cladding corrosion, embrittlement, or breaching; and the creation of a flammable environment via radiolysis of water. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents details on the quantification of residual water in the Advanced Drying Cycle Simulator (ADCS), an apparatus built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. The ADCS is outfitted with thermocouples to measure the thermal response of the ADCS to simulated decay heats and internal helium fill pressures relevant to commercial drying procedures. The ADCS is also instrumented with pressure transducers to measure the pressures and vacuum levels observed during simulated commercial drying. The most unique instrumentation used for quantifying residual water in the ADCS is a Hiden Analytical HPR-30 mass spectrometer (MS), which measures gas compositions of the ADCS internal free volume, based on partial pressures calculated from relative proportions of gas molecules detected by the MS. This report details the methodology used to implement MS measurements in quantifying residual water in the ADCS. This methodology includes the calibration of the HPR-30 MS to a Buck Research Instruments CR-4 chilled mirror hygrometer, which itself is calibrated to a NIST-traceable standard. Data collected by both the MS and the chilled mirror hygrometer from water/helium mixtures ranging from 150 to 500,000 ppmv water in helium were used to generate calibration curves, establishing a source of verification of MS measured water contents. Details regarding water content measurement uncertainties are included in this report, defining the accuracy and verifiability of the HPR-30 MS in measuring residual water content in simulated dry storage canister environments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

On the nonlinear development of the most unstable Goertler vortex mode

The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

On the nonlinear development of the most unstable Goertler vortex mode TLSP - Final report

The nonlinear development of the most unstable Goertler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Goertler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Goertler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Goertler equations valid for O(1) Goertler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Kepler Light Curves of V1504 Cygni and V344 Lyrae: A Study of the Outburst Properties

We examine the Kepler light curves of V1504 Cyg and V344 Lyr, encompassing approximately 460 d at 1 min cadence. During this span each system exhibited approximately 40 outbursts, including four superoutbursts. We find that, in both systems, the normal outbursts lying between two superoutbursts increase in duration by a factor approximately 1.2 - 1.7, and then reset to a small value after the following superoutburst. In V344 Lyr the trend of quiescent intervals between normal outbursts is to increase to a local maximum about half way through the supercycle the interval from one superoutburst to the next - and then to decrease back to a small value by the time of the next superoutburst, whereas for V1504 Cyg the quiescent intervals are relatively constant during the supercycle. Both of these trends are inconsistent with the Osaki's thermal-tidal model, which robustly predicts a secular increase in the quiescent intervals between normal outbursts during a supercycle. Also, most of the normal outbursts have an asymmetric, fast-rise/slower-decline shape, which would be consistent with outbursts triggered at large radii. The exponential rate of decay of the plateau phase of the superoutbursts is 8 d mag(sup -1) for approximately 1504 Cyg and 12 d mag(sup -1) for V344 Lyr. This time scale gives a direct measure of the VISCOUS time scale III the outer accretion disk given the expectation that the entire disk is in the hot, viscous state during superoutburst. The resulting constraint on the Shakura-Sunyaev parameter, alpha(sub hot) approximately equal to 0.1, is consistent with the value inferred from the fast dwarf nova decays. By looking at the slow decay rate for superoutbursts, which occur in systems below the period gap, in combination with the slow decay rate in one long outburst above the period gap (in U Gem), we infer a steep dependence of the decay rate on orbital period for long outbursts. We argue that this relation implies a steep dependence of alpha(sub cold) on orbital period, which may be consistent with recent findings of Patterson, and is consistent with tidal torquing as being the dominant angular momentum transport mechanism in quiescent disks in interacting binary systems.

Cannizzo, John K.↗

Preliminary Design of a New Pulsed-Neutron Die--Away Experiments with Absorbers (IER-552 CED-1)

This report presents the preliminary design (CED-1) of IER-552, pulsed-neutron die-away (PNDA) experiments with absorbers, to be conducted by Lawrence Livermore National Laboratory (LLNL). The goal of the PNDA experiments with absorbing materials is to produce high-quality, low-cost integral benchmarks that can be used to validate cross sections that are high priority nuclear data for the Department of Energy’s Nuclear Criticality Safety Program (NCSP). PNDA experiments can be used to generate nuclear data validation benchmarks quickly and cheaply, as they do not require fissile material or a nuclear facility to conduct them. The original PNDA experiments were set up at LLNL under IER-501. Extending the PNDA experiments to include absorbing materials leverages existing facilities and resources at LLNL to validate cross sections, with a basic experiment design that uses a neutron generator to impinge a short, mono-energetic neutron pulse on a poisoned, moderating target sample. The neutron pulse lasts 100 µs, after which the neutron population in the sample reaches thermal equilibrium and later spatial equilibrium. In this fully equilibrized state, the neutron population has a characteristic decay-time eigenvalue. The eigenvalue can be used as an integral parameter to validate nuclear data involved with neutron scattering and absorption. Two observables can be measured: the neutron population and the population of secondary γ-rays produced by inelastic and capture reactions. This report explores the feasibility of both approaches. For γ-ray measurements, this allows for direct measurement of reactions of interest for validation. During the neutron burst time, inelastic (prompt) γ-rays dominate, and these data can be isolated for validation of neutron inelastic scattering data. After the pulse ends, neutron capture (delayed) γ-rays dominate and are produced by many reaction channels as neutrons slow down to thermal energies. Two sets of γ-ray spectra (inelastic, capture) can be produced for validation. Additionally, the time-dependent decay of the capture γ-ray population when the neutron population is thermalized and in its fundamental mode can be used for validation. Two experiments can be imagined where pulsed neutrons are used to validate the nuclear data of absorbing materials. The first integrates well into the existing PNDA testbed at LLNL. The absorber material would be introduced into a moderating target and the same neutron generator, neutron detectors, electronics, and shielding box used for IER-501 would be employed. From the same experimental setup, many absorbing materials can be validated by simply changing the neutron poison in the moderator. The only significant modification would be the introduction of a γ-ray detector. Such an experiment would primarily serve to validate nuclear data in the thermal energy range. The second kind of experiment would validate the fast and epithermal energy range. It would involve pulsing neutron into a bare assembly of the absorber material and measuring the die-away of neutrons. The fast leakage of such an assembly requires nanosecond neutron-generator pulse widths and fast detectors/electronics that could measure ~200 ns die-away curves. This kind of experiment would ideally be performed at a cyclotron, like that available at Lawrence Berkeley National Laboratory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advances in phosphor two-color ratio method thermography for full-field surface temperature measurements

We report thermographic phosphors can be employed for optical sensing of surface, gas phase, and bulk material temperatures through different strategies including the time-decay method, time-integrated method, and frequency-domain method. We focus on the time-integrated method, also known as the ratio method, as it can be more practical in many situations. This work advances the ratio method using two machine vision cameras with CMOS detectors for full-field temperature measurements of a solid surface. A phosphor calibration coupon is fabricated using aerosol deposition and employed for in situ determination of the temperature-versus-intensity ratio relationship. Algorithms from digital image correlation are employed to determine the stereoscopic imaging system intrinsic and extrinsic parameters, and accurately register material points on the sample to subpixel locations in each image with 0.07 px or better accuracy. Detector nonlinearity is carefully characterized and corrected. Temperature-dependent, spatial non-uniformity of the full-field intensity ratio—posited to be caused by a blue-shift effect of the bandpass filter for non-collimated light and/or a wavelength-dependent transmission efficiency of the lens—is assessed and treated for cases where a standard flat-field correction fails to correct the non-uniformity. In sum, pixel-wise calibration curves relating the computed intensity ratio to temperature in the range of T = 300–430 K are generated, with an embedded error of less than 3 K. This work offers a full calibration methodology and several improvements on two-color phosphor thermography, opening the door for full-field temperature measurements in dynamic tests with deforming test articles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗