Search NASA⌕ Search

SEARCH · Search NASA

Results for “Activity coefficient”

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,657 records · Page 92

Biomechanics of the Treadmill Locomotion on the International Space Station

Exercise prescriptions completed by International Space Station (ISS) crewmembers are typically based upon evidence obtained during ground-based investigations, with the assumption that the results of long-term training in weightlessness will be similar to that attained in normal gravity. Coupled with this supposition are the assumptions that exercise motions and external loading are also similar between gravitational environments. Normal control of locomotion is dependent upon learning patterns of muscular activation and requires continual monitoring of internal and external sensory input [1]. Internal sensory input includes signals that may be dependent on or independent of gravity. Bernstein hypothesized that movement strategy planning and execution must include the consideration of segmental weights and inertia [2]. Studies of arm movements in microgravity showed that individuals tend to make errors but that compensation strategies result in adaptations, suggesting that control mechanisms must include peripheral information [3-5]. To date, however, there have been no studies examining a gross motor activity such as running in weightlessness other than using microgravity analogs [6-8]. The objective of this evaluation was to collect biomechanical data from crewmembers during treadmill exercise before and during flight. The goal was to determine locomotive biomechanics similarities and differences between normal and weightless environments. The data will be used to optimize future exercise prescriptions. This project addresses the Critical Path Roadmap risks 1 (Accelerated Bone Loss and Fracture Risk) and 11 (Reduced Muscle Mass, Strength, and Endurance). Data were collected from 7 crewmembers before flight and during their ISS missions. Before launch, crewmembers performed a single data collection session at the NASA Johnson Space Center. Three-dimensional motion capture data were collected for 30 s at speeds ranging from 1.5 to 9.5 mph in 0.5 mph increments with a 12-camera system. During flight, each crewmember completed up to 6 data collection sessions spread across their missions, performing their normal exercise prescription for the test day, resulting in varying data collection protocols between sessions. Motion data were collected by a single HD video camera positioned to view the crewmembers' left side, and tape markers were placed on their feet, legs, and neck on specific landmarks. Before data collection, the crewmembers calibrated the video camera. Video data were collected during the entire exercise session at 30 Hz. Kinematic data were used to determine left leg hip, knee, and ankle range of motion and contact time, flight time, and stride time for each stride. 129 trials in weightlessness were analyzed. Mean time-normalized strides were found for each trial, and cross-correlation procedures were used to examine the strength and direction of relationships between segment movement pattern timing in each gravitational condition. Cross-correlation analyses between gravitational conditions revealed highly consistent movement patterns at each joint. Peak correlation coefficients occurred at 0% phase, indicating there were no lags in movement timing. Joint ranges of motion were similar between gravitational conditions, with some slight differences between subjects. Motion patterns in weightlessness were highly consistent at a given speed with those occurring in 1G, indicating that despite differing sensory input, subjects maintain running kinematics. The data suggest that individuals are capable of compensating for loss of limb weight when creating movement strategies. These results have important implications for creating training programs for use in weightlessness as practitioners can have greater confidence in running motions transferring across gravitational environments. Furthermore, these results have implications for use by researchers investigating motor control mechanisms and investigating hypotheses related to movement strategies when using sensory input that is dependent upon gravity.

DeWitt, John↗

Technology Development and Field Monitoring in nZEB - US Country report IEA HPT Annex 49 Task 3

The International Energy Agency (IEA) Heat Pumping Technologies (HPT) Annex 49, “Design and Integration of Heat Pumps for Nearly Zero Energy Buildings,” deals with the application of heat pumps (HPs) as a core component of the HVAC system for nearly or net-zero energy buildings. This report covers the Task 3 activities of the US team. Three institutions are involved on the US team and have worked on the following projects. 1) Oak Ridge National Laboratory (ORNL) summarized development activities since the conclusion of IEA HPT Annex 40 for several integrated HP (IHP) systems - electric ground-source IHP and air-source IHP versions and engine-driven AS-IHP version. 2) The University of Maryland partnered with ORNL and Blue Bear Management to develop a personal cooling device called RoCo that can provide personalized conditioned air to occupants in inadequately or unconditioned environments. With RoCo, building facility management can elevate the HVAC thermostat settings without compromising occupants’ thermal comfort. Researchers have found that a 4°F increase in thermostat settings can save 12–30% energy savings. Therefore, RoCo is a promising technology that can help reduce building energy consumption and facilitate achievement of net-zero energy building performance. 3) The National Institute of Standards and Technology (NIST) is working on a field study effort on the NIST Net-Zero Energy Residential Test Facility. Two air-source split-system HPs were installed in a residential, net-zero energy home that was constructed as a laboratory on the NIST campus in Gaithersburg, Maryland. The first HP was a two-stage, 7 kW (2 ton), 15.8 seasonal energy efficiency ratio (SEER), 9.05 heating seasonal performance factor (HSPF) conventionally ducted system; the second HP was a variable-speed, 10.6 kW (3 ton), 14 SEER, 8.35 HSPF, high-velocity ducted system. These two systems operated side by side, using separate supply ducts and a common return duct, on a weekly alternating schedule to condition the home that was operated with very consistent simulated thermal loads. The team wanted to determine whether the high-velocity system could provide comparable energy-use efficiency to the conventional system. The results of this study showed that it did meet the required loads and had slightly greater efficiency; the average cooling coefficient of performance (COP) was (0.40 ± 0.11) higher, and the average heating COP was statistically equal. A new firmware was provided at the end of the heating season that greatly improved the performance of the high-velocity system; its average heating COP went from (1.8 ± 0.9) to (2.5 ± 1.1) at a 95% confidence level. The new firmware heating COP averaged (1.05 ± 0.23) higher than the old firmware over the same outdoor temperatures. The defrost performance is very different for these two systems, yet they consumed equivalent energy per HDD. The conventional system uses a timed-initiate, temperature-terminate algorithm with auxiliary electric resistive heating, whereas the high-velocity system uses calculated evaporator parameters with a hot-gas bypass before a full reverse-cycle defrost with no supplementary resistive heat.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A hybrid analytical and numerical model for cross-over and performance decay in a unit cell vanadium redox flow battery

