Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interval analysis”

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,495 records · Page 83

End-To-End Uncertainty Quantification with Analytical Derivatives for Design Under Uncertainty

Uncertainty quantification (UQ) is a rapidly growing and evolving discipline, especially within the aerospace community. Performing analysis with UQ can provide decision makers with a wealth of information about a candidate design. However, the value of UQ is fully realized when the information gained during UQ analysis is leveraged in a feedback loop of a design optimization process, often referred to as design under uncertainty. Although design under uncertainty can be a powerful risk mitigation technique, there are a number of roadblocks that prevent its implementation. Two primary factors are computational costs and added complexity of the analysis. High fidelity simulations on the order tens of uncertain variables quickly become computationally infeasible. Also, implementing UQ into an existing multidisciplinary design and optimization (MDO) process often requires extensive knowledge of the UQ methods and careful treatment of the problem formulation. The objective of this work is to address these two primary roadblocks and enable practitioners to efficiently perform design under uncertainty with limited knowledge of the UQ discipline. Methods outlined in this paper demonstrate MDO incorporating UQ into the design process, leveraging an analytic derivative tool chain through the entire optimization. The proposed approach leverages machine learning techniques to generate a differentiable confidence interval output from polynomial chaos models. This technique, coupled with the incorporation of analytical derivatives through the Polynomial Chaos Expansion (PCE) process, eliminates the need to estimate derivatives which are usually obtained from finite difference, complex step, or similar methods. Developing a differentiable confidence interval allows mixed uncertainty problems (both epistemic and aleatory) to be modeled. Without such modeling, these problems cannot accurately predict objective functions containing statistical quantities such as mean and variance. The addition of analytic derivatives to a polynomial chaos-based UQ method decreases the computational costs of performing design under uncertainty by orders of magnitude in comparison with methods such as complex step. The method and codes developed are modular in nature and are a drop-in solution for design under uncertainty within existing MDO problems. A low-fidelity analytical multidisciplinary optimization under uncertainty for a wing design in OpenMDAO is detailed in this paper. This demonstration case will include both objective functions and constraints which are influenced by uncertain parameters.

Ben D Phillips↗

End-To-End Uncertainty Quantification with Analytical Derivatives for Design Under Uncertainty

Uncertainty quantification (UQ) is a rapidly growing and evolving discipline, especially within the aerospace community. Performing analysis with UQ can provide decision makers with a wealth of information about a candidate design. However, the value of UQ is fully realized when the information gained during UQ analysis is leveraged in a feedback loop of a design optimization process, often referred to as design under uncertainty. Although design under uncertainty can be a powerful risk mitigation technique, there are a number of roadblocks that prevent its implementation. Two primary factors are computational costs and added complexity of the analysis. High fidelity simulations on the order tens of uncertain variables quickly become computationally infeasible. Also, implementing UQ into an existing multidisciplinary design and optimization (MDO) process often requires extensive knowledge of the UQ methods and careful treatment of the problem formulation. The objective of this work is to address these two primary roadblocks and enable practitioners to efficiently perform design under uncertainty with limited knowledge of the UQ discipline. Methods outlined in this paper demonstrate MDO incorporating UQ into the design process, leveraging an analytic derivative tool chain through the entire optimization. The proposed approach leverages machine learning techniques to generate a differentiable confidence interval output from polynomial chaos models. This technique, coupled with the incorporation of analytical derivatives through the Polynomial Chaos Expansion (PCE) process, eliminates the need to estimate derivatives which are usually obtained from finite difference, complex step, or similar methods. Developing a differentiable confidence interval allows mixed uncertainty problems (both epistemic and aleatory) to be modeled. Without such modeling, these problems cannot accurately predict objective functions containing statistical quantities such as mean and variance. The addition of analytic derivatives to a polynomial chaos-based UQ method decreases the computational costs of performing design under uncertainty by orders of magnitude in comparison with methods such as complex step. The method and codes developed are modular in nature and are a drop-in solution for design under uncertainty within existing MDO problems. A low-fidelity analytical multidisciplinary optimization under uncertainty for a wing design in OpenMDAO is detailed in this paper. This demonstration case will include both objective functions and constraints which are influenced by uncertain parameters.

