Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observational 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.

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

At least 1,099 records · Page 61

Variance Estimation in Monte Carlo Eigenvalue Simulations Using Spectral Analysis Method

Monte Carlo (MC) simulation is widely used to solve the eigenvalue form of the Boltzmann transport equation that mathematically represents the neutron transport process through complex multiplying (fissionable) systems. Monte Carlo eigenvalue simulation starts with an assumed fission source distribution and uses the fission sites from the previous iteration (cycle) as the starting source in the current iteration. Important system parameters (MC tallies) such as fuel pin-power distribution are estimated over several cycles after the convergence of the fission source distribution to a stationary distribution. However, the MC fission source iteration algorithm that uses fission source sites from the previous cycle introduces a cycle-to-cycle correlation. Monte Carlo simulations that do not account for the cycle-to-cycle correlation systematically underestimate the variance of the estimated system parameters (sample mean). This paper presents the relationship between the spectral density in the frequency domain at frequency zero and the variance of the sample mean. This paper introduces a novel method in the frequency domain for the MC variance estimation. For the three test problems used in this paper, researchers have observed that the new method results in an improvement of more than one order of magnitude to the standard deviation of the sample mean. The new method also compares favorably with the previously introduced batch, bootstrap, and covariance-adjusted methods when applied to the three test problems investigated in this paper. This new method does not require modification of the MC eigenvalue algorithm (power iteration), is code agnostic, and is therefore easy to use when implementing in any existing MC code. In conclusion, the new estimate can be calculated without saving tally results of all active/stationary cycles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Broad absorption line quasars in the Dark Energy Spectroscopic Instrument Early Data Release

Broad absorption line (BAL) quasars are characterized by gas clouds that absorb flux at the wavelength of common quasar spectral features, although blueshifted by velocities that can exceed $0.1c$. BAL features are interesting as signatures of significant feedback, yet they can also compromise cosmological studies with quasars by distorting the shape of the most prominent quasar emission lines, impacting redshift accuracy and measurements of the matter density distribution traced by the Lyman $\alpha$ forest. We present a catalogue of BAL quasars discovered in the Dark Energy Spectroscopic Instrument (DESI) survey Early Data Release, which were observed as part of DESI Survey Validation, as well as the first two months of the main survey. We describe our method to automatically identify BAL quasars in DESI data, the quantities we measure for each BAL, and investigate the completeness and purity of this method with mock DESI observations. We mask the wavelengths of the BAL features and re-evaluate each BAL quasar redshift, finding new redshifts which are $243\, {\rm km}\, {\rm s}^{-1}$ smaller on average for the BAL quasar sample. These new, more accurate redshifts are important to obtain the best measurements of quasar clustering, especially at small scales. Finally, we present some spectra of rarer classes of BALs that illustrate the potential of DESI data to identify such populations for further study.

79 ASTRONOMY AND ASTROPHYSICS↗

Scatter and crosstalk corrections for {sup 99m}Tc/{sup 123}I dual-radionuclide imaging using a CZT SPECT system with pinhole collimators