Developing an accurate and efficient model for cross-over is critical for improving the long-term performance of redox flow batteries (RFBs). Here, this work presents a hybrid analytical and numerical model that combines a two-dimensional analytical solution to the active species, a one-dimensional analytical model for cross-over mechanisms, and a zero-dimensional numerical model for outlet concentrations of reactants. By comparing with experiment over 41 cycles (ca. 144 h), the model reported a mean voltage difference of 0.0089 V, mean time difference of 14 s per cycle (of 3.5 h), maximum relative difference for capacity and energy of 1.34% and 1.63%, respectively. The predicted mean concentrations for V 2+ , V 3+ , and VO$_{2}^{+}$ in membrane are 65%~77% of measured values from the literature. Upon validation, the model reproduced behaviors in electrolyte imbalance similar to those observed in experiments and numerical models, and revealed the control of cross-over, self-discharge, stoichiometry of side reactions, and Coulombic efficiency on electrolyte imbalance. The model also demonstrates excellent computational efficiency for simulating 41 cycles (around 37,000 points) within 3~4 s. The demonstrated efficiency and accuracy for predicting cross-over and its impacts on voltage, capacity&energy decay, membrane concentrations, and electrolyte imbalance makes it a reliable tool for optimizing RFBs’ long-term performance.

25 ENERGY STORAGE↗

Mechanisms underpinning the ultrahigh piezoelectricity in Sm-doped 0.705Pb(Mg 1/3 Nb 2/3 )O 3 -0.295PbTiO 3 : Temperature-induced metastable local structure and field-induced polarization rotation

The solid solution of (100 – x )%Pb(Mg 1/3 Nb 2/3 )O 3 - x %PbTiO 3 (PMN- x PT) exhibits ultrahigh piezoelectric and dielectric properties near the morphotropic phase boundary compositions and, as such, has been extensively studied in recent years. Recently, 2.5 mol. % Sm-doped PMN-29PT polycrystalline ceramics were reported to possess the highest piezoelectric coefficients (~1500 pC/N) among all reported piezoceramics, but the atomic-scale mechanisms for such high piezoelectric properties are not yet clear. Here, in situ X-ray diffraction and X-ray total scattering measurements during the application of an electric field, together with in situ total scattering measurement at different temperatures, were conducted for 2.5 mol. % Sm-doped PMN-29.5PT (2.5Sm-PMN-29.5PT). Both the largest field-induced strain and the piezoelectric response were found in the crystallites oriented with their &#12296 100 &#12297 PC directions parallel to the applied field. The local Pb displacement was analyzed using the reverse Monte Carlo method based on the pair distribution functions at different temperatures, where a temperature-induced directional change of Pb displacement was observed. Based on the experimental observations, a field-induced polarization rotation is suggested to be the dominant mechanism for the ultrahigh piezoelectricity of the 2.5Sm-PMN-29.5PT ceramic, while the ease of polarization rotation is possibly attributed to the temperature-induced metastable local monoclinic symmetries with their polar axes close to &#12296 111 &#12297 PC .

36 MATERIALS SCIENCE↗

Thermal Protection System to Enable Ice Giant Aerocapture Mission for Delivering Both an Orbiter and an In Situ Probe