Ben Phillips↗

Prompt D 0 , D + , and D* + production in Pb–Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 5.02 TeV

The production of prompt D 0 , D + , and D* + mesons was measured at midrapidity (|y| < 0.5) in Pb–Pb collisions at the centre-of-mass energy per nucleon–nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV with the ALICE detector at the LHC. The D mesons were reconstructed via their hadronic decay channels and their production yields were measured in central (0–10%) and semicentral (30–50%) collisions. The measurement was performed up to a transverse momentum (p T ) of 36 or 50 GeV/c depending on the D meson species and the centrality interval. For the first time in Pb–Pb collisions at the LHC, the yield of D 0 mesons was measured down to p T = 0, which allowed a model-independent determination of the p T -integrated yield per unit of rapidity (dN/dy). A maximum suppression by a factor 5 and 2.5 was observed with the nuclear modification factor (R AA ) of prompt D mesons at p T = 6–8 GeV/c for the 0–10% and 30–50% centrality classes, respectively. The D-meson RAA is compared with that of charged pions, charged hadrons, and J/ψ mesons as well as with theoretical predictions. The analysis of the agreement between the measured R AA , elliptic (v 2 ) and triangular (v 3 ) flow, and the model predictions allowed us to constrain the charm spatial diffusion coefficient D s . Furthermore the comparison of R AA and v 2 with different implementations of the same models provides an important insight into the role of radiative energy loss as well as charm quark recombination in the hadronisation mechanisms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A multi-item maintenance center inventory model for low-demand reparable items

In many military and commercial contexts, complex equipment undergoes scheduled maintenance overhauls at regular intervals during which all failed components are replaced. Failure to have replacements on hand for all failed parts requires emergency measures at premium cost. When reparable parts are highly reliable and expensive, both holding and shortage costs are high. This model determines the reparable parts inventory for a maintenance center under three alternative criteria: (1) maximizing job-completion rate subject to constraint on total holding costs, (2) minimizing total holding costs plus expected job noncompletion costs, and (3) minimizing total holding costs subject to a required minimum job-completion rate. Exact solutions may be obtained using dynamic programming. Approximate solutions, found easily by marginal analysis, have readily computed bounds on possible error. The solution methods for the three formulations are illustrated in a simple example.

Schaefer, M. K.↗

Analysis of AIS radiometry with emphasis on determination of atmospheric properties and surface reflectance

Airborne Imaging Spectrometer data from Mono Lake, California, were studied in order to establish spectral radiance of test areas under solar illumination. The objective is to provide a method of atmospheric correction for major absorbers from the spectrometer data themselves. Crucial to the analysis is radiometric calibration of the instrument. Good agreement is found between calculated and measured radiances for uniform surface targets (beaches), but simulations of atmospheric properties with LOWTRAN lead to unreasonably low values of atmospheric precipitable water. Absorptions from carbon dioxide are not detected in the AIS data, but are strongly present in the LOWTRAN model. The apparent low contrast of all atmospheric absorption bands leads to a study of contamination from overlapping spectral orders in the AIS data. The suspected contamination is shown unambiguously to be present beyond approximately 1500 nm and consists of an extra radiance term including atmospheric bands from the delta/2 wavelength interval. A rigorous removal of the unwanted spectral contamination does not seem possible for any data taken in the rock mode. Rough estimates for tree-mode observations might be pieced together is a suitable after-the-fact radiometric claibration of the instrument can be formulated.

Conel, J. E.↗

Flight motor set 360L008 (STS-32R). Volume 1: System overview