Purpose: The energy spectrum for a cadmium zinc telluride (CZT) detector has a low energy tail due to incomplete charge collection and intercrystal scattering. Due to these solid-state detector effects, scatter would be overestimated if the conventional triple-energy window (TEW) method is used for scatter and crosstalk corrections in CZT-based imaging systems. The objective of this work is to develop a scatter and crosstalk correction method for {sup 99m}Tc/{sup 123}I dual-radionuclide imaging for a CZT-based dedicated cardiac SPECT system with pinhole collimators (GE Discovery NM 530c/570c). Methods: A tailing model was developed to account for the low energy tail effects of the CZT detector. The parameters of the model were obtained using {sup 99m}Tc and {sup 123}I point source measurements. A scatter model was defined to characterize the relationship between down-scatter and self-scatter projections. The parameters for this model were obtained from Monte Carlo simulation using SIMIND. The tailing and scatter models were further incorporated into a projection count model, and the primary and self-scatter projections of each radionuclide were determined with a maximum likelihood expectation maximization (MLEM) iterative estimation approach. The extracted scatter and crosstalk projections were then incorporated into MLEM image reconstruction as an additive term in forward projection to obtain scatter- and crosstalk-corrected images. The proposed method was validated using Monte Carlo simulation, line source experiment, anthropomorphic torso phantom studies, and patient studies. The performance of the proposed method was also compared to that obtained with the conventional TEW method. Results: Monte Carlo simulations and line source experiment demonstrated that the TEW method overestimated scatter while their proposed method provided more accurate scatter estimation by considering the low energy tail effect. In the phantom study, improved defect contrasts were observed with both correction methods compared to no correction, especially for the images of {sup 99m}Tc in dual-radionuclide imaging where there is heavy contamination from {sup 123}I. In this case, the nontransmural defect contrast was improved from 0.39 to 0.47 with the TEW method and to 0.51 with their proposed method and the transmural defect contrast was improved from 0.62 to 0.74 with the TEW method and to 0.73 with their proposed method. In the patient study, the proposed method provided higher myocardium-to-blood pool contrast than that of the TEW method. Similar to the phantom experiment, the improvement was the most substantial for the images of {sup 99m}Tc in dual-radionuclide imaging. In this case, the myocardium-to-blood pool ratio was improved from 7.0 to 38.3 with the TEW method and to 63.6 with their proposed method. Compared to the TEW method, the proposed method also provided higher count levels in the reconstructed images in both phantom and patient studies, indicating reduced overestimation of scatter. Using the proposed method, consistent reconstruction results were obtained for both single-radionuclide data with scatter correction and dual-radionuclide data with scatter and crosstalk corrections, in both phantom and human studies. Conclusions: The authors demonstrate that the TEW method leads to overestimation in scatter and crosstalk for the CZT-based imaging system while the proposed scatter and crosstalk correction method can provide more accurate self-scatter and down-scatter estimations for quantitative single-radionuclide and dual-radionuclide imaging.

60 APPLIED LIFE SCIENCES↗

GenSpec: A Genetic Algorithm for Neutron Energy Spectrum Adjustment

Presented in this report is the description of a new method for neutron energy spectrum adjustment which uses a genetic algorithm to minimize the difference between calculated and measured reaction probabilities. The measured reaction probabilities are found using neutron activation analysis. The method adjusts a trial spectrum provided by the user which is typically calculated using a neutron transport code such as MCNP. Observed benefits of this method over currently existing methods include the reduction in unrealistic artefacts in the spectral shape as well as a reduced sensitivity to increases in the energy resolution of the derived spectrum. This report presents the adjustment results for various spectrum altering bucket environments in the central cavity of the Annular Core Research Reactor, as well as the adjustment results for the spectrum in the Sandia Pulse Reactor III. In each case, the results are compared to those generated using LSL-M2, which is a code commonly used for the purpose of spectrum adjustment. The genetic algorithm produces spectrum-averaged reaction probabilities with agreement to measured values, and comparable to those resulting from LSL-M2. The true benefit to this method, the reduction of shape artefacts in the spectrum, is difficult to quantify but can be clearly seen in the comparison of the final adjustments. Beyond these preliminary results, this report also gives a thorough description of the genetic algorithm and presents instructions for running the code using the graphical user interface. In its present state, the code does not provide uncertainties or correlations for the adjusted spectrum. This capability is currently being added, and will be presented in future work.

42 ENGINEERING↗

The cosmic Mach number - Direct comparisons of observations and models

A robust method is used to derive the cosmic Mach number, M, from existing peculiar velocity data sets. For the IR Tully-Fisher distances of spirals in the Local Supercluster of Aaronson et al. (1979, 1982), M of 1.03 and a characteristic distance from the Local Group of 1464 km/s are found. For the elliptical galaxy sample of Faber et al. (1987, 1988), M of 0.57, with a characteristic distance of 2572 km/s is found. These results are compared with Monte Carlo simulations of the observational realizations drawn from numerical simulations of the universe based on various scenarios, including the standard cold dark matter (CDM) scenario. The effect of velocity bias on the derived Mach number is found to be small. Only 5 percent of the CDM simulations have M as large as or larger than that derived for the Aaronson et al. sample.

