Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel estimation”

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 199 records · Page 11

Experimental evaluation of damping of graphite fiber composites

An investigation is described that was conducted to measure the damping of graphite fiber composites such as are to be used on the Hubble Space Telescope. The measurements were made on a composite cylinder under normal atmospheric conditions and in a vacuum chamber with varying thermal conditions. Such composites can outgas under vacuum conditions and this emission of trapped gases may have some influence on the material damping. Some results are presented for the damping coefficient measured under normal atmospheric conditions for a graphite fiber composite. The experimental damping data thus obtained were compared with the least square estimate obtained through a statistical analysis. The statistical analysis was based on a linear mathematical model consisting of a spring in parallel with a dashpot. It is observed that the discrepancy in the damping ratio between logarithmic decrement result and the least square fit is less than 6 percent for all the different tests made on the composite cylinder. Signal-to-noise ratio was maintained to be greater than 10.

Rao, M. D.↗

The growth of Goertler vortices in compressible boundary layers

The linear instability of Goertler vortices in compressible boundary layers is considered. Using asymptotic methods in the high wavenumber regime, it is shown that a growth rate estimate can be found by solving a sequence of linear equations. The growth rate obtained in this way takes non-parallel effects into account and can be found much more easily than by ordinary differential equation eigenvalue calculations associated with parallel flow theories.

Hall, Philip↗

Absorption of isotropic radiation

Estimates of the absorption of isotropic radiation are given for an absorbing medium whose geometrical shape is either a plane-parallel sided slab or a cylinder. By including the directional variation of absorption path length, the results show appreciable differences to the 'unidirectional radiation' formulas.

Harries, Wynford L.↗

Progress in the Simulation of Steady and Time-Dependent Flows with 3D Parallel Unstructured Cartesian Methods

The proposed paper will present recent extensions in the development of an efficient Euler solver for adaptively-refined Cartesian meshes with embedded boundaries. The paper will focus on extensions of the basic method to include solution adaptation, time-dependent flow simulation, and arbitrary rigid domain motion. The parallel multilevel method makes use of on-the-fly parallel domain decomposition to achieve extremely good scalability on large numbers of processors, and is coupled with an automatic coarse mesh generation algorithm for efficient processing by a multigrid smoother. Numerical results are presented demonstrating parallel speed-ups of up to 435 on 512 processors. Solution-based adaptation may be keyed off truncation error estimates using tau-extrapolation or a variety of feature detection based refinement parameters. The multigrid method is extended to for time-dependent flows through the use of a dual-time approach. The extension to rigid domain motion uses an Arbitrary Lagrangian-Eulerlarian (ALE) formulation, and results will be presented for a variety of two- and three-dimensional example problems with both simple and complex geometry.

Aftosmis, M. J.↗

Studying Functions of All Yeast Genes Simultaneously

A method of studying the functions of all the genes of a given species of microorganism simultaneously has been developed in experiments on Saccharomyces cerevisiae (commonly known as baker's or brewer's yeast). It is already known that many yeast genes perform functions similar to those of corresponding human genes; therefore, by facilitating understanding of yeast genes, the method may ultimately also contribute to the knowledge needed to treat some diseases in humans. Because of the complexity of the method and the highly specialized nature of the underlying knowledge, it is possible to give only a brief and sketchy summary here. The method involves the use of unique synthetic deoxyribonucleic acid (DNA) sequences that are denoted as DNA bar codes because of their utility as molecular labels. The method also involves the disruption of gene functions through deletion of genes. Saccharomyces cerevisiae is a particularly powerful experimental system in that multiple deletion strains easily can be pooled for parallel growth assays. Individual deletion strains recently have been created for 5,918 open reading frames, representing nearly all of the estimated 6,000 genetic loci of Saccharomyces cerevisiae. Tagging of each deletion strain with one or two unique 20-nucleotide sequences enables identification of genes affected by specific growth conditions, without prior knowledge of gene functions. Hybridization of bar-code DNA to oligonucleotide arrays can be used to measure the growth rate of each strain over several cell-division generations. The growth rate thus measured serves as an index of the fitness of the strain.

Stolc, Viktor↗

Near Real-Time Probabilistic Damage Diagnosis Using Surrogate Modeling and High Performance Computing

