Theory and methods related to the singular-function expansion and Landweber's iteration for integral equations of the first kind
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
Large space structures (LSSs) and other dynamical systems of current interest are often extremely complex assemblies of rigid and flexible bodies subjected to kinematical constraints. A formulation is presented for the governing equations of constrained multibody systems via the application of singular value decomposition (SVD). The resulting equations of motion are shown to be of minimum dimension.
Review of recent work in recasting the energy type of approximation to aircraft flight in terms of singular perturbation theory and in extending this theory to three-dimensional maneuvers. Singular perturbations for differential equations arising in optimal control are first examined for a system of fairly general form but low order. The attitude dynamics for optimal flight of a rocket in vacuum is then studied as an introductory, and fairly transparent, example. The question of the choice of variables is then discussed. Finally, optimal aircraft flight in various reduced-order approximations is investigated. In particular, the problem of three-dimensional aircraft flight is formulated for singular perturbation treatment, and possibilities for decoupling into several lower-order problems are illustrated. The use of reduced-order approximation facilitates numerical computations by reducing the number of multiplier initial values that must be determined simultaneously and by improving the conditioning of the differential equations.
We present a computational methodology (a novel Nystrom approach based on use of a non-overlapping patch technique and Chebyshev discretizations) for efficient solution of problems of acoustic and electromagnetic scattering by open surfaces. Our integral equation formulations (1) Incorporate, as ansatz, the singular nature of open-surface integral-equation solutions, and (2) For the Electric Field Integral Equation (EFIE), use analytical regularizes that effectively reduce the number of iterations required by iterative linear-algebra solution based on Krylov-subspace iterative solvers.
Using a time domain method based on the Ffowcs Williams-Hawkings equation, a reliable explanation is provided for the origin of singularities observed in the numerical prediction of supersonic propeller noise. In the last few years Tam and, more recently, Amiet have analyzed the phenomenon from different points of view. The method proposed here offers a clear interpretation of the singularities based on a new description of sources, relating to the behavior of lines where the propeller blade surface exhibit slope discontinuity.
In designing a robot control system for dual arm configurations, the control engineer is faced with two challenges: to derive the equations of motion for a given situation, and to meet certain desired control requirements (for instance, minimum energy). The former may involve closed kinematic chains, such as the case when the two arms are grasping a common object. The latter usually involves nonlinear optimization. These issues are considered in the context of transporting an inertial load using two planar three-link arms. A generalized 'reduction transformation' is applied to the dynamics to remove the singularity in the system equations. A suboptimal minimum energy method is presented to reduce a difficult 12-state, six-control nonlinear optimization to two independent, nonconflicting suboptimizations. A simulation example is provided to illustrate the degree of energy reduction possible using the optimal arm torque distribution that was developed.
The problem of plane stagnation point flow with freestream turbulence is examined from a basic theoretical standpoint. It is argued that the singularity which arises from the standard kappa-epsilon model is not due to a defect in the model but results from the use of an inconsistent freestream boundary condition. The inconsistency lies in the implementation of a production equals dissipation equilibrium hypothesis in conjunction with a freestream mean velocity field that corresponds to homogeneous plane strain - a turbulent flow which does not reach such a simple equilibrium. Consequently, the adjustment that has been made in the constants of the epsilon-transport equation to eliminate this singularity is not self-consistent since it is tantamount to artificially imposing an equilibrium structure on a turbulent flow which is known not to have one.
The problem of plane stagnation point flow with freestream turbulence is examined from a basic theoretical standpoint. It is argued that the singularity which arises from the standard kappa-epsilon model is not due to a defect in the model but results from the use of an inconsistent freestream boundary condition. The inconsistency lies in the implementation of a production equals dissipation equilibrium hypothesis in conjunction with a freestream mean velocity field that corresponds to homogeneous plane strain - a turbulent flow which does not reach such a simple equilibrium. Consequently, the adjustment that has been made in the constants of the epsilon-transport equation to eliminate this singularity is not self-consistent since it is tantamount to artificially imposing an equilibrium structure on a turbulent flow which is known not to have one.
We implemented a two- and three-dimensional finite difference/front tracking technique to solve liquid-vapor phase change problems. The mathematical and the numerical features of the method were explained in great detail in our previous reports, Briefly, we used a single formula representation which incorporated jump conditions into the governing equations. The interfacial terms were distributed as singular terms using delta functions so that the governing equations would be the same as conventional conservation equations away from the interface and in the vicinity of the interface they would provide correct jump conditions. We used a fixed staggered grid to discretize these equations and an unstructured grid to explicitly track the front. While in two dimensions the front was simply a connection of small line segments, in three dimensions it was represented by a connection of small triangular elements. The equations were written in conservative forms and during the course of computations we used regriding to control the size of the elements of the unstructured grid. Moreover, we implemented a coalescence in two dimensions which allowed the merging of different fronts or two segments of the same front when they were sufficiently close. We used our code to study thermocapillary migration of bubbles, burst of bubbles at a free surface, buoyancy-driven interactions of bubbles, evaporation of drops, rapid evaporation of an interface, planar solidification of an undercooled melt, dendritic solidification, and a host of other problems cited in the reference.
Equations for inversion of singular matrix - special form of generalized inverse of arbitrary complex matrix
The problem of an elastic solid containing a semi-infinite plane crack subjected to concentrated shears parallel to the edge of the crack is considered. A closed form solution using four distinct harmonic functions (none of which can be taken arbitrarily) is found to satisfy the finite displacement and inverse square root stress singularity at the edge of the crack. Explicit expressions in terms of elementary functions are given for the distribution of stress and displacement in the solid. These are obtained by employing Fourier and Kontorovich-Lebedev integral transforms and certain singular solutions of Laplace equations in three dimensions. The variations of the intensity of the local stress field along the crack border are shown graphically. An example is presented, which is in contrast with the conclusion established in the literature that one of the four Papkovich-Neuber functions in three-dimensional elasticity may be arbitrarily set to zero.
The problem of an elastic solid containing a semi-infinite plane crack subjected to concentrated shears parallel to the edge of the crack is considered in this paper. A closed form solution using four harmonic functions is found to satisfy the finite displacement and inverse square root stress singularity at the edge of the crack. Explicit expressions in terms of elementary functions are given for the distribution of stress and displacement in the solid. These are obtained by employing Fourier and Kontorovich-Lebedev integral transforms and certain singular solutions of Laplace equations in three dimensions. The variations of the intensity of the local stress field along the crack border are shown graphically.
The aerodynamic interaction between the wing and an inviscid jet with Mach number nonuniformity is formulated by using a two-vortex-sheet model for the jet. One of the vortex sheets accounts for the induced jet flow and the other the induced outer flow. No additional source distribution is needed for the jet at an angle of attack. The above problem is solved by satisfying the jet and wing tangency and the jet pressure-continuity conditions and using a quasi vortex lattice method for computing the induced flow field. The latter method is derived through theoretical consideration by properly accounting for singularities present in the equations and possesses the same simplicity and generality as the conventional vortex lattice method but has a better rate of numerical convergence. The resulting system of algebraic equations is solved by Purcell's vector method. The numerical formulation is first applied to the wing-slipstream interaction problem. Results for one centered-jet configuration are compared with those predicted by some existing theories.
The evolution of long, finite amplitude Rossby waves in a horizontally sheared zonal current is studied. The wave evolution is described by the Korteweg-de Vries equation or the modified Korteweg-de Vries equation depending on the atmospheric stratification. In either case, the cross-stream modal structure of these waves is given by the long-wave limit of the neutral eigensolutions of the barotropic stability equation. Both non-singular and singular eigensolutions are considered and the appropriate analysis is developed to yield a uniformly valid description of the motion in the critical-layer region where the wave speed matches the flow velocity. The analysis demonstrates that coherent, propagating, eddy structures can exist in stable shear flows and that these eddies have peculiar interaction properties quite distinct from the traditional views of turbulent motion.
The Sequential Unconstrained Minimization Technique (SUMT) offers an easy way of solving nonlinearly constrained problems. However, this algorithm frequently suffers from the need to minimize an ill-conditioned penalty function. An ill-conditioned minimization problem can be solved very effectively by posing the problem as one of integrating a system of stiff differential equations utilizing concepts from singular perturbation theory. This paper evaluates the robustness and the reliability of such a singular perturbation based SUMT algorithm on two different problems of structural optimization of widely separated scales. The report concludes that whereas conventional SUMT can be bogged down by frequent ill-conditioning, especially in large scale problems, the singular perturbation SUMT has no such difficulty in converging to very accurate solutions.
Mickens (1981, 1984) has considered analytic techniques for obtaining approximate solutions to one-dimensional nonlinear oscillatory systems x(double-dot) + x = lambda f(x, x/dot/, lambda) where lambda is a small positive parameter and f is a nonlinear polynomial function of its arguments. However, in certain cases there is interest in the analysis of physical systems for which the nonlinear function f(x, x/dot/, lambda) is singular for finite values of x or x(dot). The present investigation is concerned with the use of existing approximate analytic schemes to obtain solutions to singular nonlinear oscillatory differential equations.
A new method for solving problems in relative astrometry is proposed. In it, the relationship between the measured quantities and the components of the position vector of a star is modeled using quaternions, in effect replacing the plate constants of a standard four-plate-constant solution with the four components of a quaternion. The method allows a direct solution for the position vectors of the stars, and hence for the equatorial coordinates. Distortions, magnitude, and color effects are readily incorporated into the formalism, and the method is directly applicable to overlapping-plate problems. The advantages of the method include the simplicity of the resulting equations, their freedom from singularities, and the fact that trigonometric functions and tangential point transformations are not needed to model the plate material. A global solution over the entire sky is possible.
The solution of bifurcation problems with closely-spaced critical points is achieved by first separating the singular part of the equation system encountered during a Newton iteration and carrying the Taylor expansion of the reduced system out to higher order. This separation is accomplished by transforming the equation system into an equivalent system in which some of the original unknowns are replaced with an equal number of modal amplitude coefficients. This method was used to continue the analysis of two significant example problems well past multiple fiburcation points, allowing a detailed examination of postbuckling behavior in the presence of modal interaction.