SEARCH · Engineering Papers
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
A technique to construct a uniformly valid perturbation series solution to a particular class of nonlinear difference equations is shown. The method allows the determination of approximations to the periodic solutions to these equations. An example illustrating the technique is presented.
A uniform geometrical optics (UGO) and an extended uniform geometrical theory of diffraction (EUTD) solution is developed for the scattering and diffraction from perfectly conducting cubic polynomial strips. The new solution overcomes the difficulties of the classic GO/UTD solution near caustics and composite shadow boundaries. The approach for constructing the UGO/EUTD solution is based on a spatial domain physical optics (PO) radiation integral representation for the scattered field which is then reduced using a uniform asymptotic procedure. New uniform reflection, zero-curvature diffraction, and edge diffraction coefficients are derived and involve the ordinary and incomplete Airy integrals as canonical functions. Higher order effects such as double edge diffraction, edge-excited creeping waves, and whispering gallery modes are not examined in this work. The UGO/EUTD solution is very efficient and provides useful physical insight into the various scattering and diffraction processes. It is also universal in nature and can be used to effectively describe the scattered fields from flat, strictly concave or convex, and concave-convex boundaries containing edges. Its accuracy is confirmed via comparison with some reference moment method (MM) results.
Perturbation theories for equations of celestial mechanics, obtaining analytical solutions in terms of Chebyshev polynomial series
We present a Backlund transformation operating on solutions of the Korteweg-de Vries equation. It provides a unifying approach to the theory of solitary waves (solutions). Recursive application of the transformation to any solution of the Korteweg-de Vries equation generates a hierarchy of solutions for which an algebraic recursion relation is obtained. Every such hierarchy (one of which is the known family of pure multisoliton solutions) thus admits a nonlinear superposition principle.
Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.
A problem which has now come to a head in general relativity, primarily as a result of the many recent detections of gravitational waves from the coalescence of compact binaries, is that there are very few viable theories against which to “test” it. To perform realistic tests of theories of gravity, we need to be able to look beyond general relativity and evaluate the consistency of a parametrized, physically acceptable, family of black hole metric alternatives with observational data from, especially, gravitational wave detections using, for example, an agnostic Bayesian approach. In this paper we further examine properties of one class of such metrics, which in fact arise as solutions of ungauged supergravity. In particular, we examine the massless, neutral, minimally coupled scalar wave equation in a general stationary, axisymmetric background metric such as that of a charged rotating black hole, when the scalar field is either time independent or in the low-frequency, near-zone limit, with a view to calculating the Love numbers of tidal perturbations, and of obtaining harmonic coordinates for the background metric. For a four-parameter family of charged asymptotically flat rotating black hole solutions of ungauged supergravity theory known as STU black holes, which includes Kaluza-Klein black holes and the Kerr-Sen black hole as special cases, we find that all time-independent solutions, and hence the harmonic coordinates of the metrics, are identical to those of the Kerr solution. In the low-frequency limit we find the scalar fields exhibit the same SL(2,R) symmetry as holds in the case of the Kerr solution. We point out extensions of our results to a wider class of metrics, which includes solutions of Einstein-Maxwell-dilaton theory.
The vortices forming in flowing water behind solid bodies are not represented correctly by the solution of the potential theory nor by Helmholtz's jets. Potential theory is unable to satisfy the condition that the water adheres at the wetted bodies, and its solutions of the fundamental hydrodynamic equations are at variance with the observation that the flow separates from the body at a certain point and sends forth a highly turbulent boundary layer into the free flow. Helmholtz's theory attempts to imitate the latter effect in such a way that it joins two potential flows, jet and still water, nonanalytical along a stream curve. The admissibility of this method is based on the fact that, at zero pressure, which is to prevail at the cited stream curve, the connection of the fluid, and with it the effect of adjacent parts on each other, is canceled. In reality, however, the pressure at these boundaries is definitely not zero, but can even be varied arbitrarily. Besides, Helmholtz's theory with its potential flows does not satisfy the condition of adherence nor explain the origin of the vortices, for in all of these problems, the friction must be taken into account on principle, according to the vortex theorem.
A method of integrating the functions defined by Scheifele and Graf (1974), arising in the elimination of time from the arguments of satellite theory using Delaunay elements in the extended phase space, is proposed. By repeated applications of an identity of Stiefel and Scheifele (1971), neglecting certain terms, and assuming a solution in the form of a product of Bessel functions, a solution expanded in the eccentricity is obtained.
Not Available
The diffraction of an arbitrary cylindrical wave by a half-plane has been treated by Rahmat-Samii and Mittra who used a spectral domain approach. In this paper, their exact solution for the total field is expressed in terms of a new integral representation. For large wave number k, two rigorous procedures are described for the exact uniform asymptotic expansion of the total field solution. The uniform expansions obtained are valid in the entire space, including transition regions around the shadow boundaries. The final results are compared with the formulations of two leading uniform theories of edge diffraction, namely, the uniform asymptotic theory and the uniform theory of diffraction. Some unique observations and conclusions are made in relating the two theories.
Simultaneous effects of temperature and concentration on tensile viscosity are investigated using a semi-empirical scaling law established to relate the concentration, temperature and deformation rate by a reduced deformation rate variable. Results show the onset of high viscosity occurring in the vicinity of the onset of non-Newtonian viscosity observed in a simple shear flow experiment. High viscosity is also compared with steady and dynamic shear data. Finally, the asymptotic solution of Hassager's kinetic theory of a bead-rod model for macromolecular solutions is compared with the data, indicating that tensile viscosity depends strongly on the molecular weight, the higher the molecular weight, the more sensitive the tensile viscosity in terms of the deformation rate.
Smart composite structures with distributed sensors and actuators have the capability to actively respond to a changing environment while offering significant weight savings and additional passive controllability through ply tailoring. Piezoelectric sensing and actuation of composite laminates is the most promising concept due to the static and dynamic control capabilities. Essential to the implementation of these smart composites are the development of accurate and efficient modeling techniques and experimental validation. This research addresses each of these important topics. A refined higher order theory is developed to model composite structures with surface bonded or embedded piezoelectric transducers. These transducers are used as both sensors and actuators for closed loop control. The theory accurately captures the transverse shear deformation through the thickness of the smart composite laminate while satisfying stress free boundary conditions on the free surfaces. The theory is extended to include the effect of debonding at the actuator-laminate interface. The developed analytical model is implemented using the finite element method utilizing an induced strain approach for computational efficiency. This allows general laminate geometries and boundary conditions to be analyzed. The state space control equations are developed to allow flexibility in the design of the control system. Circuit concepts are also discussed. Static and dynamic results of smart composite structures, obtained using the higher order theory, are correlated with available analytical data. Comparisons, including debonded laminates, are also made with a general purpose finite element code and available experimental data. Overall, very good agreement is observed. Convergence of the finite element implementation of the higher order theory is shown with exact solutions. Additional results demonstrate the utility of the developed theory to study piezoelectric actuation of composite laminates with pre-existing debonding. Significant changes in the modes shapes and reductions in the control authority result due to partially debonded actuators. An experimental investigation addresses practical issues, such as circuit design and implementation, associated with piezoelectric sensing and actuation of composite laminates. Composite specimens with piezoelectric transducers were designed, constructed and tested to validate the higher order theory. These specimens were tested with various stacking sequences, debonding lengths and gains for both open and closed loop cases. Frequency changes of 15% and damping on the order of more than 20% of critical damping, via closed loop control, was achieved. Correlation with the higher order theory is very good. Debonding is shown to adversely affect the open and closed loop frequencies, damping ratios, settling time and control authority.
The problem of acoustic radiation generated by instability waves of a compressible plane turbulent shear layer is solved. The solution provided is valid up to the acoustic far-field region. It represents a significant improvement over the solution obtained by classical hydrodynamic-stability theory which is essentially a local solution with the acoustic radiation suppressed. The basic instability-wave solution which is valid in the shear layer and the near-field region is constructed in terms of an asymptotic expansion using the method of multiple scales. This solution accounts for the effects of the slightly divergent mean flow. It is shown that the multiple-scales asymptotic expansion is not uniformly valid far from the shear layer. Continuation of this solution into the entire upper half-plane is described. The extended solution enables the near- and far-field pressure fluctuations associated with the instability wave to be determined. Numerical results show that the directivity pattern of acoustic radiation into the stationary medium peaks at 20 degrees to the axis of the shear layer in the downstream direction for supersonic flows. This agrees qualitatively with the observed noise-directivity patterns of supersonic jets.
Solution by iteration estimates performance of twisted and cambered wings of arbitrary planform. SUBAER based on linearized theory lifting-surface solution provides spanwise distribution of theoretical leading edge thrust, in addition to surface distribution of perturbation velocities. Program written in FORTRAN IV.
Simple state-space formulas are derived for all controllers solving the standard H(infinity) problem of finding, for a given number gamma greater than 0, all controllers such that the H(infinity) norm of the closed-loop transfer function is (strictly) less than gamma. It is known that a controller exists if and only if the unique stabilizing solutions to two algebraic Riccati equations are positive definite and the spectral radius of their product is less than gamma squared. Under these conditions, a parameterization of all controllers solving the problem is given as a linear fractional transformation (LFT) on a contractive stable free parameter. The state dimension of the coefficient matrix for the LFT, constructed using the two Riccati solutions, equals that of the plant and has a separation structure reminiscent of classical LQG (i.e., H2) theory. A standard H2 solution is developed in parallel.
A method developed for obtaining solutions to the aerodynamic braking problem, using impulses in the exoatmospheric phases is discussed. The solution combines primer vector theory and the results of a suboptimal atmospheric guidance program. For a specified initial and final orbit, the solution determines: (1) the minimum impulsive cost using a maximum of four impulses, (2) the optimal atmospheric entry and exit-state vectors subject to equality and inequality constraints, and (3) the optimal coast times. Numerical solutions which illustrate the characteristics of the solution are presented.
General, closed-form, analytical solutions are determined for certain classes of C-minimax control problems, several alternative mathematical theories are derived, and a controller design theory is developed to give optimal control in the presence of unmeasureable external disturbances.