Search NASA⌕ Search

SEARCH · Search NASA

Results for “dependable computing”

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 235 records · Page 13

Efficient matrix partitioning for optical computing

Techniques for partitioning optical linear algebra problems to make them amenable to solution using optical processors programmed with simple algorithms are explored. Generalized methods for splitting a linear algebra matrix into a series of submatrices are reviewed, showing that simple forms can be pipelined smoothly and that parallel accumulation can be achieved by beam combining on detectors or by summing electronically. The techniques offer simplified bookkeeping, algorithmic independence, and high efficiency. The computational speed will depend on the number of multiplier-accumulators devoted to the task.

Caulfield, H. J.↗

Simulation Of Unsteady, Viscous, Incompressible Flow

Method for numerical solution of Navier-Stokes equations of viscous, incompressible flow developed based on use of fractional-step procedure. Accurate to second order in both space and time. Attempt made to minimize Poisson-equation difficulties by choosing pressures at centers and volume fluxes across faces of computational cells as dependent variables instead of familiar Cartesian components of velocity. Choice ensures satisfaction of discrete equation of conservation of mass to within round-off errors in any coordinate system and has favorable effects on convergence properties.

Rosenfeld, Moshe↗

Mass transfer from a sphere in an oscillating flow with zero mean velocity

A pseudospectral numerical method is used for the solution of the Navier-Stokes and mass transport equations for a sphere in a sinusoidally oscillating flow with zero mean velocity. The flow is assumed laminar and axisymmetric about the sphere's polar axis. Oscillating flow results were obtained for Reynolds numbers (based on the free-stream oscillatory flow amplitude) between 1 and 150, and Strouhal numbers between 1 and 1000. Sherwood numbers were computed and their dependency on the flow frequency and amplitude discussed. An assessment of the validity of the quasi-steady assumption for mass transfer is based on these results.

Drummond, Colin K.↗

A novel gas-droplet numerical method for spray combustion

This paper presents a non-iterative numerical technique for computing time-dependent gas-droplet flows. The method is a fully-interacting combination of Eulerian fluid and Lagrangian particle calculation. The interaction calculations between the two phases are formulated on a pressure-velocity coupling procedure based on the operator-splitting technique. This procedure eliminates the global iterations required in the conventional particle-source-in-cell (PSIC) procedure. Turbulent dispersion calculations are treated by a stochastic procedure. Numerical calculations and comparisons with available experimental data, as well as efficiency assessments are given for some sprays typical of spray combustion applications.

Chen, C. P.↗

System reliability approaches for advanced propulsion system structures

This paper identifies significant issues that pertain to the estimation and use of system reliability in the design of advanced propulsion system structures. Linkages between the reliabilities of individual components and their effect on system design issues such as performance, cost, availability, and certification are examined. The need for system reliability computation to address the continuum nature of propulsion system structures and synergistic progressive damage modes has been highlighted. Available system reliability models are observed to apply only to discrete systems. Therefore a sequential structural reanalysis procedure is formulated to rigorously compute the conditional dependencies between various failure modes. The method is developed in a manner that supports both top-down and bottom-up analyses in system reliability.

Cruse, T. A.↗

Computational study of duct and pipe flows using the method of pseudocompressibility

A viscous, three-dimensional, incompressible, Navier-Stokes Computational Fluid Dynamics code employing pseudocompressibility is used for the prediction of laminar primary and secondary flows in two 90-degree bends of constant cross section. Under study are a square cross section duct bend with 2.3 radius ratio and a round cross section pipe bend with 2.8 radius ratio. Sensitivity of predicted primary and secondary flow to inlet boundary conditions, grid resolution, and code convergence is investigated. Contour and velocity versus spanwise coordinate plots comparing prediction to experimental data flow components are shown at several streamwise stations before, within, and after the duct and pipe bends. Discussion includes secondary flow physics, computational method, computational requirements, grid dependence, and convergence rates.

Williams, Robert W.↗

Investigation of advancing front method for generating unstructured grid

The advancing front technique is used to generate an unstructured grid about simple aerodynamic geometries. Unstructured grids are generated using VGRID2D and VGRID3D software. Specific problems considered are a NACA 0012 airfoil, a bi-plane consisting of two NACA 0012 airfoil, a four element airfoil in its landing configuration, and an ONERA M6 wing. Inviscid time dependent solutions are computed on these geometries using USM3D and the results are compared with standard test results obtained by other investigators. A grid convergence study is conducted for the NACA 0012 airfoil and compared with a structured grid. A structured grid is generated using GRIDGEN software and inviscid solutions computed using CFL3D flow solver. The results obtained by unstructured grid for NACA 0012 airfoil showed an asymmetric distribution of flow quantities, and a fine distribution of grid was required to remove this asymmetry. On the other hand, the structured grid predicted a very symmetric distribution, but when the total number of points were compared to obtain the same results it was seen that structured grid required more grid points.

