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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Computational Model Prediction and Biological Validation Using Simplified Mixed Field Exposures for the Development of a GCR Reference Field

The yield of chromosomal aberrations has been shown to increase in the lymphocytes of astronauts after long-duration missions of several months in space. Chromosome exchanges, especially translocations, are positively correlated with many cancers and are therefore a potential biomarker of cancer risk associated with radiation exposure. Although extensive studies have been carried out on the induction of chromosomal aberrations by low- and high-LET radiation in human lymphocytes, fibroblasts, and epithelial cells exposed in vitro, there is a lack of data on chromosome aberrations induced by low dose-rate chronic exposure and mixed field beams such as those expected in space. Chromosome aberration studies at NSRL will provide the biological validation needed to extend the computational models over a broader range of experimental conditions (more complicated mixed fields leading up to the galactic cosmic rays (GCR) simulator), helping to reduce uncertainties in radiation quality effects and dose-rate dependence in cancer risk models. These models can then be used to answer some of the open questions regarding requirements for a full GCR reference field, including particle type and number, energy, dose rate, and delivery order. In this study, we designed a simplified mixed field beam with a combination of proton, helium, oxygen, and iron ions with shielding or proton, helium, oxygen, and titanium without shielding. Human fibroblasts cells were irradiated with these mixed field beam as well as each single beam with acute and chronic dose rate, and chromosome aberrations (CA) were measured with 3-color fluorescent in situ hybridization (FISH) chromosome painting methods. Frequency and type of CA induced with acute dose rate and chronic dose rates with single and mixed field beam will be discussed. A computational chromosome and radiation-induced DNA damage model, BDSTRACKS (Biological Damage by Stochastic Tracks), was updated to simulate various types of CA induced by acute exposures of the mixed field beams used for the experiments. The chromosomes were simulated by a polymer random walk algorithm with restrictions to their respective domains in the nucleus [1]. The stochastic dose to the nucleus was calculated with the code RITRACKS [2]. Irradiation of a target volume by a mixed field of ions was implemented within RITRACKs, and the fields of ions can be delivered over specific periods of time, allowing the simulation of dose-rate effects. Similarly, particles of various types and energies extracted from a pre-calculated spectra of galactic cosmic rays (GCR) can be used in RITRACKS. The number and spatial location of DSBs (DNA double-strand breaks) were calculated in BDSTRACKS using the simulated chromosomes and local (voxel) dose. Assuming that DSBs led to chromosome breaks, and simulating the rejoining of damaged chromosomes occurring during repair, BDSTRACKS produces the yield of various types of chromosome aberrations as a function of time (only final yields are presented). A comparison between experimental and simulation results will be shown.

Hada, M.↗
DOE OSTI2014

Analytical solution of the Poisson-Nernst-Planck equations for an electrochemical system close to electroneutrality

Single charge densities and the potential are used to describe models of electrochemical systems. These quantities can be calculated by solving a system of time dependent nonlinear coupled partial differential equations, the Poisson-Nernst-Planck equations. Assuming small deviations from the electroneutral equilibrium, the linearized and decoupled equations are solved for a radial symmetric geometry, which represents the interface between a cell and a sensor device. The densities and the potential are expressed by Fourier-Bessels series. The system considered has a ratio between the Debye-length and its geometric dimension on the order of 10{sup −4} so the Fourier-Bessel series can be approximated by elementary functions. The time development of the system is characterized by two time constants, τ{sub c} and τ{sub g}. The constant τ{sub c} describes the approach to the stationary state of the total charge and the potential. τ{sub c} is several orders of magnitude smaller than the geometry-dependent constant τ{sub g}, which is on the order of 10 ms characterizing the transition to the stationary state of the single ion densities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Detection and Interpretation of Long-lived X-Ray Quasi-periodic Pulsations in the X-class Solar Flare on 2013 May 14

