Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic parameter estimation”

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 883 records · Page 49

Investigating the Relativistic Motion of the Stars Near the Supermassive Black Hole in the Galactic Center

The S-star cluster in the Galactic center allows us to study the physics close to a supermassive black hole, including distinctive dynamical tests of general relativity. Our best estimates for the mass of and the distance to Sgr A* using the three stars with the shortest period (S2, S38, and S55/S0-102) and Newtonian models are M {sub BH} = (4.15 ± 0.13 ± 0.57) × 10{sup 6} M {sub ⊙} and R {sub 0} = 8.19 ± 0.11 ± 0.34 kpc. Additionally, we aim at a new and practical method to investigate the relativistic orbits of stars in the gravitational field near Sgr A*. We use a first-order post-Newtonian approximation to calculate the stellar orbits with a broad range of periapse distance r {sub p} . We present a method that employs the changes in orbital elements derived from elliptical fits to different sections of the orbit. These changes are correlated with the relativistic parameter defined as ϒ ≡ r {sub s} / r {sub p} (with r {sub s} being the Schwarzschild radius) and can be used to derive ϒ from observational data. For S2 we find a value of ϒ = 0.00088 ± 0.00080, which is consistent, within the uncertainty, with the expected value of ϒ = 0.00065 derived from M {sub BH} and the orbit of S2. We argue that the derived quantity is unlikely to be dominated by perturbing influences such as noise on the derived stellar positions, field rotation, and drifts in black hole mass.

79 ASTRONOMY AND ASTROPHYSICS↗

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗

CONSTRAINT-INDEPENDENT CONSTANT CTOA DETERMINATION FOR DUCTILE STABLE CRACK GROWTH

Crack tip opening angle (CTOA) has been used as a reliable fracture toughness parameter for decades to characterize stable ductile crack growth for thin-walled aerospace structures in the low-constraint conditions. Recently, the CTOA parameter was also applied to the pipeline industry, and a CTOA test standard ASTM E3039 was thus developed for testing a critical constant CTOA. Research showed that the constant CTOA can reasonably describe fracture toughness required to arrest a dynamic crack propagation for a modern gas pipeline. However, the CTOA fracture criterion requires constraint-independent CTOA toughness against stable ductile crack growth. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen for CTOA testing. Since a shallow crack is used, DWTT measured CTOA may depend on constraint level at the crack tip. To understand if it is the case, this paper evaluates the critical CTOA for a set of fracture toughness tests on single edge notched bend (SENB) specimens with shallow and deep cracks based on four CTOA estimation models. In which, the Ln(P)-LLD linear fit model is similar to that used by ASTM E3039 in the CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel are considered in the CTOA evaluation. The results show that the four CTOA models can determine a crack size-independent constant CTOA over stable ductile crack growth for the SENB specimens. As a result, CTOA determined by ASTM E3039 is constraintindependent and transferable to use for an actual crack propagating in a gas pipeline.

Zhu, Xian-Kui↗

REFLECTION OF PROPAGATING SLOW MAGNETO-ACOUSTIC WAVES IN HOT CORONAL LOOPS: MULTI-INSTRUMENT OBSERVATIONS AND NUMERICAL MODELING

Slow MHD waves are important tools for understanding coronal structures and dynamics. In this paper, we report a number of observations from the X-Ray Telescope (XRT) on board HINODE and Solar Dynamic Observatory /Atmospheric Imaging Assembly (AIA) of reflecting longitudinal waves in hot coronal loops. To our knowledge, this is the first report of this kind as seen from the XRT and simultaneously with the AIA. The wave appears after a micro-flare occurs at one of the footpoints. We estimate the density and temperature of the loop plasma by performing differential emission measure (DEM) analysis on the AIA image sequence. The estimated speed of propagation is comparable to or lower than the local sound speed, suggesting it to be a propagating slow wave. The intensity perturbation amplitude, in every case, falls very rapidly as the perturbation moves along the loop and eventually vanishes after one or more reflections. To check the consistency of such reflection signatures with the obtained loop parameters, we perform a 2.5D MHD simulation, which uses the parameters obtained from our observation as inputs, and perform forward modeling to synthesize AIA 94 Å images. Analyzing the synthesized images, we obtain the same properties of the observables as for the real observation. From the analysis we conclude that a footpoint heating can generate a slow wave which then reflects back and forth in the coronal loop before fading. Our analysis of the simulated data shows that the main agent for this damping is anisotropic thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗

