Partial Integrated Guidance and Control of Interceptors for High-Speed Ballistic Targets

A new partial integrated guidance and control design approach is proposed in this paper, which combines the benefits of both integrated guidance and control as well as the conventional guidance and control design philosophies. The proposed technique essentially operates in a two-loop structure. In the outer loop, an optimal guidance problem is formulated considering the nonlinear six degrees-of-freedom equation of motion of the interceptor. From this loop, the required pitch and yaw rates are generated by solving a nonlinear suboptimal guidance formulation in a computationally efficient manner while simultaneously assuring roll stabilization. Next, the inner loop tracks these outer loop body rate commands. This manipulation of the six degrees-of-freedom dynamics in both loops preserves the inherent time scale separation property between the translational and rotational dynamics, while retaining the philosophy of integrated guidance and control design as well. Because of this, the tuning process is quite straightforward and nontedious as well. Extensive six degrees-of-freedom simulations studies have been carried out, considering three-dimensional engagement geometry, to demonstrate the effectiveness of the proposed new design approach engaging high-speed ballistic targets. A variety of comparison studies have also been carried out to demonstrate the effectiveness of the proposed approach.

Multi-site downscaling of maximum and minimum daily temperature using support vector machine

Climate change impact assessment studies involve downscaling large-scale atmospheric predictor variables (LSAPVs) simulated by general circulation models (GCMs) to site-scale meteorological variables. This article presents a least-square support vector machine (LS-SVM)-based methodology for multi-site downscaling of maximum and minimum daily temperature series. The methodology involves (1) delineation of sites in the study area into clusters based on correlation structure of predictands, (2) downscaling LSAPVs to monthly time series of predictands at a representative site identified in each of the clusters, (3) translation of the downscaled information in each cluster from the representative site to that at other sites using LS-SVM inter-site regression relationships, and (4) disaggregation of the information at each site from monthly to daily time scale using k-nearest neighbour disaggregation methodology. Effectiveness of the methodology is demonstrated by application to data pertaining to four sites in the catchment of Beas river basin, India. Simulations of Canadian coupled global climate model (CGCM3.1/T63) for four IPCC SRES scenarios namely A1B, A2, B1 and COMMIT were downscaled to future projections of the predictands in the study area. Comparison of results with those based on recently proposed multivariate multiple linear regression (MMLR) based downscaling method and multi-site multivariate statistical downscaling (MMSD) method indicate that the proposed method is promising and it can be considered as a feasible choice in statistical downscaling studies. The performance of the method in downscaling daily minimum temperature was found to be better when compared with that in downscaling daily maximum temperature. Results indicate an increase in annual average maximum and minimum temperatures at all the sites for A1B, A2 and B1 scenarios. The projected increment is high for A2 scenario, and it is followed by that for A1B, B1 and COMMIT scenarios. Projections, in general, indicated an increase in mean monthly maximum and minimum temperatures during January to February and October to December.

Multimodal shape oscillations of droplets excited by an air stream

The shape dynamics of droplets exposed to an air jet at intermediate droplet Reynolds numbers is investigated. High speed imaging and hot-wire anemometry are employed to examine the mechanism of droplet oscillation. The theory that the vortex shedding behind the droplet induces oscillation is examined. In these experiments, no particular dominant frequency is found in the wake region of the droplet. Hence the inherent free-stream disturbances prove to be driving the droplet oscillations. The modes of droplet oscillation show a band of dominant frequencies near the corresponding natural frequency, further proving that there is no particular forcing frequency involved. In the frequency spectrum of the lowest mode of oscillation for glycerol at the highest Reynolds number, no response is observed below the threshold frequency corresponding to the viscous dissipation time scale. This selective suppression of lower frequencies in the case of glycerol is corroborated by scaling arguments. The influence of surface tension on the droplet oscillation is studied using ethanol as a test fluid. Since a lower surface tension reduces the natural frequency, ethanol shows lower excited frequencies. The oscillation levels of different fluids are quantified using the droplet aspect ratio and correlated in terms of Weber number and Ohnesorge number. (C) 2014 Elsevier Ltd. All rights reserved.

Ising model on a random network with annealed or quenched disorder

We study the equilibrium properties of an Ising model on a disordered random network where the disorder can be quenched or annealed. The network consists of fourfold coordinated sites connected via variable length one-dimensional chains. Our emphasis is on nonuniversal properties and we consider the transition temperature and other equilibrium thermodynamic properties, including those associated with one-dimensional fluctuations arising from the chains. We use analytic methods in the annealed case, and a Monte Carlo simulation for the quenched disorder. Our objective is to study the difference between quenched and annealed results with a broad random distribution of interaction parameters. The former represents a situation where the time scale associated with the randomness is very long and the corresponding degrees of freedom can be viewed as frozen, while the annealed case models the situation where this is not so. We find that the transition temperature and the entropy associated with one-dimensional fluctuations are always higher for quenched disorder than in the annealed case. These differences increase with the strength of the disorder up to a saturating value. We discuss our results in connection to physical systems where a broad distribution of interaction strengths is present.