Quasi-periodic pulsations (QPP) seen in the time derivative of the GOES soft X-ray light curves are analyzed for the X3.2 event on 2013 May 14. The pulsations are apparent for a total of at least two hours from the impulsive phase to well into the decay phase, with a total of 163 distinct pulses evident to the naked eye. A wavelet analysis shows that the characteristic timescale of these pulsations increases systematically from ∼25 s at 01:10 UT, the time of the GOES peak, to ∼100 s at 02:00 UT. A second “ridge” in the wavelet power spectrum, most likely associated with flaring emission from a different active region, shows an increase from ∼40 s at 01:40 UT to ∼100 s at 03:10 UT. We assume that the QPP that produced the first ridge result from vertical kink-mode oscillations of the newly formed loops following magnetic reconnection in the coronal current sheet. This allows us to estimate the magnetic field strength as a function of altitude given the density, loop length, and QPP timescale as functions of time determined from the GOES light curves and Ramaty High Energy Solar Spectroscopic Imager ( RHESSI ) images. The calculated magnetic field strength of the newly formed loops ranges from ∼500 G at an altitude of 24 Mm to a low value of ∼10 G at 60 Mm, in general agreement with the expected values at these altitudes. Fast sausage-mode oscillations are also discussed and cannot be ruled out as an alternate mechanism for producing the QPP.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Effective mie-scattering and CO2 absorption in the dust-laden Martian atmosphere and its impact on radiative-convective temperature changes in the lower scale heights

A time dependent computer model of radiative-convective-conductive heat transfer in the Martian ground-atmosphere system was refined by incorporating an intermediate line strength CO2 band absorption which together with the strong-and weak-line approximation closely simulated the radiative transmission through a vertically inhomogeneous stratification. About 33,000 CO2 lines were processed to cover the spectral range of solar and planetary radiation. Absorption by silicate dust particulates, was taken into consideration to study its impact on the ground-atmosphere temperature field as a function of time. This model was subsequently attuned to IRIS, IR-radiometric and S-band occultation data. Satisfactory simulations of the measured IRIS spectra were accomplished for the dust-free condition. In the case of variable dust loads, the simulations were sufficiently fair so that some inferences into the effect of dust on temperature were justified.

Pallmann, A. J.↗
NASA NTRS1997

Light and Gravity Effects on Circadian Rhythms of Rhesus Macaques

Temporal integration of a biological organism's physiological, behavioral and biochemical systems depends upon its circadian timing system. The endogenous period of this timing system is typically synchronized to the 24- hour day by environmental cues. The daily alternation of light and dark has long been known as one of the most potent environmental synchronizers influencing the circadian timing system. Alterations in the lighting environment (length or intensity of light exposure) can also affect the homeostatic state of the organism. A series of experiments was performed using rhesus monkeys with the objective of defining the fundamental properties of the circadian rhythm of body temperature. Three major experiments were performed in addition to several preliminary studies. These experiments explored 1.) the response of the rhesus body temperature rhythm to varying day length and light intensity; 2.) the response of the body temperature rhythm to light exposure as a function of time of day; and 3.) the characteristics of the metabolic heat production rhythm which is responsible for the daily cycle in body temperature. Results of these three completed experiments will be reported here. In addition, preliminary experiments were also performed in social entrainment of rhesus circadian rhythms and the properties of rhesus body temperature rhythms in constant conditions, where no external time cues were provided. Four adult male rhesus monkeys served as subjects in all experiments. All experiments were performed at the California Regional Primate Research Center. Each animal was implanted with a biotelemetry unit that measured deep body temperature. All surgeries were performed by a board certified veterinary surgeon under sterile conditions. The biotelemetry implants also provided an index of activity level in each animal. For metabolic heat production measurements, oxygen consumption and carbon dioxide production were measured and the caloric equivalent of these was calculated. Specific methodologies are described in detail.

Fuller, Charles↗
NASA NTRS2001

H-alpha Proxies for EIT Crinkles: Further Evidence for Pre-Flare "Breakout"-Type Activity in an Ejective Solar Eruption

