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 721 records · Page 40

Comparison of Two Acoustic Waveguide Methods for Determining Liner Impedance

Acoustic measurements taken in a flow impedance tube are used to assess the relative accuracy of two waveguide methods for impedance eduction in the presence of grazing flow. The aeroacoustic environment is assumed to contain forward and backward-traveling acoustic waves, consisting of multiple modes, and uniform mean flow. Both methods require a measurement of the complex acoustic pressure profile over the length of the test liner. The Single Mode Method assumes that the sound pressure level and phase decay-rates of a single progressive mode can be extracted from this measured complex acoustic pressure profile. No a priori assumptions are made in the Finite Element. Method regarding the modal or reflection content in the measured acoustic pressure profile. The integrity of each method is initially demonstrated by how well their no-flow impedances match those acquired in a normal incidence impedance tube. These tests were conducted using ceramic tubular and conventional perforate liners. Ceramic tubular liners were included because of their impedance insensitivity to mean flow effects. Conversely, the conventional perforate liner was included because its impedance is known to be sensitive to mean flow velocity effects. Excellent comparisons between impedance values educed with the two waveguide methods in the absence of mean flow and the corresponding values educed with the normal incident impedance tube were observed. The two methods are then compared for mean flow Mach numbers up to 0.5, and are shown to give consistent results for both types of test liners. The quality of the results indicates that the Single Mode Method should be used when the measured acoustic pressure profile is clearly dominated by a single progressive mode, and the Finite Element Method should be used for all other cases.

Jones, Michael G.↗

Determination of variations of the solar radius from solar eclipse observations

This paper describes the method to determine the solar radius and its variations from observations made during total solar eclipses. In particular, the procedure to correct the spherical moon predictions for the effects of lunar mountains and valleys on the width and location of the path of totality is addressed in detail. The errors affecting this technique are addressed, a summary of the results of its application to three solar eclipses are presented, and the implications of the results on the constancy of the solar constant are described.

Sofia, S.↗

Observational constraints on Tachyon and DBI inflation

We present a systematic method for evaluation of perturbation observables in non-canonical single-field inflation models within the slow-roll approximation, which allied with field redefinitions enables predictions to be established for a wide range of models. We use this to investigate various non-canonical inflation models, including Tachyon inflation and DBI inflation. The Lambert W function will be used extensively in our method for the evaluation of observables. In the Tachyon case, in the slow-roll approximation the model can be approximated by a canonical field with a redefined potential, which yields predictions in better agreement with observations than the canonical equivalents. For DBI inflation models we consider contributions from both the scalar potential and the warp geometry. In the case of a quartic potential, we find a formula for the observables under both non-relativistic (sound speed c{sub s}{sup 2} ∼ 1) and relativistic behaviour (c{sub s}{sup 2} || 1) of the scalar DBI inflaton. For a quadratic potential we find two branches in the non-relativistic c{sub s}{sup 2} ∼ 1 case, determined by the competition of model parameters, while for the relativistic case c{sub s}{sup 2} → 0, we find consistency with results already in the literature. We present a comparison to the latest Planck satellite observations. Most of the non-canonical models we investigate, including the Tachyon, are better fits to data than canonical models with the same potential, but we find that DBI models in the slow-roll regime have difficulty in matching the data.

79 ASTRONOMY AND ASTROPHYSICS↗

High speed spectroscopy using temporal positioned optical fibers with an optical scanner mirror

Systems and methods for examining spectral data over the course of a high speed event are described. The systems and methods can enable observation of the spectral evolution of a transient phenomenon into segment intervals on the order of, milliseconds or microseconds. The methods include reflecting light from an event off of a rotating mirror and sequentially delivering light from the mirror to a series of optical waveguides for sequential spectral analysis. The systems and methods can be useful in a wide variety of applications such as, LIBS applications; examination of high energy devices such as explosions or simulations of explosions; examination of deposition processes, e.g., coating formations; examination of chemical reactions; etc.

47 OTHER INSTRUMENTATION↗

The nature and energetics of AGN-driven perturbations in the hot gas in the Perseus Cluster