Strauss, Michael A.↗

Nano-engineered Multiwall Carbon Nanotube-copper Composite Thermal Interface Material for Efficient Heat Conduction

Efforts in integrated circuit (IC) packaging technologies have recently been focused on management of increasing heat density associated with high frequency and high density circuit designs. While current flip-chip package designs can accommodate relatively high amounts of heat density, new materials need to be developed to manage thermal effects of next-generation integrated circuits. Multiwall carbon nanotubes (MWNT) have been shown to significantly enhance thermal conduction in the axial direction and thus can be considered to be a candidate for future thermal interface materials by facilitating efficient thermal transport. This work focuses on fabrication and characterization of a robust MWNT-copper composite material as an element in IC package designs. We show that using vertically aligned MWNT arrays reduces interfacial thermal resistance by increasing conduction surface area, and furthermore, the embedded copper acts as a lateral heat spreader to efficiently disperse heat, a necessary function for packaging materials. In addition, we demonstrate reusability of the material, and the absence of residue on the contacting material, both novel features of the MWNT-copper composite that are not found in most state-of-the-art thermal interface materials. Electrochemical methods such as metal deposition and etch are discussed for the creation of the MWNT-Cu composite, detailing issues and observations with using such methods. We show that precise engineering of the composite surface affects the ability of this material to act as an efficient thermal interface material. A thermal contact resistance measurement has been designed to obtain a value of thermal contact resistance for a variety of different thermal contact materials.

Ngo, Quoc↗

Black Hole Mergers as Probes of Structure Formation

Intense structure formation and reionization occur at high redshift, yet there is currently little observational information about this very important epoch. Observations of gravitational waves from massive black hole (MBH) mergers can provide us with important clues about the formation of structures in the early universe. Past efforts have been limited to calculating merger rates using different models in which many assumptions are made about the specific values of physical parameters of the mergers, resulting in merger rate estimates that span a very wide range (0.1 - 104 mergers/year). Here we develop a semi-analytical, phenomenological model of MBH mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z 5 - 30 structure formation era. We do this by generating synthetic LISA observable data (total BH mass, BH mass ratio, redshift, merger rates), which are then analyzed using a Markov Chain Monte Carlo method. This allows us to constrain the physical parameters of the mergers. We find that our methodology works well at estimating merger parameters, consistently giving results within 1- of the input parameter values. We also discover that the number of merger events is a key discriminant among models. This helps our method be robust against observational uncertainties. Our approach, which at this stage constitutes a proof of principle, can be readily extended to physical models and to more general problems in cosmology and gravitational wave astrophysics.

Alicea-Munoz, E.↗

Diagnosis of Magnetic Structures and Intermittency in Space Plasma Turbulence using the Method of Surrogate Data

Several observations in space plasmas have reported the presence of coherent structures at different plasma scales. Structure formation is believed to be a direct consequence of nonlinear interactions between the plasma modes, which depend strongly on phase synchronization of those modes. Despite this important role of the phases in turbulence, very limited work has been however devoted to study the phases as a potential tracers of nonlinearities in comparison with the wealth of literature on power spectra of turbulence where phases are totally missed. We present a method based on surrogate data to systematically detect coherent structures in turbulent signals. The new method has been applied successfully to magnetosheath turbulence (Sahraoui, Phys. Rev. E, 2008, in press), where the relationship between the identified phase coherence and intermittency (classically identified as non Gaussian tails of the PDFs) as well as the energy cascade has been studied. Here we review the main results obtained in that study and show further applications to small scale solar wind turbulence. Implications of the results on theoretical modelling of space turbulence (applicability of weak/wave turbulence, its validity limits and its connection to intermittency) will be discussed.

Sahraoui, Fouad↗

Radiometric calibration of the Earth observing system's imaging sensors

Philosophy, requirements, and methods of calibration of multispectral space sensor systems as applicable to the Earth Observing System (EOS) are discussed. Vicarious methods for calibration of low spatial resolution systems, with respect to the Advanced Very High Resolution Radiometer (AVHRR), are then summarized. Finally, a theoretical introduction is given to a new vicarious method of calibration using the ratio of diffuse-to-global irradiance at the Earth's surfaces as the key input. This may provide an additional independent method for in-flight calibration.