Flight motor set 360L008 was launched as part of NASA space shuttle mission STS-32R. As with all previous redesigned solid rocket motor launches, overall motor performance was excellent. All ballistic contract end item specification parameters were verified with the exception of ignition interval and rise rates, which could not be verified due to elimination of developmental flight instrumentation. But the available low sample rate data showed nominal propulsion performance. All ballistic and mass property parameters closely matched the predicted values and were well within the required contract end item specification levels that could be assessed. All field joint heaters and igniter joint heaters performed without anomalies. Redesigned field joint heaters and the redesigned left-hand igniter heater were used on this flight. The changes to the heaters were primarily to improve durability and reducing handling damage. Evaluation of the ground environment instrumentation measurements again verified thermal mode analysis data and showed agreement with predicted environmental effects. No launch commit criteria violation occurred. Postflight inspection again verified superior performance of the insulation, phenolics, metal parts, and seals. Postflight evaluation indicated both nozzles performed as expected during flight. All combustion gas was contained by insulation in the field and case-to-nozzle joints. Recommendations were made concerning improved thermal modeling and measurements. The rationale for these recommendations and complete result details are presented.

Garecht, D. M.↗

Energy Finite Element Analysis for Computing the High Frequency Vibration of the Aluminum Testbed Cylinder and Correlating the Results to Test Data

The Energy Finite Element Analysis (EFEA) is a finite element based computational method for high frequency vibration and acoustic analysis. The EFEA solves with finite elements governing differential equations for energy variables. These equations are developed from wave equations. Recently, an EFEA method for computing high frequency vibration of structures either in vacuum or in contact with a dense fluid has been presented. The presence of fluid loading has been considered through added mass and radiation damping. The EFEA developments were validated by comparing EFEA results to solutions obtained by very dense conventional finite element models and solutions from classical techniques such as statistical energy analysis (SEA) and the modal decomposition method for bodies of revolution. EFEA results have also been compared favorably with test data for the vibration and the radiated noise generated by a large scale submersible vehicle. The primary variable in EFEA is defined as the time averaged over a period and space averaged over a wavelength energy density. A joint matrix computed from the power transmission coefficients is utilized for coupling the energy density variables across any discontinuities, such as change of plate thickness, plate/stiffener junctions etc. When considering the high frequency vibration of a periodically stiffened plate or cylinder, the flexural wavelength is smaller than the interval length between two periodic stiffeners, therefore the stiffener stiffness can not be smeared by computing an equivalent rigidity for the plate or cylinder. The periodic stiffeners must be regarded as coupling components between periodic units. In this paper, Periodic Structure (PS) theory is utilized for computing the coupling joint matrix and for accounting for the periodicity characteristics.

Vlahopoulos, Nickolas↗

The Effect of Information Feedback Upon Psychophysical Judgments

An analysis was made of the role of presentation schedules and information feedback on performance in a forced-choice signal detection task. The experimental results indicate that information feedback facilitates performance, but only for certain presentation schedules. The present study was designed to assess performance in a signal detection task under two conditions of information feedback. In the I-condition, S was told on each trial whether his detection response was correct or incorrect; in the !-condition S was given no feedback regarding the correctness of his response. The task involved a 2-response, forced-choice auditory detection problem. On each trial 2 temporal intervals were defined and S was required to report which interval he believed contained the signal; i. e., in one interval a tone burst in a background of white noise was presented, while the other interval contained only white noise. A trial will be denoted as s1 or s2, depending on whether the signal was embedded in the 1st or 2nd interval; the S's response will be denoted A1 or A2 to indicate which interval he reported contained the signal. The probability of an s1 trial will be denoted as y. In this study two values of y were used (.50 and.75) and, as indicated above, two conditions of information feedback. Thus there were 4 experimental conditions (501, · 50I, 751, 75I); each S was run under all 4 conditions. Method Gaussian noise was presented binaurally in S's headphones throughout a test session and the signal was a 1000-cps sinusoid tone; the tone was presented for 100 msec. including equal fall and rise times of 20 msec. The ratio of signal energy to noise power in a unit bandwidth was 2.9, and was constant throughout the study. The. S was seated before a stimulus display board. On each trial a red warning light was flashed for 100 msec. Two amber lights then came on successively each for 1 sec.; these lights defined the 2 observation intervals. The onset of the signal occurred 500 msec. after the onset of one of the observation intervals. After the second amber light went off, S indicated his response by pressing 1 of 2 wand switches under cards reading "1st interval" and "2nd interval." For the !-condition a green light flashed on above the correct response key after S's response; the green light was omitted in the !-condition. Each trial lasted 6 sec. The S's were 12 male college students with normal hearing. They were run for two practice sessions followed by 20 test sessions. Test sessions were run on consecutive days, 350 trials/day. Each day S ran on 1 of the 4 experimental conditions; in successive 4-day blocks S ran one day on each of the 4 experimental conditions in a random order. Thus, over 20 days each of the experimental conditions was repeated 5 times.

