Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “KERNELS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

The Determination of the Large-Scale Circulation of the Pacific Ocean from Satellite Altimetry using Model Green's Functions

A Green's function method for obtaining an estimate of the ocean circulation using both a general circulation model and altimetric data is demonstrated. The fundamental assumption is that the model is so accurate that the differences between the observations and the model-estimated fields obey a linear dynamics. In the present case, the calculations are demonstrated for model/data differences occurring on very a large scale, where the linearization hypothesis appears to be a good one. A semi-automatic linearization of the Bryan/Cox general circulation model is effected by calculating the model response to a series of isolated (in both space and time) geostrophically balanced vortices. These resulting impulse responses or 'Green's functions' then provide the kernels for a linear inverse problem. The method is first demonstrated with a set of 'twin experiments' and then with real data spanning the entire model domain and a year of TOPEX/POSEIDON observations. Our present focus is on the estimate of the time-mean and annual cycle of the model. Residuals of the inversion/assimilation are largest in the western tropical Pacific, and are believed to reflect primarily geoid error. Vertical resolution diminishes with depth with 1 year of data. The model mean is modified such that the subtropical gyre is weakened by about 1 cm/s and the center of the gyre shifted southward by about 10 deg. Corrections to the flow field at the annual cycle suggest that the dynamical response is weak except in the tropics, where the estimated seasonal cycle of the low-latitude current system is of the order of 2 cm/s. The underestimation of observed fluctuations can be related to the inversion on the coarse spatial grid, which does not permit full resolution of the tropical physics. The methodology is easily extended to higher resolution, to use of spatially correlated errors, and to other data types.

Stammer, Detlef↗
NASA NTRS1997

Measurement and Analysis of Failures in Computer Systems

This thesis presents a study of software failures spanning several different releases of Tandem's NonStop-UX operating system running on Tandem Integrity S2(TMR) systems. NonStop-UX is based on UNIX System V and is fully compliant with industry standards, such as the X/Open Portability Guide, the IEEE POSIX standards, and the System V Interface Definition (SVID) extensions. In addition to providing a general UNIX interface to the hardware, the operating system has built-in recovery mechanisms and audit routines that check the consistency of the kernel data structures. The analysis is based on data on software failures and repairs collected from Tandem's product report (TPR) logs for a period exceeding three years. A TPR log is created when a customer or an internal developer observes a failure in a Tandem Integrity system. This study concentrates primarily on those TPRs that report a UNIX panic that subsequently crashes the system. Approximately 200 of the TPRs fall into this category. Approximately 50% of the failures reported are from field systems, and the rest are from the testing and development sites. It has been observed by Tandem developers that fewer cases are encountered from the field than from the test centers. Thus, the data selection mechanism has introduced a slight skew.

Thakur, Anshuman↗
NASA NTRS1997

A General Reversible Hereditary Constitutive Model: Theoretical Developments - Part 1

Using an internal-variable formalism as a starting point, we describe the viscoelastic extension of a previously-developed viscoplasticity formulation of the complete potential structure type. It is mainly motivated by experimental evidence for the presence of rate/time effects in the so-called quasilinear, reversible, material response range. Several possible generalizations are described, in the general format of hereditary-integral representations for non-equilibrium, stress-type, state variables, both for isotropic as well as anisotropic materials. In particular, thorough discussions are given on the important issues of thermodynamic admissibility requirements for such general descriptions, resulting in a set of explicit mathematical constraints on the associated kernel (relaxation and creep compliance) functions. In addition, a number of explicit, integrated forms are derived, under stress and strain control to facilitate the parametric and qualitative response characteristic studies reported here, as well as to help identify critical factors in the actual experimental characterizations from test data that will be reported in Part II.

Saleeb, A. F.↗
NASA NTRS1997

Frictionless Contact of Multilayered Composite Half Planes Containing Layers With Complex Eigenvalues

A previously developed local-global stiffness matrix methodology for the response of a composite half plane, arbitrarily layered with isotropic, orthotropic or monoclinic plies, to indentation by a rigid parabolic punch is further extended to accommodate the presence of layers with complex eigenvalues (e.g., honeycomb or piezoelectric layers). First, a generalized plane deformation solution for the displacement field in an orthotropic layer or half plane characterized by complex eigenvalues is obtained using Fourier transforms. A local stiffness matrix in the transform domain is subsequently constructed for this class of layers and half planes, which is then assembled into a global stiffness matrix for the entire multilayered half plane by enforcing continuity conditions along the interfaces. Application of the mixed boundary condition on the top surface of the half plane indented by a rigid punch results in an integral equation for the unknown pressure in the contact region. The integral possesses a divergent kernel which is decomposed into Cauchy-type and regular parts using the asymptotic properties of the local stiffness matrix and a relationship between Fourier and finite Hilbert transform of the contact pressure. The solution of the resulting singular integral equation is obtained using a collocation technique based on the properties of orthogonal polynomials developed by Erdogan and Gupta. Examples are presented that illustrate the important influence of low transverse properties of layers with complex eigenvalues, such as those exhibited by honeycomb, on the load versus contact length response and contact pressure distributions for half planes containing typical composite materials.