Slater, P. N.↗

On determining magnetospheric diffusion coefficients from the observed effects of Jupiter's satellite Io

A method is derived for determining the radial diffusion coefficient from observed satellite effects of the inner Jovian satellites on the energetic particle fluxes. The method is based on data from L values which are significantly removed from the actual sweeping region. With regard to the large losses to the protons at Io's L shell, it is suggested that in addition to satellite sweepup, the losses may be associated with an enhanced precipitation due to resonant interaction with ion cyclotron waves near Io's orbit. It is noted that such additional loss mechanisms may also apply to electrons, and that such losses may significantly affect the estimated diffusion coefficient.

Thomsen, M. F.↗

Estimation of time averages from irregularly spaced observations - With application to coastal zone color scanner estimates of chlorophyll concentration

The sampling error of an arbitrary linear estimate of a time-averaged quantity constructed from a time series of irregularly spaced observations at a fixed located is quantified through a formalism. The method is applied to satellite observations of chlorophyll from the coastal zone color scanner. The two specific linear estimates under consideration are the composite average formed from the simple average of all observations within the averaging period and the optimal estimate formed by minimizing the mean squared error of the temporal average based on all the observations in the time series. The resulting suboptimal estimates are shown to be more accurate than composite averages. Suboptimal estimates are also found to be nearly as accurate as optimal estimates using the correct signal and measurement error variances and correlation functions for realistic ranges of these parameters, which makes it a viable practical alternative to the composite average method generally employed at present.

Chelton, Dudley B.↗

MicroPulse DIAL (MPD) – a diode-laser-based lidar architecture for quantitative atmospheric profiling

Abstract. Continuous water vapor and temperature profiles are critically needed for improved understanding of the lower atmosphere and potential advances in weather forecasting skill. Ground-based, national-scale profiling networks are part of a suite of instruments to provide such observations; however, the technological method must be cost-effective and quantitative. We have been developing an active remote sensing technology based on a diode-laser-based lidar technology to address this observational need. Narrowband, high-spectral-fidelity diode lasers enable accurate and calibration-free measurements requiring a minimal set of assumptions based on direct absorption (Beer–Lambert law) and a ratio of two signals. These well-proven quantitative methods are known as differential absorption lidar (DIAL) and high-spectral-resolution lidar (HSRL). This diode-laser-based architecture, characterized by less powerful laser transmitters than those historically used for atmospheric studies, can be made eye-safe and robust. Nevertheless, it also requires solar background suppression techniques such as narrow-field-of-view receivers with an ultra-narrow bandpass to observe individual photons backscattered from the atmosphere. We discuss this diode-laser-based lidar architecture's latest generation and analyze how it addresses a national-scale profiling network's need to provide continuous thermodynamic observations. The work presented focuses on general architecture changes that pertain to both the water vapor and the temperature profiling capabilities of the MicroPulse DIAL (MPD). However, the specific subcomponent testing and instrument validation presented are for the water vapor measurements only. A fiber-coupled seed laser transmitter optimization is performed and shown to meet all of the requirements for the DIAL technique. Further improvements – such as a fiber-coupled near-range receiver, the ability to perform quality control via automatic receiver scanning, advanced multi-channel scalar capabilities, and advanced processing techniques – are discussed. These new developments increase narrowband DIAL technology readiness and are shown to allow higher-quality water vapor measurements closer to the surface via preliminary intercomparisons within the MPD network itself and with radiosondes.

Spuler, Scott M.↗

Magnetic Field Intermittency in the Solar Wind: Parker Solar Probe and SolO Observations Ranging from the Alfvén Region up to 1AU