Atkinson, Richard C.↗

Cryo-Temp Data Logger

Biological experiments on the ISS are frequently frozen at ultra-low temperatures on orbit prior to return to Earth for analysis. Recording temperatures during the freezing process, and subsequent handling to the point of post-flight processing, can provide important information of sample integrity to the scientist. The MadgeTech Cryo-Temp Data Logger has been modified for space flight applications by the Payload Development Team at Kennedy Space Center. This ultra-low temperature data logger is rated for temperatures from -86 C to +35 C with 0.1 C resolution and +/-1 C accuracy, but has been successfully tested at Kennedy Space Centers Cryogenic Lab down to -100 C. The logger provides date and time stamped temperature recordings and can be set with logging intervals from 5 seconds up to 30 minutes using the MadgeTech Data Logger Software. Data can be recorded, depending on logging interval, for greater than 90 days. The Cryo-Temp Data logger will first be flown on Space X-12 in the BRIC-Petri Dish Fixation Unit (PDFU) canister hardware, recording temperature data for a Space Life and Physical Sciences (SLPS) flight experiment. The modified Cryo-Temp Data Logger is small in size (Figure 1) and is the ideal temperature monitoring solution for any application involving independent temperature data logging.

Logger↗

Interactions between CO2 chemoreflexes and arterial baroreflexes

We studied interactions between CO2 chemoreflexes and arterial baroreflexes in 10 supine healthy young men and women. We measured vagal carotid baroreceptor-cardiac reflexes and steady-state fast Fourier transform R-R interval and photoplethysmographic arterial pressure power spectra at three arterial pressure levels (nitroprusside, saline, and phenylephrine infusions) and three end-tidal CO2 levels (3, 4, and 5%, fixed-frequency, large-tidal-volume breathing, CO2 plus O2). Our study supports three principal conclusions. First, although low levels of CO2 chemoreceptor stimulation reduce R-R intervals and R-R interval variability, statistical modeling suggests that this effect is indirect rather than direct and is mediated by reductions of arterial pressure. Second, reductions of R-R intervals during hypocapnia reflect simple shifting of vagally mediated carotid baroreflex responses on the R-R interval axis rather than changes of baroreflex gain, range, or operational point. Third, the influence of CO2 chemoreceptor stimulation on arterial pressure (and, derivatively, on R-R intervals and R-R interval variability) depends critically on baseline arterial pressure levels: chemoreceptor effects are smaller when pressure is low and larger when arterial pressure is high.

NASA Discipline Cardiopulmonary↗

Observation and Numerical Simulation of Cold Ions Energized by EMIC Waves

This is the first report of significant energization (up to 7,000 eV) of low-energy He + ions, which occurred simultaneously with H-band electromagnetic ion cyclotron (EMIC) wave activity, in a direction mostly perpendicular to the ambient magnetic field. The event was detected by the Arase satellite in the dayside plasmatrough region off the magnetic equator on 15 May 2019. The peak energy of the He + flux enhancements is mostly above 1,000 eV. At some interval, the He + ions are energized up to ~7,000 eV. The H-band waves are excited in a frequency band between the local crossover and helium gyrofrequencies and are close to a linear polarization state with weakly left-handed or right-handed polarization. The normal angle of the waves exhibits significant variation between 0° and 80°, indicating a non-parallel propagation. Here, we run a hybrid code with parameters estimated from the Arase observations to examine the He + energization. The simulations show that cold He + ions are energized up to more than 1,000 eV, similar to the spacecraft observations. From the analysis of the simulated wave fields and cold plasma motions, we found that the ratio of the wave frequency to He + gyrofrequency is a primary factor for transverse energization of cold He + ions. As a consequence of the numerical analysis, we suggest that the significant transverse energization of He + ions observed by Arase is attributed to H-band EMIC waves excited near the local helium gyrofrequency.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic field and particle pressure in the plasma sheet of Jupiter

