State estimation of spiralling target using UKF with generalized sinusoidal target acceleration model

Hit-to-kill interception of high velocity spiraling target requires accurate state estimation of relative kinematic parameters describing spiralling motion. In this pa- per, spiraling target motion is captured by representing target acceleration through sinusoidal function in inertial frame. A nine state unscented Kalman filter (UKF) formulation is presented here with three relative positions, three relative velocities, spiraling frequency of target, inverse of ballistic coefficient and maneuvering coef-ficient. A key advantage of the target model presented here is that it is of generic nature and can capture spiraling as well as pure ballistic motions without any change of tuning parameters. Extensive Six-DOF simulation experiments, which includes a modified PN guidance and dynamic inversion based autopilot, show that near Hit-to-Kill performance can be obtained with noisy RF seeker measurements of gimbal angles, gimbal angle rates, range and range rate.

Guided-wave-based damage detection in a composite T-joint using 3D scanning laser Doppler vibrometer

Composite T-joints are commonly used in modern composite airframe, pressure vessels and piping structures, mainly to increase the bending strength of the joint and prevents buckling of plates and shells, and in multi-cell thin-walled structures. Here we report a detailed study on the propagation of guided ultrasonic wave modes in a composite T-joint and their interactions with delamination in the co-cured co-bonded flange. A well designed guiding path is employed wherein the waves undergo a two step mode conversion process, one is due to the web and joint filler on the back face of the flange and the other is due to the delamination edges close to underneath the accessible surface of the flange. A 3D Laser Doppler Vibrometer is used to obtain the three components of surface displacements/velocities of the accessible face of the flange of the T-joint. The waves are launched by a piezo ceramic wafer bonded on to the back surface of the flange. What is novel in the proposed method is that the location of any change in material/geometric properties can be traced by computing a frequency domain power flow along a scan line. The scan line can be chosen over a grid either during scan or during post-processing of the scan data off-line. The proposed technique eliminates the necessity of baseline data and disassembly of structure for structural interrogation.

Transmission of elastic waves through a frictional boundary

Incident energy gets transmitted, reflected and absorbed across an interface in jointed rock mass leading to energy dissipation and alteration of waves. Wave velocities get attenuated during their propagation across joints and this behavior is studied using bender/extender element tests. The velocity attenuation and modulus reduction observed in experimental tests are modeled with three dimensional distinct element code and results are validated. Normal propagation of an incident shear wave through a jointed rock mass cause slip of the rock blocks if shear stress of wave exceeds the shear strength of the joint. As the properties of joint determine the transmission of energy across an interface, a parametric study is then conducted with the validated numerical model by varying the parameters that may determine the energy transmission across a joint using modified Miller's method. Results of the parametric study are analyzed and presented in the paper. (C) 2014 Elsevier Ltd. All rights reserved.

Solving Axisymmetric Stability Problems by Using Upper Bound Finite Elements, Limit Analysis, and Linear Optimization

A numerical formulation has been proposed for solving an axisymmetric stability problem in geomechanics with upper bound limit analysis, finite elements, and linear optimization. The Drucker-Prager yield criterion is linearized by simulating a sphere with a circumscribed truncated icosahedron. The analysis considers only the velocities and plastic multiplier rates, not the stresses, as the basic unknowns. The formulation is simple to implement, and it has been employed for finding the collapse loads of a circular footing placed over the surface of a cohesive-frictional material. The formulation can be used to solve any general axisymmetric geomechanics stability problem.

Current kinematics and dynamics of Africa and the East African Rift System