Conceptual design of hollow electron lenses for beam halo control in the Large Hadron Collider

Collimation with hollow electron beams is a technique for halo control in high-power hadron beams. It is based on an electron beam (possibly pulsed or modulated in intensity) guided by strong axial magnetic fields which overlaps with the circulating beam in a short section of the ring. The concept was tested experimentally at the Fermilab Tevatron collider using a hollow electron gun installed in one of the Tevatron electron lenses. We are proposing a conceptual design for applying this technique to the Large Hadron Collider at CERN. A prototype hollow electron gun for the LHC was built and tested. The expected performance of the hollow electron beam collimator was based on Tevatron experiments and on numerical tracking simulations. Halo removal rates and enhancements of halo diffusivity were estimated as a function of beam and lattice parameters. Proton beam core lifetimes and emittance growth rates were checked to ensure that undesired effects were suppressed. Hardware specifications were based on the Tevatron devices and on preliminary engineering integration studies in the LHC machine. Required resources and a possible timeline were also outlined, together with a brief discussion of alternative halo-removal schemes and of other possible uses of electron lenses to improve the performance of the LHC.

43 PARTICLE ACCELERATORS↗

A Dynamic Simulation of the Indirect Land Use Implications of Recent Biofuel Production and Use in the United States.

The global indirect land use change (ILUC) implications of biofuel use in the United States of America (USA) from 2001 to 2010 are evaluated with a dynamic general equilibrium model. The effects of biofuels production on agricultural land area vary by year; from a net expansion of 0.17 ha per 1000 gallons produced (2002) to a net contraction of -0.13 ha per 1000 gallons (2018) in Case 1 of our simulation. In accordance with the general narrative about the implications of biofuel policy, agricultural land area increased in many regions of the world. However, oil-export dependent economies experienced agricultural land contraction because of reductions in their revenues. Reducing crude oil imports is a major goal of biofuel policy, but the land use change implications have received little attention in the literature. Simulations evaluating the effects of doubling supply elasticities for land and fossil resources show that these parameters can significantly influence the land use change estimates. Therefore, research that provides empirically-based and spatially-detailed agricultural land-supply curves and capability to project future fossil energy prices is critical for improving estimates of the effects of biofuel policy on land use.

Oladosu, Gbadebo A.↗

Constraining Bedrock Groundwater Residence Times in a Mountain System With Environmental Tracer Observations and Bayesian Uncertainty Quantification

Groundwater residence time distributions provide fundamental insights on the hydrological processes within watersheds. Yet, observations that can constrain groundwater residence times over broad timescales remain scarce in mountain catchment studies. We use environmental tracers (CFC-12, SF 6 , 3 H, and 4 He) to investigate groundwater residence times along a hillslope in the East River Watershed, Colorado, USA. We develop a Bayesian inference framework that applies a Markov-chain Monte Carlo (MCMC) approach to estimate noble gas recharge temperature, elevation, and excess-air parameters and the resulting environmental tracer concentrations. MCMC is then used to propagate the environmental tracer uncertainties to estimates of groundwater mean residence times inferred with lumped parameter models. All samples contain 3 H, CFC-12, and SF 6 in addition to terrigenic 4 He, suggesting a mixture of water characterized by modern and premodern residence times. 4He exponential mean residence times range from hundreds of years at the upslope well to thousands of years at the toe-slope well assuming average crustal production rates. We find that binary mixing residence time distributions with separate young and old mixing fractions are needed to predict the 4 He, CFC-12, SF 6 , and 3 H observations, supporting the importance of flow path mixing in this bedrock system. Our findings that the fractured bedrock hosts groundwater with a mixture of residence times ranging from decades to millennia suggest variable recharge dynamics and flow path mixing along the hillslope and highlight the importance of characterizing groundwater systems with observations that are sensitive to transport over a broad range of residence times.