In this paper, cores of relaxed galaxy clusters are often disturbed by AGN. Their Chandra observations revealed a wealth of structures induced by shocks, subsonic gas motions, bubbles of relativistic plasma, etc. In this paper, we determine the nature and energy content of gas fluctuations in the Perseus core by probing statistical properties of emissivity fluctuations imprinted in the soft- and hard-band X-ray images. About 80 per cent of the total variance of perturbations on ~8–70 kpc scales in the core have an isobaric nature, i.e. are consistent with subsonic displacements of the gas in pressure equilibrium with the ambient medium. The observed variance translates to the ratio of energy in perturbations to thermal energy of ~13 per cent. In the region dominated by weak ‘ripples’, about half of the total variance is associated with isobaric perturbations on scales of a few tens of kpc. If these isobaric perturbations are induced by buoyantly rising bubbles, then these results suggest that most of the AGN-injected energy should first go into bubbles rather than into shocks. Using simulations of a shock propagating through the Perseus atmosphere, we found that models reproducing the observed features of a central shock have more than 50 per cent of the AGN-injected energy associated with the bubble enthalpy and only about 20 per cent is carried away with the shock. Such energy partition is consistent with the AGN-feedback model, mediated by bubbles of relativistic plasma, and supports the importance of turbulence in the cooling–heating balance.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Linear Estimation of Particle Bulk Parameters from Multi-Wavelength Lidar Measurements

An algorithm for linear estimation of aerosol bulk properties such as particle volume, effective radius and complex refractive index from multiwavelength lidar measurements is presented. The approach uses the fact that the total aerosol concentration can well be approximated as a linear combination of aerosol characteristics measured by multiwavelength lidar. Therefore, the aerosol concentration can be estimated from lidar measurements without the need to derive the size distribution, which entails more sophisticated procedures. The definition of the coefficients required for the linear estimates is based on an expansion of the particle size distribution in terms of the measurement kernels. Once the coefficients are established, the approach permits fast retrieval of aerosol bulk properties when compared with the full regularization technique. In addition, the straightforward estimation of bulk properties stabilizes the inversion making it more resistant to noise in the optical data. Numerical tests demonstrate that for data sets containing three aerosol backscattering and two extinction coefficients (so called 3 + 2 ) the uncertainties in the retrieval of particle volume and surface area are below 45% when input data random uncertainties are below 20 %. Moreover, using linear estimates allows reliable retrievals even when the number of input data is reduced. To evaluate the approach, the results obtained using this technique are compared with those based on the previously developed full inversion scheme that relies on the regularization procedure. Both techniques were applied to the data measured by multiwavelength lidar at NASA/GSFC. The results obtained with both methods using the same observations are in good agreement. At the same time, the high speed of the retrieval using linear estimates makes the method preferable for generating aerosol information from extended lidar observations. To demonstrate the efficiency of the method, an extended time series of observations acquired in Turkey in May 2010 was processed using the linear estimates technique permitting, for what we believe to be the first time, temporal-height distributions of particle parameters.

linear estimation↗

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 1 0A4 mergers/year). We develop a semi-analytical, phenomenological model of massive black hole 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 approx. 5 - 30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology.

Alicea-Munoz, E.↗

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 10(exp 4) mergers/year). We develop a semi-analytical, phenomenological model of massive black hole 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 approximately equal to 5-30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology. This work is supported in part by the Cooperative Education Program at NASA/GSFC.

Alicea-Munoz, E.↗

Uncertainty of Observation Impact Estimation in an Adjoint Model Investigated with an Observing System Simulation Experiment