Although the East African Rift System (EARS) is an archetype continental rift, the forces driving its evolution remain debated. Some contend buoyancy forces arising from gravitational potential energy (GPE) gradients within the lithosphere drive rifting. Others argue for a major role of the diverging mantle flow associated with the African Superplume. Here we quantify the forces driving present-day continental rifting in East Africa by (1) solving the depth averaged 3-D force balance equations for 3-D deviatoric stress associated with GPE, (2) inverting for a stress field boundary condition that we interpret as originating from large-scale mantle tractions, (3) calculating dynamic velocities due to lithospheric buoyancy forces, lateral viscosity variations, and velocity boundary conditions, and (4) calculating dynamic velocities that result from the stress response of horizontal mantle tractions acting on a viscous lithosphere in Africa and surroundings. We find deviatoric stress associated with lithospheric GPE gradients are similar to 8-20 MPa in EARS, and the minimum deviatoric stress resulting from basal shear is similar to 1.6 MPa along the EARS. Our dynamic velocity calculations confirm that a force contribution from GPE gradients alone is sufficient to drive Nubia-Somalia divergence and that additional forcing from horizontal mantle tractions overestimates surface kinematics. Stresses from GPE gradients appear sufficient to sustain present-day rifting in East Africa; however, they are lower than the vertically integrated strength of the lithosphere along most of the EARS. This indicates additional processes are required to initiate rupture of continental lithosphere, but once it is initiated, lithospheric buoyancy forces are enough to maintain rifting.

Turbulence and cooling in galaxy cluster cores

We study the interplay between turbulent heating, mixing, and radiative cooling in an idealized model of cool cluster cores. Active galactic nuclei (AGN) jets are expected to drive turbulence and heat cluster cores. Cooling of the intracluster medium (ICM) and stirring by AGN jets are tightly coupled in a feedback loop. We impose the feedback loop by balancing radiative cooling with turbulent heating. In addition to heating the plasma, turbulence also mixes it, suppressing the formation of cold gas at small scales. In this regard, the effect of turbulence is analogous to thermal conduction. For uniform plasma in thermal balance (turbulent heating balancing radiative cooling), cold gas condenses only if the cooling time is shorter than the mixing time. This condition requires the turbulent kinetic energy to be a parts per thousand(3) the plasma internal energy; such high velocities in cool cores are ruled out by observations. The results with realistic magnetic fields and thermal conduction are qualitatively similar to the hydrodynamic simulations. Simulations where the runaway cooling of the cool core is prevented due to mixing with the hot ICM show cold gas even with subsonic turbulence, consistent with observations. Thus, turbulent mixing is the likely mechanism via which AGN jets heat cluster cores. The thermal instability growth rates observed in simulations with turbulence are consistent with the local thermal instability interpretation of cold gas in cluster cores.

Flocking at a distance in active granular matter

The self-organized motion of vast numbers of creatures in a single direction is a spectacular example of emergent order. Here, we recreate this phenomenon using actuated nonliving components. We report here that millimetre-sized tapered rods, rendered motile by contact with an underlying vibrated surface and interacting through a medium of spherical beads, undergo a phase transition to a state of spontaneous alignment of velocities and orientations above a threshold bead area fraction. Guided by a detailed simulation model, we construct an analytical theory of this flocking transition, with two ingredients: a moving rod drags beads; neighbouring rods reorient in the resulting flow like a weathercock in the wind. Theory and experiment agree on the structure of our phase diagram in the plane of rod and bead concentrations and power-law spatial correlations near the phase boundary. Our discovery suggests possible new mechanisms for the collective transport of particulate or cellular matter.

Phase Diagram Study of Succinonitrile-Vanillin Organic Alloy System

Phase diagram studies of succinonitrile-vanillin system show the formation of 2:1 congruent melting type compound. Crystallization velocities of pure components, succinonitrile-vanillin complex, and two eutectics have been determined at different undercoolings. On the basis of heat of fusion measurements, excess thermodynamic functions have been calculated. Microstructural studies revealed that impurities modify the morphology. FTIR spectral studies and computer simulation have shown the existence of hydrogen bonding in the eutectics and the congruent melting compound. On the basis of experimental results, the mechanism of formation of eutectics and its solidification behavior are discussed.

Airblast spray in crossflow - Structure, trajectory and droplet sizing