The Ice Giants have been identified as high priority science destinations in the last Decadal Survey [1] and could benefit from aerocapture as the primary method for orbit insertion [2]. A mass-efficient aerocapture system will enable the delivery of an orbiter along with an atmospheric probe (for in situ measurements to anchor global data collected by the orbiter) and possibly a lander at Triton [3]. Aerocapture could be executed either using low L/D rigid aeroshell with lift modulation (LMA) [4] or using deployable aeroshell using drag modulation (DMA) [3]. Nearly two decades ago, a NASA-funded team performed Neptune-Triton aerocapture studies with a mid-L/D lifting configuration [5] for achieving orbit using LMA. This study showed aerocapture challenges. Due to very high peak entry conditions combined with very high heat-load, a suite of TPS materials was required and this suite was deemed problematic from a qualification perspective, due test facility limitations. In the past 20 years, progress made in GN&C for lift-guided entry missions such as MSL, Orion EFT1, Mars 2020 and the upcoming Artemis missions, and the expertise in blunt body aerodynamics at large scale (~ 5m) has led the EDL community to conclude that aerocapture is a “go do” engineering activity and most technologies are in hand to propose missions with aerocapture [6] [7]. Aerocapture using DMA, currently in development, is an option for Ice Giant Missions. While DMA is simpler in some sense, due to ballistic entry and no need for lift-guided maneuvering, it has challenges and it’s maturity is lower. LMA and DMA both require one or more ablative Thermal Protection System (TPS) materials for the rigid aeroshell element. The ablative TPS needs to be robust and mass efficient due to the high heat loads and size of the rigid aeroshell. Currently, there are capable ablative thermal protection materials, e.g., Heatshield for Extreme Entry Environments Technology (HEEET), 3-D woven Mid-Density Carbon- Phenolic (3MDCP), and PICA (Phenolic-Impregnated Carbon Ablator) that are mature, i.e., at TRL 6 or higher. NASA also invested in Conformal PICA that was matured to TRL 5. Our goal is to evaluate the applicability of high TRL TPS and consider other design options. We first establish bounding aerocapture trajectories for a wide range of arrival conditions and the associated aerothermal environment. Based on the environments, we then determine the predicted TPS mass for the aeroshell [4]. In this presentation, we will outline the process by which we establish bounding aerocapture trajectories for hyperbolic excess velocities ranging from 27 km/s to 35 km/s, which are shown to be a range of velocities that can reduce the trip time from ~14 years to 8 years. The above velocity range translates to ~12 km/s to ~24 km/s arrival velocities at the planetary entry interface [2]. The velocity reduction required to achieve orbit ranges between ~2.5 km/s to 9.5 km/s for both Neptune and Triton. Propulsive insertion alone, due to the amount of fuel required to achieve the required velocity reduction, limits the science returned [2]. We establish the bounding aerocapture trajectories for a low L/D (~ 0.4) configuration for three different ballistic coefficients. The ballistic coefficient range is determined from three different aeroshell diameters of 3m, 4m and 5m and with an entry system mass of 2200 kg. With the above range of design parameters, we then determine conservative/bounding estimates of aerothermal environments by using a combination of CFD simulations and stagnation point heating estimates [7]. This engineering approach allows us to first assess the TPS need vs. TPS capability and determine the applicability of existing TPS. Once an applicable suite of TPS is determined, the TPS thickness and mass are computed. We show that the TPS mass fraction can be as low as 5% to as high as 20%, depending on total trip time reduction and other design parameters for a range of TPS. This is a large range for TPS mass fraction. We show PICA and HEEET can indeed enable aerocapture missions, but the missions incur a mass penalty. TPS mass savings, can be further reduced with the use of conformal PICA. Advancing the development of Conformal PICA to make it robust across the entire aerothermal environment (peak heat-flux, pressure and shear) range will result in TPS mass fractions of < 10% for Ice Giant aerocapture missions such as the Neptune-Triton mission. Aerocapture allows for not only shortening the trip time but enables larger mass to be placed in orbit. Furthermore, probes deployed from orbit will benefit in reduced entry environments allowing for a lower risk TPS implementation as compared to mission designs where the probe is released prior to orbit insertion. One of the challenges for the Ice Giant community is to ensure mission designs that maximize science and allow flexibility in the placement of the entry probe. The traditional approach to release the probe ahead of the orbiter may not optimize returned science. In this presentation, we will make the case for mature TPS such as HEEET and PICA. While these materials can enable aerocapture missions, completing the development of conformal PICA and extending Conformal PICA to be more robust, will have significant impact to TPS mass efficiency and significantly enhance science return for future Gas- and Ice-Giant missions.

E Venkatapathy↗

Requirements Description of DASSH-F

This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A robust bi-orthogonal/dynamically-orthogonal method using the covariance pseudo-inverse with application to stochastic flow problems

We develop a new robust methodology for the stochastic Navier–Stokes equations based on the dynamically-orthogonal (DO) and bi-orthogonal (BO) methods . Both approaches are variants of a generalized Karhunen–Loève (KL) expansion in which both the stochastic coefficients and the spatial basis evolve according to system dynamics, hence, capturing the low-dimensional structure of the solution. The DO and BO formulations are mathematically equivalent , but they exhibit computationally complimentary properties. Specifically, the BO formulation may fail due to crossing of the eigenvalues of the covariance matrix, while both BO and DO become unstable when there is a high condition number of the covariance matrix or zero eigenvalues. To this end, we combine the two methods into a robust hybrid framework and in addition we employ a pseudo-inverse technique to invert the covariance matrix. The robustness of the proposed method stems from addressing the following issues in the DO/BO formulation: (i) eigenvalue crossing: we resolve the issue of eigenvalue crossing in the BO formulation by switching to the DO near eigenvalue crossing using the equivalence theorem and switching back to BO when the distance between eigenvalues is larger than a threshold value; (ii) ill-conditioned covariance matrix: we utilize a pseudo-inverse strategy to invert the covariance matrix; (iii) adaptivity: we utilize an adaptive strategy to add/remove modes to resolve the covariance matrix up to a threshold value. In particular, we introduce a soft-threshold criterion to allow the system to adapt to the newly added/removed mode and therefore avoid repetitive and unnecessary mode addition/removal. When the total variance approaches zero, we show that the DO/BO formulation becomes equivalent to the evolution equation of the Optimally Time-Dependent modes . We demonstrate the capability of the proposed methodology with several numerical examples, namely (i) stochastic Burgers equation: we analyze the performance of the method in the presence of eigenvalue crossing and zero eigenvalues; (ii) stochastic Kovasznay flow: we examine the method in the presence of a singular covariance matrix; and (iii) we examine the adaptivity of the method for an incompressible flow over a cylinder where for large stochastic forcing thirteen DO/BO modes are active.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cesium Ion Exchange Using Crystalline Silicotitanate with 5.6 M Sodium Simulant

