SEARCH · Search NASA
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.
Explore the source record for details and available documents.
The small disturbance potential flow theory is applied to determine the lift of an airfoil in a nonuniform parallel stream. The given stream is replaced by an equivalent stream with a certain number of velocity discontinuities, and the influence of these discontinuities is obtained by the method of images. Next, this method is extended to the problem of an airfoil in a nonuniform stream of smooth velocity profile. This model allows perturbation velocity potential in a rotational undisturbed stream. A comparison of these results with numerical solutions of Euler equations indicates that, although approximate, the present method provides useful information about the interaction problem while avoiding the need to solve the Euler equations.
In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.
To determine the feasibility of coupling the output of an optical fiber to a rib waveguide in a temperature environment ranging from 20 C to 300 C, a theoretical calculation of the coupling efficiency between the two was investigated. This is a significant problem which needs to be addressed to determine whether an integrated optic device can function in a harsh temperature environment. Because the behavior of the integrated-optic device is polarization sensitive, a polarization-preserving optic fiber, via its elliptical core, was used to couple light with a known polarization into the device. To couple light energy efficiently from an optical fiber into a channel waveguide, the design of both components should provide for well-matched electric field profiles. The rib waveguide analyzed was the light input channel of an integrated-optic pressure sensor. Due to the complex geometry of the rib waveguide, there is no analytical solution to the wave equation for the guided modes. Approximation or numerical techniques must be utilized to determine the propagation constants and field patterns of the guide. In this study, three solution methods were used to determine the field profiles of both the fiber and guide: the effective-index method (EIM), Marcatili's approximation, and a Fourier method. These methods were utilized independently to calculate the electric field profile of a rib channel waveguide and elliptical fiber at two temperatures, 20 C and 300 C. These temperatures were chosen to represent a nominal and a high temperature that the device would experience. Using the electric field profile calculated from each method, the theoretical coupling efficiency between the single-mode optical fiber and rib waveguide was calculated using the overlap integral and results of the techniques compared. Initially, perfect alignment was assumed and the coupling efficiency calculated. Then, the coupling efficiency calculation was repeated for a range of transverse offsets at both temperatures. Results of the calculation indicate a high coupling efficiency can be achieved when the two components were properly aligned. The coupling efficiency was more sensitive to alignment offsets in the y direction than the x, due to the elliptical modal profile of both components. Changes in the coupling efficiency over temperature were found to be minimal.
A least squares method is presented for computing approximate solutions of indefinite partial differential equations of the mixed type such as those that arise in connection with transonic flutter analysis. The method retains the advantages of finite difference schemes namely simplicity and sparsity of the resulting matrix system. However, it offers some great advantages over finite difference schemes. First, the method is insensitive to the value of the forcing frequency, i.e., the resulting matrix system is always symmetric and positive definite. As a result, iterative methods may be successfully employed to solve the matrix system, thus taking full advantage of the sparsity. Furthermore, the method is insensitive to the type of the partial differential equation, i.e., the computational algorithm is the same in elliptic and hyperbolic regions. In this work the method is formulated and numerical results for model problems are presented. Some theoretical aspects of least squares approximations are also discussed.
The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.
A simple method for milling accurate off-axis parabolic mirrors with a computer-controlled milling machine is discussed. For machines with a built-in circle-cutting routine, an exact paraboloid can be milled with few computer commands and without the use of the spherical or linear approximations. The proposed method can be adapted easily to cut off-axis sections of elliptical or spherical mirrors.
The conservation and convergence properties of spectral Fourier methods for the numerical approximation of the Korteweg-de Vries equation are analyzed. It is proved that the (aliased) collocation pseudospectral method enjoys the same convergence properties as the spectral Galerkin method, which is less effective from the computational point of view. This result provides a precise mathematical answer to a question raised by several authors in recent years.
The NASA Engine Performance Program (NEPP) can configure and analyze almost any type of gas turbine engine that can be generated through the interconnection of a set of standard physical components. In addition, the code can optimize engine performance by changing adjustable variables under a set of constraints. However, for engine cycle problems at certain operating points, the NEPP code can encounter difficulties: nonconvergence in the currently implemented Powell's optimization algorithm and deficiencies in the Newton-Raphson solver during engine balancing. A project was undertaken to correct these deficiencies. Nonconvergence was avoided through a cascade optimization strategy, and deficiencies associated with engine balancing were eliminated through neural network and linear regression methods. An approximation-interspersed cascade strategy was used to optimize the engine's operation over its flight envelope. Replacement of Powell's algorithm by the cascade strategy improved the optimization segment of the NEPP code. The performance of the linear regression and neural network methods as alternative engine analyzers was found to be satisfactory. This report considers two examples-a supersonic mixed-flow turbofan engine and a subsonic waverotor-topped engine-to illustrate the results, and it discusses insights gained from the improved version of the NEPP code.
The expectation of the solution process in a stochastic operator equation can be obtained from averaged equations only under very special circumstances. Conditions for validity are given and the significance and validity of the approximation in widely used hierarchy methods and the ?self-consistent field' approximation in nonequilibrium statistical mechanics are clarified. The error at any level of the hierarchy can be given and can be avoided by the use of the iterative method.
Proposed method of obtaining approximately desired radiation or reception pattern from antenna that includes reflector based on concept of superposition of electromagnetic fields generated by multiple feedhorns or feed antenna elements arrayed at various positions near reflector and excited at electronically adjustable magnitudes and phases. In intended application, reflector nominally paraboloidal, feed elements N feedhorns in hexagonal array, and method used to compensate for deviations of real reflector surface from nominal paraboloidal shape. Method and concept also applicable to electronic beam steering and electronic antenna compensation in other situations.
This study will demonstrate the importance of application of computer numerically controlled (CNC) machines in generation of gear tooth surfaces with new topology. This topology decreases gear vibration and will extend the gear capacity and service life. A preliminary investigation by a tooth contact analysis (TCA) program has shown that gear tooth surfaces in line contact (for instance, involute helical gears with parallel axes, worm gear drives with cylindrical worms, etc.) are very sensitive to angular errors of misalignment that cause edge contact and an unfavorable shape of transmission errors and vibration. The new topology of gear tooth surfaces is based on the localization of bearing contact, and the synthesis of a predesigned parabolic function of transmission errors that is able to absorb a piecewise linear function of transmission errors caused by gear misalignment. The report will describe the following topics: description of kinematics of CNC machines with six degrees of freedom that can be applied for generation of gear tooth surfaces with new topology. A new method for grinding of gear tooth surfaces by a cone surface or surface of revolution based on application of CNC machines is described. This method provides an optimal approximation of the ground surface to the given one. This method is especially beneficial when undeveloped ruled surfaces are to be ground. Execution of motions of the CNC machine is also described. The solution to this problem can be applied as well for the transfer of machine tool settings from a conventional generator to the CNC machine. The developed theory required the derivation of a modified equation of meshing based on application of the concept of space curves, space curves represented on surfaces, geodesic curvature, surface torsion, etc. Condensed information on these topics of differential geometry is provided as well.
The commonly used minimum mass power law representation of the pre-solar nebula is reanalyzed using a new cumulative-mass-model. This model predicts a smoother surface density approximation compared with methods based on direct computation of surface density. The density is quantified using two independent analytical formulations. First, a best-fit transcendental function is applied directly to the basic planetary data. Next a solution to the time-dependent disk evolution equation is parametrically adapted to the solar nebula data. The latter model is shown to be a good approximation to the finite-size early Solar Nebula, and by extension to other extra solar protoplanetary disks.
Current high-order methods such as discontinuous Galerkin and/or flux reconstruction can provide effective discretization for the spatial derivatives. Together with a time discretization, such methods result in either too small a time step size in the case of an explicit scheme or a very large system in the case of an implicit one. To tackle these problems, two new high-order space-time schemes for conservation laws are introduced: the first is explicit and the second, implicit. The explicit method here, also called the moment scheme, achieves a Courant-Friedrichs-Lewy (CFL) condition of 1 for the case of one-spatial dimension regardless of the degree of the polynomial approximation. (For standard explicit methods, if the spatial approximation is of degree p, then the time step sizes are typically proportional to 1/p(exp 2)). Fourier analyses for the one and two-dimensional cases are carried out. The property of super accuracy (or super convergence) is discussed. The implicit method is a simplified but optimal version of the discontinuous Galerkin scheme applied to time. It reduces to a collocation implicit Runge-Kutta (RK) method for ordinary differential equations (ODE) called Radau IIA. The explicit and implicit schemes are closely related since they employ the same intermediate time levels, and the former can serve as a key building block in an iterative procedure for the latter. A limiting technique for the piecewise linear scheme is also discussed. The technique can suppress oscillations near a discontinuity while preserving accuracy near extrema. Preliminary numerical results are shown
The problem of the quantum treatment of nonreactive collisions of two simple molecules or of an atom and a molecule is considered mathematically. Kinematical structure theorems analogous to the Wigner-Eckart theorem for scalar operators are derived for the transition operator (T-operator) of a two-molecule system, where, in addition to the usual conservation laws, certain additional conservation laws are presumed satisfied in a collision. The additional conservation laws were examined in order to develop a testable approximation scheme (CS method) for the specification of the kinematics of a system of two rigid diatomic molecules. The rand p-helicity basis states for a two-rigid rotor system are defined and unitary transformations established between these and conventional bases. It is found that only questionable criteria exist for testing the on-the-energy-shell matrix of a T(energy) operator for the presence of properties of two of the additional conservation laws.
A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.
An abstract approximation theory and computational methods are developed for the determination of optimal linear-quadratic feedback control, observers and compensators for infinite dimensional discrete-time systems. Particular attention is paid to systems whose open-loop dynamics are described by semigroups of operators on Hilbert spaces. The approach taken is based on the finite dimensional approximation of the infinite dimensional operator Riccati equations which characterize the optimal feedback control and observer gains. Theoretical convergence results are presented and discussed. Numerical results for an example involving a heat equation with boundary control are presented and used to demonstrate the feasibility of the method.
The surface finish of parts produced by additive manufacturing processes is much rougher than the surface finish generated by machining processes, and a rougher surface can reduce the fatigue strength of a part. This paper discusses an effort to quantify that reduction of strength in high-cycle fatigue for selective laser melt (SLM) coupons. A high-cycle fatigue (HCF) knockdown factor was estimated for Inconel 718, manufactured with the SLM process. This factor is the percentage reduction from the maximum stress in fatigue for low-stress ground (LSG) specimens to the maximum stress of those left with the original surface condition at the same fatigue life. Specimens were provided by a number of vendors, free to use their "best practice"; only one heat treat condition was considered; and several test temperatures were characterized, including room temperature, 800F, 1000F, and 1200F. The 1000F data had a large variance, and was omitted from consideration in this document. A first method used linear approximations extracted from the graphs, and only where data was available for both. A recommended knockdown factor of the as-built surface condition (average roughness of approximately 245 micro-inches/inch) versus low-stress ground condition (roughness no more than 4 micro-inches/inch) was established at approximately 1/3 or 33%. This is to say that for the as-built surface condition, a maximum stress of 2/3 of the stress for LSG can be expected to produce a similar life in the as-built surface condition. In this first evaluation, the knockdown factor did not appear to be a function of temperature. A second approach, the "KP method", incorporated the surface finish measure into a new parameter termed the pseudo-stress intensity factor, Kp, which was formulated to be similar to the fracture mechanics stress intensity factor. Using Kp, the variance seemed to be reduced across all sources, and knockdown factors were estimated using Kp over the range where data occurred. A plot of the results suggests that the knockdown factor is a function of temperature, and that for low lives the knockdown might be lower than the knockdown observed above about one million cycles, where it tended to stabilize. This was not universal for all temperatures tested. The higher temperature tests are thought to be influenced by the test temperature, which perhaps continued the aging process. Further evaluation of the method is suggested.