An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations is given, specifically for the steady Euler equations. Solutions of the equations were obtained by Newton's linearization procedure, commonly used to solve the roots of nonlinear algebraic equations. In application of the same procedure for solving a set of differential equations we give a theorem showing that a quadratic convergence rate can be achieved. While the domain of quadratic convergence depends on the problems studied and is unknown a priori, we show that firstand second-order derivatives of flux vectors determine whether the condition for quadratic convergence is satisfied. The first derivatives enter as an implicit operator for yielding new iterates and the second derivatives indicates smoothness of the flows considered. Consequently flows involving shocks are expected to require larger number of iterations. First-order upwind discretization in conjunction with the Steger-Warming flux-vector splitting is employed on the implicit operator and a diagonal dominant matrix results. However the explicit operator is represented by first- and seond-order upwind differencings, using both Steger-Warming's and van Leer's splittings. We discuss treatment of boundary conditions and solution procedures for solving the resulting block matrix system. With a set of test problems for one- and two-dimensional flows, we show detailed study as to the efficiency, accuracy, and convergence of the present method

Comment on ``Validity of certain soft-photon amplitudes''

The criteria suggested by Welsh and Fearing (nucl-th/9606040) to judge the validity of certain soft-photon amplitudes are examined. We comment on aspects of their analysis which lead to incorrect conclusions about published amplitudes and point out important criteria which were omitted from their analysis.Comment: 6 pages plus 1 postscript figure, Revte

Interaction between a normal shock wave and a turbulent boundary layer at high transonic speeds. Part 2: Wall shear stress

An analysis is presented of the flow in the two inner layers, the Reynolds stress sublayer and the wall layer. Included is the calculation of the shear stress at the wall in the interaction region. The limit processes considered are those used for an inviscid flow

Navier-Stokes analysis of transonic cascade flow

A new kind of C-type grid is proposed, this grid is non-periodic on the wake and allows minimum skewness for cascades with high turning and large camber. Reynolds-averaged Navier-Stokes equations are solved on this type of grid using a finite volume discretization and a full multigrid method which uses Runge-Kutta stepping as the driving scheme. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. A detailed numerical study is proposed for a highly loaded transonic blade. A grid independence analysis is presented in terms of pressure distribution, exit flow angles, and loss coefficient. Comparison with experiments clearly demonstrates the capability of the proposed procedure

Development of new flux splitting schemes

Maximizing both accuracy and efficiency has been the primary objective in designing a numerical algorithm for CFD. This is especially important for solution of complex three-dimensional systems of Navier-Stokes equations which often include turbulence modeling and chemistry effects. Recently, upwind schemes have been well received for both their capability of resolving discontinuities and their sound theoretical basis in characteristic theory for hyperbolic systems. With this in mind, two new flux splitting techniques are presented for upwind differencing

Migration of Interplanetary Dust

We numerically investigate the migration of dust particles with initial orbits close to those of the numbered asteroids, observed trans-Neptunian objects, and Comet Encke. The fraction of silicate asteroidal particles that collided with the Earth during their lifetime varied from 1.1% for 100 micron particles to 0.008% for 1 micron particles. Almost all asteroidal particles with diameter d>4 microns collided with the Sun. The peaks in the migrating asteroidal dust particles' semi-major axis distribution at the n:(n+1) resonances with Earth and Venus and the gaps associated with the 1:1 resonances with these planets are more pronounced for larger particles. The probability of collisions of cometary particles with the Earth is smaller than for asteroidal particles, and this difference is greater for larger particles.Comment: Annals of the New York Academy of Sciences, 15 pages, 8 Figures, submitte

Multigrid calculation of three-dimensional viscous cascade flows

A 3-D code for viscous cascade flow prediction was developed. The space discretization uses a cell-centered scheme with eigenvalue scaling to weigh the artificial dissipation terms. Computational efficiency of a four stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full multigrid method. The Baldwin-Lomax eddy viscosity model is used for turbulence closure. A zonal, nonperiodic grid is used to minimize mesh distortion in and downstream of the throat region. Applications are presented for an annular vane with and without end wall contouring, and for a large scale linear cascade. The calculation is validated by comparing with experiments and by studying grid dependency

Interaction between a normal shock wave and a turbulent boundary layer at high transonic speeds. Part 1: Pressure distribution. Part 2: Wall shear stress. Part 3: Simplified formulas for the prediction of surface pressures and skin friction

An asymptotic description is derived for the interaction between a shock wave and a turbulent boundary layer in transonic flow, for a particular limiting case. The dimensionless difference between the external flow velocity and critical sound speed is taken to be much smaller than one, but large in comparison with the dimensionless friction velocity. The basic results are derived for a flat plate, and corrections for longitudinal wall curvature and for flow in a circular pipe are also shown. Solutions are given for the wall pressure distribution and the shape of the shock wave. Solutions for the wall shear stress are obtained, and a criterion for incipient separation is derived. Simplified solutions for both the wall pressure and skin friction distributions in the interaction region are given. These results are presented in a form suitable for use in computer programs