54 ENVIRONMENTAL SCIENCES↗

KINEMATICAL AND CHEMICAL VERTICAL STRUCTURE OF THE GALACTIC THICK DISK. II. A LACK OF DARK MATTER IN THE SOLAR NEIGHBORHOOD

We estimated the dynamical surface mass density {Sigma} at the solar position between Z = 1.5 and 4 kpc from the Galactic plane, as inferred from the kinematics of thick disk stars. The formulation is exact within the limit of validity of a few basic assumptions. The resulting trend of {Sigma}(Z) matches the expectations of visible mass alone, and no dark component is required to account for the observations. We extrapolate a dark matter (DM) density in the solar neighborhood of 0 {+-} 1 mM{sub Sun} pc{sup -3}, and all the current models of a spherical DM halo are excluded at a confidence level higher than 4{sigma}. A detailed analysis reveals that a small amount of DM is allowed in the volume under study by the change of some input parameter or hypothesis, but not enough to match the expectations of the models, except under an exotic combination of non-standard assumptions. Identical results are obtained when repeating the calculation with kinematical measurements available in the literature. We demonstrate that a DM halo would be detected by our method, and therefore the results have no straightforward interpretation. Only the presence of a highly prolate (flattening q > 2) DM halo can be reconciled with the observations, but this is highly unlikely in {Lambda}CDM models. The results challenge the current understanding of the spatial distribution and nature of the Galactic DM. In particular, our results may indicate that any direct DM detection experiment is doomed to fail if the local density of the target particles is negligible.

79 ASTRONOMY AND ASTROPHYSICS↗

APOSTLE: LONGTERM TRANSIT MONITORING AND STABILITY ANALYSIS OF XO-2b

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 10 transits of XO-2b over a period of 3 yr. We present measurements that confirm previous estimates of system parameters like the normalized semi-major axis (a/R{sub *}), stellar density ({rho}{sub *}), impact parameter (b), and orbital inclination (i{sub orb}). Our errors on system parameters like a/R{sub *} and {rho}{sub *} have improved by {approx}40% compared to previous best ground-based measurements. Our study of the transit times show no evidence for transit timing variations (TTVs) and we are able to rule out co-planar companions with masses {>=}0.20 M{sub Circled-Plus} in low order mean motion resonance with XO-2b. We also explored the stability of the XO-2 system given various orbital configurations of a hypothetical planet near the 2:1 mean motion resonance. We find that a wide range of orbits (including Earth-mass perturbers) are both dynamically stable and produce observable TTVs. We find that up to 51% of our stable simulations show TTVs that are smaller than the typical transit timing errors ({approx}20 s) measured for XO-2b, and hence remain undetectable.

79 ASTRONOMY AND ASTROPHYSICS↗

Scoping review of Culex mosquito life history trait heterogeneity in response to temperature

Abstract Background Mosquitoes in the genus Culex are primary vectors in the US for West Nile virus (WNV) and other arboviruses. Climatic drivers such as temperature have differential effects on species-specific changes in mosquito range, distribution, and abundance, posing challenges for population modeling, disease forecasting, and subsequent public health decisions. Understanding these differences in underlying biological dynamics is crucial in the face of climate change. Methods We collected empirical data on thermal response for immature development rate, egg viability, oviposition, survival to adulthood, and adult lifespan for Culex pipiens, Cx. quinquefasciatus, Cx. tarsalis , and Cx. restuans from existing literature according to the PRISMA scoping review guidelines. Results We observed linear relationships with temperature for development rate and lifespan, and nonlinear relationships for survival and egg viability, with underlying variation between species. Optimal ranges and critical minima and maxima also appeared varied. To illustrate how model output can change with experimental input data from individual Culex species, we applied a modified equation for temperature-dependent mosquito type reproduction number for endemic spread of WNV among mosquitoes and observed different effects. Conclusions Current models often input theoretical parameters estimated from a single vector species; we show the need to implement the real-world heterogeneity in thermal response between species and present a useful data resource for researchers working toward that goal. Graphical Abstract