We present Halpha observations from Big Bear Solar Observatory of an eruptive flare in NOAA AR 8210, occurring near 22:30 UT on 1998 May 1. Previously, using the EUV Imaging Telescope (EIT) on the SOHO spacecraft, we found that a pattern of transient, localized brightenings, which we call "EIT crinkles," appears in the neighborhood of the eruption near the time of flare onset. These EIT crinkles occur at a location in the active region well separated from the sheared core magnetic fields, which is where the most intense features of the eruption are concentrated. We also previously found that high-cadence images from the Soft X-ray Telescope (SXT) on Yohkoh indicate that soft X-ray intensity enhancements in the core begin after the start of the EIT crinkles. With the Halpha data, we find remote flare brightening counterparts to the EIT crinkles. Lightcurves as functions of time of various areas of the active region show that several of the remote flare brightenings undergo intensity increases prior to onset of principle brightenings in the core region, consistent with our earlier findings from EIT and SXT data. These timing relationships are consistent with the eruption onset mechanism known as the breakout model, introduced by Antiochos and colleagues, which proposes that eruptions begin with reconnection at a magnetic null high above the core region. Our observations are also consistent with other proposed mechanisms which do not involve early reconnection in the core region. As a corollary, our observations are not consistent with the so-called tether cutting models, which say that the eruption begins with reconnection in the core. The Halpha data further show that a filament in the core region becomes activated near the time of EIT crinkle onset, but little if any of the filament actually erupts, despite the presence of a halo Coronal Mass Ejection (CME) associated with this event.

Sterling, Alphonse C.↗
NASA NTRS2009

Microdosimeter Instrument (MIDN-II) for Personnel Dosimetry

This group of collaborators has been funded, by the National Space Biomedical Research Institute (NSBRI) and the U. S. National Aeronautics and Space Administration (NASA) as of February, 2009 to develop a prototypic dosimeter/microdosimeter instrument that measures in real time the radiation risk to personnel in a time varying radiation field of possibly unknown composition. The microdosimetric detectors will be solid-state devices more compact, more rugged than conventional microdosimetric proportional counters used previously in space and will consume less power. The proposed detector will be based on the heritage of the MIDN-MidSTAR microdosimeter launched on the MidSTAR spacecraft on March, 2007, the only solid-state microdosimeter ever flown in space. The system will be suitable for measurements in spacesuits, spacecrafts, remote rovers, or other dynamic or static environments. Measurements as a function of time will not only provide the instantaneous and average absorbed physical dose but also corresponding microdosimetric spectra, dose rates, the dose equivalents, and the dose equivalent rates. Values of the dose equivalents are used to establish relative risks for humans exposed to radiation and to determine regulatory limits. Because these parameters will be available to the astronauts and mission control in real time, these systems can be used not only to quantify exposures and limits, but also to allow appropriate actions to be taken to reduce radiation exposures and their consequences. Earlier prototypes have been used successfully to characterize beams of energetic protons and heavier ions including carbon, oxygen, silicon, titanium, and iron at the NASA Space Radiation Laboratory at Brookhaven National Laboratory. The initial design will be presented and discussed and spectra will be presented.

Pisacane, V. L.↗
DOE OSTI2015

Intersystem-crossing and phosphorescence rates in fac-Ir{sup III}(ppy){sub 3}: A theoretical study involving multi-reference configuration interaction wavefunctions