This work investigates novel approaches to probabilistic damage diagnosis that utilize surrogate modeling and high performance computing (HPC) to achieve substantial computational speedup. Motivated by Digital Twin, a structural health management (SHM) paradigm that integrates vehicle-specific characteristics with continual in-situ damage diagnosis and prognosis, the methods studied herein yield near real-time damage assessments that could enable monitoring of a vehicle's health while it is operating (i.e. online SHM). High-fidelity modeling and uncertainty quantification (UQ), both critical to Digital Twin, are incorporated using finite element method simulations and Bayesian inference, respectively. The crux of the proposed Bayesian diagnosis methods, however, is the reformulation of the numerical sampling algorithms (e.g. Markov chain Monte Carlo) used to generate the resulting probabilistic damage estimates. To this end, three distinct methods are demonstrated for rapid sampling that utilize surrogate modeling and exploit various degrees of parallelism for leveraging HPC. The accuracy and computational efficiency of the methods are compared on the problem of strain-based crack identification in thin plates. While each approach has inherent problem-specific strengths and weaknesses, all approaches are shown to provide accurate probabilistic damage diagnoses and several orders of magnitude computational speedup relative to a baseline Bayesian diagnosis implementation.

Warner, James E.↗

Determination of relative site motions in the western United States using Mark III very long baseline interferometry

Four years of very long baseline interferometry (VLBI) measurements involving fixed and mobile stations in the western U.S. are analyzed using weighted least squares. Baseline solution time series are used to estimate rates of change of the vector baseline between pairs of VLBI stations. Horizontal vector motions are determined for 18 sites with respect to a site in the Mojave desert. Six sites west of and two sites just to the east of the San Andreas fault show northwestward motion closely parallel to the San Andreas fault. Their motions range from 18 to 43 mm/yr and are significant at confidence levels greater than 99 percent. Motions in California relative to interior North American plate sites are used to estimate an integrated contribution from Basin and Range spreading to the total relative Pacific plate motion as 9-10 mm/yr. A maximum relative velocity of 50-51 mm/yr is found between the interior North American sites and the western California sites, in close agreement with the NUVEL-1 plate motion model.

Clark, T. A.↗

Unveiling the Galaxy Population at 1.3 < z < 4: the HUDF05 NICMOS Parallel Fields

Using the Hubble Ultra Deep Field Near Infrared Camera and Multi-Object Spectrometer (HUDF-NICMOS) UDF05 parallel fields, we cross-matched 301 out of 630 galaxies with the ACS filters V606 and z850, NICMOS filters J110 and H160, and Spitzer IRAC filters at 3.6, 4.5, 5.8 , and 8.0 (mu)m. We modeled the spectral energy distributions (SEDs) to estimate: photometric redshifts, dust extinction, stellar mass, bolometric luminosity, starburst age and metallicity. To validate the photometric redshifts, comparisons with 16 spectroscopic redshifts give 75% within Delta < 0.2, which agrees with the sensitivities expected from the Balmer-break in our dataset. Five parallel fields observed by NICMOS have sensitivities in the H160-band of 80% at mAB = 25.4 and 50% at mAB = 26.7. Because the sample is H160-band selected, it is sensitive to stellar mass rather than UV luminosities. We also use Monte Carlo simulations to determine that the parameters from the best-fit SEDs are robust for the redshift ranges z > or approx. 1.3. Based on the robustness of the photometric redshifts, we analyze a subsample of the 301 galaxies at 1.3 < or = z < or = 2 (35 objects) and 3 < or = z < or = 4 (31 objects) and determine that L(BoI) and the star formation rate increase significantly from z approx. 1.5 to 4. The Balmer decrement is indicative of more evolved galaxies, and at high redshifts, they serve as records of some of the first galaxies. Therefore, the galaxies in this sample are great candidates for future surveys with the James Webb Space Telescope and Atacama Large Millimeter Array.

Petty, Sara M.↗

A Posteriori Error Estimation for Finite Volume and Finite Element Approximations Using Broken Space Approximation

We consider a posteriori error estimates for finite volume and finite element methods on arbitrary meshes subject to prescribed error functionals. Error estimates of this type are useful in a number of computational settings: (1) quantitative prediction of the numerical solution error, (2) adaptive meshing, and (3) load balancing of work on parallel computing architectures. Our analysis recasts the class of Godunov finite volumes schemes as a particular form of discontinuous Galerkin method utilizing broken space approximation obtained via reconstruction of cell-averaged data. In this general framework, weighted residual error bounds are readily obtained using duality arguments and Galerkin orthogonality. Additional consideration is given to issues such as nonlinearity, efficiency, and the relationship to other existing methods. Numerical examples are given throughout the talk to demonstrate the sharpness of the estimates and efficiency of the techniques. Additional information is contained in the original.

Barth, Timothy J.↗

Satellite Estimation of Spectral Surface UV Irradiance: Effect of Horizontally Homogeneous Clouds - 2

