SEARCH · Engineering Papers

Results for “level set methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2026

North Slope of Alaska XSAPR b1 Data Processing Report: May 2025-April 2026

The U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) User Facility operates an X-band Scanning ARM Precipitation Radar (XSAPR) at the North Slope of Alaska (NSA) atmospheric observatory in Utqiaġvik. This radar provides year-round dual-polarization measurements of snow and ice in an arctic environment. Over the past few years, there has been greater focus on the NSA XSAPR, and recent efforts by ARM staff have resulted in quality-controlled, b-level data (Rocque et al. 2026b). This report outlines the processing of the second set of b-level data released for this radar from May 2025 through April 2026 following similar methods from Rocque et al. (2026b).

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1990

Compiling global name-space programs for distributed execution

Distributed memory machines do not provide hardware support for a global address space. Thus programmers are forced to partition the data across the memories of the architecture and use explicit message passing to communicate data between processors. The compiler support required to allow programmers to express their algorithms using a global name-space is examined. A general method is presented for analysis of a high level source program and along with its translation to a set of independently executing tasks communicating via messages. If the compiler has enough information, this translation can be carried out at compile-time. Otherwise run-time code is generated to implement the required data movement. The analysis required in both situations is described and the performance of the generated code on the Intel iPSC/2 is presented.

Koelbel, Charles↗

NASA NTRS2010

Method of Obtaining High Resolution Intrinsic Wire Boom Damping Parameters for Multi-Body Dynamics Simulations

The goal of NASA's Magnetospheric MultiScale (MMS) mission is to understand magnetic reconnection with sensor measurements from four spinning satellites flown in a tight tetrahedron formation. Four of the six electric field sensors on each satellite are located at the end of 60- meter wire booms to increase measurement sensitivity in the spin plane and to minimize motion coupling from perturbations on the main body. A propulsion burn however, might induce boom oscillations that could impact science measurements if oscillations do not damp to values on the order of 0.1 degree in a timely fashion. Large damping time constants could also adversely affect flight dynamics and attitude control performance. In this paper, we will discuss the implementation of a high resolution method for calculating the boom's intrinsic damping, which was used in multi-body dynamics simulations. In summary, experimental data was obtained with a scaled-down boom, which was suspended as a pendulum in vacuum. Optical techniques were designed to accurately measure the natural decay of angular position and subsequently, data processing algorithms resulted in excellent spatial and temporal resolutions. This method was repeated in a parametric study for various lengths, root tensions and vacuum levels. For all data sets, regression models for damping were applied, including: nonlinear viscous, frequency-independent hysteretic, coulomb and some combination of them. Our data analysis and dynamics models have shown that the intrinsic damping for the baseline boom is insufficient, thereby forcing project management to explore mitigation strategies.

Yew, Alvin G.↗

DOE OSTI2016

A Solid of Revolution Time Study using COG11.1 and MCNP6.1

The performance of two continuous energy Monte Carlo transport codes, COG11.1 and MCNP6.1, has been compared for the {sup 239}Pu Jezebel Benchmark using cross sections based on the ENDF/B-VII.1 evaluated nuclear data library. The purpose of this study is to determine the length of time it takes each code to complete a series of criticality calculations, using two different methods to model a solid of revolution. A solid of revolution is essential in many fields of engineering and mathematics, and has been a standard feature in computer-aided design and manufacturing since the 1970's. In the field of criticality safety, a solid of revolution can be used to reproduce and simplify complex engineering drawings and systems that would otherwise require excessive use of geometric primitives and computational resources to achieve a comparable level of fidelity. There are two basic methods to model a solid of revolution: Method 1 (COG) specifies a set of x-y pairs or polar coordinates and rotate the curve about the x-axis. Method 2 (COG and MCNP) creates a set of end-to-end conical frustums and/or cylinders. Method 1 is only available in COG and can be implemented using one surface and sector specification. Method 2 can be implemented in both COG and MCNP. This time study was initially conceived as a result of criticality safety analysis performed for a complex system with a high degree of rotational symmetry. Each COG input deck typically took 1-2 days to run in serial mode. COG input decks were then converted to MCNP format and run on the same computer cluster. Instead of expected 1-2 day runtimes, it took MCNP6.1 more than one week to complete each run. Since MCNP typically runs faster than many codes, these lengthy runtimes were considered unusual. This study attempts to investigate the runtime discrepancy between the two codes. The methods and benchmark in this time study were selected to provide a consistent, generalized approach to modeling an arbitrary solid of revolution. While there are many different ways to model the {sup 239}Pu Jezebel Benchmark, the focus of this study is on investigating the practical difference between the two basic methods to model a surface of revolution. Throughout this study, the only code that maintains reasonably low runtimes is COG11.1 (Method 1). Based on the observed performance of each code and method, it is likely that the divergent runtimes are the result of tracking particles across planes that lie between conical frustums, as opposed to defining the object boundary as a single surface. Slight fluctuations in graphed results were also observed, due to random sampling and CPU loading. This was confirmed by re-running input decks and comparing average runtimes to a preliminary curve fit. From the perspective of a criticality safety analyst, the results of this study are expected to contribute towards understanding one of the practical trade-offs between codes and methods as they apply to real-world analysis. Ultimately, it's about choosing the right tool for the job.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

