Generalized model predictive static programming and its application to 3D impact angle constrained guidance of air-to-surface missiles

A new `generalized model predictive static programming (G-MPSP)' technique is presented in this paper in the continuous time framework for rapidly solving a class of finite-horizon nonlinear optimal control problems with hard terminal constraints. A key feature of the technique is backward propagation of a small-dimensional weight matrix dynamics, using which the control history gets updated. This feature, as well as the fact that it leads to a static optimization problem, are the reasons for its high computational efficiency. It has been shown that under Euler integration, it is equivalent to the existing model predictive static programming technique, which operates on a discrete-time approximation of the problem. Performance of the proposed technique is demonstrated by solving a challenging three-dimensional impact angle constrained missile guidance problem. The problem demands that the missile must meet constraints on both azimuth and elevation angles in addition to achieving near zero miss distance, while minimizing the lateral acceleration demand throughout its flight path. Both stationary and maneuvering ground targets are considered in the simulation studies. Effectiveness of the proposed guidance has been verified by considering first order autopilot lag as well as various target maneuvers.

Effects of elastase and collagenase on the nonlinearity and anisotropy of porcine aorta

We use enzymatic manipulation methods to investigate the individual and combined roles of elastin and collagen on arterial mechanics. Porcine aortic tissues were treated for differing amounts of time using enzymes elastase and collagenase to cause degradation in substrate proteins elastin and collagen and obtain variable tissue architecture. We use equibiaxial mechanical tests to quantify the material properties of control and enzyme treated tissues and histological methods to visualize the underlying tissue microstructure in arterial tissues. Our results show that collagenase treated tissues were more compliant in the longitudinal direction as compared to control tissues. Collagenase treatment also caused a decrease in the tissue nonlinearity as compared to the control samples in the study. A one hour collagenase treatment was sufficient to cause fragmentation and degradation of the adventitial collagen. In contrast, elastase treatment leads to significantly stiffer tissue response associated with fragmented and incomplete elastin networks in the tissue. Thus, elastin in arterial walls distributes tensile stresses whereas collagen serves to reinforce the vessel wall in the circumferential direction and also contributes to tissue anisotropy. A microstructurally motivated strain energy function based on circumferentially oriented medial fibers and helically oriented collagen fibers in the adventitia is useful in describing these experimental results.

On generalized gravitational entropy, squashed cones and holography

We consider generalized gravitational entropy in various higher derivative theories of gravity dual to four dimensional CFTs using the recently proposed regularization of squashed cones. We derive the universal terms in the entanglement entropy for spherical and cylindrical surfaces. This is achieved by constructing the Fefferman-Graham expansion for the leading order metrics for the bulk geometry and evaluating the generalized gravitational entropy. We further show that the Wald entropy evaluated in the bulk geometry constructed for the regularized squashed cones leads to the correct universal parts of the entanglement entropy for both spherical and cylindrical entangling surfaces. We comment on the relation with the Iyer-Wald formula for dynamical horizons relating entropy to a Noether charge. Finally we show how to derive the entangling surface equation in Gauss-Bonnet holography.

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.

Cubic Sieve Congruence of the Discrete Logarithm Problem, and fractional part sequences

The Cubic Sieve Method for solving the Discrete Logarithm Problem in prime fields requires a nontrivial solution to the Cubic Sieve Congruence (CSC) x(3) equivalent to y(2)z (mod p), where p is a given prime number. A nontrivial solution must also satisfy x(3) not equal y(2)z and 1 <= x, y, z < p(alpha), where alpha is a given real number such that 1/3 < alpha <= 1/2. The CSC problem is to find an efficient algorithm to obtain a nontrivial solution to CSC. CSC can be parametrized as x equivalent to v(2)z (mod p) and y equivalent to v(3)z (mod p). In this paper, we give a deterministic polynomial-time (O(ln(3) p) bit-operations) algorithm to determine, for a given v, a nontrivial solution to CSC, if one exists. Previously it took (O) over tilde (p(alpha)) time in the worst case to determine this. We relate the CSC problem to the gap problem of fractional part sequences, where we need to determine the non-negative integers N satisfying the fractional part inequality {theta N} < phi (theta and phi are given real numbers). The correspondence between the CSC problem and the gap problem is that determining the parameter z in the former problem corresponds to determining N in the latter problem. We also show in the alpha = 1/2 case of CSC that for a certain class of primes the CSC problem can be solved deterministically in <(O)over tilde>(p(1/3)) time compared to the previous best of (O) over tilde (p(1/2)). It is empirically observed that about one out of three primes is covered by the above class. (C) 2013 Elsevier B.V. All rights reserved.

Generalized State Estimation and Model Predictive Guidance for Spiraling and Ballistic Targets

An extended Kalman filter based generalized state estimation approach is presented in this paper for accurately estimating the states of incoming high-speed targets such as ballistic missiles. A key advantage of this nine-state problem formulation is that it is very much generic and can capture spiraling as well as pure ballistic motion of targets without any change of the target model and the tuning parameters. A new nonlinear model predictive zero-effort-miss based guidance algorithm is also presented in this paper, in which both the zero-effort-miss as well as the time-to-go are predicted more accurately by first propagating the nonlinear target model (with estimated states) and zero-effort interceptor model simultaneously. This information is then used for computing the necessary lateral acceleration. Extensive six-degrees-of-freedom simulation experiments, which include noisy seeker measurements, a nonlinear dynamic inversion based autopilot for the interceptor along with appropriate actuator and sensor models and magnitude and rate saturation limits for the fin deflections, show that near-zero miss distance (i.e., hit-to-kill level performance) can be obtained when these two new techniques are applied together. Comparison studies with an augmented proportional navigation based guidance shows that the proposed model predictive guidance leads to a substantial amount of conservation in the control energy as well.

A fractional order proportional integral controller for path following and trajectory tracking of miniature air vehicles

In this paper, a fractional order proportional-integral controller is developed for a miniature air vehicle for rectilinear path following and trajectory tracking. The controller is implemented by constructing a vector field surrounding the path to be followed, which is then used to generate course commands for the miniature air vehicle. The fractional order proportional-integral controller is simulated using the fundamentals of fractional calculus, and the results for this controller are compared with those obtained for a proportional controller and a proportional integral controller. In order to analyze the performance of the controllers, four performance metrics, namely (maximum) overshoot, control effort, settling time and integral of the timed absolute error cost, have been selected. A comparison of the nominal as well as the robust performances of these controllers indicates that the fractional order proportional-integral controller exhibits the best performance in terms of ITAE while showing comparable performances in all other aspects.

Generalized scaling of misorientation angle distributions at meso-scale in deformed materials

Scaling behaviour has been observed at mesoscopic level irrespective of crystal structure, type of boundary and operative micro-mechanisms like slip and twinning. The presence of scaling at the meso-scale accompanied with that at the nano-scale clearly demonstrates the intrinsic spanning for different deformation processes and a true universal nature of scaling. The origin of a 1/2 power law in deformation of crystalline materials in terms of misorientation proportional to square root of strain is attributed to importance of interfaces in deformation processes. It is proposed that materials existing in three dimensional Euclidean spaces accommodate plastic deformation by one dimensional dislocations and their interaction with two dimensional interfaces at different length scales. This gives rise to a 1/2 power law scaling in materials. This intrinsic relationship can be incorporated in crystal plasticity models that aim to span different length and time scales to predict the deformation response of crystalline materials accurately.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

SOLUTIONS OF A SYSTEM OF FORCED BURGERS EQUATION IN TERMS OF GENERALIZED LAGUERRE POLYNOMIALS

In this article, we obtain explicit solutions of a linear PDE subject to a class of radial square integrable functions with a monotonically increasing weight function |x|(n-1)e(beta vertical bar x vertical bar 2)/2, beta >= 0, x is an element of R-n. This linear PDE is obtained from a system of forced Burgers equation via the Cole-Hopf transformation. For any spatial dimension n > 1, the solution is expressed in terms of a family of weighted generalized Laguerre polynomials. We also discuss the large time behaviour of the solution of the system of forced Burgers equation.

Generalized Model Predictive Static Programming and Angle-Constrained Guidance of Air-to-Ground Missiles

A new generalized model predictive static programming technique is presented for rapidly solving a class of finite-horizon nonlinear optimal control problems with hard terminal constraints. Two key features for its high computational efficiency include one-time backward integration of a small-dimensional weighting matrix dynamics, followed bya static optimization formulation that requires only a static Lagrange multiplier to update the control history. It turns out that under Euler integration and rectangular approximation of finite integrals it is equivalent to the existing model predictive static programming technique. In addition to the benchmark double integrator problem, usefulness of the proposed technique is demonstrated by solving a three-dimensional angle-constrained guidance problem for an air-to-ground missile, which demands that the missile must meet constraints on both azimuth and elevation angles at the impact point in addition to achieving near-zero miss distance, while minimizing the lateral acceleration demand throughout its flight path. Simulation studies include maneuvering ground targets along with a first-order autopilot lag. Comparison studies with classical augmented proportional navigation guidance and modern general explicit guidance lead to the conclusion that the proposed guidance is superior to both and has a larger capture region as well.

Long range emission enhancement and anisotropy in coupled quantum dots induced by aligned gold nanoantenna

Quantum dot arrays have been projected as the material of choice for next generation displays and photodetectors. Extensive ongoing research aims at improving optical and electrical efficiencies of such devices. We report experimental results on non-local long range emission intensity enhancement and anisotropy in quantum dot assemblies induced by isolated and partially aligned gold nanoantennas. Spatially resolved photoluminescence clearly demonstrate that the effect is maximum, when the longitudinal surface plasmon resonance of the nanoantenna is resonant with the emission maxima of the quantum dots. We estimated the decay length of this enhancement to be similar to 2.6 mu m, which is considerably larger than the range of near field interaction of metal nanoantenna. Numerical simulations qualitatively capture the near field behavior of the nanorods but fail to match the experimentally observed non-local effects. We have suggested how strong interactions of quantum dots in the close packed assemblies, mediated by the nanoantennas, could lead to such observed behavior. (C) 2014 AIP Publishing LLC.

Effect of Anisotropy of Fibers on the Stress-Strain Response of Fiber-Reinforced Soil

Reinforcing soil with fibers is a useful method for improving the strength and settlement response of soil. The soil and fiber characteristics and their interaction are some of the major factors affecting the strength of reinforced soil. The fibers are usually randomly distributed in the soil, and their orientation has a significant effect on the behavior of the reinforced soil. In the paper, a study of the effect of anisotropic distribution of fibers on the stress-strain response is presented. Based on the concept of the modified Cam clay model, an analytical model was formulated for the fiber-reinforced soil, and the effect of fiber orientation on the stress-strain behavior of soil was studied in detail. The results show that, as the inclination of fibers with the horizontal plane increased, the contribution of fibers in improving the strength of fiber-reinforced soil decreased. The effect of fibers is maximum when they are in the direction of extension, and vice versa. (C) 2014 American Society of Civil Engineers.

Interphase anisotropy effects on lamellar eutectics: A numerical study

In directional solidification of binary eutectics, it is often observed that two-phase lamellar growth patterns grow tilted with respect to the direction z of the imposed temperature gradient. This crystallographic effect depends on the orientation of the two crystal phases alpha and beta with respect to z. Recently, an approximate theory was formulated that predicts the lamellar tilt angle as a function of the anisotropy of the free energy of the solid(alpha)-solid(beta) interphase boundary. We use two different numerical methods-phase field (PF) and dynamic boundary integral (BI)-to simulate the growth of steady periodic patterns in two dimensions as a function of the angle theta(R) between z and a reference crystallographic axis for a fixed relative orientation of alpha and beta crystals, that is, for a given anisotropy function (Wulff plot) of the interphase boundary. For Wulff plots without unstable interphase-boundary orientations, the two simulation methods are in excellent agreement with each other and confirm the general validity of the previously proposed theory. In addition, a crystallographic ``locking'' of the lamellae onto a facet plane is well reproduced in the simulations. When unstable orientations are present in the Wulff plot, it is expected that two distinct values of the tilt angle can appear for the same crystal orientation over a finite theta(R) range. This bistable behavior, which has been observed experimentally, is well reproduced by BI simulations but not by the PF model. Possible reasons for this discrepancy are discussed.

Tailoring local density of optical states to control emission intensity and anisotropy of quantum dots in hybrid photonic-plasmonic templates

We report results of controlled tuning of the local density of states (LDOS) in versatile, flexible, and hierarchical self assembled plasmonic templates. Using 5 nm diameter gold (Au) spherical nanoantenna within a polymer template randomly dispersed with quantum dots, we show how the photoluminescence intensity and lifetime anisotropy of these dots can be significantly enhanced through LDOS tuning. Finite difference time domain simulations corroborate the experimental observations and extend the regime of enhancement to a wider range of geometric and spectral parameters bringing out the versatility of these functional plasmonic templates. It is also demonstrated how the templates act as plasmonic resonators for effectively engineer giant enhancement of the scattering efficiency of these nano antenna embedded in the templates. Our work provides an alternative method to achieve spontaneous emission intensity and anisotropy enhancement with true nanoscale plasmon resonators. (C) 2015 AIP Publishing LLC.