Zhang, Wang↗
NASA NTRS1997

Land Surface Temperature Measurements from EOS MODIS Data

We made modifications to the linear kernel bidirectional reflectance distribution function (BRDF) models from Roujean et al. and Wanner et al. that extend the spectral range into the thermal infrared (TIR). With these TIR BRDF models and the IGBP land-cover product, we developed a classification-based emissivity database for the EOS/MODIS land-surface temperature (LST) algorithm and used it in version V2.0 of the MODIS LST code. Two V2.0 LST codes have been delivered to the MODIS SDST, one for the daily L2 and L3 LST products, and another for the 8-day 1km L3 LST product. New TIR thermometers (broadband radiometer with a filter in the 10-13 micron window) and an IR camera have been purchased in order to reduce the uncertainty in LST field measurements due to the temporal and spatial variations in LST. New improvements have been made to the existing TIR spectrometer in order to increase its accuracy to 0.2 C that will be required in the vicarious calibration of the MODIS TIR bands.

Wan, Zhengming↗
NASA NTRS1994

Radiative Heat Transfer in Finite Cylindrical Enclosures with Nonhomogeneous Participating Media

Results of a numerical solution for radiative heat transfer in homogeneous and nonhomogeneous participating media are presented. The geometry of interest is a finite axisymmetric cylindrical enclosure. The integral formulation for radiative transport is solved by the YIX method. A three-dimensional solution scheme is applied to two-dimensional axisymmetric geometry to simplify kernel calculations and to avoid difficulties associated with treating boundary conditions. As part of the effort to improve modeling capabilities for turbulent jet diffusion flames, predicted distributions for flame temperature and soot volume fraction are used to calculate radiative heat transfer from soot particles in such flames. It is shown that the nonhomogeneity of radiative property has very significant effects. The peak value of the divergence of radiative heat flux could be underestimated by 2 factor of 7 if a mean homogeneous radiative property is used. Since recent studies have shown that scattering by soot agglomerates is significant in flames, the effect of magnitude of scattering is also investigated and found to be nonnegligible.

Hsu, Pei-Feng↗
NASA NTRS1998

Long-Term file activity patterns in a UNIX workstation environment

As mass storage technology becomes more affordable for sites smaller than supercomputer centers, understanding their file access patterns becomes crucial for developing systems to store rarely used data on tertiary storage devices such as tapes and optical disks. This paper presents a new way to collect and analyze file system statistics for UNIX-based file systems. The collection system runs in user-space and requires no modification of the operating system kernel. The statistics package provides details about file system operations at the file level: creations, deletions, modifications, etc. The paper analyzes four months of file system activity on a university file system. The results confirm previously published results gathered from supercomputer file systems, but differ in several important areas. Files in this study were considerably smaller than those at supercomputer centers, and they were accessed less frequently. Additionally, the long-term creation rate on workstation file systems is sufficiently low so that all data more than a day old could be cheaply saved on a mass storage device, allowing the integration of time travel into every file system.

Gibson, Timothy J.↗
NASA NTRS1960

Free-Spinning-Tunnel Investigation of a 1/20-Scale Model of an Unswept-Wing Jet-Propelled Trainer Airplane

A flutter analysis employing the kernel function for three- dimensional, subsonic, compressible flow is applied to a flutter-tested tail surface which has an aspect ratio of 3.5, a taper ratio of 0.15, and a leading-edge sweep of 30 deg. Theoretical and experimental results are compared at Mach numbers from 0.75 to 0.98. Good agreement between theoretical and experimental flutter dynamic pressures and frequencies is achieved at Mach numbers to 0.92. At Mach numbers from 0.92 to 0.98, however, a second solution to the flutter determinant results in a spurious theoretical flutter boundary which is at a much lower dynamic pressure and at a much higher frequency than the experimental boundary.

Bowman, James S., Jr.↗
NASA NTRS1998

Robust Multigrid Smoothers for Three Dimensional Elliptic Equations with Strong Anisotropies