Relativistic and Correlation Effects in CuH, AgH and AuH: Comparison of Various Relativistic Methods

The effects of relativity on the bond lengths, dissociation energies, and harmonic vibrational frequencies of the 1Epsilon(+) electronic ground states of the group IB hydrides CuH, AgH and AuH have been evaluated with a variety of ab initio methods. These properties were investigated with moderately-sized basis sets at the self-consistent field Hartree Fock (SCF HF) level and with second-order Moller-Plesset (MP2) perturbation theory for electron correlation. Comparisons were made between all-electron results using the nonrelativistic Hamiltonian, perturbation theory (PT) at first-order with only the one-electron non-fine structure terms of the Breit-Pauli Hamiltonian, the spin-free Douglas-Kroll (DK) transformed Dirac Hamiltonian and the untransformed Dirac Hamiltonian, and results using two sets of relativistic effective core potentials (RECPs). The expected trends of bond length decrease, dissociation energy increase and harmonic frequency increase with both relativity and correlation are found. Both sets of RECPs are shown to give good results, if accompanied by a reasonable basis set. The DK method is demonstrated to be an inexpensive, reliable approximation to the DHF method.

Collins, Charlene L.↗

NASA NTRS1994

Accuracy assessment of the large-scale dynamic ocean topography from TOPEX/POSEIDON altimetry

The quality of TOPEX/POSEIDON determinations of the global scale dynamic ocean topography have been assessed by determining mean topography solutions for successive 10-day repeat cycles and by examining the temporal changes in the sea surface topography to identify known features. The assessment is based on the analysis of TOPEX altimeter data cycles 1 through 36. Important errors in the tide model used to correct the altimeter data have been identified. The errors were reduced significantly by use of a new tide model derived with the TOPEX/POSEIDON measurements. Maps of the global 1-year mean topography, produced using four of the most accurate of the marine geoid, show that the largest error in the dynamic ocean topography show expected features, such as the known annual hemispherical sea surface rise and fall and the seasonal variability due to monsoon influence in the Indian Ocean. Changes in the sequence of 10-day topography maps show the development and propagation of an equatorial Kelvin wave in the Pacific beginning in December 1992 with a propagation velocity of approximately 3 m/s. The observations are consistent with observed changes in the equatorial trade winds, and with tide gauge and other in situ observations of the strengthening of the El Nino. Comparison of TOPEX-determine sea surface height at points near oceanic tide gauges shows agreement at the 4 cm root-mean-square (RMS) level over the tropical Pacific. The results show that the TOPEX altimeter data set can be used to map the ocean surface with a temporal resolution of 10 days and an accuracy which is insonsistent with traditional in situ methods for the determination of sea level variations.

Tapley, B. D.↗

DOE OSTI2019

Metallic Low-Level Waste (LLW) Handling Options: An Optimal Life cycle-Based Approach for North American Operations - 19413

