Determination of optimal trajectory under design constraints for a satellite launch vehicle

An optimal pitch steering programme of a solid-fuel satellite launch vehicle to maximize either (1) the injection velocity at a given altitude, or (2) the size of circular orbit, for a given payload is presented. The two-dimensional model includes the rotation of atmosphere with the Earth, the vehicle's lift and drag, variation of thrust with time and altitude, inverse-square gravitational field, and the specified initial vertical take-off. The inequality constraints on the aerodynamic load, control force, and turning rates are also imposed. Using the properties of the central force motion the terminal constraint conditions at coast apogee are transferred to the penultimate stage burnout. Such a transformation converts a time-free problem into a time-fixed one, reduces the number of terminal constraints, improves accuracy, besides demanding less computer memory and time. The adjoint equations are developed in a compact matrix form. The problem is solved on an IBM 360/44 computer using a steepest ascent algorithm. An illustrative analysis of a typical launch vehicle establishes the speed of convergence, and accuracy and applicability of the algorithm.

Design and Evaluation of a Robust Optical Beam-Interruption-Based Vehicle Classifier System

This paper presents the design and development of a novel optical vehicle classifier system, which is based on interruption of laser beams, that is suitable for use in places with poor transportation infrastructure. The system can estimate the speed, axle count, wheelbase, tire diameter, and the lane of motion of a vehicle. The design of the system eliminates the need for careful optical alignment, whereas the proposed estimation strategies render the estimates insensitive to angular mounting errors and to unevenness of the road. Strategies to estimate vehicular parameters are described along with the optimization of the geometry of the system to minimize estimation errors due to quantization. The system is subsequently fabricated, and the proposed features of the system are experimentally demonstrated. The relative errors in the estimation of velocity and tire diameter are shown to be within 0.5% and to change by less than 17% for angular mounting errors up to 30 degrees. In the field, the classifier demonstrates accuracy better than 97.5% and 94%, respectively, in the estimation of the wheelbase and lane of motion and can classify vehicles with an average accuracy of over 89.5%.

Robust Reentry Guidance of a Reusable Launch Vehicle Using Model Predictive Static Programming

A robust suboptimal reentry guidance scheme is presented for a reusable launch vehicle using the recently developed, computationally efficient model predictive static programming. The formulation uses the nonlinear vehicle dynamics with a spherical and rotating Earth, hard constraints for desired terminal conditions, and an innovative cost function having several components with associated weighting factors that can account for path and control constraints in a soft constraint manner, thereby leading to smooth solutions of the guidance parameters. The proposed guidance essentially shapes the trajectory of the vehicle by computing the necessary angle of attack and bank angle that the vehicle should execute. The path constraints are the structural load constraint, thermal load constraint, bounds on the angle of attack, and bounds on the bank angle. In addition, the terminal constraints include the three-dimensional position and velocity vector components at the end of the reentry. Whereas the angle-of-attack command is generated directly, the bank angle command is generated by first generating the required heading angle history and then using it in a dynamic inversion loop considering the heading angle dynamics. Such a two-loop synthesis of bank angle leads to better management of the vehicle trajectory and avoids mathematical complexity as well. Moreover, all bank angle maneuvers have been confined to the middle of the trajectory and the vehicle ends the reentry segment with near-zero bank angle, which is quite desirable. It has also been demonstrated that the proposed guidance has sufficient robustness for state perturbations as well as parametric uncertainties in the model.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.

Fast decay of the velocity autocorrelation function in dense shear flow of inelastic hard spheres

We find in complementary experiments and event-driven simulations of sheared inelastic hard spheres that the velocity autocorrelation function psi(t) decays much faster than t(-3/2) obtained for a fluid of elastic spheres at equilibrium. Particle displacements are measured in experiments inside a gravity-driven flow sheared by a rough wall. The average packing fraction obtained in the experiments is 0.59, and the packing fraction in the simulations is varied between 0.5 and 0.59. The motion is observed to be diffusive over long times except in experiments where there is layering of particles parallel to boundaries, and diffusion is inhibited between layers. Regardless, a rapid decay of psi(t) is observed, indicating that this is a feature of the sheared dissipative fluid, and is independent of the details of the relative particle arrangements. An important implication of our study is that the non-analytic contribution to the shear stress may not be present in a sheared inelastic fluid, leading to a wider range of applicability of kinetic theory approaches to dense granular matter.

On the relationship between finger width, velocity, and fluxes in thermohaline convection

