Implicit and multigrid procedures for steady-state computations with upwind algorithms

A computational study for the convergence acceleration of Euler and Navier-Stokes computations with upwind schemes has been conducted in a unified framework. It involves the flux-vector splitting algorithms due to Steger-Warming and Van Leer, the flux-difference splitting algorithms due to Roe and Osher and the hybrid algorithms, AUSM (Advection Upstream Splitting Method) and HUS (Hybrid Upwind Splitting). Implicit time integration with line Gauss-Seidel relaxation and multigrid are among the procedures which have been systematically investigated on an individual as well as cumulative basis. The upwind schemes have been tested in various implicit-explicit operator combinations such that the optimal among them can be determined based on extensive computations for two-dimensional flows in subsonic, transonic, supersonic and hypersonic flow regimes. In this study, the performance of these implicit time-integration procedures has been systematically compared with those corresponding to a multigrid accelerated explicit Runge-Kutta method. It has been demonstrated that a multigrid method employed in conjunction with an implicit time-integration scheme yields distinctly superior convergence as compared to those associated with either of the acceleration procedures provided that effective smoothers, which have been identified in this investigation, are prescribed in the implicit operator.

Electrical Characterization of Airborne Vehicle Exhaust Plume

he induced current and voltage on the skin of an airborne vehicle due to the coupling of external electromagnetic field could be altered in the presence of ionized exhaust plume. So in the present work, a theoretical analysis is done to estimate the electrical parameters such as electrical conductivity and permittivity and their distribution in the axial and radial directions of the exhaust plume of an airborne vehicle. The electrical conductivity depends on the distribution of the major ionic species produced from the propellant combustion. In addition it also depends on temperature and pressure distribution of the exhaust plume as well as the generated shock wave. The chemically reactive rocket exhaust flow is modeled in two stages. The first part is simulated from the combustion chamber to the throat of the supersonic nozzle by using NASA Chemical Equilibrium with Application (CEA) package and the second part is simulated from the nozzle throat to the downstream of the plume by using a commercial Computational Fluid Dynamics (CFD) solver. The contour plots of the exhaust parameters are presented. Eight barrel shocks which influence the distribution of the vehicle exhaust parameters are obtained in this simulation. The computed peak value of the electrical conductivity of the plume is 0.123 S/m and the relative permittivity varies from 0.89 to 0.99. The attenuation of the microwave when it is passing through the conducting exhaust plume has also been presented.

Force Measurements on a Reflex Cambered Delta Wing

AREFLEX spanwise cambered delta wing with a conical camber designed for M= 1.4, using the method of Ref. 1, was tested at the design Mach number as well as off-design Mach number M=0.15 and 2.3, respectively. The test results are compared with those of a plane wing and also with the available theoretical results at the design condition. At subsonic speed, the cambered wing has less lift at a given incidence and higher lift-to-drag ratio at a given lift than the plane wing, while at supersonic speeds, both of these quantities were less on the cambered wing. At supersonic speed, at the design incidence and Mach number, there is good agreement between results from theory and experiment. The center of pressure on the cambered wing is ahead of that on the plane wing at subsonic speed, while the reverse is true at supersonic speeds. Finally, it is found that over a useful range of lift the cambered wing is aerodynamically more efficient at subsonic speeds, and less so at supersonic speeds, than the plane wing.

Approximation of the perturbation equations of a quasi-linear hyperbolic system in the neighborhood of a bicharacteristic

In 1956 Whitham gave a nonlinear theory for computing the intensity of an acoustic pulse of an arbitrary shape. The theory has been used very successfully in computing the intensity of the sonic bang produced by a supersonic plane. [4.] derived an approximate quasi-linear equation for the propagation of a short wave in a compressible medium. These two methods are essentially nonlinear approximations of the perturbation equations of the system of gas-dynamic equations in the neighborhood of a bicharacteristic curve (or rays) for weak unsteady disturbances superimposed on a given steady solution. In this paper we have derived an approximate quasi-linear equation which is an approximation of perturbation equations in the neighborhood of a bicharacteristic curve for a weak pulse governed by a general system of first order quasi-linear partial differential equations in m + 1 independent variables (t, x1,…, xm) and derived Gubkin's result as a particular case when the system of equations consists of the equations of an unsteady motion of a compressible gas. We have also discussed the form of the approximate equation describing the waves propagating upsteam in an arbitrary multidimensional transonic flow.