Daily activities associated with nuclear power plant operations produce large volumes of low-level contaminated waste. During operation, a small percentage of this low-level waste is metallic (< 1% by volume). However, this figure increases materially when a facility is finally decommissioned, with approximately 1000 m{sup 3} of contaminated metal waste generated for every 1 GW decommissioned (13% by volume). With global nuclear decommissioning efforts expected to increase until 2050, the amount of low-level metallic waste generated is expected to increase the global inventory by over an order of magnitude. The high cost of conventional surface storage and burial of radioactive metal, coupled with the projected increase in low-level metallic waste volumes provides strong incentive to identify cost-effective waste management solutions. Surface decontamination and thermal methods are proven options for decontaminating low-level metallic waste. Surface decontamination processes include chemical, electrochemical, or mechanical methods, which are used to remove contamination from exposed layers of metal. Surface decontamination is most effective when metal has low levels of activity localized on the outer layers. When the activity is more uniformly distributed in the metal (or the geometry is complex), alternative methods such as thermal decontamination are required. In the thermal decontamination process, contaminated scrap metal is melted in a furnace and held in a molten state for a period of time while the radioactive material is separated out through gravity. Radioactive material reports to slag which is handled separately, or to the exhaust which is handled using conventional gas handling technologies. Thermal methods of decontamination involving metal melting present an opportunity to reduce long-term storage cost with a reduction in waste volume, along with the opportunity to reuse the decontaminated metal. Methods for dealing with metallic low-level waste from nuclear facilities vary extensively across the globe. The approach employed at any given site may be set by a host of factors ranging from regulatory influence, public perception, local economics, available geology for underground storage and often, precedent established long ago by early operations. This paper examines several low-level metallic waste handling options including techniques for surface decontamination and metal melting, as well as a combination of both. The various methods are categorized according to their suitability for waste characterization and local/global installations. The effectiveness of each option is benchmarked against a range of criteria, including capital and operating costs, environmental impact, overall sustainability and risk management issues. The overall objective of this paper is to develop a methodology for determining the most beneficial low-level metallic waste management strategy for a facility once the quantity, type and nature of the waste is established on a long-term basis. It is the authors' hope that the methodology can be used as a tool to test a facility's current plans (which are primarily focused on long-term storage) and identify areas for optimization with an emphasis on overall cost and environmental footprint reduction. The trend towards long term planning to handle large amounts of contaminated metal, along with obtaining a social license to implement a sustainable waste management solution, are key decision-making variables that will be discussed in relation to the LLW handling options presented in this paper. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1998

A Study of Fixed-Order Mixed Norm Designs for a Benchmark Problem in Structural Control

This study investigates the use of H2, p-synthesis, and mixed H2/mu methods to construct full-order controllers and optimized controllers of fixed dimensions. The benchmark problem definition is first extended to include uncertainty within the controller bandwidth in the form of parametric uncertainty representative of uncertainty in the natural frequencies of the design model. The sensitivity of H2 design to unmodelled dynamics and parametric uncertainty is evaluated for a range of controller levels of authority. Next, mu-synthesis methods are applied to design full-order compensators that are robust to both unmodelled dynamics and to parametric uncertainty. Finally, a set of mixed H2/mu compensators are designed which are optimized for a fixed compensator dimension. These mixed norm designs recover the H, design performance levels while providing the same levels of robust stability as the u designs. It is shown that designing with the mixed norm approach permits higher levels of controller authority for which the H, designs are destabilizing. The benchmark problem is that of an active tendon system. The controller designs are all based on the use of acceleration feedback.

Whorton, Mark S.↗

NASA NTRS1992

An approach to constrained aerodynamic design with application to airfoils

An approach was developed for incorporating flow and geometric constraints into the Direct Iterative Surface Curvature (DISC) design method. In this approach, an initial target pressure distribution is developed using a set of control points. The chordwise locations and pressure levels of these points are initially estimated either from empirical relationships and observed characteristics of pressure distributions for a given class of airfoils or by fitting the points to an existing pressure distribution. These values are then automatically adjusted during the design process to satisfy the flow and geometric constraints. The flow constraints currently available are lift, wave drag, pitching moment, pressure gradient, and local pressure levels. The geometric constraint options include maximum thickness, local thickness, leading-edge radius, and a 'glove' constraint involving inner and outer bounding surfaces. This design method was also extended to include the successive constraint release (SCR) approach to constrained minimization.

Campbell, Richard L.↗

DOE OSTI2013

Estimating nonrigid motion from inconsistent intensity with robust shape features

