Boundary-layer effects on internal flow choking in dual-thrust solid rocket motors

Theoretical studies have been carried out to examine internal flow choking in the inert simulators of a dual-thrust motor. Using a two-dimensional k-omega turbulence model, detailed parametric studies have been carried out to examine aerodynamic choking and the existence of a fluid throat at the transition region during the startup transient of dual-thrust motors. This code solves standard k-omega turbulence equations with shear flow corrections using a coupled second-order-implicit unsteady formulation. In the numerical study, a fully implicit finite volume scheme of the compressible, Reynolds-averaged, Navier-Stokes equations is employed. It was observed that, at the subsonic inflow conditions, there is a possibility of the occurrence of internal flow choking in dual-thrust motors due to the formation of a fluid throat at the beginning of the transition region induced by area blockage caused by boundary-layer-displacement thickness. It has been observed that a 55% increase in the upstream port area of the dual-thrust motor contributes to a 25% reduction in blockage factor at the transition region, which could negate the internal How choking and supplement with an early choking of the dual-thrust motor nozzle. If the height of the upstream port relative to the motor length is too small, the developing boundary layers from either side of the port can interact, leading to a choked,flow. On the other hand, if the developing boundary layers are far enough apart, then choking does not occur. The blockage factor is greater in magnitude for the choked case than for the unchoked case. More tangible explanations are presented in this paper for the boundary-layer blockage and the internal flow choking in dual-thrust motors, which hitherto has been unexplored.

Energy identities in water wave theory for free-surface boundary condition with higher-order derivatives

A modified form of Green's integral theorem is employed to derive the energy identity in any water wave diffraction problem in a single-layer fluid for free-surface boundary condition with higher-order derivatives. For a two-layer fluid with free-surface boundary condition involving higher-order derivatives, two forms of energy identities involving transmission and reflection coefficients for any wave diffraction problem are also derived here by the same method. Based on this modified Green's theorem, hydrodynamic relations such as the energy-conservation principle and modified Haskind–Hanaoka relation are derived for radiation and diffraction problems in a single as well as two-layer fluid.

Grain boundary sliding during diffusion and dislocation creep in a Mg-0.7 pct Al alloy

Early studies on grain boundary sliding (GBS) in Mg alloys have suggested frequently that the contribution of GBS to creep is high even under conditions corresponding to dislocation creep. The role of creep strain and grain size in influencing the experimental measurements has not been clearly identified. Grain boundary sliding measurements were conducted in detail over experimental conditions corresponding to diffusion creep as well as dislocation creep in a single-phase Mg-0.7 wt pet Al alloy. The results indicated clearly that the GBS contribution to creep was Very high during,, diffusion creep at low stresses (similar to 75 pct) and substantially reduced during dislocation creep at high stresses (similar to 15 pct). These measurements were consistent with the observation of significant intragranular slip band activity observed in most grains at high stresses and very little slip band activity at low stresses. The experimental measurements and analysis indicated also that the GBS contribution to creep was high during the initial stages of creep and decreased to a steady-state value at large strains.

Evolution of isolated streamwise vortices in the late stages of boundary layer transition

Experiments were conducted with two, smooth hills, lying well within the boundary layer over a flat plate mounted in a wind tunnel. One hill was shallow, with peak height 1.5 mm and width 50 mm; the other, steep, 3 mm high and 30 mm wide. Since the hills occupied one-half of the tunnel span, streamwise vorticity formed near the hills' edge. At a freestream speed of 3.5 m/s, streaks formed with inflectional wall-normal and spanwise velocity profiles but without effecting transition. Transition, observed at 7.5 m/s, took different routes with the two hills. With the steep hill, streamwise velocity signals exhibited the passage of a wave packet which intensified before breakdown to turbulence. With the shallow hill there was a broad range of frequencies present immediately downstream of the hill. These fluctuations grew continuously and transition occurred within a shorter distance. Since the size of the streamwise vorticity generated at the hill edge is of the order of the hill height, the shallow hill generates vorticity closer to the wall and supports an earlier transition, whereas the steep hill creates a thicker vortex and associated streaks which exhibit oscillations due to their own instability as an additional precursor stage before transition.

Quantitative photoacoustic tomography from boundary pressure measurements: noniterative recovery of optical absorption coefficient from the reconstructed absorbed energy map