Nonlinear Control Design for an Air-breathing Engine with State Estimation

A nonlinear control design approach is presented in this paper for a challenging application problem of ensuring robust performance of an air-breathing engine operating at supersonic speed. The primary objective of control design is to ensure that the engine produces the required thrust that tracks the commanded thrust as closely as possible by appropriate regulation of the fuel flow rate. However, since the engine operates in the supersonic range, an important secondary objective is to ensure an optimal location of the shock in the intake for maximum pressure recovery with a sufficient margin. This is manipulated by varying the throat area of the nozzle. The nonlinear dynamic inversion technique has been successfully used to achieve both of the above objectives. In this problem, since the process is faster than the actuators, independent control designs have also been carried out for the actuators as well to assure the satisfactory performance of the system. Moreover, an extended Kalman Filter based state estimation design has been carried out both to filter out the process and sensor noises as well as to make the control design operate based on output feedback. Promising simulation results indicate that the proposed control design approach is quite successful in obtaining robust performance of the air-breathing system.

Theoretical gain optimization studies in 10.6 µm CO2-N2 gasdynamic lasers. IV. Further results of parametric study

Based on a method proposed by Reddy and Shanmugasundaram, similar solutions have been obtained for the steady inviscid quasi-one-dimensional nonreacting flow in the supersonic nozzle of CO2-N2-H2O and CO2-N2-He gasdynamic laser systems. Instead of using the correlations of a nonsimilar function NS for pure N2 gas, as is done in previous publications, the NS correlations are computed here for the actual gas mixtures used in the gasdynamic lasers. Optimum small-signal optical gain and the corresponding optimum values of the operating parameters like reservoir pressure and temperature and nozzle area ratio are computed using these correlations. The present results are compared with the previous results and the main differences are discussed. Journal of Applied Physics is copyrighted by The American Institute of Physics.

The hydrogen bond: a molecular beam microwave spectroscopist's view with a universal appeal

In this manuscript, we propose a criterion for a weakly bound complex formed in a supersonic beam to be characterized as a `hydrogen bonded complex'. For a `hydrogen bonded complex', the zero point energy along any large amplitude vibrational coordinate that destroys the orientational preference for the hydrogen bond should be significantly below the barrier along that coordinate so that there is at least one bound level. These are vibrational modes that do not lead to the breakdown of the complex as a whole. If the zero point level is higher than the barrier, the `hydrogen bond' would not be able to stabilize the orientation which favors it and it is no longer sensible to characterize a complex as hydrogen bonded. Four complexes, Ar-2-H2O, Ar-2-H2S, C2H4-H2O and C2H4-H2S, were chosen for investigations. Zero point energies and barriers for large amplitude motions were calculated at a reasonable level of calculation, MP2(full)/aug-cc-pVTZ, for all these complexes. Atoms in molecules (AIM) theoretical analyses of these complexes were carried out as well. All these complexes would be considered hydrogen bonded according to the AIM theoretical criteria suggested by Koch and Popelier for C-H center dot center dot center dot O hydrogen bonds (U. Koch and P. L. A. Popelier, J. Phys. Chem., 1995, 99, 9747), which has been widely and, at times, incorrectly used for all types of contacts involving H. It is shown that, according to the criterion proposed here, the Ar-2-H2O/H2S complexes are not hydrogen bonded even at zero kelvin and C2H4-H2O/H2S complexes are. This analysis can naturally be extended to all temperatures. It can explain the recent experimental observations on crystal structures of H2S at various conditions and the crossed beam scattering studies on rare gases with H2O and H2S.