Purpose: To develop a nonrigid motion estimation method that is robust to heterogeneous intensity inconsistencies amongst the image pairs or image sequence. Methods: Intensity and contrast variations, as in dynamic contrast enhanced magnetic resonance imaging, present a considerable challenge to registration methods based on general discrepancy metrics. In this study, the authors propose and validate a novel method that is robust to such variations by utilizing shape features. The geometry of interest (GOI) is represented with a flexible zero level set, segmented via well-behaved regularized optimization. The optimization energy drives the zero level set to high image gradient regions, and regularizes it with area and curvature priors. The resulting shape exhibits high consistency even in the presence of intensity or contrast variations. Subsequently, a multiscale nonrigid registration is performed to seek a regular deformation field that minimizes shape discrepancy in the vicinity of GOIs. Results: To establish the working principle, realistic 2D and 3D images were subject to simulated nonrigid motion and synthetic intensity variations, so as to enable quantitative evaluation of registration performance. The proposed method was benchmarked against three alternative registration approaches, specifically, optical flow, B-spline based mutual information, and multimodality demons. When intensity consistency was satisfied, all methods had comparable registration accuracy for the GOIs. When intensities among registration pairs were inconsistent, however, the proposed method yielded pronounced improvement in registration accuracy, with an approximate fivefold reduction in mean absolute error (MAE = 2.25 mm, SD = 0.98 mm), compared to optical flow (MAE = 9.23 mm, SD = 5.36 mm), B-spline based mutual information (MAE = 9.57 mm, SD = 8.74 mm) and mutimodality demons (MAE = 10.07 mm, SD = 4.03 mm). Applying the proposed method on a real MR image sequence also provided qualitatively appealing results, demonstrating good feasibility and applicability of the proposed method. Conclusions: The authors have developed a novel method to estimate the nonrigid motion of GOIs in the presence of spatial intensity and contrast variations, taking advantage of robust shape features. Quantitative analysis and qualitative evaluation demonstrated good promise of the proposed method. Further clinical assessment and validation is being performed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2024

Developing a Control Strategy for Minimum Airflow Setting Considering CO2 Level and Energy Consumption in a Variable Air Volume System

In an office building equipped with a Variable Air Volume (VAV) system, this paper introduces a novel method for controlling the minimum supply airflow fraction in each zone’s VAV box, having a capability to consider indoor CO2 level and energy consumption. The EnergyPlus simulation using the medium office prototype model was employed, which evaluated the performance of the energy and CO2 concentration for five VAV box airflow control strategies. The paper focuses on CO2 concentration-based airflow control method and compares it with other four methods including conventional single-max, reduced minimum single-max, demand-controlled ventilation(DCV), and dualmax control methods according to guidelines and common practices. The newly proposed control strategy directly correlates the minimum airflow fraction to CO2 concentration. A general trend emerged when comparing CO2 concentrations—lower minimum airflow fractions were associated with higher concentrations. The proposed control method effectively maintained low CO2 concentrations and enabled a lower airflow fraction contributing to energy consumption reduction. It was confirmed that heating energy consumption in climate zone 4A, 5B, and 6A showed a maximum saving of approximately 30% compared to the conventional single-max and dual max control strategies. It was found that cooling energy consumption in climate zone 4A and 6A can achieve a maximum saving of approximately 10% compared to the conventional control strategies. The proposed CO2 concentration-based control logic is promising as it not only improves the indoor air quality lowering the CO2 concentration in the occupied spaces, but also contributes to HVAC energy savings.

Lee, Jong Man↗

NASA NTRS2011

Femoral Blood Flow and Cardiac Output During Blood Flow Restricted Leg Press Exercise