We discuss the behavior of several plane relaxation methods as multigrid smoothers for the solution of a discrete anisotropic elliptic model problem on cell-centered grids. The methods compared are plane Jacobi with damping, plane Jacobi with partial damping, plane Gauss-Seidel, plane zebra Gauss-Seidel, and line Gauss-Seidel. Based on numerical experiments and local mode analysis, we compare the smoothing factor of the different methods in the presence of strong anisotropies. A four-color Gauss-Seidel method is found to have the best numerical and architectural properties of the methods considered in the present work. Although alternating direction plane relaxation schemes are simpler and more robust than other approaches, they are not currently used in industrial and production codes because they require the solution of a two-dimensional problem for each plane in each direction. We verify the theoretical predictions of Thole and Trottenberg that an exact solution of each plane is not necessary and that a single two-dimensional multigrid cycle gives the same result as an exact solution, in much less execution time. Parallelization of the two-dimensional multigrid cycles, the kernel of the three-dimensional implicit solver, is also discussed. Alternating-plane smoothers are found to be highly efficient multigrid smoothers for anisotropic elliptic problems.

Llorente, Ignacio M.↗
NASA NTRS1998

CAPRI (Computational Analysis PRogramming Interface): A Solid Modeling Based Infra-Structure for Engineering Analysis and Design Simulations

CAPRI is a CAD-vendor neutral application programming interface designed for the construction of analysis and design systems. By allowing access to the geometry from within all modules (grid generators, solvers and post-processors) such tasks as meshing on the actual surfaces, node enrichment by solvers and defining which mesh faces are boundaries (for the solver and visualization system) become simpler. The overall reliance on file 'standards' is minimized. This 'Geometry Centric' approach makes multi-physics (multi-disciplinary) analysis codes much easier to build. By using the shared (coupled) surface as the foundation, CAPRI provides a single call to interpolate grid-node based data from the surface discretization in one volume to another. Finally, design systems are possible where the results can be brought back into the CAD system (and therefore manufactured) because all geometry construction and modification are performed using the CAD system's geometry kernel.

Haimes, Robert↗
NASA NTRS1997

Nonlinear Rescaling and Proximal-Like Methods in Convex Optimization

The nonlinear rescaling principle (NRP) consists of transforming the objective function and/or the constraints of a given constrained optimization problem into another problem which is equivalent to the original one in the sense that their optimal set of solutions coincides. A nonlinear transformation parameterized by a positive scalar parameter and based on a smooth scaling function is used to transform the constraints. The methods based on NRP consist of sequential unconstrained minimization of the classical Lagrangian for the equivalent problem, followed by an explicit formula updating the Lagrange multipliers. We first show that the NRP leads naturally to proximal methods with an entropy-like kernel, which is defined by the conjugate of the scaling function, and establish that the two methods are dually equivalent for convex constrained minimization problems. We then study the convergence properties of the nonlinear rescaling algorithm and the corresponding entropy-like proximal methods for convex constrained optimization problems. Special cases of the nonlinear resealing algorithm are presented. In particular a new class of exponential penalty-modified barrier functions methods is introduced.

Polyak, Roman↗
NASA NTRS1999

Object-Oriented Design for Sparse Direct Solvers

We discuss the object-oriented design of a software package for solving sparse, symmetric systems of equations (positive definite and indefinite) by direct methods. At the highest layers, we decouple data structure classes from algorithmic classes for flexibility. We describe the important structural and algorithmic classes in our design, and discuss the trade-offs we made for high performance. The kernels at the lower layers were optimized by hand. Our results show no performance loss from our object-oriented design, while providing flexibility, case of use, and extensibility over solvers using procedural design.

Dobrian, Florin↗
NASA NTRS1999

Analysis of an Interface Crack for a Functionally Graded Strip Sandwiched between Two Homogeneous Layers of Finite Thickness

The interface crack problem for a composite layer that consists of a homogeneous substrate, coating and a non-homogeneous interface was formulated for singular integral equations with Cauchy kernels and integrated using the Lobatto-Chebyshev collocation technique. Mixed-mode Stress Intensity Factors and Strain Energy Release Rates were calculated. The Stress Intensity Factors were compared for accuracy with relevant results previously published. The parametric studies were conducted for the various thickness of each layer and for various non-homogeneity ratios. Particular application to the Zirconia thermal barrier on steel substrate is demonstrated.

Shbeeh, N. I.↗
NASA NTRS1999

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations That Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + definite integral of g(x, t)F(t,y(t))dt with limits between 0 and 1,0 less than or equal to x les than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integral equations arise, e.g., when one applied Green's function techniques to nonlinear two-point boundary value problems of the form y "(x) =f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and y(l) = y(sub l), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trepezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal rule, thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram↗
NASA NTRS1998

Northern and Southern Hemisphere Ground-Based Infrared Spectroscopic Measurements of Tropospheric Carbon Monoxide and Ethane