Force-free equilibria of magnetized jets

Force-free equilibrium configurations of magnetic-pressure-dominated magnetized supersonic jets confined by slowly varying external pressure are investigated analytically. For the case where internal dissipation mechanisms are active, the lowest-energy field configuration is found to be the superposition of an axisymmetric mode and a helical mode with a wavelength equal to 5 times the jet radius, and the pressure below which the nonaxisymmetric mode becomes energetically favorable is given as 2700 times the product of the 4th power of the magnetic helicity per unit length and the -6th power of the magnetic flux. A model of the total and polarized emission of such a configuration is developed and applied to the extended well-collimated astronomically resolved jet NGC 6251. The model is shown to reproduce significant features such as transverse oscillations of the ridge line, width oscillations and emission knots, the projected magnetic-field configuration, oscillations of the degree of polarization, and the distribution of the Faraday rotation measure.

Magnetic energy dissipation in force-free jets

It is shown that a magnetic-pressure-dominated, supersonic jet which expands (or contracts) in response to variations in the confining external pressure can dissipate magnetic energy through field-line reconnection as it relaxes to a minimum-energy configuration. In order for a continuous dissipation to take place, the effective reconnection time must be a fraction ɛ ⪉ 1 of the expansion time. The amount of energy dissipation is calculated, and it is concluded that magnetic energy dissipation could, in principle, power the observed synchrotron emission in extragalactic radio jets such as NGC 6251. However, this mechanism is only viable if the reconnection time is substantially shorter than the nominal resistive tearing time in the jet.

Nonlinear Sliding Mode Control Design for an Air-breathing Engine with State Estimation

A modern system theory based nonlinear control design is discussed in this paper for successful operation of an air-breathing engine operating at supersonic speed. The primary objective of the control design of such an air-breathing engine is to ensure that the engine dynamically produces the thrust that tracks a commanded value of thrust as closely as possible by regulating the fuel flow to the combustion system. However, since the engine operates in the supersonic range, an important secondary objective is to manage the shock wave configuration in the intake section of the engine which is manipulated by varying the throat area of the nozzle. A nonlinear sliding mode control technique has been successfully used to achieve both of the above objectives. In this problem, since the process is faster than the actuators, independent control designs are also carried out for the actuators as well to assure the satisfactory performance of the system. Moreover, to filter out the sensor and process noises and to estimate the states for making the control design operate based on output feedback, an Extended Kalman Filter based state estimation design is also carried out. The promising simulation results suggest that the proposed control design approach is quite successful in obtaining robust performance of the air-breathing engine.

Accelerometer-Based Force Balance for High Enthalpy Facilities

A single component accelerometer-based force balance is developed, calibrated, and used for high enthalpy applications. Functionality of this force balance, for such applications, is demonstrated for the first time during high enthalpy tests in a newly established free piston-driven shock tunnel, HST3, using a 60 degrees apex angle blunt cone model at 0 degrees angle of incidence. Usefulness of this force balance is also confirmed, for much complicated high enthalpy flow situations, during the drag reduction studies with counterflow supersonic jet from the stagnation point.

Theoretical gain‐optimization studies in CO2–N2 gasdynamic lasers. II. Results of parametric study

Based on a method presented in detail in a previous work by the authors, similar solutions have been obtained for the steady inviscid quasi‐one‐dimensional nonreacting flow in the supersonic nozzle of a CO2–N2 gasdynamic laser system, with either H2O or He as the catalyst. It has been demonstrated how these solutions could be used to optimize the small‐signal gain coefficient on a specified vibrational‐rotational transition. Results presented for a wide range of mixture compositions include optimum values for the small‐signal gain, area ratio, reservoir temperature, and a binary scaling parameter, which is the product of reservoir pressure and nozzle shape factor.

Theoretical gain‐optimization studies in CO2–N2 gasdynamic lasers. I. Theory