We have employed combined density functional theory and multi-reference configuration interaction methods including spin–orbit coupling (SOC) effects to investigate the photophysics of the green phosphorescent emitter fac-tris-(2-phenylpyridine)iridium (fac-Ir(ppy){sub 3}). A critical evaluation of our quantum chemical approaches shows that a perturbational treatment of SOC is the method of choice for computing the UV/Vis spectrum of this heavy transition metal complex while multi-reference spin–orbit configuration interaction is preferable for calculating the phosphorescence rates. The particular choice of the spin–orbit interaction operator is found to be of minor importance. Intersystem crossing (ISC) rates have been determined by Fourier transformation of the time correlation function of the transition including Dushinsky rotations. In the electronic ground state, fac-Ir(ppy){sub 3} is C{sub 3} symmetric. The calculated UV/Vis spectrum is in excellent agreement with experiment. The effect of SOC is particularly pronounced for the metal-to-ligand charge transfer (MLCT) band in the visible region of the absorption spectrum which does not only extend its spectral onset towards longer wavelengths but also experiences a blue shift of its maximum. Pseudo-Jahn-Teller interaction leads to asymmetric coordinate displacements in the lowest MLCT states. Substantial electronic SOC and a small energy gap make ISC an ultrafast process in fac-Ir(ppy){sub 3}. For the S{sub 1}↝T{sub 1} non-radiative transition, we compute a rate constant of k{sub ISC} = 6.9 × 10{sup 12} s{sup −1} which exceeds the rate constant of radiative decay to the electronic ground state by more than six orders of magnitude, in agreement with the experimental observation of a subpicosecond ISC process and a triplet quantum yield close to unity. As a consequence of the geometric distortion in the T{sub 1} state, the T{sub 1} → S{sub 0} transition densities are localized on one of the phenylpyridyl moieties. In our best quantum chemical model, we obtain phosphorescence decay times of 264 μs, 13 μs, and 0.9 μs, respectively, for the T{sub 1,I}, T{sub 1,II}, and T{sub 1,III} fine-structure levels in dichloromethane (DCM) solution. In addition to reproducing the correct orders of magnitude for the individual phosphorescence emission probabilities, our theoretical study gives insight into the underlying mechanisms. In terms of intensity borrowing from spin-allowed transitions, the low emission probability of the T{sub 1,I} substate is caused by the mutual cancellation of contributions from several singlet states to the total transition dipole moment. Their contributions do not cancel but add up in case of the much faster T{sub 1,III} → S{sub 0} emission while the T{sub 1,II} → S{sub 0} emission is dominated by intensity borrowing from a single spin-allowed process, i.e., the S{sub 2} → S{sub 0} transition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Structure of Asphaltenes during Precipitation Investigated by Ultra-Small-Angle X-ray Scattering

Time-resolved size and structure measurements of asphaltenes while in the process of precipitating were monitored for the first time using ultra-small-angle X-ray scattering. The results revealed that asphaltenes precipitating from a heptane–toluene mixture demonstrate a hierarchical structure of an asphaltene-rich phase (e.g., droplet) that further agglomerates into fractal flocs. The fractal flocs that form by the agglomeration of the asphaltene-rich phase are what is commonly detected by optical microscopy above the precipitation detection point. The surface of the asphaltene-rich phase is initially rough and transitions to a smooth interface, as would be expected for a highly viscous liquid. Simultaneous small-angle X-ray scattering measurements were also performed to investigate the structure of soluble asphaltenes, providing comprehensive structural characterization from the nanometer- to micrometer-length scales as a function of time. Further, the results demonstrate that the size and concentration of asphaltenes remaining in solution (e.g., soluble asphaltenes) do not change during precipitation, whereas the structure of insoluble asphaltenes varies. The ability to measure the properties of asphaltenes as they undergo precipitation opens new opportunities for understanding the fundamental mechanisms of asphaltene deposition and aggregation and the impact of chemical inhibitors to alter these processes. Here, the universality of these conclusions and how specific properties vary as a function of asphaltene source and solution properties can provide valuable insight into asphaltene behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Machine learning in materials research: Developments over the last decade and challenges for the future

The number of studies that apply machine learning (ML) to materials science has been growing at a rate of approximately 1.67 times per year over the past decade. In this review, I examine this growth in various contexts. First, I present an analysis of the most commonly used tools (software, databases, materials science methods, and ML methods) used within papers that apply ML to materials science. The analysis demonstrates that despite the growth of deep learning techniques, the use of classical machine learning is still dominant as a whole. It also demonstrates how new research can effectively build upon past research, particular in the domain of ML models trained on density functional theory calculation data. Next, I present the progression of best scores as a function of time on the matbench materials science benchmark for formation enthalpy prediction. In particular, a dramatic improvement of 7 times reduction in error is obtained when progressing from feature-based methods that use conventional ML (random forest, support vector regression, etc.) to the use of graph neural network techniques. Finally, I provide views on future challenges and opportunities, focusing on data size and complexity, extrapolation, interpretation, access, and relevance.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A STRONGLY COUPLED REACTOR CORE ISOLATION COOLING SYSTEM MODEL FOR EXTENDED STATION BLACK-OUT ANALYSES