This study reports results of an experimental investigation of airblast spray of water and ethanol in crossflow. Laser shadowgraphy and Particle/Droplet Imaging Analysis (PDIA) are used to derive spray trajectory and drop size information while Particle Tracking Velocimetry (PTV) is used to measure droplet velocities. A new phenomenon of spray bifurcation is observed for low Gas to Liquid Ratio (GLR) cases. The reasons for the spatial bifurcation can be attributed to a combination of reasons. These are (a) presence of large ligaments and droplets in the near-nozzle region for low GLRs (b) secondary breakup experienced by ligaments/droplets leading to formation of a large number of small droplets, and (c) the crossflow causing differential dispersion of the small and large droplets. A novel correlation for spray trajectory is proposed incorporating the momentum ratio and liquid surface tension. This correlation is shown to be effective in predicting the non-linear spray trajectory over a large range of conditions for not only water but ethanol and Jet-A also. It is observed that the larger droplets penetrate further into the crossflow, in the direction of injection. Thus, with increase in height of the measurement location from the injection plane, the droplet Sauter Mean Diameter (SMD) is found to increase. Moreover, as the droplets travel downstream in the crossflow direction, the droplet SMD is observed to decrease. The effect of drag is assessed by comparing velocity of different sizes of droplets at various locations. Smaller droplets are entrained into the crossflow at much lower elevations, whereas larger droplets tend to penetrate further into the crossflow. (C) 2015 Elsevier Ltd. All rights reserved.

EXPERIMENTAL INVESTIGATION OF UNSTEADY FLOW IN A TRANSONIC UNI-STAGE AXIAL COMPRESSOR

Detailed steady and unsteady experimental measurements and analysis were performed on a Single stage Transonic Axial Compressor with asymmetric rotor tip clearance for studying the compressor stall phenomena. The installed compressor had asymmetric tip clearance around the rotor casing varying from about 0.65mm to 1.25mm. A calibrated 5-hole aerodynamic probe was traversed radially at exit of rotor and showed the characteristics of increased flow angle at lower mass flow rates for all the speeds. Mach number distribution and boundary layer effects were also clearly captured. Unsteady measurements for velocity were carried out to study the stall cell behavior using a single component calibrated hotwire probe oriented in axial and tangential directions for choke/free flow and near stall conditions. The hotwire probe was traversed radially across the annulus at inlet to the compressor and showed that the velocity fluctuations were dissimilar when probe was aligned axial and tangential to the flow. Averaged velocities across the annulus showed the reduction in velocity as stall was approached. Axial mean flow velocity decreased across the annulus for all the speeds investigated. Tangential velocity at free flow condition was higher at the tip region due to larger radius. At stall condition, the tangential velocity showed decreased velocities at the tip and slightly increased velocities at the hub section indicating that the flow has breakdown at the tip region of the blade and fluid is accelerated below the blockage zone. The averaged turbulent intensity in axial and tangential flow directions increased from free flow to stall condition for all compressor rated speeds. Fast Fourier Transform (FFT) of the raw signals at stall flow condition showed stall cell and its corresponding frequency of occurrence. The stalling frequency of about half of rotational speed of the rotor along with large tip clearance suggests that modal type stall inception was occurring.

Homogeneous isotropic superfluid turbulence in two dimensions: Inverse and forward cascades in the Hall-Vinen-Bekharevich-Khalatnikov model

We present the first direct-numerical-simulation study of the statistical properties of two-dimensional superfluid turbulence in the simplified, Hall-Vinen-Bekharevich-Khalatnikov two-fluid model. We show that both normalfluid and superfluid energy spectra can exhibit two power-law regimes, the first associated with an inverse cascade of energy and the second with the forward cascade of enstrophy. We quantify the mutual-friction-induced alignment of normal and superfluid velocities by obtaining probability distribution functions of the angle between them and the ratio of their moduli.

QUANTITATIVE OH PLANAR LASER INDUCED FLUORESCENCE DIAGNOSTICS OF SYNGAS AND METHANE COMBUSTION IN A CAVITY COMBUSTOR