Based on a method proposed by Reddy and Daum, the equations governing the steady inviscid nonreacting gasdynamic laser (GDL) flow in a supersonic nozzle are reduced to a universal form so that the solutions depend on a single parameter which combines all the other parameters of the problem. Solutions are obtained for a sample case of available data and compared with existing results to validate the present approach. Also, similar solutions for a sample case are presented.

Linear stability, transient growth and the role of viscosity stratification in compressible Couette flow

Linear stability and the nonmodal transient energy growth in compressible plane Couette flow are investigated for two prototype mean flows: (a) the uniform shear flow with constant viscosity, and (b) the nonuniform shear flow with stratified viscosity. Both mean flows are linearly unstable for a range of supersonic Mach numbers (M). For a given M, the critical Reynolds number (Re) is significantly smaller for the uniform shear flow than its nonuniform shear counterpart; for a given Re, the dominant instability (over all streamwise wave numbers, α) of each mean flow belongs to different modes for a range of supersonic M. An analysis of perturbation energy reveals that the instability is primarily caused by an excess transfer of energy from mean flow to perturbations. It is shown that the energy transfer from mean flow occurs close to the moving top wall for “mode I” instability, whereas it occurs in the bulk of the flow domain for “mode II.” For the nonmodal transient growth analysis, it is shown that the maximum temporal amplification of perturbation energy, Gmax, and the corresponding time scale are significantly larger for the uniform shear case compared to those for its nonuniform counterpart. For α=0, the linear stability operator can be partitioned into L∼L̅ +Re2 Lp, and the Re-dependent operator Lp is shown to have a negligibly small contribution to perturbation energy which is responsible for the validity of the well-known quadratic-scaling law in uniform shear flow: G(t∕Re)∼Re2. In contrast, the dominance of Lp is responsible for the invalidity of this scaling law in nonuniform shear flow. An inviscid reduced model, based on Ellingsen-Palm-type solution, has been shown to capture all salient features of transient energy growth of full viscous problem. For both modal and nonmodal instability, it is shown that the viscosity stratification of the underlying mean flow would lead to a delayed transition in compressible Couette flow.

CFD Based Design Optimizations of the Diffuser of a Gas Dynamically Driver Laser Cavity

Based on the an earlier CFD analysis of the performance of the gas-dynamically controlled laser cavity [1]it was found that there is possibility of optimizing the geometry of the diffuser that can bring about reductions in both size and cost of the system by examining the critical dimensional requirements of the diffuser. Consequently,an extensive CFD analysis has been carried out for a range of diffuser configurations by simulating the supersonic flow through the arrangement including the laser cavity driven by a bank of converging – diverging nozzles and the diffuser. The numerical investigations with 3D-RANS code are carried out to capture the flow patterns through diffusers past the cavity that has multiple supersonic jet interactions with shocks leading to complex flow pattern. Varying length of the diffuser plates is made to be the basic parameter of the study. The analysis reveals that the pressure recovery pattern during the flow through the diffuser from the simulation, being critical for the performance of the laser device shows its dependence on the diffuser length is weaker beyond a critical lower limit and this evaluation of this limit would provide a design guideline for a more efficient system configuration.The observation based on the parametric study shows that the pressure recovery transients in the near vicinity of the cavity is not affected for the reduction in the length of the diffuser plates up to its 10% of the initial size, indicating the design in the first configuration that was tested experimentally has a large factor of margin. The flow stability in the laser cavity is found to be unaffected since a strong and stable shock is located at the leading edge of the diffuser plates while the downstream shock and flow patterns are changed, as one would expect. Results of the study for the different lengths of diffusers in the range of 10% to its full length are presented, keeping the experimentally tested configuration used in the earlier study [1] as the reference length. The conclusions drawn from the analysis is found to be of significance since it provides new design considerations based on the understanding of the intricacies of the flow, allowing for a hardware optimization that can lead to substantial size reduction of the device with no loss of performance.