Parker Solar Probe (PSP) and SolO data are utilized to investigate magnetic field intermittency in the solar wind (SW). Small-scale intermittency (20−100 di) is observed to radially strengthen when methods relying on higher-order moments are considered (SFq; SDK), but no clear trend is observed at larger scales. However, lower-order moment-based methods (e.g., partial variance of increments; PVI) are deemed more appropriate for examining the evolution of the bulk of coherent structures (CSs), PVI ≥ 3. Using PVI, we observe a scale-dependent evolution in the fraction of the data set occupied by CSs, fPVI≥3. Specifically, regardless of the SW speed, a subtle increase is found in fPVI≥3 for ℓ = 20 di, in contrast to a more pronounced radial increase in CSs observed at larger scales. Intermittency is investigated in relation to plasma parameters. Though, slower SW speed intervals exhibit higher fPVI≥6 and higher kurtosis maxima, no statistical differences are observed for fPVI≥3. Highly Alfvénic intervals display lower levels of intermittency. The anisotropy with respect to the angle between the magnetic field and SW flow, ΘVB is investigated. Intermittency is weaker at ΘVB ≈ 0° and is strengthened at larger angles. Considering the evolution at a constant alignment angle, a weakening of intermittency is observed with increasing advection time of the SW. Our results indicate that the strengthening of intermittency in the inner heliosphere is driven by the increase in comparatively highly intermittent perpendicular intervals sampled by the probes with increasing distance, an effect related directly to the evolution of the Parker spiral.

Nikos Sioulas↗

Magnetic force microscopy study of domain walls in Co{sub 2}Z ferrite

Graphical abstract: - Highlights: • Hexaferrite Co{sub 2}Z is synthesized through the modified Pechini method. • Magnetic domains are observed in anisotropic Co{sub 2}Z single grain using MFM. • Observed single grain domain thickness is in good agreement with Dotsh model. - Abstract: Hexaferrite Co{sub 2}Z was synthesized through the modified Pechini method. Partially oriented samples were obtained after consolidation with uniaxial pressing and calcination/sintering at 1300 °C/1330 °C. The sample composition and morphology was identified with X-ray diffractometry (XRD) and scanning electron microscopy (SEM) with energy-dispersive X-ray spectrometry (EDS). MFM studies of the single grains revealed a domain structure with 0.7 μm wide. The Co{sub 2}Z static magnetization was measured with a vibrating sample magnetometer (VSM), and was used to calculate a single grain domain with a thickness of 4.8 μm. This result is in good agreement with SEM observations of the single grain thickness.

36 MATERIALS SCIENCE↗

Moving beyond the total sea ice extent in gauging model biases

Here, reproducing characteristics of observed sea ice extent remains an important climate modeling challenge. This study describes several approaches to improve how model biases in total sea ice distribution are quantified, and applies them to historically forced simulations contributed to phase 5 of the Coupled Model Intercomparison Project (CMIP5). The quantity of hemispheric total sea ice area, or some measure of its equatorward extent, is often used to evaluate model performance. A new approach is introduced that investigates additional details about the structure of model errors, with an aim to reduce the potential impact of compensating errors when gauging differences between simulated and observed sea ice. Using multiple observational datasets, several new methods are applied to evaluate the climatological spatial distribution and the annual cycle of sea ice cover in 41 CMIP5 models. It is shown that in some models, error compensation can be substantial, for example resulting from too much sea ice in one region and too little in another. Error compensation tends to be larger in models that agree more closely with the observed total sea ice area, which may result from model tuning. The results herein suggest that consideration of only the total hemispheric sea ice area or extent can be misleading when quantitatively comparing how well models agree with observations. Further work is needed to fully develop robust methods to holistically evaluate the ability of models to capture the finescale structure of sea ice characteristics; however, the “sector scale” metric used here aids in reducing the impact of compensating errors in hemispheric integrals.

54 ENVIRONMENTAL SCIENCES↗

Minimum-size design of regulation systems and the application to Space Station

A method is proposed to determine the sizes of the regulator and observer and to design the feedback and Kalman gain matrices for the reduced system. The sizes determined are minimum in the sense of satisfying the specified eigenvalue limits when the interaction between the regulator and the observer can be neglected. The resulting sizes will depend on the properties of the original system and the desired response. If the interaction between the regulator and the observer becomes substantial, more modes in the observer should be included, or other methods should be used for reducing the observation spillover or control spillover.

Reddy, A. S. S. R.↗

ICESat Observations of Inland Surface Water Stage, Slope, and Extent: a New Method for Hydrologic Monitoring