The local variability of UV irradiance at the Earth's surface is mostly caused by clouds in addition to the seasonal variability. Parametric representations of radiative transfer RT calculations are presented for the convenient solution of the transmission T of ultraviolet radiation through plane parallel clouds over a surface with reflectivity R(sub s). The calculations are intended for use with the Total Ozone Mapping Spectrometer (TOMS) measured radiances to obtain the calculated Lambert equivalent scene reflectivity R for scenes with and without clouds. The purpose is to extend the theoretical analysis of the estimation of UV irradiance from satellite data for a cloudy atmosphere. Results are presented for a range of cloud optical depths and solar zenith angles for the cases of clouds over a low reflectivity surface R(sub s) less than 0.1, over a snow or ice surface R(sub s) greater than 0.3, and for transmission through a non-conservative scattering cloud with single scattering albedo omega(sub 0) = 0.999. The key finding for conservative scattering is that the cloud-transmission function C(sub T), the ratio of cloudy-to clear-sky transmission, is roughly C(sub T) = 1 - R(sub c) with an error of less than 20% for nearly overhead sun and snow-free surfaces. For TOMS estimates of UV irradiance in the presence of both snow and clouds, independent information about snow albedo is needed for conservative cloud scattering. For non-conservative scattering with R(sub s) greater than 0.5 (snow) the satellite measured scene reflectance cannot be used to estimate surface irradiance. The cloud transmission function has been applied to the calculation of UV irradiance at the Earth's surface and compared with ground-based measurements.

Krothov, N.↗

Estimation of stiffening effect of shaft and housing material outside projected area of a rolling element bearing

In the analysis of distortions occurring in rolling-element bearings, it is common to neglect the stiffening effect of shafting outside the bearing region. The magnitude of such an effect will be dependent primarily on the bearing width-to-bore ratio, the shaft geometry, and the location of the bearing on the shaft. An estimate is given of the stiffening effect for a wide range of these variables. In addition, brief consideration is given to the parallel situation existing at the outer ring housing.

Taylor, C. M.↗

Modal identification using single-mode projection filters and comparison with ERA and MLE results

The Single-Mode Projection Filter (SPF) is a newly developed algorithm for eigensystem parameter identification from both analytical results and test data. The SPF is formulated with a single mode only and practical for parallel processing implementation. Explicit formulations of SPF are derived for the multi-input multi-output (MIMO) system by using the orthogonal matrices of the controllability and observability matrices in the general sense. The modal parameters of SPF are initially obtained from an analytical model in modal space. The experimental data are then processed through SPF to update its modal parameters and to minimize a cost function defined by the norm of an error matrix. The updated modal parameters represent the characteristics of the test data. A two-dimensional global minimum optimization algorithm is developed and applied for the filter update by using the interval analysis method. The SPF is developed based on a single-mode subsystem and identifies only one modal frequency and one modal damping within a specified region. For an n-modes structure, n SPF can be implemented for parallel processing to reduce the computational burden. The SPF is applied to analyze the simulated data for the MAST beam structure. The estimated modal parameters are comparable to those from the Eigensystem Realization Algorithm (ERA) and repeated modal frequencies are identified. The modal analysis of the Spacecraft Control Laboratory Experiment (SCOLE) data is also performed by using the ERA and the Maximum Likelihood Estimate (MLE). The result shows that the first five modal frequencies are very close from ERA and MLE. However, there are slight disparities in the damping rates and the computational burdens are quite different among these two algorithms.

Huang, Jen-Kuang↗

Using parallel banded linear system solvers in generalized eigenvalue problems

Subspace iteration is a reliable and cost effective method for solving positive definite banded symmetric generalized eigenproblems, especially in the case of large scale problems. This paper discusses an algorithm that makes use of two parallel banded solvers in subspace iteration. A shift is introduced to decompose the banded linear systems into relatively independent subsystems and to accelerate the iterations. With this shift, an eigenproblem is mapped efficiently into the memories of a multiprocessor and a high speed-up is obtained for parallel implementations. An optimal shift is a shift that balances total computation and communication costs. Under certain conditions, we show how to estimate an optimal shift analytically using the decay rate for the inverse of a banded matrix, and how to improve this estimate. Computational results on iPSC/2 and iPSC/860 multiprocessors are presented.

Zhang, Hong↗

A comparison of energetic ions in the plasma depletion layer and the quasi-parallel magnetosheath