The Low-Activity Waste Pretreatment System (LAWPS) facility is planned to pretreat Hanford tank waste supernate by filtering to remove solids and processing through ion exchange columns to remove cesium. The ion exchange media is currently (since 2018) targeted to be crystalline silicotitanate (CST). This material has been produced as an engineered form (generally spherical beads) by Honeywell UOP LLC (Des Plaines, IL). Previous forms of this media have been tested in a wide array of simulants and process scales, but column breakthrough and batch contact testing to date on actual tank waste has been limited to a few tank wastes. Washington River Protection Solutions requested a series of studies to demonstrate Cs load behavior in a small column format and determine if a simple batch contact test with CST can be used to qualify tank waste supernate prior to processing it in the LAWPS. The goal of the batch contact testing is to determine if there are any issues with the tank waste that would preclude Cs removal by the ion exchanger. Testing was conducted with IONSIV TM R9140-B CST provided by Honeywell UOP in 2018 (Batch 2081000057). The CST was received in the sodium (Na) form, so no pretreatment conditioning to convert to the Na form was required. Before use in testing, the CST was passed through a #25 sieve (710-micron openings) and collected on a #60 sieve (250-micron openings) with the intent of removing any large and fine particles. Physical properties were measured on the CST including particle size distribution (d50 = 570 microns) and scanning electron micrographic examination. The CST was rinsed with water until the contact solution was visually nearly clear. The CST bed density (1.00 g/mL) and CST bed void fraction (0.656) were measured on the sieved and rinsed CST. The batch contact testing varied Cs concentration in a 5.6 M Na simple simulant solution. Testing was conducted in a liquid volume to solids mass ratio of 200. The distribution coefficient (K d ) was determined to be 770 mL/g at the Cs equilibrium condition of 8.0 μg/mL; with a bed density of 1.00 g/mL, this corresponds to a predicted 50% Cs breakthrough of 770 bed volumes (BVs). The Cs load capacity at 8.0 μg/mL equilibrium condition was determined to be 6.16 mg Cs/g dry CST. The column testing was prototypic to the intended LAWPS operations in a lead-lag column format, albeit on a small-scale basis with 10-mL CST beds. The feed was processed downflow through the lead column and then through the lag column. Three flowrates were tested: 1.19, 1.99, and 4.56 BV/h. The feed was displaced with 0.1 M NaOH, and then the columns were rinsed with water sequentially through the lead then lag columns. One column system was evaluated top to bottom with a gamma scanner to assess the axial Cs loading. Table ES.1 summarizes the relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of an Assessment Method for Radiation Exposure from NORM in Construction Materials - 19032

Once radiological impact is usually considered to be related to nuclear energy accidents and atomic bombs, most of the research on radiation protection has been focused on artificial radionuclides and radioactive waste. Less attention had been given to radiation exposure due to the presence of naturally occurring radioactive materials in raw materials usually considered as not radioactive and used in manufacturing of many building materials. At national level, radiological data from building materials is not available as the control of building materials regarding the content of naturally occurring radionuclides is addressed for the first time in the EU Directive 2013/59/EURATOM. At international level the research developed is largely incomplete: data on the activity concentration measurements of Ra-226, Th-232 and K-40 in building materials used in the construction industry in the EU is non-statistically representative. The data collected are not enough to have an idea of the European situation, but they confirm high concentration of Ra-226, Th-232 and K-40, in some cases. Moreover, most of the research studies carried out in the EU are based on the assessment of the activity concentration index, which was derived for practical monitoring purposes to identify whether a dose criterion is met; to be used only as a screening tool for identifying materials which might be of concern. The dose criterion considers only effective dose due to γ-radiation, but population at home is also exposed to Rn-222 (and Rn-220) from soil, water and building material itself, sometimes to β radiation from superficial materials, and even to Cs-137 from some types of wood. The indoors exposure can be actually much higher than the value assessed with the dose criterion method. Also, in most cases toxic and radioactive components are inseparably. Therefore, all technical questions related to the toxic part of these materials should be addressed: an integrated approach is needed to assess the impacts of using these materials (new materials, recyclates and by-products). Until recently, radiation from natural sources seemed both unremarkable and unalterable - a background phenomenon. It is now known that doses from natural sources and in particular, decay products of radon can be remarkably high within houses in some areas. Therefore, concerning exposure to natural radiation sources, special provisions have been laid down in the EU Directive 1996/28 EURATOM. Still, the radiological evaluation of building materials was not mandatory in this Directive. The most recent EU Directive 2013/59/EURATOM brings some changes concerning the exposure to natural occurring radioactive materials but the previous screening method was kept. An effective dose of 1 mSv/year was set, as a basis for the exposure to radionuclides of natural origin (except for radon). This work seeks to develop an integrated radiological and chemical-mineralogical assessment methodology for realistic exposure scenarios concerning different types of raw materials, by-products and recyclates and used as building materials. The methodology includes the characterization of selected materials accounting for radiological exposure - calculation of effective dose; toxicity assessment - determination of chemical composition and leaching tests, study of changes in the microstructure and consequently in the radon emanation coefficient. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Global Top-Down Smoke-Aerosol Emissions Estimation Using Satellite Fire Radiative Power Measurements