The reactor core isolation cooling (RCIC) system in a boiling water reactor (BWR) provides makeup cooling water to the reactor pressure vessel (RPV) when the main steam lines are isolated and the normal supply of water to the reactor vessel is lost. The RCIC system operates independently of AC power, service air, or external cooling water systems. The only required external energy source is from the battery to maintain the logic circuits to control the opening and/or closure of valves in the RCIC systems in order to control the RPV water level by shutting down the RCIC pump to avoid overfilling the RPV and flooding the steam line to the RCIC turbine. It is generally considered in almost all the existing station black-out accidents (SBO) analyses that loss of the DC power would result in overfilling the steam line and allowing liquid water to flow into the RCIC turbine, where it is assumed that the turbine would then be disabled. This behavior, however, was not observed in the Fukushima Daiichi accidents, where the Unit 2 RCIC functioned without DC power for nearly three days. Therefore, more detailed mechanistic models for RCIC system components are needed to understand the extended SBO for BWRs. As part of the effort to develop the next generation reactor system safety analysis code RELAP-7, we have developed a strongly coupled RCIC system model, which consists of a turbine model, a pump model, a check valve model, a wet well model, and their coupling models. Unlike the traditional SBO simulations where mass flow rates are typically given in the input file through time dependent functions, the real mass flow rates through the turbine and the pump loops in our model are dynamically calculated according to conservation laws and turbine/pump operation curves. A simplified SBO demonstration RELAP-7 model with this RCIC model has been successfully developed. The demonstration model includes the major components for the primary system of a BWR, as well as the safety system components such as the safety relief valve (SRV), the RCIC system, the wet well, and the dry well. The results show reasonable system behaviors while exhibiting rich dynamics such as variable flow rates through RCIC turbine and pump during the SBO transient. The model has the potential to resolve the Fukushima RCIC mystery after adding the off-design two-phase turbine operation model and other additional improvements.

RCIC↗
NASA NTRS2001

H(alpha) Proxies for EIT Crinkles: Further Evidence for Preflare "Breakout"-Type Activity in an Ejective Solar Eruption

We present H(alpha) observations from Big Bear Solar Observatory of an eruptive flare in NOAA Active Region 8210, occurring near 22:30 UT on 1998 May 1. Previously, using the Extreme Ultraviolet Imaging Telescope (EIT) on the Solar and Heliospheric Observatory (SOHO) spacecraft, we found that a pattern of transient, localized brightenings, which we call 'EIT crinkles,' appears in the neighborhood of the eruption near the time of flare onset. These EIT crinkles occur at a location in the active region well separated from the sheared core magnetic fields, which is where the most intense features of the eruption are concentrated. We also previously found that high-cadence images from the Soft X-ray Telescope (SXT) on Yohkoh indicate that soft X-ray intensity enhancements in the core begin after the start of the EIT crinkles. With the H(alpha) data, we find remote flare brightening counterparts to the EIT crinkles. Light curves as functions of time of various areas of the active region show that several of the remote flare brightenings undergo intensity increases prior to the onset of principal brightenings in the core region, consistent with our earlier findings from EIT and SXT data. These timing relationships are consistent with the eruption onset mechanism known as the breakout model, introduced by Antiochos and colleagues, which proposes that eruptions begin with reconnection at a magnetic null high above the core region. Our observations are also consistent with other proposed mechanisms that do not involve early reconnection in the core region. As a corollary, our observations are not consistent with the so-called tether-cutting models, which say that the eruption begins with reconnection in the core. The H(alpha) data further show that a filament in the core region becomes activated near the time of EIT crinkle onset, but little if any of the filament actually erupts, despite the presence of a halo coronal mass ejection (CME) associated with this event.