Time series of CO and C2H, measurements have been derived from high-resolution infrared solar spectra recorded in Lauder, New Zealand (45.0 degrees S, 169.7 degrees E, altitude 0.37 km), and at the U.S. National Solar Observatory (31.9 degrees N, 11, 1.6 degrees W, altitude 2.09 km) on Kitt Peak. Lauder observations were obtained between July 1993 and November 1997, while the Kitt Peak measurements were recorded between May 1977 and December 1997. Both databases were analyzed with spectroscopic parameters that included significant improvements for C2H6 relative to previous studies. Target CO and C2H6 lines were selected to achieve similar vertical samplings based on averaging kernels. These calculations show that partial columns from layers extending from the surface to the mean tropopause and from the mean tropopause to 100 km are nearly independent. Retrievals based on a semiempirical application of the Rodgers optimal estimation technique are reported for the lower layer, which has a broad maximum in sensitivity in the upper troposphere. The Lauder CO and C2H, partial columns exhibit highly asymmetrical seasonal cycles with minima in austral autumn and sharp peaks in austral spring. The spring maxima are the result of tropical biomass burning emissions followed by deep convective vertical transport to the upper troposphere and long-range horizontal transport. Significant year-to-year variations are observed for both CO and C2H6, but the measured trends, (+0.37 +/- 0.57)% yr(exp -1) and (-0.64 +/- 0.79)% yr(exp -1), I sigma, respectively, indicate no significant long-term changes. The Kitt Peak data also exhibit CO and C2H6, seasonal variations in the lower layer with trends equal to (-0.27 +/- 0.17)% yr(exp -1) and (-1.20 +/- 0.35)% yr(exp -1), 1 sigma, respectively. Hence a decrease in the Kitt Peak tropospheric C2H6 column has been detected, though the CO trend is not significant. Both measurement sets are compared with previous observations, reported trends, and three-dimensional model calculations.

Rinsland, Curtis P.↗
NASA NTRS1999

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations that Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + integral(0 to 1) g(x,t) F(t, y(t)) dt, 0 less than or equal to x less than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integrals equations arise, e.g., when one applies Green's function techniques to nonlinear two-point boundary value problems of the form U''(x) = f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and g(l) = y(sub 1), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trapezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram↗
NASA NTRS1998

Northern and Southern Hemisphere Ground-Based Infrared Spectroscopic Measurements of Tropospheric Carbon Monoxide and Ethane

Time series of CO and C2H6 measurements have been derived from high resolution infrared solar spectra recorded in Lauder, New Zealand (45.0 deg S, 169.7 deg E, altitude 0.37 km) and at the U. S. National Solar Observatory (31.90 deg N, 111.6 deg W, altitude 2.09 km) on Kitt Peak. Lauder observations were obtained between July 1993 and November 1997 while the Kitt Peak measurements were recorded between May 1977 and December 1997. Both databases were analyzed with spectroscopic parameters that included significant improvements for C2H6 relative to previous studies. Target CO and C2H6 lines were selected to achieve similar vertical samplings based on averaging kernels. These calculations show that partial columns from layers extending from the surface to the mean tropopause and from the mean tropopause to 100 km are nearly independent. Retrievals based on a semiempirical application of the Rodgers optimal estimation technique are reported for the lower layer, which has a broad maximum in sensitivity in the upper troposphere. The Lauder CO and C2H6 partial columns exhibit highly asymmetrical seasonal cycles with minima in austral autumn and sharp peaks in austral spring. The spring maxima are the result of tropical biomass burning emissions followed by deep convective vertical transport to the upper troposphere and long-range horizontal transport. Significant year-to-year variations are observed for both CO and C2H6, but the measured trends, (+0.37 +/- 0.57)%/ yr and (-0.64 +/- 0.79)%/ yr, 1 sigma, respectively, indicate no significant long-term changes. The Kitt Peak data also exhibit CO and C2H6 seasonal variations in the lower layer with trends equal to (-0.27 +/- 0.17)%/ yr and (-1.20 +/- 0.35)%/ yr, 1 sigma, respectively. Hence, a decrease in the Kitt Peak tropospheric C2H6 column has been detected, though the CO trend is not significant. Both measurement sets are compared with previous observations, reported trends, and three-dimensional model calculations.

Rinsland, Curtis P.↗
NASA NTRS1999

Analysis of a Generally Oriented Crack in a Functionally Graded Strip Sandwiched Between Two Homogeneous Half Planes

The driving forces for a generally oriented crack embedded in a Functionally Graded strip sandwiched between two half planes are analyzed using singular integral equations with Cauchy kernels, and integrated using Lobatto-Chebyshev collocation. Mixed-mode Stress Intensity Factors (SIF) and Strain Energy Release Rates (SERR) are calculated. The Stress Intensity Factors are compared for accuracy with previously published results. Parametric studies are conducted for various nonhomogeneity ratios, crack lengths. crack orientation and thickness of the strip. It is shown that the SERR is more complete and should be used for crack propagation analysis.

Shbeeb, N.↗