Fire emissions estimates have long been based on bottom-up approaches that are not only complex, but also fraught with compounding uncertainties. We present the development of a global gridded (1 deg ×1 deg) emission coefficients (Ce) product for smoke total particulate matter (TPM) based on a top-down approach using coincident measurements of fire radiative power (FRP) and aerosol optical thickness (AOT) from the Moderate-resolution Imaging Spectroradiometer (MODIS) sensors aboard the Terra and Aqua satellites. This new Fire Energetics and Emissions Research version 1.0 (FEER.v1) Ce product has now been released to the community and can be obtained from http://feer.gsfc. nasa.gov/, along with the corresponding 1-to-1 mapping of their quality assurance (QA) flags that will enable the Ce values to be filtered by quality for use in various applications. The regional averages of Ce values for different ecosystem types were found to be in the ranges of 16-21/gMJ−1 for savanna and grasslands, 15-32/gMJ−1 for tropical forest, 9-12/gMJ−1 for North American boreal forest, and 18- 26/MJ−1 for Russian boreal forest, croplands and natural vegetation. The FEER.v1 Ce product was multiplied by time-integrated FRP data to calculate regional smoke TPM emissions, which were compared with equivalent emissions products from three existing inventories. FEER.v1 showed higher and more reasonable smoke TPM estimates than two other emissions inventories that are based on bottom-up approaches and already reported in the literature to be too low, but portrayed an overall reasonable agreement with another top-down approach. This suggests that top-down approaches may hold better promise and need to be further developed to accelerate the reduction of uncertainty associated with fire emissions estimation in air-quality and climate research and applications. Results of the analysis of FEER.v1 data for 2004-2011 show that 65-85 Tg yr−1 of TPM is emitted globally from open biomass burning, with a generally decreasing trend over this short time period. The FEER.v1 Ce product is the first global gridded product in the family of "emission factors", that is based essentially on satellite measurements, and requires only direct satellite FRP measurements of an actively burning fire anywhere to evaluate its emission rate in near-real time, which is essential for operational activities, such as the monitoring and forecasting of smoke emission impacts on air quality.

biomass burning↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

Assimilation of Wave Imaging Radar Observations for Real-time Wave-by-Wave Forecasting

This project addressed Topic 3: “Wave Measurement Instrumentation for Feed Forward Controls” under the FOA number DE-FOA-0000971. The overall goal of the program was to develop a phase-resolving wave forecasting technique for application to the active control of Wave Energy Conversion (WEC) devices. We have developed an approach that couples a wave imaging marine radar with a phase-resolving linear wave model for real-time wave field reconstruction and forward propagation of the wave field in space and time. The scope of the project was to develop and assess the performance of this novel forecasting system. Specific project goals were as follows: Develop and verify a fast, GPU-based (Graphical Processing Unit) wave propagation model suitable for phase-resolved computation of nearshore wave transformation over variable bathymetry; Compare the accuracy and speed of performance of the wave model against a deep water model in their ability to predict wave field transformation in the intermediate water depths (50 to 70 m) typical of planned WEC sites; Develop and implement a variational assimilation algorithm that can ingest wave imaging radar observations and estimate the time-varying wave conditions offshore of the domain of interest such that the observed wave field is best reconstructed throughout the domain and then use this to produce model forecasts for a given WEC location; Collect wave-resolving marine radar data, along with relevant in situ wave data, at a suitable wave energy test site, apply the algorithm to the field data, assess performance, and identify any necessary improvements; and Develop a production cost estimate that addresses the affordability of the wave forecasting technology and include in the Final Report. The developed forecasting algorithm (“Wavecast”) was evaluated for both speed and accuracy against a substantial synthetic dataset. Early in the project, performance tests definitively demonstrated that the system was capable of forecasting in real-time, as the GPU-based wave model backbone was very computationally efficient. The data assimilation algorithm was developed on a polar grid domain in order to match the sampling characteristics of the observation system (wave imaging marine radar). For verification purposes, a substantial set of synthetic wave data (i.e. forward runs of the wave model) were generated to be used as ground truth for comparison to the reconstructions and forecasts produced by Wavecast. For these synthetic cases, Wavecast demonstrated very good accuracy, for example, typical forecast correlation coefficients were between 0.84-0.95 when compared to the input data. Dependencies on shadowing, observational noise, and forecast horizon were also identified. During the second year of the project, a short field deployment was conducted in order to assess forecast accuracy under field conditions. For this, a radar was installed on a fishing vessel and observations were collected at the South Energy Test Site (SETS) off the coast of Newport, OR. At the SETS site, simultaneous in situ wave observations were also available owing to an ongoing field project funded separately. Unfortunately, the position and heading information that was available for the fishing vessel were not of sufficient accuracy in order to validate the forecast in a phase-resolving sense. Instead, a spectral comparison was made between the Wavecast forecast and the data from the in situ wave buoy. Although the wave and wind conditions during the field test were complex, the comparison showed a promising reconstruction of the wave spectral shape, where both peaks in the bimodal spectrum were represented. However, the total reconstructed spectral energy (across all directions and frequencies) was limited to 44% of the observed spectrum. Overall, wave-by-wave forecasting using a data assimilation approach based on wave imaging radar observations and a physics-based wave model shows promise for short-term phase-resolved predictions. Two recommendations for future work are as follows: first, we would recommend additional focused field campaigns for algorithm validation. The field campaign should be long enough to capture a range of wave conditions relevant to the target application and WEC site. In addition, it will be crucial to make sure the vessel of choice has high accuracy position and heading instrumentation (this instrumentation is commercially available but not standard on commercial fishing vessels). The second recommendation is to expand the model physics in the wave model backbone to include some nonlinear effects. Specifically, the third-order correction to the wave speed due to amplitude dispersion would be the next step in order to more accurately represent the phase speeds of large amplitude waves.

16 TIDAL AND WAVE POWER↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

Variation in Global Chemical Composition of PM2.5: Emerging Results from SPARTAN