Sterling, Alphonse C.↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Probability of Loss of Assured Safety in Systems with Multiple Time-Dependent Failure Modes: Incorporation of Delayed Link Failure in the Presence of Aleatory Uncertainty

Probability of loss of assured safety (PLOAS) is modeled for weak link (WL)/strong link (SL) systems in which one or more WLs or SLs could potentially degrade into a precursor condition to link failure that will be followed by an actual link failure after some amount of elapsed time. The descriptor loss of assured safety (LOAS) is used because failure of the WL system places the entire system in an inoperable configuration while failure of the SL system before failure of the WL system, although undesirable, does not necessarily result in an unintended operation of the entire system. Thus, safety is “assured” by failure of the WL system before failure of the SL system. Here, the following topics are considered: (i) Definition of precursor occurrence time cumulative distribution functions (CDFs) for individual WLs and SLs, (ii) Formal representation, approximation and illustration of PLOAS with (a) constant delay times, (b) aleatory uncertainty in delay times, and (c) delay times defined by functions of link properties at occurrence times for link failure precursors, and (iii) Procedures for the verification of PLOAS calculations for the three indicated definitions of delayed link failure.

58 GEOSCIENCES↗
DOE OSTI2017

Shielding Analysis Capability of UNF-ST&DARDS

The Used Nuclear Fuel Storage, Transportation and Disposal Analysis Resource and Data System (UNF-ST&DARDS) is used to perform dose rate calculations for spent nuclear fuel (SNF) transportation packages based on the actual physical and nuclear characteristics (i.e., assembly design, burnup, initial enrichment, and cooling time) of the as-loaded SNF. Nuclear fuel data, transportation package model templates, and SNF canister loading map information residing within the tool facilitate automated generation of SCALE input files for radiation source term and dose rate calculations. Transportation package specific models developed for UNF-ST&DARDS dose rate analyses are described in detail. UNF-ST&DARDS dose rate analyses were performed in this paper for over 400 SNF canisters from 16 sites in their designated transportation casks. For simplicity, representative dose rate calculation results are presented as a function of time (i.e., selected calendar years between 2020 and 2100) for 73 SNF canisters in dry storage at four sites. For these canisters, the projected maximum dose rate values at 2 m from the lateral surfaces of the vehicle under normal conditions of transport (NCT) would vary between 1.9 and 11.5 mrem/h in 2020. Five SNF canisters will exceed the regulatory dose rate limit of 10 mrem/h at 2 m in 2020, and the analyzed SNF canisters will comply with regulatory dose rate limits by 2030. An analysis of the impact on the dose rate of fuel failure and reconfiguration during transportation indicated that the maximum dose rate for hypothetical accident conditions will be unaffected, and the NCT maximum dose rate at 2 m would have a maximum increase by a factor of 1.7 for a representative pressurized water reactor package and by a factor of 2.6 for a representative boiling water reactor package relative to intact fuel. Finally, analysis of the actual heat loading and the dose rate at 2 m from the lateral surface of the vehicle for the five SNF canisters exceeding the NCT regulatory dose rate limit of 10 mrem/h in 2020 showed that the dose rate could be more limiting with respect to regulatory requirements than the heat loading; i.e., the canister transportability date may be evaluated based on the transportation package’s external dose rate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Development of Adaptive Model Refinement (AMoR) for Multiphysics and Multifidelity Problems