The current work reports quantitative OH species concentration in the cavity of a trapped vortex combustor (TVC) in the context of mixing and flame stabilization studies using both syngas and methane fuels. Planar laser induced fluorescence (PLIF) measurements of OH radical obtained using a Nd: YAG pumped dye laser are quantified using a flat flame McKenna burner. The momentum flux ratio (MFR), defined as the ratio of the cavity fuel jet momentum to that of the guide vane air stream, is observed to be a key governing parameter. At high MFRs similar to 4.5, the flame front is observed to form at the interface of the fuel jet and the air jet stream. This is substantiated by velocity vector field measurements. For syngas, as the MFR is lowered to similar to 0.3, the fuel-air mixing increases and a flame front is formed at the bottom and downstream edge of the cavity where a stratified charge is present. This trend is observed for different velocities at similar equivalence ratios. In case of methane combustion in the cavity, where the MFRs employed are extremely low at similar to 0.01, a different mechanism is observed. A fuel-rich mixture is now observed at the center of the cavity and this mixture undergoes combustion. On further increase of the cavity equivalence ratio, the rich mixture exceeds the flammability limit and forms a thin reaction zone at the interface with air stream. As a consequence, a shear layer flame at the top of the cavity interface with the mainstream is also observed. The equivalence ratio in the cavity also determines the combustion characteristics in the case of fuel-air mixtures that are formed as a result of the mixing. Overall, flame stabilization mechanisms have been proposed, which account for the wide range of MFRs and premixing in the mainstream as well.

Kinematic flow patterns in slow deformation of a dense granular material

The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using in situ imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale non-smooth contact dynamics simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.

Experimental studies on evaporation of fuel droplets under forced convection using spray in crossflow methodology

Experimental data on evaporation of droplets of decane, Jet-A1, and Jet-A1 surrogate are generated using a spray in crossflow configuration. The advantage of a crossflow configuration is that it enables us to study droplet evaporation under forced convective conditions involving droplet diameters of size relevant in practical combustors. Specifically, spray from an airblast atomizer is injected into a preheated crossflow of air and the resulting spray is characterized in terms of spray structure along with droplet size and velocity. An existing correlation for the spray trajectory is modified to incorporate the effect of elevated temperature, and is found to be in good agreement with the experimental data. Droplet sizes and velocities are measured at different locations along the crossflow direction to assess droplet evaporation. Specifically, droplets having size less than 25-mu m are selected for further analysis since these droplets are observed to exhibit velocities which are aligned with the crossflow. By comparing the droplet diameter profiles at upstream and downstream locations, the evaporation constant k for the d(2)-law is obtained iteratively. To assess the efficacy of the values of k obtained, the calculated droplet size distribution using the proposed k values at the downstream location is compared with the measured droplet size distribution at that location. A reasonably good match is found for all the three liquids confirming the validity of the analysis. (C) 2015 Elsevier Ltd. All rights reserved.

Intertwined vorticity and elastodynamics in flapping wing propulsion

We performed numerical experiments on a one-dimensional elastic solid oscillating in a two-dimensional viscous incompressible fluid with the intent of discerning the interplay of vorticity and elastodynamics in flapping wing propulsion. Perhaps for the first time, we have established the role of foil deflection topology and its influence on vorticity generation, through spatially and temporally evolving foil slope and curvature. Though the frequency of oscillation of the foil has a definite role, it is the phase relation between foil slope and pressure that determines thrust or drag. Similarly, the phase difference between flapping velocity, and pressure and inertial forces, determine the power input to the foil, and in turn drives propulsive efficiency. At low frequencies of oscillation, the sympathetic slope and curvature of deformation of the foil allow generation of leading-edge vortices that do not separate; they cause substantial rise in pressure between the leading edge and mid-chord. The circulatory component of pressure is determined primarily by the leading-edge vortex and therefore thrust too is predominantly circulatory in origin at low frequencies. In the intermediate and high-frequency range, thrust and drag on the foil spatially alternate and non-circulatory forces dominate over circulatory and viscous forces. For the mass ratios we simulated, thrust due to flapping varies quadratically as a function of Strouhal number or trailing-edge flapping velocity; further, the trailing edge flapping velocities peak at the same set of frequencies where the thrust is also a maximum. Propulsive efficiency, on the other hand, is roughly a mirror image of the thrust variation with respect to Strouhal number. Given that most instances of flapping propulsion in nature are primarily through distributed muscular actuation that enables precise control of deformation shape, leading to high thrust and efficiency, the results presented here are pointers towards understanding some of the mechanisms that drive thrust and propulsive efficiency.