Thomas, A. M.↗

An efficient pressure-velocity procedure for gas-droplet two-phase flow calculations

This paper presents a noniterative numerical technique for computing time-dependent gas-droplet flows. The method is a fully interacting combination of Eulerian fluid and Lagrangian particle calculations. The interaction calculations between the two phases are formulated on a pressure-velocity-coupling procedure based on the operator-splitting technique. This procedure eliminates the global iterations required in the conventional particle-source-in-cell procedure. Turbulent dispersion calculations are treated by a stochastic procedure. Numerical calculations and comparisons with available experimental data as well as efficiency assessments are given for some sprays typical of spray combustion applications.

Chen, C. P.↗

Verification of a two-dimensional infiltration model for the resin transfer molding process

A two-dimensional finite element model for the infiltration of a dry textile preform by an injected resin was verified. The model, which is based on the finite element/control volume technique, determines the total infiltration time and the pressure increase at the mold inlet associated with the RTM process. Important input data for the model are the compaction and permeability behavior of the preform along with the kinetic and rheological behavior of the resin. The compaction behavior for several textile preforms was determined by experimental methods. A power law regression model was used to relate fiber volume fraction to the applied compaction pressure. Results showed a large increase in fiber volume fraction with the initial application of pressure. However, as the maximum fiber volume fraction was approached, the amount of compaction pressure required to decrease the porosity of the preform rapidly increased. Similarly, a power law regression model was used to relate permeability to the fiber volume fraction of the preform. Two methods were used to measure the permeability of the textile preform. The first, known as the steady state method, measures the permeability of a saturated preform under constant flow rate conditions. The second, denoted the advancing front method, determines the permeability of a dry preform to an infiltrating fluid. Water, corn oil, and an epoxy resin, Epon 815, were used to determine the effect of fluid type and viscosity on the steady state permeability behavior of the preform. Permeability values measured with the different fluids showed that fluid viscosity had no influence on the permeability behavior of 162 E-glass and TTI IM7/8HS preforms. Permeabilities measured from steady state and advancing front experiments for the warp direction of 162 E-glass fabric were similar. This behavior was noticed for tests conducted with corn oil and Epon 815. Comparable behavior was observed for the warp direction of the TTI IM7/8HS preform and corn oil. Mold filling and flow visualization experiments were performed to verify the analytical computer model. Frequency dependent electromagnetic sensors were used to monitor the resin flow front as a function of time. For the flow visualization tests, a video camera and high resolution tape recorder were used to record the experimental flow fronts. Comparisons between experimental and model predicted flow fronts agreed well for all tests. For the mold filling tests conducted at constant flow rate injection, the model was able to accurately predict the pressure increase at the mold inlet during the infiltration process. A kinetics model developed to predict the degree of cure as a function of time for the injected resin accurately calculated the increase in the degree of cure during the subsequent cure cycle.

Hammond, Vincent H.↗

Cirrus microphysics observations made during FIRE 2: Small particles, high concentrations, and probe comparisons

Aircraft observations of cirrus cloud microphysics were made near Coffeyville, Kansas during Nov. and Dec. 1991. Cloud microphysics measurements were made using both a PMS 2DC probe and an ice particle replicator, both were mounted on the UND Citation aircraft. Intercomparison is made of the size, area, and ice mass spectra determined from these probes. The PMS 2DC undercounts particles for D less than 70 microns and the replicator oversizes particles for D greater than 150 microns, at least when column rosettes are encountered. High concentrations of particles with D less than 50 microns are noted in selected portions of the 22 Nov. 91 replicator data set. Relations between the maximum dimension of a crystal and its shadow area (known as area dimensional relationships) are computed from the PMS data. Area and mass dimensional relationships are used to give a simple analytical expression for computing the wavelength dependent absorption coefficient averaged over a size bin. Calculations based upon the replicator data show that crystals with D less than 50 microns contribute significantly to the solar extinction and infrared absorption coefficients during some time intervals.

Arnott, W. P.↗

An analysis of whistler waves at interplanetary shocks