This project investigated the development and utilization of Adaptive Model Refinement (AMoR) for nuclear systems simulation applications. AMoR refers to utilization of several models of physical phenomena which differ in prediction fidelity. If the highest fidelity model is judged to always provide or exceeded the desired fidelity, than if one can determine the difference in a Quantity of Interest (QoI) between the highest fidelity model and lower fidelity models, one could utilize the fidelity model that would just provide the magnitude of the QoI desired. Assuming lower fidelity models require less computational resources, in this manner computational efficiency can be realized provided the QoI value can be accurately and efficiently evaluated. This work utilized Generalized Perturbation Theory (GPT) to evaluate the QoI, by convoluting the GPT solution with the residual of the highest fidelity model determined using the solution from lower fidelity models. Specifically, a reactor core neutronics problem and thermal-hydraulics problem were studied to develop and utilize AMoR. The highest fidelity neutronics model was based upon the 3D space-time, two-group, nodal diffusion equations as solved in the NESTLE computer code. Added to the NESTLE code was the ability to determine the time-dependent GPT neutron flux. The lower fidelity neutronics model was based upon the point kinetics equations along with utilization of a prolongation operator to determine the 3D space-time, two-group flux. The highest fidelity thermal-hydraulics model was based upon the space-time equations governing fluid flow in a closed channel around a heat generating fuel rod. The Homogenous Equilibrium Mixture (HEM) model was used for the fluid and Finite Difference Method was applied to both the coolant and fuel pin energy conservation equations. The lower fidelity thermal-hydraulic model was based upon the same equations as used for the highest fidelity model but now with coarse spatial meshing, corrected somewhat by employing effective fuel heat conduction values. The effectiveness of switching between the highest fidelity model and lower fidelity model as a function of time was assessed using the neutronics problem. Based upon work completed to date, one concludes that the time switching is effective in annealing out differences between the highest and lower fidelity solutions. The effectiveness of using a lower fidelity GPT solution, along with a prolongation operator, to estimate the QoI was also assessed. The utilization of a lower fidelity GPT solution was done in an attempt to avoid the high computational burden associated with solving for the highest fidelity GPT solution. Based upon work completed to date, one concludes that the lower fidelity adjoint solution is not sufficiently accurate with regard to estimating the QoI; however, a formulation has been revealed that may provide a path for addressing this shortcoming.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1982

A system for measuring methane fluxes from inland and coastal wetland environments

A CH4 measuring system with a sensitivity of 0.01 ppmv and a response time of 1 sec is described. Pulses are separable from moving and stationary readings, the device floats, and air velocity in the chamber can be controlled in the range 0-4.4 m/sec. The system is set on the water so that controlled speed air from beneath can enter an inner chamber and be pumped through a gas filter correlation system (GFC) before being returned to the chamber. The GFC mechanism yields an absorption signal in the form of a dc voltage, which a processor orders printed on a strip chart as a function of time. Detection of the CH4 amounts is through IR absorption, and the sensitivity of the instrument has been measured at .002-.051 g/sq m per day with air flows in the range 0.90-4.40 m/sec. The unit can be modified to detect CO and CO2, and further applications for low plant canopies and land surfaces are indicated.

Sebacher, D. I.↗
NASA NTRS2017

Detection and Interpretation of Long-lived X-Ray Quasi-Periodic Pulsations in the X-Class Solar Flare on 2013 May 14

Quasi-periodic pulsations (QPP) seen in the time derivative of the GOES soft X-ray light curves are analyzed for the X3.2 event on 2013 May 14. The pulsations are apparent for a total of at least two hours from the impulsive phase to well into the decay phase, with a total of 163 distinct pulses evident to the naked eye. A wavelet analysis shows that the characteristic timescale of these pulsations increases systematically from approx. 25 s at 01:10 UT, the time of the GOES peak, to approx. 100 s at 02:00 UT. A second ridge in the wavelet power spectrum, most likely associated with flaring emission from a different active region, shows an increase from approx. 40 s at 01:40 UT to approx. 100 s at 03:10 UT. We assume that the QPP that produced the first ridge result from vertical kink-mode oscillations of the newly formed loops following magnetic reconnection in the coronal current sheet. This allows us to estimate the magnetic field strength as a function of altitude given the density, loop length, and QPP timescale as functions of time determined from the GOES light curves and Ramaty High Energy Solar Spectroscopic Imager (RHESSI) images. The calculated magnetic field strength of the newly formed loops ranges from approx. 500 G at an altitude of 24Mm to a low value of approx. 10 G at 60 Mm, in general agreement with the expected values at these altitudes. Fast sausage-mode oscillations are also discussed and cannot be ruled out as an alternate mechanism for producing the QPP.

Dennis, Brian R.↗