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 217 records · Page 12

Exploring Discretization Error in Simulation-Based Aerodynamic Databases

This work examines the level of discretization error in simulation-based aerodynamic databases and introduces strategies for error control. Simulations are performed using a parallel, multi-level Euler solver on embedded-boundary Cartesian meshes. Discretization errors in user-selected outputs are estimated using the method of adjoint-weighted residuals and we use adaptive mesh refinement to reduce these errors to specified tolerances. Using this framework, we examine the behavior of discretization error throughout a token database computed for a NACA 0012 airfoil consisting of 120 cases. We compare the cost and accuracy of two approaches for aerodynamic database generation. In the first approach, mesh adaptation is used to compute all cases in the database to a prescribed level of accuracy. The second approach conducts all simulations using the same computational mesh without adaptation. We quantitatively assess the error landscape and computational costs in both databases. This investigation highlights sensitivities of the database under a variety of conditions. The presence of transonic shocks or the stiffness in the governing equations near the incompressible limit are shown to dramatically increase discretization error requiring additional mesh resolution to control. Results show that such pathologies lead to error levels that vary by over factor of 40 when using a fixed mesh throughout the database. Alternatively, controlling this sensitivity through mesh adaptation leads to mesh sizes which span two orders of magnitude. We propose strategies to minimize simulation cost in sensitive regions and discuss the role of error-estimation in database quality.

Aftosmis, Michael J.↗

Passive fetal heart rate monitoring apparatus and method with enhanced fetal heart beat discrimination

An apparatus for acquiring signals emitted by a fetus, identifying fetal heart beats and determining a fetal heart rate is presented. Multiple sensor signals are outputted by a passive fetal heart rate monitoring sensor. Multiple parallel nonlinear filters filter these multiple sensor signals to identify fetal heart beats in the signal data. A processor determines a fetal heart rate based on these identified fetal heart beats. The processor includes the use of a figure of merit weighting of heart rate estimates based on the identified heart beats from each filter for each signal. The fetal heart rate thus determined is outputted to a display, storage, or communications channel. A method for enhanced fetal heart beat discrimination includes acquiring signals from a fetus, identifying fetal heart beats from the signals by multiple parallel nonlinear filtering, and determining a fetal heart rate based on the identified fetal heart beats. A figure of merit operation in this method provides for weighting a plurality of fetal heart rate estimates based on the identified fetal heart beats and selecting the highest ranking fetal heart rate estimate.

Zahorian, Stephen A.↗

Passive fetal heart rate monitoring apparatus and method with enhanced fetal heart beat discrimination

An apparatus for acquiring signals emitted by a fetus, identifying fetal heart beats and determining a fetal heart rate. Multiple sensor signals are outputted by a passive fetal heart rate monitoring sensor. Multiple parallel nonlinear filters filter these multiple sensor signals to identify fetal heart beats in the signal data. A processor determines a fetal heart rate based on these identified fetal heart beats. The processor includes the use of a figure of merit weighting of heart rate estimates based on the identified heart beats from each filter for each signal. The fetal heart rate thus determined is outputted to a display, storage, or communications channel. A method for enhanced fetal heart beat discrimination includes acquiring signals from a fetus, identifying fetal heart beats from the signals by multiple parallel nonlinear filtering, and determining a fetal heart rate based on the identified fetal heart beats. A figure of merit operation in this method provides for weighting a plurality of fetal heart rate estimates based on the identified fetal heart beats and selecting the highest ranking fetal heart rate estimate.

Zahorian, Stephen A.↗

Toroidal transformer design program with application to inverter circuitry

Estimates of temperature, weight, efficiency, regulation, and final dimensions are included in the output of the computer program for the design of transformers for use in the basic parallel inverter. The program, written in FORTRAN 4, selects a tape wound toroidal magnetic core and, taking temperature, materials, core geometry, skin depth, and ohmic losses into account, chooses the appropriate wire sizes and number of turns for the center tapped primary and single secondary coils. Using the program, 2- and 4-kilovolt-ampere transformers are designed for frequencies from 200 to 3200 Hz and the efficiency of a basic transistor inverter is estimated.

Dayton, J. A., Jr.↗

Massively Parallel Assimilation of TOGA/TAO and Topex/Poseidon Measurements into a Quasi Isopycnal Ocean General Circulation Model Using an Ensemble Kalman Filter

A massively parallel ensemble Kalman filter (EnKF)is used to assimilate temperature data from the TOGA/TAO array and altimetry from TOPEX/POSEIDON into a Pacific basin version of the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. The EnKF is an approximate Kalman filter in which the error-covariance propagation step is modeled by the integration of multiple instances of a numerical model. An estimate of the true error covariances is then inferred from the distribution of the ensemble of model state vectors. This inplementation of the filter takes advantage of the inherent parallelism in the EnKF algorithm by running all the model instances concurrently. The Kalman filter update step also occurs in parallel by having each processor process the observations that occur in the region of physical space for which it is responsible. The massively parallel data assimilation system is validated by withholding some of the data and then quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The distributions of the forecast and analysis error covariances predicted by the ENKF are also examined.

Keppenne, Christian L.↗

High density constraint on the entropy instability

The entropy instability squared is a nonisothermal effect which is eliminated by parallel ion pressure at high densities (k sub z lambda sub e 1/2 sq root of m/M), reducing previous growth rate estimates and the range of unstable parameters.

Hudson, M. K.↗

Multispacecraft observations of microscale fluctuations in the solar wind

Data from the Explorer 33 and 35 satellites were used to study Alfvenic fluctuations (i.e., fluctuations with a high correlation between perturbation in velocity and magnetic field) in the solar wind. It was found that the average minimum variance method does not give a good estimate of the phase front direction of these fluctuations. There is a distinct class of Alfvenic fluctuations with phase planes nearly parallel to the mean field B sub 0. The characteristic length of the phase fronts of the fluctuations is probably significantly greater than 50 R-E. It is concluded that since tangential discontinuities and possibly other types of static structures appear in Alfvenic fluctuations, these fluctuations are not solely Alfven waves.

Denskat, K. U.↗

Structure and Growth of the Marine Boundary Layer

LANDSAT visible imagery and a one-dimensional Lagrangian boundary layer model were used to hypothesize the nature and the development of the marine boundary layer during a winter episode of strong seaward cold air advection. Over-water heating and moistening of the cold, dry continental air is estimable from linear relations involving horizontal gradients of the near-surface air temperature and humidity. A line of enhanced convection paralleling the Atlantic U.S. coast from south of New York Bay to the vicinity of Virginia Beach, VA was attributed to stronger convergence at low levels. This feature was characterized as a mesoscale front. With the assistance of a three-dimensional mesoscale boundary layer model, initialized with data obtained from the MASEX, the marine boundary layer can be mapped over the entire Atlantic coastal domain and the evolution of the boundary layer can be studied as a function of different characteristics of important surface level forcings. The effects on boundary layer growth due to the magnitude and pattern of sea surface temperature, to the shape of the coastline, and to atmospheric conditions, such as the orientation of the prevailing wind are examined.

Mccumber, M.↗

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.↗