We present an analysis of whistler wave magnetic and electric field amplitude ratios from which we compute wave propagation angles and energies of electrons in resonance with the waves. To do this analysis, we compute the theoretical dependence of ratios of wave components on the whistler wave propagation angle Theta for various combinations of orthogonal wave components. Ratios of wave components that would be observed by a spinning spacecraft are determined, and the effects of arbitrary inclinations of the spacecraft to the ambient magnetic field and to the whistler wave vector are studied. This analysis clearly demonstrates that B/E, the ratio of magnetic to electric field amplitudes, cannot be assumed to be the wave index of refraction, contrary to assumptions of some earlier studies. Therefore previous interpretations of whistler wave observations based on this assumption must be reinvestigated. B/E ratios derived using three orthogonal wave components can be used to unambiguously determine Theta. Using spin plane observations alone, a significant uncertainty occurs in the determination of Theta. Nevertheless, for whistler waves observed downstream of several interplanetary shocks by the Ulysses plasma wave experiment we find that Theta is highly oblique. We suggest that the analysis of wave amplitude ratios used in conjunction with traditional stability analyses provide a promising tool for determining which particle distributions and resonances are likely to be dominant contributors to wave growth.

Lengyel-Frey, D.↗

EUV/soft x-ray spectra for low B neutron stars

Recent ROSAT and EUVE detections of spin-powered neutron stars suggest that many emit 'thermal' radiation, peaking in the EUV/soft X-ray band. These data constrain the neutron stars' thermal history, but interpretation requires comparison with model atmosphere computations, since emergent spectra depend strongly on the surface composition and magnetic field. As recent opacity computations show substantial change to absorption cross sections at neutron star photospheric conditions, we report here on new model atmosphere computations employing such data. The results are compared with magnetic atmosphere models and applied to PSR J0437-4715, a low field neutron star.

Romani, Roger W.↗

An Evaluation of Attitude-Independent Magnetometer-Bias Determination Methods

Although several algorithms now exist for determining three-axis magnetometer (TAM) biases without the use of attitude data, there are few studies on the effectiveness of these methods, especially in comparison with attitude dependent methods. This paper presents the results of a comparison of three attitude independent methods and an attitude dependent method for computing TAM biases. The comparisons are based on in-flight data from the Extreme Ultraviolet Explorer (EUVE), the Upper Atmosphere Research Satellite (UARS), and the Compton Gamma Ray Observatory (GRO). The effectiveness of an algorithm is measured by the accuracy of attitudes computed using biases determined with that algorithm. The attitude accuracies are determined by comparison with known, extremely accurate, star-tracker-based attitudes. In addition, the effect of knowledge of calibration parameters other than the biases on the effectiveness of all bias determination methods is examined.

Hashmall, J. A.↗

Global Asymptotic Behavior of Iterative Implicit Schemes

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing three models of 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed using the theory of dynamical systems. The iterative procedures include simple iteration and full and modified Newton iterations. The results are compared with standard Runge-Kutta explicit methods, a noniterative implicit procedure, and the Newton method of solving the steady part of the ODEs. Studies showed that aside from exhibiting spurious asymptotes, all of the four implicit LMMs can change the type and stability of the steady states of the differential equations (DEs). They also exhibit a drastic distortion but less shrinkage of the basin of attraction of the true solution than standard nonLMM explicit methods. The simple iteration procedure exhibits behavior which is similar to standard nonLMM explicit methods except that spurious steady-state numerical solutions cannot occur. The numerical basins of attraction of the noniterative implicit procedure mimic more closely the basins of attraction of the DEs and are more efficient than the three iterative implicit procedures for the four implicit LMMs. Contrary to popular belief, the initial data using the Newton method of solving the steady part of the DEs may not have to be close to the exact steady state for convergence. These results can be used as an explanation for possible causes and cures of slow convergence and nonconvergence of steady-state numerical solutions when using an implicit LMM time-dependent approach in computational fluid dynamics.

Yee, H. C.↗

MHD Streamer Structure, Slow Solar Wind, and the Streamer Brightness Boundary

Flow tubes adjacent to closed magnetic field lines on the boundaries of streamers can have extremely large geometric spreading factors. Numerical models in this thin layer are subject to grid definition uncertainties. Therefore, we compute flow tube geometry using the analytic model of streamer structure described by Pneuman. This model has been found to be more widely applicable than commonly believed as a consequence of observations made with SOHO/UVCS and YOKHOH/SXT. We use the model to compute the radial dependence of flow tube geometry (the "spreading factors") for several different streamer models. The results are used to analyze the hypothesis that extremely slow flows in these open flow tubes may cause high densities relative to adjacent coronal hole flow. Such high density could mean that the streamer brightness boundary is defined by the open flow tubes adjacent to streamers rather than closed field lines.

