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.
Iterative method for optimal control problems for Markov chains, discussing degenerate nonlinear elliptic equations solution in optimal stochastic control theory
Abstract. We consider robust optimal experimental design (ROED) for nonlinear Bayesian inverse problems governed by partial differential equations (PDEs). An optimal design is one that maximizes some utility quantifying the quality of the solution of an inverse problem. However, the optimal design is dependent on elements of the inverse problem such as the simulation model, the prior, or the measurement error model. ROED aims to produce an optimal design that is aware of the additional uncertainties encoded in the inverse problem and remains optimal even after variations in them. We follow a worst-case scenario approach to develop a new framework for robust optimal design of nonlinear Bayesian inverse problems. The proposed framework (a) is scalable and designed for infinite-dimensional Bayesian nonlinear inverse problems constrained by PDEs; (b) develops efficient approximations of the utility, namely the expected information gain; (c) employs eigenvalue sensitivity techniques to develop analytical forms and efficient evaluation methods of the gradient of the utility with respect to the uncertainties against which we wish to be robust; and (d) employs a probabilistic optimization paradigm that properly defines and efficiently solves the resulting combinatorial max-min optimization problem. The effectiveness of the proposed approach is illustrated for optimal sensor placement problem in an inverse problem governed by an elliptic PDE.
Earth-moon and other trajectories calculation in planar elliptic restricted three body problem by slowly varying Jacobi function
By the use of a Hamiltonian formulation, a basic group velocity is defined as the derivative of frequency with respect to wavenumber keeping action density constant, and is shown to represent an incremental action velocity in the general nonlinear case. The stability treatment of Whitham and Lighthill is extended to several dimensions. The water-wave analysis of Whitham (1967) is extended to two space dimensions, and is shown to predict oblique-mode instabilities for kh smaller than 1.36. A treatment of Lighthill's (1965) solution in the one-dimensional elliptic case resolves the problem of the energy distribution in the solution past the critical time.
In the analysis of rolling-sliding concentrated contacts, such as gears, bearings and traction drives, the traction characteristics of the lubricant are of prime importance. The elastic shear modulus and limiting shear stress properties of the lubricant dictate the traction/slip characteristics and power loss associated with an EHD contact undergoing slip and/or spin. These properties can be deducted directly from the initial slope m and maximum traction coefficient micron of an experimental traction curve. In this investigation, correlation equations are presented to predict m and micron for two modern traction fluids based on the regression analysis of 334 separate traction disk machine experiments. The effects of contact pressure, temperature, surface velocity, ellipticity ratio are examined. Problems in deducing lubricant shear moduli from disk machine tests are discussed.
Origins of spectral methods, especially their relation to the Method of Weighted Residuals, are surveyed. Basic Fourier, Chebyshev, and Legendre spectral concepts are reviewed, and demonstrated through application to simple model problems. Both collocation and tau methods are considered. These techniques are then applied to a number of difficult, nonlinear problems of hyperbolic, parabolic, elliptic, and mixed type. Fluid dynamical applications are emphasized.
Origins of spectral methods, especially their relation to the Method of Weighted Residuals, are surveyed. Basic Fourier, Chebyshev, and Legendre spectral concepts are reviewed, and demonstrated through application to simple model problems. Both collocation and tau methods are considered. These techniques are then applied to a number of difficult, nonlinear problems of hyperbolic, parabolic, elliptic, and mixed type. Fluid-dynamical applications are emphasized.
In the analysis of rolling-sliding concentrated contacts, such as gears, bearings and traction drives, the traction characteristics of the lubricant are of prime importance. The elastic shear modulus and limiting shear stress properties of the lubricant dictate the traction/slip characteristics and power loss associated with an EHD contact undergoing slip and/or spin. These properties can be deducted directly from the initial slope m and maximum traction coefficient micron of an experimental traction curve. In this investigation, correlation equations are presented to predict m and micron for two modern traction fluids based on the regression analysis of 334 separate traction disk machine experiments. The effects of contact pressure, temperature, surface velocity, ellipticity ratio are examined. Problems in deducing lubricant shear moduli from disk machine tests are discussed. Previously announced in STAR as N83-20116
Origins of spectral methods, especially their relation to the method of weighted residuals, are surveyed. Basic Fourier and Chebyshev spectral concepts are reviewed and demonstrated through application to simple model problems. Both collocation and tau methods are considered. These techniques are then applied to a number of difficult, nonlinear problems of hyperbolic, parabolic, elliptic and mixzed type. Fluid dynamical applications are emphasized.
The shock cell structure and screech tone frequencies of nonaxisymmetric supersonic jets are examined. A linear shock cell model that utilizes a vortex sheet as the mixing layer of the jet is developed in order to estimate the characteristics of the shock structures and tone frequencies. Eigenvalue problems for rectangular and elliptical jets are solved, and shock cell structure solutions for the jets are derived. The shock cell spacings calculated with the model are compared to the data of Powell (1953) and Hammitt (1961); good correlation is obtained over a wide range of Mach numbers. The shock cell spacing formula is applied to the computation of screech tone frequencies of rectangular jets. It is observed that the calculated frequencies agree with the measurements of Powell and Krothapalli et. al. (1986) for the Mach number range 1.15-1.80.
The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.
Spectral element methods are high-order weighted-residual techniques for partial differential equations that combine the geometric flexibility of finite element techniques with the rapid convergence rate of spectral schemes. The theoretical foundations and numerical implementation of spectral element methods for the incompressible Navier-Stokes equations are presented, considering the construction and analysis of optimal-order spectral element discretizations for elliptic and saddle (Stokes) problems, as well as the efficient solution of the resulting discrete equations by rapidly convergent tensor-product-based iterative procedures. Several examples of spectral element simulation of moderate Reynolds number unsteady flow in complex geometry are presented.
The present demonstration of the usefulness of the integral series expansion of the stream function as a far-field computational boundary condition shows the method to require only a 10-percent/time-step increase in computational effort over alternative boundary conditions, in the case of implementation of unsteady problems using a direct elliptic solver. So long as the vorticity was encompassed within the computational domain, the method proved sufficiently accurate to yield virtually identical results for two widely different domains. While the integral-series condition yielded the best results for periodic flow, the Neumann condition gave comparable accuracy with less computation time for the steady-flow case despite its inability to treat periodic flow with vortex shedding.
Considerable effort was spent investigating the applicability of a Hamiltonian averaged over high frequency terms, where long period and secular terms up to second order in eccentricity were kept. The Hamiltonian that is given from the planar, elliptic, restricted three body problem applied to Titan-Hyperion, when the Kepler terms are also expanded to second order in small quantities and several conical transformations are carried out, is presented and discussed.
MPERRE computer program accelerates convergence of sequence of vectors by use of minimum-polynomial extrapolation (MPE) and reduced-rank extrapolation (RRE). Effective in accelerating convergences of such sequences of vectors as those obtained from iterative solution of systems of linear and nonlinear equations. In conjunction with various iterative techniques, successfully employed in finite-difference solution of large-scale elliptic boundary-value problems and problems in computational-fluid-dynamics. Only input required is sequence of vectors, convergence of which is accelerated. Program economical and easy to use. Written in FORTRAN 77.
The state of the art in algebraic multgrid (AMG) methods is discussed. The interaction between the relaxation process and the coarse grid correction necessary for proper behavior of the solution probes is discussed in detail. Sufficient conditions on relaxation and interpolation for the convergence of the V-cycle are given. The relaxation used in AMG, what smoothing means in an algebraic setting, and how it relates to the existing theory are considered. Some properties of the coarse grid operator are discussed, and results on the convergence of two-level and multilevel convergence are given. Details of an algorithm particularly studied for problems obtained by discretizing a single elliptic, second order partial differential equation are given. Results of experiments with such problems using both finite difference and finite element discretizations are presented.
Method is modification of generalized Newton-Ralphson algorithm for analyzing two-point boundary problems. It constructs sequence of solutions that converge to precise dynamic solution in the sequence limit. Program calculates periodic orbits in either circular or elliptical restricted three-body problems.
A complete numerical solution is presented to the problem of isothermal elastohydrodynamic lubrication of elliptical contacts for low-elastic-modulus materials operating under fully loaded conditions. No assumption is made for the pressure or film thickness within the contact, and compressibility and viscous effects are taken into account. Because of the dimensionless representation of the coordinates, the actual Hertzian contact ellipse becomes a circle regardless of the value of the ellipticity parameter. A minimum-film-thickness relation and a central-film-thickness relation are derived from examining 17 different cases. Contour plots for detailed illustration of the pressure distribution and film thickness in the conjunction are provided.