The Surface PARTiculate mAtter Network (SPARTAN) is a long-term project that includes characterization of chemical and physical attributes of aerosols from filter samples collected worldwide. This paper discusses the ongoing efforts of SPARTAN to define and quantify major ions and trace metals found in fine particulate matter (PM (sub 2.5). Our methods infer the spatial and temporal variability of PM (sub 2.5) in a cost-effective manner. Gravimetrically weighed filters represent multi-day averages of PM (sub 2.5), with a collocated nephelometer sampling air continuously. SPARTAN instruments are paired with AErosol RObotic NETwork (AERONET) sun photometers to better understand the relationship between ground-level PM (sub 2.5) and columnar aerosol optical depth (AOD). We have examined the chemical composition of PM (sub 2.5) at 12 globally dispersed, densely populated urban locations and a site at Mammoth Cave (US) National Park used as a background comparison. So far, each SPARTAN location has been active between the years 2013 and 2016 over periods of 2-26 months, with an average period of 12 months per site. These sites have collectively gathered over 10 years of quality aerosol data. The major PM (sub 2.5) constituents across all sites (relative contribution plus or minus Standard Deviation) are ammoniated sulfate (20 percent plus or minus 11 percent), crustal material (13.4 percent plus or minus 9.9 percent), equivalent black carbon (11.9 percent plus or minus 8.4 percent), ammonium nitrate (4.7 percent plus or minus 3.0 percent), sea salt (2.3 percent plus or minus 1.6 percent), trace element oxides (1.0 percent plus or minus 1.1 percent), water (7.2 percent plus or minus 3.3 percent) at 35 percent relative humidity, and residual matter (40 percent plus or minus 24 percent). Analysis of filter samples reveals that several PM (sub 2.5) chemical components varied by more than an order of magnitude between sites. Ammoniated sulfate ranges from 1.1 microns per cubic meter (Buenos Aires, Argentina) to 17 microns per cubic meter (Kanpur, India in the dry season). Ammonium nitrate ranged from 0.2 microns per cubic meter (Mammoth Cave, in summer) to 6.8 microns per cubic meter (Kanpur, dry season). Equivalent black carbon ranged from 0.7 microns per cubic meter (Mammoth Cave) to over 8 microns per cubic meter (Dhaka, Bangladesh and Kanpur, India). Comparison of SPARTAN vs. coincident measurements from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network at Mammoth Cave yielded a high degree of consistency for daily PM (sub 2.5) (r squared equals 0.76, slope equals 1.12), daily sulfate (r squared equals 0.86, slope equals 1.03), and mean fractions of all major PM (sub 2.5) components (within 6 percent). Major ions generally agree well with previous studies at the same urban locations (e.g. sulfate fractions agree within 4 percent for 8 out of 11 collocation comparisons). Enhanced anthropogenic dust fractions in large urban areas (e.g. Singapore, Kanpur, Hanoi, and Dhaka) are apparent from high Zn to Al ratios. The expected water contribution to aerosols is calculated via the hygroscopicity parameter kappa (sub v (volume)) for each filter. Mean aggregate values ranged from 0.15 (Ilorin) to 0.28 (Rehovot). The all-site parameter mean is 0.20 plus or minus 0.04. Chemical composition and water retention in each filter measurement allows inference of hourly PM (sub 2.5) at 35 percent relative humidity by merging with nephelometer measurements. These hourly PM (sub 2.5) estimates compare favourably with a beta attenuation monitor (MetOne) at the nearby US embassy in Beijing, with a coefficient of variation r squared equals 0.67 (number equals 3167), compared to r squared equals 0.62 when v (volume) was not considered. SPARTAN continues to provide an open-access database of PM (sub 2.5) compositional filter information and hourly mass collected from a global federation of instruments.

Snider, Graydon↗

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Influence of Nucleation Mechanisms on the Radiative Properties of Deep Convective Clouds and Subvisible Cirrus in CRYSTAL/FACE

During the past few years we have conducted work on several different topics, as reflected by our publications. As one of the Co-Project scientists for The Cirrus Regional Study of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL FACE) we worked to help design the mission and then conduct it in the field. Another major activity during the past two years has been to pull together various groups to formulate plans for follow on missions to CRYSTAL FACE. We organized a workshop at the University of Colorado during the summer of 2003 to assess the best locations for future missions. Working with a group of about 10 scientists from around the country we prepared a science-planning document (Tropical Composition, Cloud and Climate Coupling Experiment (TC(sup 4)) that outlined the rationale, locations, strategy to accomplish the goals, and possible payloads for a set of three tropical missions. We also prepared background materials for various NRAs being prepared at NASA Headquarters for missions in Costa Rica, Darwin and Guam. In conjunction with the group at NASA Ames we have helped build a new numerical model for deep convection and have applied that model to simulate the CRYSTAL data. Our goal in particular has been to better understand how convection distributes water vapor isotopes. CRYSTAL observations of water isotopes are very different from those suggested by previous workers who assumed the isotopes would obey Rayleigh fractionation. The water isotope study has several implications. First it is a check on the realism of the deep convection model. Second, the isotopes are a measure of the precipitation removal in the atmosphere. Hence they provide a constraint on a parameter that is difficult to otherwise measure. Finally it has been suggested that isotopes may be the key to unraveling the water transport into the stratosphere and upper troposphere. Such transport is critical both for the radiation balance and for stratospheric chemistry. Ours is the first model that is able to treat this transport. Our initial results have just been submitted to Geophys. Res. Lett (Smith et al., 2005). Essentially we are able to explain the vertical profiles of isotopes in the tropical tropopause transition layer. We are also able to account for stratospheric humidity and isotope abundances with this model. We have also been heavily involved in trying to improve our understanding of nitric acid condensation on ice. Gao et al (2004) have shown that water supersaturations above ice occur when the atmosphere is supersaturated with respect to nitric acid trihydrate. As one of the co-authors of that work, we suggested the mechanism that may explain why this is occurring. Essentially, ice does not like to grow near unit supersaturation, but does so because the water molecules can find sites on the ice surface to attach themselves to before they fly off the ice surface. This phenomena was well known in the 1960s when it was a source of debate about whether condensation and evaporation coefficients for ice would be the same. Evaporation does not require any molecular orientation, while condensation does, so it was possible that the coefficients would differ. They don't differ because the water molecules rapidly move across the surface and find places to attach. Nitric acid may be occupying these preferred sites and therefore the water molecules can't find a desirable place to attach. We anticipate that this research will be the subject of laboratory work during the coming few years. Another possibility that has been suggested is that cubic ice is forming in clouds. We have measured the vapor pressure of cubic ice, and plan to publish that result in the next few months. We have also been working on additional aspects of the condensation of nitric acid on ice. With Y. Kondo we studied the condensation of NOy on ice using the SOLVE data. Gamblin et al. have continued this work. The CRYSTAL NOy and HNO, groups have shown th their data can be fit using standard Langmuir isotherms as suggested in some, but not all, laboratory studies. We have found in the SOLVE data set that this is not the case. Moreover some laboratory studies show there are important kinetic effects that may be occurring in the atmosphere limiting the transfer of nitric acid to the ice. The SOLVE data seem consistent with these studies. We are currently re-analyzing the CRYSTAL data to look for these kinetic effects. There are a number of implications of these studies. One of the more interesting is that the nitric acid coating on ice can be used as a cloud clock to determine how long the cloud parcel has been in existence. We have also been involved with several laboratory studies. We have worked to improve the database on ice optical constants, which are critical for remote sensing. We have also studied the ways in which ice nucleates on clays. We suspect now that the standard theories used for depositional ice nucleation are completely incorrect. Further work will be needed to develop a new theory.

Toon, Owen B.↗

Assessment of Impurities in Neptunium Metal Using Double-Focusing SF-ICP-MS

In the transmutation process, minor actinides such as neptunium, americium and curium obtained from spent nuclear fuel are incorporated into new fuel that can be placed into a fast reactor, where these elements fission into products with shorter half-lives which reduces the health and safety hazards. In order to better understand the all the involved processes, fundamental physical properties of actinides should be investigated, and for that reason, high-purity materials need to be produced. One of the main challenges in characterizing these materials is the lack of pure metal available for study, and lack of analytical methods to determine the purity. Trans-uranic elements are challenging to study because of the need to be handled in specialized facilities, safeguards concerns and they are not easily purchased and/or produced. High-purity Np metal is being produced at the Materials and Fuels Complex (MFC) of Idaho National Laboratory (INL), in a process of direct oxide reduction developed by Squires and Lessing, with a target purity of 99.999%. The goal of this research is to develop methods to determine the impurities in the Np metal, with detection limits low enough to achieve this purity level. In this work, we have used several ICP-MS and ICP OES techniques, but the focus was the method development using a recently acquired Double-Focusing Sector Field Inductively Coupled Plasma Mass Spectrometer (SF-ICP-MS). To produce the Np metal, approximately 20 g of the neptunium oxide powder was calcined at 800 oC for 6 h in a magnesium oxide crucible, which is placed inside a stainless steel crucible, that goes into the Hot Uniaxial Press (HUP) furnace insert. After the calcining step, the Np oxide was mixed with approximately 9 g of calcium metal and 100 g of ultra-dry calcium chloride. The HUP furnace was heated in a ramp rate of 50 oC min -1 up to 1000 oC. When the reaction mixture reached 850 oC, a tantalum stirrer was activated at 100 rpm. After the tantalum stirrer was removed and the mixture cooled overnight, the MgO crucible was opened with a hammer and chisel to retrieve the Np metal button at the bottom of the reaction mixture. Prior to the dissolution process, the surface contaminants of the sample were removed. The Np button was washed in 10 mL of ultra-pure water for 1 min, with a faint reaction observed on the surface, indicating the presence of Ca metal. The water-washed sample was placed in a tube with 2 M HCl for only 3 s, as a vigorous reaction in the button surface was observed. The leachates were saved for further analysis. The sample was then rinsed with ultra-pure water and air-dried prior to dissolution. The cleaned sample was placed in a PTFE beaker and 25 mL of 6 M HCl and 150 µL of 24 M HF were added. After 10 minutes of reaction, the beaker was placed on a hot plate for 1 h and 30 min, followed by the addition of 1 mL of 16 M HNO3 and allowed to react for a further 1 h and 20 min. The sample was removed from heat and 10 mL of HNO3 16 M was added to react with the HCl. The beaker was covered and left for overnight reaction, producing a solution with approximately 20 mL of final volume. A Quadrupole Inductively Coupled Plasma Mass Spectrometer (Q-ICP-MS, Perkin Elmer ELAN DRC II Q-ICP-MS) was used to obtain qualitative mass spectra for the dissolution blank, each of the leachates and the dissolved sample. In the water leachate spectrum mainly Ca and Mg were observed, while the acid wash had a similar peak pattern to the one observed for the Np metal sample, indicating that the process not only removed surface contamination, but also reacted with the sample itself. In the Np material spectrum, Ta and actinides were identified as candidates for major contaminants, as the most intense peaks corresponded to the mass-to-charge ratios of these elements. Minor contaminants included Al, Au, Ca, Cr, Fe, I, Na, Nb, Pb, Sn, Ti, Tl, V, Y. The presence of impurities in the material could have originated from several sources, such as the materials used for the casting process, the use of impure NpO2 or contamination of reagents and materials used for the analysis. Some of the isotopes listed as potential contaminants are susceptible to spectral interferences, so Double-Focusing SF-ICP-MS (Nu Instruments Attom ES HR-ICP-MS) was used for the quantitative analysis. For each isotope determined, several items were evaluated: raw signal of a blank solution, sensitivity and determination coefficient of calibration curves, possibility of interferences and comparison of results obtained for at least two resolving powers (RP): 300 and 4000. The evaluation of the signal in a blank solution can indicate if a spectral interference from the plasma gas and/or the solvent was successfully resolved (commonly observed for 52Cr+/40Ar12C+, for example). It is important to highlight that a signal reduction is expected when a higher RP is used, as the ion beam transmission is decreased by the use of narrower slits

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ASSESMENT OF IMPURITIES IN NP METAL USING DOUBLE-FOCUSING SF-ICP-MS

In the transmutation process, minor actinides such as neptunium, americium and curium obtained from spent nuclear fuel are incorporated into new fuel that can be placed into a fast reactor, where these elements fission into products with shorter half-lives which reduces the health and safety hazards [1]. In order to better understand the all the involved processes, fundamental physical properties of actinides should be investigated, and for that reason, high-purity materials need to be produced. One of the main challenges in characterizing these materials is the lack of pure metal available for study, and lack of analytical methods to determine the purity. Trans-uranic elements are challenging to study because of the need to be handled in specialized facilities, safeguards concerns and they are not easily purchased and/or produced. High-purity Np metal is being produced at the Materials and Fuels Complex (MFC) of Idaho National Laboratory (INL), in a process of direct oxide reduction developed by Squires and Lessing [2], with a target purity of 99.999%. The goal of this research is to develop methods to determine the impurities in the Np metal, with detection limits low enough to achieve this purity level. In this work, we have used several ICP-MS and ICP OES techniques, but the focus was the method development using a recently acquired Double-Focusing Sector Field Inductively Coupled Plasma Mass Spectrometer (SF-ICP-MS). To produce the Np metal, approximately 20 g of the neptunium oxide powder was calcined at 800 oC for 6 h in a magnesium oxide crucible, which is placed inside a stainless steel crucible, that goes into the Hot Uniaxial Press (HUP) furnace insert. After the calcining step, the Np oxide was mixed with approximately 9 g of calcium metal and 100 g of ultra-dry calcium chloride. The HUP furnace was heated in a ramp rate of 50 oC min -1 up to 1000 oC. When the reaction mixture reached 850 oC, a tantalum stirrer was activated at 100 rpm. After the tantalum stirrer was removed and the mixture cooled overnight, the MgO crucible was opened with a hammer and chisel to retrieve the Np metal button at the bottom of the reaction mixture [2]. Prior to the dissolution process, the surface contaminants of the sample were removed. The Np button was washed in 10 mL of ultra-pure water for 1 min, with a faint reaction observed on the surface, indicating the presence of Ca metal. The water-washed sample was placed in a tube with 2 M HCl for only 3 s, as a vigorous reaction in the button surface was observed. The leachates were saved for further analysis. The sample was then rinsed with ultra-pure water and air-dried prior to dissolution. The cleaned sample was placed in a PTFE beaker and 25 mL of 6 M HCl and 150 µL of 24 M HF were added. After 10 minutes of reaction, the beaker was placed on a hot plate for 1 h and 30 min, followed by the addition of 1 mL of 16 M HNO3 and allowed to react for a further 1 h and 20 min. The sample was removed from heat and 10 mL of HNO3 16 M was added to react with the HCl. The beaker was covered and left for overnight reaction, producing a solution with approximately 20 mL of final volume. A Quadrupole Inductively Coupled Plasma Mass Spectrometer (Q-ICP-MS, Perkin Elmer ELAN DRC II Q-ICP-MS) was used to obtain qualitative mass spectra for the dissolution blank, each of the leachates and the dissolved sample. In the water leachate spectrum mainly Ca and Mg were observed, while the acid wash had a similar peak pattern to the one observed for the Np metal sample, indicating that the process not only removed surface contamination, but also reacted with the sample itself. In the Np material spectrum, Ta and actinides were identified as candidates for major contaminants, as the most intense peaks corresponded to the mass-to-charge ratios of these elements. Minor contaminants included Al, Au, Ca, Cr, Fe, I, Na, Nb, Pb, Sn, Ti, Tl, V, Y. The presence of impurities in the material could have originated from several sources, such as the materials used for the casting process, the use of impure NpO2 or contamination of reagents and materials used for the analysis. Some of the isotopes listed as potential contaminants are susceptible to spectral interferences, so Double-Focusing SF-ICP-MS (Nu Instruments Attom ES HR-ICP-MS) was used for the quantitative analysis. For each isotope determined, several items were evaluated: raw signal of a blank solution, sensitivity and determination coefficient of calibration curves, possibility of interferences and comparison of results obtained for at least two resolving powers (RP): 300 and 4000. The evaluation of the signal in a blank solution can indicate if a spectral interference from the plasma gas and/or the solvent was successfully resolved (commonly observed for 52Cr+/40Ar12C+, for example). It is important to highlight that a signal reduction is expected when a higher RP is used, as the ion beam transmission is decreased by the use of na

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