Low load blood flow restricted resistance exercise (LBFR) causes muscle hypertrophy that may be stimulated by the local ischemic environment created by the cuff pressure. However, local blood flow (BF) during such exercise is not well understood. PURPOSE: To characterize femoral artery BF and cardiac output (CO) during leg press exercise (LP) performed at a high load (HL) and low load (LL) with different levels of cuff pressure. METHODS: Eleven subjects (men/women 4/7, age 31.4+/-12.8 y, weight 68.9+/-13.2 kg, mean+/-SD) performed 3 sets of supine left LP to fatigue with 90 s of rest in 4 conditions: HL (%1-RM/cuff pressure: 80%/0); LL (20%/0); LBFR(sub DBP) (20%/1.3 x diastolic blood pressure, BP); LBFR(sub SBP) (20%/1.3 x supine systolic BP). The cuff remained inflated throughout the LBFR exercise sessions. Artery diameter, velocity time integral (VTI), and stroke volume (SV) were measured using Doppler ultrasound at rest and immediately after each set of exercise. Heart rate (HR) was monitored using a 3-lead ECG. BF was calculated as VTI x vessel cross-sectional area. CO was calculated as HR x SV. The data obtained after each set of exercise were averaged and used for analyses. Multi-level modeling was used to determine the effect of exercise condition on dependent variables. Statistical significance was set a priori at p< 0.05. RESULTS: Artery diameter did not change from baseline. BF increased (p<0.05) after exercise in each condition except LBFR(sub SBP) in the order of HL (12.73+/-1.42 cm3,mean+/-SE) > LL (9.92+/-0.82 cm3) > LBFR(sub dBP)(6.47+/-0.79 cm3) > LBFR(sub SBP) (3.51+/-0.59 cm3). Blunted exercise induced increases occurred in HR, SV, and CO after LBFR compared to HL and LL. HR increased 45% after HL and LL and 28% after LBFR (p<0.05), but SV increased (p<0.05) only after HL. Consequently, the increase (p<0.05) in CO was greater in HL and LL (approximately 3 L/min) than in LBFR (approximately 1 L/min). CONCLUSION: BF during LBFR(sub SBP) was 1/3 of that observed in LL, which supports the hypothesis that local ischemia stimulates the LBFR hypertrophic response. As the cuff did not compress the artery, the ischemia may have occurred because of the blunted rise in CO or because arterial BP cannot overcome the cuff pressure. As LBFR(sub DBP) effectively reduced BF and CO with cuff pressures less than systolic BP, future studies should investigate the hypertrophic potential of LBFR at even lower cuff pressures.

Everett, M. E.↗

DOE OSTI2020

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Optimizing Sensor Count and Placement to Detect Bond Wire Lift-Offs and Surface Defects in High-Power IGBT Modules Using Low-Cost Piezo-Electric Resonators

This manuscript presents the most recent results and findings to identify bond wire lift-offs and surface defects in high-power isolated gate bipolar junction transistor (IGBT) modules. The authors of this manuscript formerly proposed a low-cost, piezoelectric resonator-based measurement unit to detect bond wire related degradation in larger IGBT modules. Since high-power IGBT modules are expensive, it was apparent that evaluating our proposed method by inducing controlled damage to fresh (new) IGBTs may not be cost-effective especially when multiple sets of data need to be captured by inducing damage to different levels. In order to overcome this limitation, IGBT bond wires have been mimicked using a 3D printed enclosure, a PCB, and copper wires with dimensions very closely resembling a real IGBT. Using this method, multiple test devices can be built at the cost of a real IGBT, and the proposed technique could be fine-tuned without damaging expensive, real IGBTs. Our recent findings can be used to determine real IGBT degradation and bond wire lift-offs using only two sensors, as opposed to six transducers used in the first iteration of the setup. In addition to optimizing the sensor count, we have also identified the best possible locations of these sensors by attempting multiple placements inside the IGBT casing.

condition monitoring↗

DOE OSTI2017

Constraints on Galactic Neutrino Emission with Seven Years of IceCube Data

The origins of high-energy astrophysical neutrinos remain a mystery despite extensive searches for their sources. We present constraints from seven years of IceCube Neutrino Observatory muon data on the neutrino flux coming from the Galactic plane. This flux is expected from cosmic-ray interactions with the interstellar medium or near localized sources. Two methods were developed to test for a spatially extended flux from the entire plane, both of which are maximum likelihood fits but with different signal and background modeling techniques. We consider three templates for Galactic neutrino emission based primarily on gamma-ray observations and models that cover a wide range of possibilities. Based on these templates and in the benchmark case of an unbroken E{sup −2.5} power-law energy spectrum, we set 90% confidence level upper limits, constraining the possible Galactic contribution to the diffuse neutrino flux to be relatively small, less than 14% of the flux reported in Aartsen et al. above 1 TeV. A stacking method is also used to test catalogs of known high-energy Galactic gamma-ray sources.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

MFANS 2024 - Formally Proving Characteristics of Cyber-Physical Systems