The results of an analysis of the energetic particle and magnetic field data acquired by the Voyager 2 spacecraft at distances of about 40-70 Jupiter radii on the nightside of the planet are reported. As in a previous study of similar data at distances of greater than about 80 Jupiter radii, the energy densities of ions (primarily protons) is found to be sufficient to provide the diamagnetic depressions measured in the magnetic field intensity as the spacecraft made successive encounters with the nightside plasma sheet. There is some evidence that the percent contribution of the protons to the energy balance decreases with increasing distance from the planet over this radial interval, although this conclusion is dependent upon the assumption that the proton and heavier ion (oxygen) energy spectra are similar.

Lanzerotti, L. J.↗

Observations of IMF and seasonal effects in high-latitude convection

Strong interplanetary magnetic field (IMF) and seasonal effects in the convection of nightside ionospheric plasma are described. The findings are based on a statistical analysis of observations made with the Johns Hopkins University/ Applied Physics Lab (JHU/APL) HF radar located at Goose Bay, Labrador. For positive sign of the IMF dusk-dawn component, By greater than 0 the dawn cell is more crescent shaped and the dusk cell more round while for BY less than 0 these pairings of size and shape are reversed. The more extreme crescent /round cell dichotomy is obtained for BY greater than 0. The return flows associated with the crescent-shaped cell dominate at midnight MLT (magnetic local time); the reversal in the zonal velocity in the 67 deg-69 deg lambda (magnetic latitude) interval occurs 2.5 hr earlier in summer than in winter. The maximum effects are obtained on the nightside for the pairings By greater than 0, summer and BY less than 0, winter; the first produces the more structured cell in the morning, the second in the evening, and this cell dominates the return flow at midnight. The difference in the zonal flow reversals for these pairings exceeds 4 hr in MLT.

Ruohoniemi, J. M.↗

Sustainable Campus with PEV and Microgrid

Market penetration of electric vehicles (EVs) is gaining momentum, as is the move towards increasingly distributed, clean and renewable electricity sources. EV charging shifts a significant portion of transportation energy use onto building electricity meters. Hence, integration strategies for energy-efficiency in buildings and transport sectors are of increasing importance. This paper focuses on a portion of that integration: the analysis of an optimal interaction of EVs with a building-serving transformer, and coupling it to a microgrid that includes PV, a fuel cell and a natural gas micro-turbine. The test-case is the Nanyang Technological University (NTU), Singapore campus. The system under study is the Laboratory of Clean Energy Research (LaCER) Lab that houses the award winning Microgrid Energy Management System (MG-EMS) project. The paper analyses three different case scenarios to estimate the number of EVs that can be supported by the building transformer serving LaCER. An approximation of the actual load data collected for the building into different time intervals is performed for a transformer loss of life (LOL) calculation. The additional EV loads that can be supported by the transformer with and without the microgrid are analyzed. The numbers of possible EVs that can be charged at any given time under the three scenarios are also determined. The possibility of using EV fleet at NTU campus to achieve demand response capability and intermittent PV output leveling through vehicle to grid (V2G) technology and building energy management systems is also explored.

Singh, Reshma↗

Cloud, Aerosol, and Radiative Properties Over the Western North Atlantic Ocean