Adjoint models are often used to estimate the impact of different observations on short-term forecast skill. A common difficulty with the evaluation of short term forecast quality is the choice of verification fields. The use of self-analysis fields for verification is typical but incestuous, and introduces uncertainty due to biases and errors in the analysis field. In this study, an observing system simulation experiment (OSSE) is used to explore the uncertainty in adjoint model estimations of observation impact. The availability of the true state for verification in the OSSE framework in the form of the Nature Run allows calculation of the observation impact without the uncertainties present in self-analysis verification. These impact estimates are compared to estimates calculated using self-analysis verification. The Global Earth Observing System version 5 (GEOS-5) forecast model with Gridpoint Statistical Interpolation (GSI) is used with the National Aeronautics and Space Administration Global Modeling and Assimilation Office (NASA/GMAO) OSSE capability. The adjoint model includes moist processes, with total wet energy selected as the norm for evaluation of observation impacts. The results show that there are measurable but small discrepancies in the adjoint model estimation of observation impact. In general, observations of temperature and winds tend to have overestimated impacts with self-analysis verification, while observations of humidity and moisture-affected observations tend to have underestimated impacts. The small magnitude of the differences in impact estimates supports the robustness of the adjoint method of estimating observation impacts.

N C Prive↗

Interrelationships among methods of estimating microbial biomass across multiple soil orders and biomes

Understanding the role of soil microbes is critical to ecosystem processes, and more thorough comparisons of measurement proxies for soil microbial biomass could broaden the inclusion of explicit microbial parameterization in soil carbon cycling and earth system models. We measured physical, chemical, and biological data from eight soil orders representing 11 major biomes and four climate regions. Four prominent methods to measure microbial abundance—chloroform fumigation extraction (CFE), total DNA yield, gene copy number by quantitative polymerase chain reaction (GCN), and phospholipid fatty acids (PLFA)—were compared to assess their relationships with each other and with soil characteristics. Correlations were observed when comparing methods, with CFE correlating strongly with total DNA yield, GCN, and PLFA; CFE with bacterial GCN and bacterial PLFA; and to a lesser extent, total PLFA and total DNA yield. Correlations improved with the removal of organic soils (Histosols, Gelisols). Comparisons involving extracted DNA were improved by correcting for clay content, due to DNA extraction inefficiencies in clay-rich soils. Correlations involving fungi (PLFA or GCN) were always less significant. These methods could serve as reliable, inter-relatable proxies for the estimation of total soil microbial biomass while recognizing that the proxies are less effective at parsing differences between bacteria and fungi. Here, we provide specific equations to relate measures of soil microbial biomass by these four different methods to enable microbial models to utilize a greater diversity of observed data sources in parameterizations and simulations. Caveats for the equations and their values are also discussed.

59 BASIC BIOLOGICAL SCIENCES↗

ROSAT X-Ray Observation of the Second Error Box for SGR 1900+14

The positions of the two error boxes for the soft gamma repeater (SGR) 1900+14 were determined by the "network synthesis" method, which employs observations by the Ulysses gamma-ray burst and CGRO BATSE instruments. The location of the first error box has been observed at optical, infrared, and X-ray wavelengths, resulting in the discovery of a ROSAT X-ray point source and a curious double infrared source. We have recently used the ROSAT HRI to observe the second error box to complete the counterpart search. A total of six X-ray sources were identified within the field of view. None of them falls within the network synthesis error box, and a 3 sigma upper limit to any X-ray counterpart was estimated to be 6.35 x 10(exp -14) ergs/sq cm/s. The closest source is approximately 3 min. away, and has an estimated unabsorbed flux of 1.5 x 10(exp -12) ergs/sq cm/s. Unlike the first error box, there is no supernova remnant near the second error box. The closest one, G43.9+1.6, lies approximately 2.dg6 away. For these reasons, we believe that the first error box is more likely to be the correct one.

Li, P.↗

Cosmic evolution of long gamma-ray burst luminosity

The cosmic evolution of gamma-ray burst (GRB) luminosity is essential for revealing the GRB physics and for using GRBs as cosmological probes. We investigate the luminosity evolution of long GRBs with a large sample of 258 Swift/BAT GRBs. By describing the peak luminosity evolution of individual GRBs as L{sub p}∝(1+z){sup k}, we get k=1.49±0.19 using the nonparametric τ statistics method without considering observational biases of GRB trigger and redshift measurement. By modeling these biases with the observed peak flux and characterizing the peak luminosity function of long GRBs as a smoothly broken power law with a break that evolves as L{sub b}∝(1+z){sup k{sub b}}, we obtain k{sub b}=1.14{sub −0.47}{sup +0.99} through simulations based on the assumption that the long GRB rate follows the star formation rate incorporating the cosmic metallicity history. The derived k and k{sub b} values are systematically smaller than those reported in previous papers. By removing the observational biases of the GRB trigger and redshift measurement based on our simulation analysis, we generate mock complete samples of 258 and 1000 GRBs to examine how these biases affect the τ statistics method. We get k=0.94±0.14 and k=0.80±0.09 for the two samples, indicating that these observational biases may lead to overestimating the k value. With the large uncertainty of k{sub b} derived from our simulation analysis, one cannot even convincingly argue for a robust evolution feature of the GRB luminosity.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