59 BASIC BIOLOGICAL SCIENCES↗

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗

Developing and Testing a Novel Stochastic Ice Microphysics Parameterization for Cloud and Climate Models Using ARM Field Campaign Data (Final Progress Report)

The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗

Dynamics of the chemo-optical response of supported films of nematic liquid crystals

We report a study that elucidates key physicochemical phenomena that underlie the dynamics of adsorbate-induced ordering transitions in micrometer-thick films of nematic liquid crystals (LCs) supported on chemically functionalized surfaces. By flowing a gas containing dimethyl methylphosphonate (DMMP) over a film of nematic 4-cyano-4'-pentylbiphenyl (5CB) supported on a surface decorated with aluminum perchlorate, we quantify the dynamics of the optical response (intensity of transmitted light) of the LC film as a function of key experimental parameters, including linear flow rate of the gas over the surface of the film and the concentration of the DMMP within the gas. Building from prior infrared spectroscopic analysis revealing that the optical response of the LC is due to an ordering transition triggered by competitive ligand-exchange reactions at the chemically functionalized surface (i.e., between the phosphonate and nitrile groups of the DMMP and LC, respectively, with the aluminum perchlorate-decorated surface), we interpret the measurements reported herein within the framework of a simple model of mass transport. The model reveals that the rate-limiting process underlying the response of the LC is the mass transport of the DMMP across the concentration boundary layer formed on the vapor side of the LC–vapor interface, and not the mass transport of DMMP across the LC film nor the intrinsic ligand exchange kinetics. The analysis also reveals that the response of the LC is triggered once the concentration of DMMP within the LC film rises to approximately 0.25 mM. By flowing a stream of gas free of DMMP over a LC film pre-exposed to DMMP, we established that the relaxation of the LC film to the initial state is also largely controlled by mass transport of the DMMP across the vapor side of the LC–vapor interface. Analysis of the relaxation of the system leads to an estimate of the sensitivity of the supported LC film to the concentration of DMMP in the gas phase (approximately 150 ppb for the non-optimized system reported in this paper). These results and others reported in this paper identify key physical phenomena that dictate the dynamics of adsorbate-induced ordering transitions in LC systems, and provide important guidance to the design of LCs as chemo-responsive materials for use in sensors.

dimethyl methylphosphonate↗

Thermoelectric DC conductivities in hyperscaling violating Lifshitz theories

We analytically compute the thermoelectric conductivities at zero frequency (DC) in the holographic dual of a four dimensional Einstein-Maxwell-Axion-Dilaton theory that admits a class of asymptotically hyperscaling violating Lifshitz backgrounds with a dynamical exponent z and hyperscaling violating parameter θ. We show that the heat current in the dual Lifshitz theory involves the energy flux, which is an irrelevant operator for z > 1. The linearized fluctuations relevant for computing the thermoelectric conductivities turn on a source for this irrelevant operator, leading to several novel and non-trivial aspects in the holographic renormalization procedure and the identification of the physical observables in the dual theory. Moreover, imposing Dirichlet or Neumann boundary conditions on the spatial components of one of the two Maxwell fields present leads to different thermoelectric conductivities. Thus, Dirichlet boundary conditions reproduce the thermoelectric DC conductivities obtained from the near horizon analysis of Donos and Gauntlett, while Neumann boundary conditions result in a new set of DC conductivities. We make preliminary analytical estimates for the temperature behavior of the thermoelectric matrix in appropriate regions of parameter space. In particular, at large temperatures we find that the only case which could lead to a linear resistivity ρ ~ T corresponds to z = 4/3.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SU-G-BRA-08: Diaphragm Motion Tracking Based On KV CBCT Projections with a Constrained Linear Regression Optimization