Cyber-physical systems (CPS) are engineered systems that rely on the smooth integration of computational algorithms and physical elements. This integration presents new challenges for verifying that systems will behave as expected. The goal of this presentation is to present current challenges and potential solutions for the formal verification of cyber-physical systems. For cyber systems, formal methods refer to systematically rigorous mathematical techniques employed in the specification, development, analysis, and verification of both software and hardware systems. Recent advancements in computer science have yielded sophisticated tools specifically designed to address challenges associated with formal methods in complex systems. These tools leverage various foundational concepts such as logic, formal languages, program semantics, type systems, type theory, and automata theory. A notable achievement in the application of formal methods is the seL4 microkernel, claimed to be the first general-purpose operating-system kernel to be verified. Its proof implies the absence of bugs and guarantees that the kernel meets specifications. For physical systems, dynamic and control theory has a history of using rigorous analytic techniques to prove functional correctness. Lyapunov, optimal, classical, modern, and robust control theories all provide rigorous mathematical methods both to analyze system performance and to design controller that can be guaranteed to meet certain objectives. Recent computational techniques like level set theory and reachability analysis provide assertions that a system's state will avoid unsafe regions. Even though success has been independently achieved for cyber systems and physical systems, the integration of such systems creates new challenges. In particular, there is an obvious discrepancy between finite-state machines and infinite-state systems, resulting in different approaches for modeling and analyzing these system. While it is possible to simulate hybrid systems, this provides only a demonstration of a performance and not proof. For hybrid systems, current formal methods and system analysis approaches typically require a workarounds to work on hybrid systems like CPS. This paper will outline the state of the art and limits of current practice for formally verifying CPS and will identify possible research directions that require attention.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1992

Validating a large geophysical data set: Experiences with satellite-derived cloud parameters

We are validating the global cloud parameters derived from the satellite-borne HIRS2 and MSU atmospheric sounding instrument measurements, and are using the analysis of these data as one prototype for studying large geophysical data sets in general. The HIRS2/MSU data set contains a total of 40 physical parameters, filling 25 MB/day; raw HIRS2/MSU data are available for a period exceeding 10 years. Validation involves developing a quantitative sense for the physical meaning of the derived parameters over the range of environmental conditions sampled. This is accomplished by comparing the spatial and temporal distributions of the derived quantities with similar measurements made using other techniques, and with model results. The data handling needed for this work is possible only with the help of a suite of interactive graphical and numerical analysis tools. Level 3 (gridded) data is the common form in which large data sets of this type are distributed for scientific analysis. We find that Level 3 data is inadequate for the data comparisons required for validation. Level 2 data (individual measurements in geophysical units) is needed. A sampling problem arises when individual measurements, which are not uniformly distributed in space or time, are used for the comparisons. Standard 'interpolation' methods involve fitting the measurements for each data set to surfaces, which are then compared. We are experimenting with formal criteria for selecting geographical regions, based upon the spatial frequency and variability of measurements, that allow us to quantify the uncertainty due to sampling. As part of this project, we are also dealing with ways to keep track of constraints placed on the output by assumptions made in the computer code. The need to work with Level 2 data introduces a number of other data handling issues, such as accessing data files across machine types, meeting large data storage requirements, accessing other validated data sets, processing speed and throughput for interactive graphical work, and problems relating to graphical interfaces.

Kahn, Ralph↗

DOE OSTI2013

Development of an Absolute Gas-Counting Capability for Low to Medium Activities

ABSTRACT Pacific Northwest National Laboratory (PNNL) is developing a capability to measure the absolute activity concentration of various gaseous radionuclides using length-compensated proportional-counting methods. This capability will enable the validation and use of low-level, gaseous radionuclide calibration standards for use in PNNL’s shallow underground laboratory. Two sets of unequal length proportional counters have been fabricated. These detector assemblies operate on a static gas-fill principle, in contrast to continuous, flow-through configurations. One set of three counters has been fabricated using ultra-low background (ULB) electroformed copper and low-background fabrication methods. Once fully operational, these ULB counters will be used in PNNL’s shallow underground counting laboratory for analysis of gases with low activity concentrations < 1 Bq/cc. A second set of four unequal length counters has been fabricated from Oxygen-Free High-Conductivity Copper (OFHC) using similar low-background cleaning and assembly methods. These OFHC counters will be operated above ground in the analysis of gases with activity concentrations in the range of 1-10 Bq/cc. A gas delivery system is being developed to actively mix the analyte gas with an appropriate amount of count-gas and uniformly deliver it to the counters with high accuracy and repeatability. A description of both detector assemblies and gas delivery system will be given along with a preliminary uncertainty analysis of a simulated 0.05 Bq/cm3 gas measurement.

Williams, Richard M.↗