River discharge and changes in lake, reservoir and wetland water storage are critical terms in the global surface water balance, yet they are poorly observed globally and the prospects for adequate observations from in-situ networks are poor (Alsdorf et al., 2003). The NASA-sponsored Surface Water Working Group has established a framework for advancing satellite observations of river discharge and water storage changes which focuses on obtaining measurements of water surface height (stage), slope, and extent. Satellite laser altimetry, which can achieve centimeter-level elevation precision for single, small laser footprints, provides a method to obtain these inland water parameters and contribute to global water balance monitoring. Since its launch in January, 2003 the Ice, Cloud, and land Elevation Satellite (ICESat), a NASA Earth Observing System mission, has achieved over 540 million laser pulse observations of ice sheet, ocean surface, land topography, and inland water elevations and cloud and aerosol height distributions. By recording the laser backscatter from 80 m diameter footprints spaced 175 m along track, ICESat acquires globally-distributed elevation profiles, using a 1064 nm laser altimeter channel, and cloud and aerosol profiles, using a 532 nm atmospheric lidar channel. The ICESat mission has demonstrated the following laser altimeter capabilities relevant to observations of inland water: (1) elevation measurements with a precision of 2 to 3 cm for flat surfaces, suitable for detecting river surface slopes along long river reaches or between multiple crossings of a meandering river channel, (2) from the laser backscatter waveform, detection of water surface elevations beneath vegetation canopies, suitable for measuring water stage in flooded forests, (3) single pulse absolute elevation accuracy of about 50 cm (1 sigma) for 1 degree sloped surfaces, with calibration work in progress indicating that a final accuracy of about 12 cm (1 sigma) will be achieved for clear atmosphere conditions, suitable for detection of stage changes through time, (4) ability to precisely point the spacecraft so as to position the laser profile on the Earth the surface with a cross-track accuracy of 50 m (1 sigma), enabling small water bodies and specific locations to be targeted and re-observed through time, (5) adequate signal levels from specular water surfaces up to 5 degrees off-nadir, enabling complete global access to any location on the Earth's surface from the ICESat repeat orbit by off-nadir pointing, and (6) day and night operation with successful laser ranging to the Earth's surface through thin to moderate cloud cover, enabling more frequent measurements than can be achieved by passive optical sensors. Here we illustrate these capabilities by showing ICESat observations through time for selected river and lake locations.

Harding, David J.↗

A Method to Retrieve Rainfall Rate over Land from TRMM Observations

Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR) observations over mesoscale convective systems (MCSs) reveal that there are localized maxima in the rain rate with a scale of about 10 to 20 km that represent thunderstorms (Cbs). Some of these Cbs are developing or intense, while others are decaying or weak. These Cbs constitute only about 20 % of the rain area of a given MCS. Outside of Cbs, the average rain rate is much weaker than that within Cbs. From an analysis of the PR data, we find that the spatial distribution of rain and its character, convective or stratiform, is highly inhomogeneous. This complex nature of rain exists on a scale comparable to that of a Cb. The 85 GHz brightness temperature, T85, observations of the TRMM Microwave Imager (TMI) radiometer taken over an MCS reflect closely the PR rain rate pattern over land. Local maxima in rain rate shown by PR are observed as local minima in T85. Where there are no minima in T85, PR observations indicate there is light rain. However, the TMI brightness temperature measurements (Tbs) have poor ability to discriminate convective rain from stratiform rain. For this reason, a TMI rain retrieval procedure that depends primarily on the magnitude of Tbs performs poorly. In order to retrieve rain rate from TMI data on land one has to include the spatial distribution information deduced from the T85 data in the retrieval method. Then, quantitative estimation of rain rate can be accomplished. A TMI rain retrieval method developed along these lines can yield estimates of rain rate and its frequency distribution which agree closely with that given by PR. We find the current TRMM project TMI (Version 5) rain retrieval algorithm on land could be improved with the retrieval scheme developed here. To support the conceptual frame work of the rain retrieval method developed here, a theoretical analysis of the TMI brightness temperatures in convective and stratiform regions is presented.

Prabhakara, C.↗