Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three dimensional problems

A review of the ONR/NAVAIR research option combustion instabilities in compact ramjets, 1983-1988

This paper consists of two parts summarizing two portions of the ONR/NAVAIR Research Option. The option began in 1983 and continued for five years, involving 11 organizations. Simultaneously, similar or related programs supported by other agencies or institutions were being carried out in several other places. Results of those programs have been briefly summarized in five papers collected in a document to be published by C.P.L.A. This paper contains two of the five papers in that document. Here we cover the subjects of approximate analyses and stability; and large-scale structures and passive control. The first is concerned chiefly with an analytical framework constructed on the basis of observations; it is intended to provide a means of correlating and interpreting data, and predicting the stability of motions in a combustion chamber. The second is a summary of recent experimental work directed to understanding the flows in dump combustors of the sort used in modern ramjet engines. Much relevant material is not included here, but may be found in the remaining papers of the document cited above. For completeness, we note briefly the substance of those reports. In their summary "Spray Combustion Processes in Ramjet Combustion Instability," Bowman (Stanford), Law (University of California, Davis) and Sirignano (University of California, Irvine) review several aspects of spray combustion relevant to combustion instabilities. The objectives of the works were: (1) to determine the effect of spray characteristics on the energy release pattern in a dump combustor and the subsequent effects on combustion instability; (2) to gain a fundamental understanding of the coupling of the spray vaporization process with an unsteady flow field; and (3) to investigate methods for controlling and enhancing spray vaporization rates in liquid-fueled ramjets. During the past five years considerable progress has been made in applying methods of computational fluid dynamics to the flow in a dump combustor including consequences of energy release due to combustion processes. Jou has summarized work done at Flow Research, Inc. and at the Naval Research Laboratory in his paper "A Summary Report on Large-Eddy Simulations of Pressure Oscillations in a Ramjet Combustor." The serious effects of combustion instabilities on the inlets of ramjet engines were discovered in the late 1970's in experimental work at the Aeropropulsion Laboratory, Wright Field, the Naval Weapons Center and the Marquardt Company. The most thorough laboratory work on the unsteady behavior of inlets has been accomplished at the McDonnell-Douglas Research Laboratory by Sajben who has reviewed the subject in his paper "The Role of Inlet in Ramjet Pressure Oscillations.

Hydrodynamic Flows on Curved Surfaces: Spectral Numerical Methods for Radial Manifold Shapes

We formulate hydrodynamic equations and spectrally accurate numerical methods for investigating the role of geometry in flows within two-dimensional fluid interfaces. To achieve numerical approximations having high precision and level of symmetry for radial manifold shapes, we develop spectral Galerkin methods based on hyperinterpolation with Lebedev quadratures for $L^2$-projection to spherical harmonics. We demonstrate our methods by investigating hydrodynamic responses as the surface geometry is varied. Relative to the case of a sphere, we find significant changes can occur in the observed hydrodynamic flow responses as exhibited by quantitative and topological transitions in the structure of the flow. We present numerical results based on the Rayleigh-Dissipation principle to gain further insights into these flow responses. We investigate the roles played by the geometry especially concerning the positive and negative Gaussian curvature of the interface. We provide general approaches for taking geometric effects into account for investigations of hydrodynamic phenomena within curved fluid interfaces.Comment: 14 figure

On Time Optimization of Centroidal Momentum Dynamics

Recently, the centroidal momentum dynamics has received substantial attention to plan dynamically consistent motions for robots with arms and legs in multi-contact scenarios. However, it is also non convex which renders any optimization approach difficult and timing is usually kept fixed in most trajectory optimization techniques to not introduce additional non convexities to the problem. But this can limit the versatility of the algorithms. In our previous work, we proposed a convex relaxation of the problem that allowed to efficiently compute momentum trajectories and contact forces. However, our approach could not minimize a desired angular momentum objective which seriously limited its applicability. Noticing that the non-convexity introduced by the time variables is of similar nature as the centroidal dynamics one, we propose two convex relaxations to the problem based on trust regions and soft constraints. The resulting approaches can compute time-optimized dynamically consistent trajectories sufficiently fast to make the approach realtime capable. The performance of the algorithm is demonstrated in several multi-contact scenarios for a humanoid robot. In particular, we show that the proposed convex relaxation of the original problem finds solutions that are consistent with the original non-convex problem and illustrate how timing optimization allows to find motion plans that would be difficult to plan with fixed timing.Comment: 7 pages, 4 figures, ICRA 201

On the Singular Neumann Problem in Linear Elasticity

The Neumann problem of linear elasticity is singular with a kernel formed by the rigid motions of the body. There are several tricks that are commonly used to obtain a non-singular linear system. However, they often cause reduced accuracy or lead to poor convergence of the iterative solvers. In this paper, different well-posed formulations of the problem are studied through discretization by the finite element method, and preconditioning strategies based on operator preconditioning are discussed. For each formulation we derive preconditioners that are independent of the discretization parameter. Preconditioners that are robust with respect to the first Lam\'e constant are constructed for the pure displacement formulations, while a preconditioner that is robust in both Lam\'e constants is constructed for the mixed formulation. It is shown that, for convergence in the first Sobolev norm, it is crucial to respect the orthogonality constraint derived from the continuous problem. Based on this observation a modification to the conjugate gradient method is proposed that achieves optimal error convergence of the computed solution