THERMAL DIAGNOSTICS WITH THE ATMOSPHERIC IMAGING ASSEMBLY ON BOARD THE SOLAR DYNAMICS OBSERVATORY: A VALIDATED METHOD FOR DIFFERENTIAL EMISSION MEASURE INVERSIONS

We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory. The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized Gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a nonlinear force-free field, and (3) thermodynamic models from a fully compressible, 3D MHD simulation of active region (AR) corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and Hinode X-ray Telescope data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of charge trapping in CZT detectors

A semiconductor detection material, Cadmium Zinc Telluride, was investigated for potential use as a particle detector for space flight applications. Our set of detectors was determined to be flawed for particle detection methods, likely due to charge transport mechanisms within the detector material, or due to reduced voltage required to avoid dielectric breakdown. It was observed that the method of crystal growth may play a role in the performance, but performance may also vary between manufacturers using the same growth method. It is recommended to perform further study at a larger proton beam facility.

47 OTHER INSTRUMENTATION↗

Assessing time-of-flight signal-to-noise ratio gains within the myocardium and subsequent reductions in administered activity in cardiac PET studies

Background: Time-of-flight (TOF) is known to increase signal-to-noise ratio (SNR) and facilitate reductions in administered activity. Established measures of SNR gain are derived from areas of uniform uptake, which is not applicable to the heterogeneous uptake in cardiac PET images using fluoro-deoxyglucose (FDG). This study aimed to develop a technique to quantify SNR gains within the myocardium due to TOF. Methods: Reference TOF SNR gains were measured in 88 FDG oncology patients. Phantom data were used to translate reference SNR gains and validate a method of quantifying SNR gains within the myocardium from parametric images produced from multiple replicate images. This technique was applied to 13 FDG cardiac viability patients. Results: Reference TOF SNR gains of +23% ± 8.5% were measured in oncology patients. Measurements of SNR gain from the phantom data were in agreement and showed the parametric image technique to be sufficiently robust. SNR gains within the myocardium in the viability patients were +21% ± 2.8%. Conclusion: A method to quantify SNR gains from TOF within the myocardium has been developed and evaluated. SNR gains within the myocardium are comparable to those observed by established methods. This allows guidance for protocol optimization for TOF systems in cardiac PET.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Spike-Free Adaptive Sliding Mode Control: Application to Permanent Magnet Synchronous Motors

A new methodology for adaptive sliding mode control (ASMC) has been widely used to improve the control performance in various systems. This method exhibits several advantages, including low sliding mode control (SMC) chattering, no knowledge of the system disturbance bound, and no overestimation of the control gain. Despite its advantages, this method can be hampered by the spike phenomenon, slow control gain convergence, and difficulty in achieving optimal performance under varying disturbances. Consequently, this article proposes a spike-free ASMC method with a disturbance observer (DOB) to address these problems. Previous ASMC methods have been analyzed via simulations to verify the aforementioned problems. Here, this analysis highlights the need for disturbance compensation and improvements in the SMC gain adaptation law. Therefore, a DOB is designed to mitigate the spike phenomenon by compensating for disturbances. Subsequently, an SMC gain adaptation law based on disturbance error estimation is designed to eliminate the spike phenomenon completely. The proposed adaptation law makes the SMC gain to converge to a slightly higher value than the disturbance estimation error. Consequently, the proposed method not only eliminates the spike phenomenon, but also ensures optimal performance under varying disturbances. The performance of the proposed method is experimen tally validated through a comparative study.

42 ENGINEERING↗

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