Purpose: To develop a novel strategy to extract the respiratory motion of the thoracic diaphragm from kilovoltage cone beam computed tomography (CBCT) projections by a constrained linear regression optimization technique. Methods: A parabolic function was identified as the geometric model and was employed to fit the shape of the diaphragm on the CBCT projections. The search was initialized by five manually placed seeds on a pre-selected projection image. Temporal redundancies, the enabling phenomenology in video compression and encoding techniques, inherent in the dynamic properties of the diaphragm motion together with the geometrical shape of the diaphragm boundary and the associated algebraic constraint that significantly reduced the searching space of viable parabolic parameters was integrated, which can be effectively optimized by a constrained linear regression approach on the subsequent projections. The innovative algebraic constraints stipulating the kinetic range of the motion and the spatial constraint preventing any unphysical deviations was able to obtain the optimal contour of the diaphragm with minimal initialization. The algorithm was assessed by a fluoroscopic movie acquired at anteriorposterior fixed direction and kilovoltage CBCT projection image sets from four lung and two liver patients. The automatic tracing by the proposed algorithm and manual tracking by a human operator were compared in both space and frequency domains. Results: The error between the estimated and manual detections for the fluoroscopic movie was 0.54mm with standard deviation (SD) of 0.45mm, while the average error for the CBCT projections was 0.79mm with SD of 0.64mm for all enrolled patients. The submillimeter accuracy outcome exhibits the promise of the proposed constrained linear regression approach to track the diaphragm motion on rotational projection images. Conclusion: The new algorithm will provide a potential solution to rendering diaphragm motion and ultimately improving tumor motion management for radiation therapy of cancer patients.

60 APPLIED LIFE SCIENCES↗

Quantifying the Role of Snowmelt in Stream Discharge in an Alaskan Watershed: An Analysis Using a Spatially Distributed Surface Hydrology Model

Here, this study uses a spatially distributed surface hydrology model to investigate the role of snowmelt in stream discharge for the Tanana Flats Basin in interior Alaska. The Parameter ESTimation code is used to calibrate the model with observed stream discharge data. The model was further evaluated using remote sensing-based snow cover product and in situ snowpack water equivalent (SWE) observations. A 36 year (1980–2015) U.S. Geological Survey Precipitation-Runoff Modeling System simulation shows (1) the monthly stream discharge from the Tanana Flats Basin in April decreased by 44%; (2) snow cover area at high altitudes (above 2000 m) decreased in summer, both SWE and snowmelt also decreased significantly, especially in spring; (3) the timings of snowmelt onset and ending shifted by 2 (earlier) and 5 (later) days per decade, respectively; and (4) snowmelt accounts for 40% of the annual stream discharge. This study provides a quantitative tool to investigating hydrological systems considering the impacts of snow dynamics in cold regions. Additionally, this study suggests that future warming will further decrease snow coverage, advance snow melting time, and hereafter change the stream discharge dynamics in the Arctic.

58 GEOSCIENCES↗

Simulation methodology for fuel assembly drop accident during handling

Fresh and irradiated fuel assembly handling is a routine part of Nuclear Power Plant activities. Fuel Assembly handling occurs during fresh fuel receipt, core refueling, spent fuel assemblies transfer from spent fuel pool to the storage containers, etc. Westinghouse, ENUSA and KEPCO-NF have collaborated in the technical area of Fuel Assembly Drop Accident (FADA) Simulation to update the approach for estimating the number of failed fuel rods as a consequence of a fresh or irradiated fuel assembly drop, which might occur during a handling accident in a PWR-type reactor. The proposed FADA simulation methodology starts with an assessment of fuel assembly impact velocity at the final stage of the drop in water or air. The vertical, horizontal or inclined drop directions are considered. This velocity is an input to determine a fuel assembly structural response during impact which defines the loading function for the fuel assembly components including spacer grids. The dynamic load on the fuel rod is calculated from the detailed spacer grid model and the corresponding loading function. Finally, the load from the spacer grid spring and dimples are applied to the fuel rod model to determine the cladding parameters for use in compliance to fuel rod cladding failure criteria. A combination of first principal, finite element analysis (FEA) and computational fluid dynamic (CFD) methods has been utilized to create a technically sound methodology. A number of specific and traditional tests for the fuel assembly and its components has been used to confirm that the proposed approaches which describe the physical phenomena related to the fuel assembly drop accident are reasonably conservative.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