Energetic ion spectra measured by the Active Magnetospheric Particle Tracer Explorers/Charge Composition Explorer (AMPTE/CCE) downstream from the Earth's quasi-parallel bow shock (in the quasi-parallel magnetosheath) and in the plasma depletion layer are compared. In the latter region, energetic ions are from a single source, leakage of magnetospheric ions across the magnetopause and into the plasma depletion layer. In the former region, both the magnetospheric source and shock acceleration of the thermal solar wind population at the quasi-parallel shock can contribute to the energetic ion spectra. The relative strengths of these two energetic ion sources are determined through the comparison of spectra from the two regions. It is found that magnetospheric leakage can provide an upper limit of 35% of the total energetic H(+) population in the quasi-parallel magnetosheath near the magnetopause in the energy range from approximately 10 to approximately 80 keV/e and substantially less than this limit for the energetic He(2+) population. The rest of the energetic H(+) population and nearly all of the energetic He(2+) population are accelerated out of the thermal solar wind population through shock acceleration processes. By comparing the energetic and thermal He(2+) and H(+) populations in the quasi-parallel magnetosheath, it is found that the quasi-parallel bow shock is 2 to 3 times more efficient at accelerating He(2+) than H(+). This result is consistent with previous estimates from shock acceleration theory and simulati ons.

Fuselier, Stephen A.↗

Using parallel banded linear system solvers in generalized eigenvalue problems

Subspace iteration is a reliable and cost effective method for solving positive definite banded symmetric generalized eigenproblems, especially in the case of large scale problems. This paper discusses an algorithm that makes use of two parallel banded solvers in subspace iteration. A shift is introduced to decompose the banded linear systems into relatively independent subsystems and to accelerate the iterations. With this shift, an eigenproblem is mapped efficiently into the memories of a multiprocessor and a high speedup is obtained for parallel implementations. An optimal shift is a shift that balances total computation and communication costs. Under certain conditions, we show how to estimate an optimal shift analytically using the decay rate for the inverse of a banded matrix, and how to improve this estimate. Computational results on iPSC/2 and iPSC/860 multiprocessors are presented.

DISTRIBUTED MEMORY MULTIPROCES↗

Constraints on the Evolution of the Dichotomy Boundary at 50-90E

The global dichotomy is a fundamental feature of Mars. It marks the boundary between the highly cratered, older southern highlands, and the northern plains. Recent analysis of buried craters in the northern lowlands confirms the long held suspicion that they are comparable in age to the southern highlands, but with surficial deposits of younger material. A variety of exogenic and endogenic models have been proposed for the origin of the dichotomy, including multiple impacts, plate tectonics, and degree one convection produced by core formation, a plume under the lowlands, or a plume under the highlands. New gravity and topography data from the Mars Global Surveyor (MGS) Mission favor endogenic processes. In this study we examine MGS topography, gravity and magnetic field data to constrain the tectonic history the dichotomy in the region 30-60N and 50-90E, which encompasses portions of the Ismenius Lacus quadrangle. The dichotomy formed very early the history of Mars and has undergone extensive modification by impact cratering, erosion, and faulting. This history must be carefully interpreted in order to reconstruct the original nature of the dichotomy boundary and ultimately discriminate between models of origin. In the study area boundary-parallel faults are well preserved, and may be the result of gravitational relaxation. The geologic history has been examined in detail, including estimates of volumes of material eroded. Further, it is one of the few regions where there is a correlation between the free air gravity, magnetic anomalies, and the geology. This allows to constrain subsurface faulting beneath the lowlands fill material. In addition to being an excellent location to unravel the complex history of the dichotomy, this area preserves the transition from a highly magnetized highlands crust to an unmagnetized or slightly magnetized lowlands crust.

S E Smrekar↗

Efficient GPU-Accelerated MultiSource Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational-wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analyzing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains massive black hole binaries, compact galactic binaries, and a parametrized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble Markov Chain Monte Carlo (MCMC) sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or transdimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data.We recover posterior distributions for all 15 (6) of the injected massive black hole binaries (MBHB) in the LDC2A training (hidden) dataset. We catalog ∼12000 galactic binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA↗

Parallel Eigenvalue extraction

A new numerical algorithm for the solution of large-order eigenproblems typically encountered in linear elastic finite element systems is presented. The architecture of parallel processing is utilized in the algorithm to achieve increased speed and efficiency of calculations. The algorithm is based on the frontal technique for the solution of linear simultaneous equations and the modified subspace eigenanalysis method for the solution of the eigenproblem. Assembly, elimination and back-substitution of degrees of freedom are performed concurrently, using a number of fronts. All fronts converge to and diverge from a predefined global front during elimination and back-substitution, respectively. In the meantime, reduction of the stiffness and mass matrices required by the modified subspace method can be completed during the convergence/divergence cycle and an estimate of the required eigenpairs obtained. Successive cycles of convergence and divergence are repeated until the desired accuracy of calculations is achieved. The advantages of this new algorithm in parallel computer architecture are discussed.

Akl, Fred A.↗