Evaluation of TRMM PR Sampling Error Over a Subtropical Basin Using Bootstrap Technique

Quantitative use of satellite-derived rainfall products for various scientific applications often requires them to be accompanied with an error estimate. Rainfall estimates inferred from low earth orbiting satellites like the Tropical Rainfall Measuring Mission (TRMM) will be subjected to sampling errors of nonnegligible proportions owing to the narrow swath of satellite sensors coupled with a lack of continuous coverage due to infrequent satellite visits. The authors investigate sampling uncertainty of seasonal rainfall estimates from the active sensor of TRMM, namely, Precipitation Radar (PR), based on 11 years of PR 2A25 data product over the Indian subcontinent. In this paper, a statistical bootstrap technique is investigated to estimate the relative sampling errors using the PR data themselves. Results verify power law scaling characteristics of relative sampling errors with respect to space-time scale of measurement. Sampling uncertainty estimates for mean seasonal rainfall were found to exhibit seasonal variations. To give a practical example of the implications of the bootstrap technique, PR relative sampling errors over a subtropical river basin of Mahanadi, India, are examined. Results reveal that the bootstrap technique incurs relative sampling errors < 33% (for the 2 degrees grid), < 36% (for the 1 degrees grid), < 45% (for the 0.5 degrees grid), and < 57% (for the 0.25 degrees grid). With respect to rainfall type, overall sampling uncertainty was found to be dominated by sampling uncertainty due to stratiform rainfall over the basin. The study compares resulting error estimates to those obtained from latin hypercube sampling. Based on this study, the authors conclude that the bootstrap approach can be successfully used for ascertaining relative sampling errors offered by TRMM-like satellites over gauged or ungauged basins lacking in situ validation data. This technique has wider implications for decision making before incorporating microwave orbital data products in basin-scale hydrologic modeling.

Surface Freshwater Storage Variations in the Orinoco Floodplains Using Multi-Satellite Observations

Variations in surface water extent and storage are poorly characterized from regional to global scales. In this study, a multi-satellite approach is proposed to estimate the water stored in the floodplains of the Orinoco Basin at a monthly time-scale using remotely-sensed observations of surface water from the Global Inundation Extent Multi-Satellite (GIEMS) and stages from Envisat radar altimetry. Surface water storage variations over 2003-2007 exhibit large interannual variability and a strong seasonal signal, peaking during summer, and associated with the flood pulse. The volume of surface water storage in the Orinoco Basin was highly correlated with the river discharge at Ciudad Bolivar (R = 0.95), the closest station to the mouth where discharge was estimated, although discharge lagged one month behind storage. The correlation remained high (R = 0.73) after removing seasonal effects. Mean annual variations in surface water volume represented similar to 170 km(3), contributing to similar to 45% of the Gravity Recovery and Climate Experiment (GRACE)-derived total water storage variations and representing similar to 13% of the total volume of water that flowed out of the Orinoco Basin to the Atlantic Ocean.

Instabilities in unsteady boundary layers with reverse flow

Instabilities arising in unsteady boundary layers with reverse flow have been investigated experimentally. Experiments are conducted in a piston driven unsteady water tunnel with a shallow angle diffuser placed in the test section. The ratio of temporal (Pi(t)) to spatial (Pi(x)) component of the pressure gradient can be varied by a controlled motion of the piston. In all the experiments, the piston velocity variation with time is trapezoidal consisting of three phases: constant acceleration from rest, constant velocity and constant deceleration to rest. The adverse pressure gradient (and reverse flow) are due to a combination of spatial deceleration of the free stream in the diffuser and temporal deceleration of the free stream caused by the piston deceleration. The instability is usually initiated with the formation of one or more vortices. The onset of reverse flow in the boundary layer, location and time of formation of the first vortex and the subsequent flow evolution are studied for various values of the ratio Pi(x) (Pi(x) + Pi(t)) for the bottom and the top walls. Instability is due to the inflectional velocity profiles of the unsteady boundary layer. The instability is localized and spreads to the other regions at later times. At higher Reynolds numbers growth rate of instability is higher and localized transition to turbulence is observed. Scalings have been proposed for initial vortex formation time and wavelength of the instability vortices. Initial vortex formation time scales with convective time, delta/Delta U, where S is the boundary layer thickness and Delta U is the difference of maximum and minimum velocities in the boundary layer. Non-dimensional vortex formation time based on convective time scale for the bottom and the top walls are found to be 23 and 30 respectively. Wavelength of instability vortices scales with the time averaged boundary layer thickness. (C) 2015 Elsevier Masson SAS. All rights reserved.