We describe a noniterative method for recovering optical absorption coefficient distribution from the absorbed energy map reconstructed using simulated and noisy boundary pressure measurements. The source reconstruction problem is first solved for the absorbed energy map corresponding to single- and multiple-source illuminations from the side of the imaging plane. It is shown that the absorbed energy map and the absorption coefficient distribution, recovered from the single-source illumination with a large variation in photon flux distribution, have signal-to-noise ratios comparable to those of the reconstructed parameters from a more uniform photon density distribution corresponding to multiple-source illuminations. The absorbed energy map is input as absorption coefficient times photon flux in the time-independent diffusion equation (DE) governing photon transport to recover the photon flux in a single step. The recovered photon flux is used to compute the optical absorption coefficient distribution from the absorbed energy map. In the absence of experimental data, we obtain the boundary measurements through Monte Carlo simulations, and we attempt to address the possible limitations of the DE model in the overall reconstruction procedure.

Boundary value problems for third-order nonlinear ordinary differential equations

In this paper, we describe how to analyze boundary value problems for third-order nonlinear ordinary differential equations over an infinite interval. Several physical problems of interest are governed by such systems. The seminumerical schemes described here offer some advantages over solutions obtained by using traditional methods such as finite differences, shooting method, etc. These techniques also reveal the analytic structure of the solution function. For illustrative purposes, several physical problems, mainly drawn from fluid mechanics, are considered; they clearly demonstrate the efficiency of the techniques presented here.

Boundary value problems for third-order nonlinear ordinary differential equations

In this paper, we describe how to analyze boundary value problems for third-order nonlinear ordinary differential equations over an infinite interval. Several physical problems of interest are governed by such systems. The seminumerical schemes described here offer some advantages over solutions obtained by using traditional methods such as finite differences, shooting method, etc. These techniques also reveal the analytic structure of the solution function. For illustrative purposes, several physical problems, mainly drawn from fluid mechanics, are considered; they clearly demonstrate the efficiency of the techniques presented here.

Effect of boundary condition description on convergence of solution in a boundary-value problem

It is well known that the numerical accuracy of a series solution to a boundary-value problem by the direct method depends on the technique of approximate satisfaction of the boundary conditions and on the stage of truncation of the series. On the other hand, it does not appear to be generally recognized that, when the boundary conditions can be described in alternative equivalent forms, the convergence of the solution is significantly affected by the actual form in which they are stated. The importance of the last aspect is studied for three different techniques of computing the deflections of simply supported regular polygonal plates under uniform pressure. It is also shown that it is sometimes possible to modify the technique of analysis to make the accuracy independent of the description of the boundary conditions.

The influence of Zr/Ti content on the morphotropic phase boundary in the PZT-PZN system

A quantitative structural investigation was carried out on (1-y)PbZrxTi1-xO3-yPbZn(1/3)Nb(2/3)O(3) where y=0.1 and 0.2 ((1-y)PZT-yPZN). High resolution XRD data have been used for quantitative phase analysis. The nominal compositions were prepared by a two-step low temperature calcining solid-state method. The sintered samples show an average grain size of 1-2 mu m. It is demonstrated that the increase in the concentration of PZN leads to the shift of the morphotropic phase boundary (MPB) of PZT towards the PbZrO3 end member. In the present work, an effort has been made to quantitatively determine the MPB phase contents and to regain the coexistence of tetragonal and monoclinic phases by varying the value of x(i.e. Zr/Ti ratio). The width of the MPB becomes considerably larger for y=0.10 and 0.20 as compared to pure PZT. This is attributed to the considerably lower grain size of our samples resulting from the adopted preparation method. (C) 2010 Elsevier B.V. All rights reserved.

Compressible axially symmetric laminar boundary layer flow in the stagnation region of a blunt body in the presence of magnetic field

In this paper, the steady laminar viscous hypersonic flow of an electrically conducting fluid in the region of the stagnation point of an insulating blunt body in the presence of a radial magnetic field is studied by similarity solution approach, taking into account the variation of the product of density and viscosity across the boundary layer. The two coupled non-linear ordinary differential equations are solved simultaneously using Runge-Kutta-Gill method. It has been found that the effect of the variation of the product of density and viscosity on skin friction coefficient and Nusselt number is appreciable. The skin friction coefficient increases but Nusselt number decreases as the magnetic field or the total enthalpy at the wall increases