Double-diffusive finger convection occurs in many natural processes.The theories for double-diffusive phenomena that exist at present consider systems with linear stratification in temperature and salinity. The double-diffusive systems with step change in salinity and temperature are, however, not amenable to simple stability analysis. Hence factors that control the width of the finger, velocity, and fluxes in systems that have step change in temperature and salinity have not been understood so far. In this paper we provide new physical insight regarding factors that influence finger convection in two-layer double-diffusive system through two-dimensional numerical simulations. Simulations have been carried out for density stability ratios (R-rho) from 1.5 to 10. For each density stability ratio, the thermal Rayleigh number (Ra-T) has been systematically varied from 7x10(3) to 7x10(8). Results from these simulations show how finger width, velocity, and flux ratios in finger convection are interrelated and the influence of governing parameters such as density stability ratio and the thermal Rayleigh number. The width of the incipient fingers at the time of onset of instability has been shown to vary as Ra-T-1/3. Velocity in the finger varies as Ra(T)1/3/R-rho. Results from simulation agree with the scale analysis presented in the paper. Our results demonstrate that wide fingers have lower velocities and flux ratios compared to those in narrow fingers. This result contradicts present notions about the relation between finger width and flux ratio. A counterflow heat-exchanger analogy is used in understanding the dependence of flux ratio on finger width and velocity.

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.

Practical Algorithms For Mean Velocity Estimation In Pulse Doppler Weather Radars Using A Small Number Of Samples

Doppler weather radars with fast scanning rates must estimate spectral moments based on a small number of echo samples. This paper concerns the estimation of mean Doppler velocity in a coherent radar using a short complex time series. Specific results are presented based on 16 samples. A wide range of signal-to-noise ratios are considered, and attention is given to ease of implementation. It is shown that FFT estimators fare poorly in low SNR and/or high spectrum-width situations. Several variants of a vector pulse-pair processor are postulated and an algorithm is developed for the resolution of phase angle ambiguity. This processor is found to be better than conventional processors at very low SNR values. A feasible approximation to the maximum entropy estimator is derived as well as a technique utilizing the maximization of the periodogram. It is found that a vector pulse-pair processor operating with four lags for clear air observation and a single lag (pulse-pair mode) for storm observation may be a good way to estimate Doppler velocities over the entire gamut of weather phenomena.

Bed Shear from Velocity Profiles: A New Approach

A new binary law of velocity distribution has been developed to describe the velocity profile for the entire flow region. The law is a combination of logarithmic law, valid in the wall (inner) region, and parabolic law, valid in the core (outer) region of the flow. The validity of the law has been established based on earlier data on flat plates, rough and smooth pipes and experimental data obtained from rigid-walled open channels with plane sand beds. A procedure of estimating bed shear stress from the proposed law of velocity distribution using the measured velocity profile has been evolved. Bed shear estimates made according to this procedure are in agreement with the values obtained from uniform flow analysis in the case of open channel flow over a sediment bed. The proposed method of estimating the bed shear stress from the observed velocity profiles is found to be particularly useful in cases where it is difficult to determine precisely the true bed level, such as in the case of flow over sediment beds.

The effects of sliding velocity and sliding time on nanocrystalline tribolayer development and properties in copper

The growth of the nanocrystalline tribolayer produced in oxygen free high conductivity copper after sliding against 440C stainless steel was studied. Tests were conducted on a pin-on-disk tribometer at sliding velocities of 0.05 and 1.0 m/s and sliding times of 0.1 to 10,000 s. Subsurface deformation and the growth of the tribolayer as a function of time were studied with the use of transmission electron microscopy and ion induced secondary electron microscopy. A continuous nanocrystalline tribolayer was produced after as little as 10 s of sliding at both sliding velocities. The tribolayer produced by sliding at 0.05 m/s continued to grow at sliding times up to 10,000 s and developed texture. Dynamic recrystallization of the tribolayer at a sliding velocity of 1.0 m/s inhibited the growth of a continuous anocrystalline tribolayer.

A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

The effects of slip velocity at a membrane surface on blood flow in the microcirculation

Closed-form solutions are presented for blood flow in the microcirculation by taking into account the influence of slip velocity at the membrane surface. In this study, the convective inertia force is neglected in comparison with that of blood viscosity on the basis of the smallness of the Reynolds number of the flow in microcirculation. The permeability property of the blood vessel is based on the well known Starling's hypothesis [11]. The effects of slip coefficient on the velocity and pressure fields are clearly depicted.

An Optimal g-Guidance Scheme for Satellite Launch Vehicles

A closed-loop steering logic based on an optimal (2-guidance is developed here. The guidance system drives the satellite launch vehicle along a two- or three- dimensional trajectory for placing the payload into a specified circular orbit. The modified g-guidance algorithm makes use of the optimal required velocity vector, which minimizes the total impulse needed for an equivalent two-impluse transfer from the present state to the final orbit. The required velocity vector is defined as velocity of the vehicle on the hypothetical transfer orbit immediately after the application of the first impulse. For this optimal transfer orbit, a simple and elegant expression for the Q-matrix is derived. A working principle for the guidance algorithm in terms of the major and minor cycles, and also for the generation of the steering command, is outlined.

Velocity and acceleration analysis of complex mechanisms by graphical iteration

A simple graphical method is presented for velocity and acceleration analysis of complex mechanisms possessing low or high degree of complexity. The method is iterative in character and generally yields the solution within a few iterations. Several examples have been worked out to illustrate the method.