This study examines the atmospheric properties of weather states (WSs) derived from the International Satellite Cloud Climatology Project (ISCCP) over the Western North Atlantic Ocean. In particular, radiation and aerosol data corresponding to two sites in the study domain, Pennsylvania State University and Bermuda, were examined to characterize the atmospheric properties of the various satellite-derived WSs. At both sites, the fair weather WS was most prevalent, followed by the cirrus WS. Differences in the seasonality of the various WSs were observed at the two sites. Fractional sky cover and effective shortwave cloud transmissivity derived from ground-based radiation measurements were able to capture differences among the satellite-derived WSs. Speciated aerosol optical thicknesses (AOT) from the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) were used to investigate potential differences in aerosol properties among the WSs. The clear sky WS exhibited below-average seasonal values of AOT at both sites year-round, as well as relatively high rates of occurrence with low AOT events. In addition, the clear sky WS showed above-average contributions from dust and black carbon to the total AOT year-round. Finally, transitions between various WSs were examined under low, high, and mid-range AOT conditions. The most common pathway was for the WSs to remain in the same state after a 3 hr interval. Some WSs, such as midlatitude storms, deep convection, middle top, and shallow cumulus, were more prevalent as ending states under high AOT conditions. This work motivates examining differences in aerosol properties between WSs in other regions.

54 ENVIRONMENTAL SCIENCES↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗

Vitrification of High-Cr Glass in Research-Scale Melter

This test report describes the experimental results from a small-scale test using the research-scale melter (RSM) at Pacific Northwest National Laboratory to demonstrate processing of a high-Cr simulated feed stream, designated HLW-HCr-16. The RSM is a small, joule-heated melter capable of processing melter feed continuously. The melter is equipped with Inconel ® 693 electrodes, Monofrax ® K-3 refractory, and an Inconel 690 pour spout. An electric kiln surrounds the melter body and minimizes heat loss from the melter body during operation. The RSM is equipped with an offgas treatment system that employs quenching, wet scrubbing, and high-efficiency mist elimination. The glass-discharge section is heated to facilitate pouring of the glass. The melter is fitted with a melt cavity that is ~25 cm (10 in.) in diameter with a nominal glass depth of 8.9 cm (3.5 in.). The melter was operated with a target glass temperature of 1150°C and target plenum temperature between 550°C and 700°C for this test. The air bubbling rate was 4.2 L/min. Overall, during the continuous operation of the melter for ~ 103 hours, ~ 141 kg of glass was produced. At the conclusion of the test, the melter and exhaust lines were visually inspected for particulate deposition and corrosion. Entrained material had adhered to the underside of the melter lid and to the exhaust piping. Enrichments in elements such as Cl, F, B, K, Li, P, Na, and S were measured in these deposits through inductively coupled plasma–optical emission spectroscopy and X-ray fluorescence analysis. When the melter electrodes and air bubbler tube were removed from the glass in the RSM, the electrodes appeared discolored, but no significant loss of metal was observed. The processing of a high-Cr simulant, HLW-HCr-16, in the RSM produced glass at an average rate of 1.36 kg/h, equaling a melter-surface-area normalized glass generation rate of 654 kg/day/m 2 . The resulting glass met the toxicity characteristic leaching procedure requirement. Test results of crystallinity, electrical conductivity, and viscosity showed good processing properties of this high-Cr high-level waste glass. RSM offgas was also sampled and analyzed at periodic intervals during steady-state operating conditions. The total decontamination factor averaged by four sampling periods was 134. The concentrations of CO and NO x in emissions were 237 to 422 and 69 to 94 parts per million by volume, respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

User guide for the Geospatial Risk Analysis Tool (V.1.0.1)

Developed by the Pacific Northwest National Laboratory (PNNL), the Geospatial Risk Analysis Tool (GRAT) links to Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT), a complete system for computing simple air parcel trajectories, as well as complex transport, dispersion, chemical transformation, and deposition simulations. The HYSPLIT model was developed by the National Oceanic and Atmospheric Administration (NOAA) Air Resources Laboratory and the Australian Bureau of Meteorology Research Center in 1998. GRAT is used to batch process meteorological data into ARL format. While HYSPLIT GUI limits users to convert meteorological data up to 6 date time points per operation, GRAT enables users to convert a significantly large amount of data (For example, 10 years of meteorological data measured with an interval of 5 minutes) within one operation. Incorporated with the Potential source distribution function (PSDF), GRAT can also be used to determine the areas influenced by the emission of hazardous chemicals.

54 ENVIRONMENTAL SCIENCES↗