Suess, Steven T.↗

An Improved Retrieval of Tropospheric Nitrogen Dioxide from GOME

We present a retrieval of tropospheric nitrogen dioxide (NO2) columns from the Global Ozone Monitoring Experiment (GOME) satellite instrument that improves in several ways over previous retrievals, especially in the accounting of Rayleigh and cloud scattering. Slant columns, which are directly fitted without low-pass filtering or spectral smoothing, are corrected for an artificial offset likely induced by spectral structure on the diffuser plate of the GOME instrument. The stratospheric column is determined from NO2 columns over the remote Pacific Ocean to minimize contamination from tropospheric NO2. The air mass factor (AMF) used to convert slant columns to vertical columns is calculated from the integral of the relative vertical NO2 distribution from a global 3-D model of tropospheric chemistry driven by assimilated meteorological data (Global Earth Observing System (GEOS)-CHEM), weighted by altitude dependent scattering weights computed with a radiative transfer model (Linearized Discrete Ordinate Radiative Transfer), using local surface albedos determined from GOME observations at NO2 wavelengths. The AMF calculation accounts for cloud scattering using cloud fraction, cloud top pressure, and cloud optical thickness from a cloud retrieval algorithm (GOME Cloud Retrieval Algorithm). Over continental regions with high surface emissions, clouds decrease the AMT by 20- 30% relative to clear sky. GOME is almost twice as sensitive to tropospheric NO2 columns over ocean than over land. Comparison of the retrieved tropospheric NO2 columns for July 1996 with GEOS-CHEM values tests both the retrieval and the nitrogen oxide radical (NOx) emissions inventories used in GEOS-CHEM. Retrieved tropospheric NO2 columns over the United States, where NOx emissions are particularly well known, are within 18% of GEOS-CHEM columns and are strongly spatially correlated (r = 0.78, n = 288, p less than 0.005). Retrieved columns show more NO2 than GEOS-CHEM columns over the Transvaal region of South Africa and industrial regions of the northeast United States and Europe. They are lower over Houston, India, eastern Asia, and the biomass burning region of central Africa, possibly because of biases from absorbing aerosols.

Martin, Randall V.↗

Understanding Heating in Active Region Cores through Machine Learning. I. Numerical Modeling and Predicted Observables

To adequately constrain the frequency of energy deposition in active region cores in the solar corona, systematic comparisons between detailed models and observational data are needed. In this paper, we describe a pipeline for forward modeling active region emission using magnetic field extrapolations and field-aligned hydrodynamic models. We use this pipeline to predict time-dependent emission from active region NOAA 1158 for low-, intermediate-, and high-frequency nanoflares. In each pixel of our predicted multi-wavelength, time-dependent images, we compute two commonly used diagnostics: the emission measure slope and the time lag. We find that signatures of the heating frequency persist in both of these diagnostics. In particular, our results show that the distribution of emission measure slopes narrows and the mean decreases with decreasing heating frequency and that the range of emission measure slopes is consistent with past observational and modeling work. Furthermore, we find that the time lag becomes increasingly spatially coherent with decreasing heating frequency while the distribution of time lags across the whole active region becomes more broad with increasing heating frequency. In a follow-up paper, we train a random forest classifier on these predicted diagnostics and use this model to classify real observations of NOAA 1158 in terms of the underlying heating frequency.

UV radiation↗

Comparison of EKF and UKF for Robotic Missions to Mars

The batched Kalman lter works well for a wide variety of orbit determination scenarios but requires linearized approximations to the underlying dynamical system. Furthermore, implementation depends on the computation of partial derivatives for all measurements and parameters used in the lter, a non-trivial task even when analytical derivations are possible. To address these limitations, many extensions and modi - cations to the linear Kalman lter have been proposed over the years since Kalman's initial publication. Our investigation examines two popular variants of the Kalman lter, namely the Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF), within the context of spacecraft missions to Mars. The EKF relies on the same mathematical basis as the classic version of the lter, however it sequentially updates the linearization as the ltering process is conducted. In contrast, the UKF is built upon the unscented transform which samples the nearby solution space and numerically propagates perturbed trajectories as part of the update / estimation process. We assess the performance of the ltering algorithms under a variety of mission scenarios with the goal of establishing decision criteria for which implementation to use under di ering circumstances.

Ely, Todd↗