Quasi-geostrophic dynamics in the presence of moisture gradients

The derivation of a quasi-geostrophic system from the rotating shallow-water equations on a midlatitude -plane coupled with moisture is presented. Condensation is prescribed to occur whenever the moisture at a point exceeds a prescribed saturation value. It is seen that a slow condensation time-scale is required to obtain a consistent set of equations at leading order. Further, since the advecting wind fields are geostrophic, changes in moisture (and hence precipitation) occur only via non-divergent mechanisms. Following observations, a saturation profile with gradients in the zonal and meridional directions is prescribed. A purely meridional gradient has the effect of slowing down the dry Rossby waves, through a reduction in the equivalent gradient' of the background potential vorticity. A large-scale unstable moist mode results on the inclusion of a zonal gradient by itself, or in conjunction with a meridional moisture gradient. For gradients that are are representative of the atmosphere, the most unstable moist mode propagates zonally in the direction of increasing moisture, matures over an intraseasonal time-scale and has small phase speed.

Deformation pathways and breakup modes in acoustically levitated bicomponent droplets under external heating

Controlled breakup of droplets using heat or acoustics is pivotal in applications such as pharmaceutics, nanoparticle production, and combustion. In the current work we have identified distinct thermal acoustics-induced deformation regimes (ligaments and bubbles) and breakup dynamics in externally heated acoustically levitated bicomponent (benzene-dodecane) droplets with a wide variation in volatility of the two components (benzene is significantly more volatile than dodecane). We showcase the physical mechanism and universal behavior of droplet surface caving in leading to the inception and growth of ligaments. The caving of the top surface is governed by a balance between the acoustic pressure field and the restrictive surface tension of the droplet. The universal collapse of caving profiles for different benzene concentration (<70% by volume) is shown by using an appropriate time scale obtained from force balance. Continuous caving leads to the formation of a liquid membrane-type structure which undergoes radial extension due to inertia gained during the precursor phase. The membrane subsequently closes at the rim and the kinetic energy leads to ligament formation and growth. Subsequent ligament breakup is primarily Rayleigh-Plateau type. The breakup mode shifts to diffusional entrapment-induced boiling with an increase in concentration of the volatile component (benzene >70% by volume). The findings are portable to any similar bicomponent systems with differential volatility.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.

Simulated long-term climate response to idealized solar geoengineering

Solar geoengineering has been proposed as a potential means to counteract anthropogenic climate change, yet it is unknown how such climate intervention might affect the Earth's climate on the millennial time scale. Here we use the HadCM3L model to conduct a 1000year sunshade geoengineering simulation in which solar irradiance is uniformly reduced by 4% to approximately offset global mean warming from an abrupt quadrupling of atmospheric CO2. During the 1000year period, modeled global climate, including temperature, hydrological cycle, and ocean circulation of the high-CO2 simulation departs substantially from that of the control preindustrial simulation, whereas the climate of the geoengineering simulation remains much closer to that of the preindustrial state with little drift. The results of our study do not support the hypothesis that nonlinearities in the climate system would cause substantial drift in the climate system if solar geoengineering was to be deployed on the timescale of a millennium.

Capillary Effect on Vertically Excited Surface Wave in Circular Cylindrical Vessel

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions was developed in inviscid fluids to investigate the motion of single free surface standing wave including the effect of surface tension.

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.

Burst event detection in wall turbulence by WVITA method

Wavelet Variable Interval Time Average (WVITA) is introduced as a method incorporating burst event detection in wall turbulence. Wavelet transform is performed to unfold the longitudinal fluctuating velocity time series measured in the near wall region of a turbulent boundary layer using hot-film anemometer. This unfolding is both in time and in space simultaneously. The splitted kinetic of the longitudinal fluctuating velocity time series among different scales is obtained by integrating the square of wavelet coefficient modulus over temporal space. The time scale that related to burst events in wall turbulence passing through the fixed probe is ascertained by maximum criterion of the kinetic energy evolution across scales. Wavelet transformed localized variance of the fluctuating velocity time series at the maximum kinetic scale is put forward instead of localized short time average variance in Variable Interval Time Average (VITA) scheme. The burst event detection result shows that WVITA scheme can avoid erroneous judgement and solve the grouping problem more effectively which is caused by VITA scheme itself and can not be avoided by adjusting the threshold level or changing the short time average interval.

Damping of Vertically Excited Surface Wave in Weakly Viscous Fluid

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.