Differential evolution based 3-D guidance law for a realistic interceptor model

This paper presents a novel, soft computing based solution to a complex optimal control or dynamic optimization problem that requires the solution to be available in real-time. The complexities in this problem of optimal guidance of interceptors launched with high initial heading errors include the more involved physics of a three dimensional missile-target engagement, and those posed by the assumption of a realistic dynamic model such as time-varying missile speed, thrust, drag and mass, besides gravity, and upper bound on the lateral acceleration. The classic, pure proportional navigation law is augmented with a polynomial function of the heading error, and the values of the coefficients of the polynomial are determined using differential evolution (DE). The performance of the proposed DE enhanced guidance law is compared against the existing conventional laws in the literature, on the criteria of time and energy optimality, peak lateral acceleration demanded, terminal speed and robustness to unanticipated target maneuvers, to illustrate the superiority of the proposed law. (C) 2013 Elsevier B. V. All rights reserved.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

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.

On the Parameterized Complexity of the Maximum Edge 2-Coloring Problem

We investigate the parameterized complexity of the following edge coloring problem motivated by the problem of channel assignment in wireless networks. For an integer q >= 2 and a graph G, the goal is to find a coloring of the edges of G with the maximum number of colors such that every vertex of the graph sees at most q colors. This problem is NP-hard for q >= 2, and has been well-studied from the point of view of approximation. Our main focus is the case when q = 2, which is already theoretically intricate and practically relevant. We show fixed-parameter tractable algorithms for both the standard and the dual parameter, and for the latter problem, the result is based on a linear vertex kernel.

Parameterized Algorithms for MAX COLORABLE INDUCED SUBGRAPH Problem on Perfect Graphs

We address the parameterized complexity ofMaxColorable Induced Subgraph on perfect graphs. The problem asks for a maximum sized q-colorable induced subgraph of an input graph G. Yannakakis and Gavril IPL 1987] showed that this problem is NP-complete even on split graphs if q is part of input, but gave a n(O(q)) algorithm on chordal graphs. We first observe that the problem is W2]-hard parameterized by q, even on split graphs. However, when parameterized by l, the number of vertices in the solution, we give two fixed-parameter tractable algorithms. The first algorithm runs in time 5.44(l) (n+#alpha(G))(O(1)) where #alpha(G) is the number of maximal independent sets of the input graph. The second algorithm runs in time q(l+o()l())n(O(1))T(alpha) where T-alpha is the time required to find a maximum independent set in any induced subgraph of G. The first algorithm is efficient when the input graph contains only polynomially many maximal independent sets; for example split graphs and co-chordal graphs. The running time of the second algorithm is FPT in l alone (whenever T-alpha is a polynomial in n), since q <= l for all non-trivial situations. Finally, we show that (under standard complexitytheoretic assumptions) the problem does not admit a polynomial kernel on split and perfect graphs in the following sense: (a) On split graphs, we do not expect a polynomial kernel if q is a part of the input. (b) On perfect graphs, we do not expect a polynomial kernel even for fixed values of q >= 2.

Exploring the Consequences of a Representative ``Disallowed'' Conformation of Aib on a 3(10)-Helical Fold

The structural effects of a representative disallowed conformation of Aib on the 3(10)-helical fold of an octapeptidomimetic are explored. The 1D (H-1, C-13) & 2D NMR, FT-IR and CD data reveal that the octapeptide 1, adopts a 3(10)-helical conformation in solution, as it does in its crystal structure. The C-terminal methyl carboxylate (CO2Me) of 1 was modified into an 1,3-oxazine (Oxa) functional group in the peptidomimetic 2. This modification results in the stabilization of the backbone of the C-terminal Aib (Aib*-Oxa) of 2, in a conformation (phi, =180, 0) that is natively disallowed to Aib. Consequent to the presence of this natively disallowed conformation, the 3(10)-helical fold is not disrupted in the body of the peptidomimetic 2. But the structural distortions that do occur in 2 are primarily in residues in the immediate vicinity of the natively disallowed conformation, rather than in the whole peptide body. Non-native electronic effects resulting from modifications in backbone functional groups can be at the origin of stabilizing residues in natively disallowed conformations. (c) 2014 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 104: 21-36, 2015.

An encaged sigma* orbital in 3-methoxy-2,4,10-trioxaadamantane. Preferred exo-anomeric effect as revealed by the crystal structure

In the titled tricyclic orthocarbonate derivative, the three endocyclic C-O bonds are longer than the exo-cyclic C-O bond (similar to 4.40 angstrom vs. similar to 1.37 angstrom). This indicates an anomeric-type interaction between the two electron lone pairs on the exocyclic oxygen atom and the antibonding orbitals of the two antiperiplanar endocyclic C-O bonds. The remaining endocyclic C-O bond - marginally shorter than the other two apparently adds to this effect. Intriguingly, the antibonding orbital of the exocyclic C-O bond extends into the interior of the adamantyl cage, and is stereoelectronically prevented from overlapping with any of the six adjacent lone pairs. The results also seem to indicate a preference for interaction between a single donor oxygen atom and multiple acceptor antibonding orbitals rather than vice versa. The results add insightfully to the substantial body of evidence favouring the antiperiplanar lone pair hypothesis (ALPH). (C) 2014 Elsevier B.V. All rights reserved.

On the impact of a concave nosed axisymmetric body on a free surface

We report on an experimental study of the vertical impact of a concave nosed axisymmetric body on a free surface. Previous studies have shown that bodies with a convex nose, like a sphere, produce a well defined splash with a relatively large cavity behind the model. In contrast, we find that with a concave nose, there is hardly a splash and the cavity extent is greatly reduced. This may be explained by the fact that in the concave nosed case, the initial impact is between a confined air pocket and the free surface unlike in the convex nosed case. From measurements of the unsteady pressure in the concave nose portion, we show that in this case, the maximum pressures are significantly lower than the classically expected ``water hammer'' pressures and also lower than those generally measured on other geometries. Thus, the presence of an air pocket in the case of a concave nosed body adds an interesting dimension to the classical problem of impact of solid bodies on to a free surface. (C) 2015 AIP Publishing LLC.

COMPUTATIONAL MODELING OF A MULTI-LAYERED PIEZO-COMPOSITE BEAM MADE UP OF MFC

The Variational Asymptotic Method (VAM) is used for modeling a coupled non-linear electromechanical problem finding applications in aircrafts and Micro Aerial Vehicle (MAV) development. VAM coupled with geometrically exact kinematics forms a powerful tool for analyzing a complex nonlinear phenomena as shown previously by many in the literature 3 - 7] for various challenging problems like modeling of an initially twisted helicopter rotor blades, matrix crack propagation in a composite, modeling of hyper elastic plates and various multi-physics problems. The problem consists of design and analysis of a piezocomposite laminate applied with electrical voltage(s) which can induce direct and planar distributed shear stresses and strains in the structure. The deformations are large and conventional beam theories are inappropriate for the analysis. The behavior of an elastic body is completely understood by its energy. This energy must be integrated over the cross-sectional area to obtain the 1-D behavior as is typical in a beam analysis. VAM can be used efficiently to approximate 3-D strain energy as closely as possible. To perform this simplification, VAM makes use of thickness to width, width to length, width multiplied by initial twist and strain as small parameters embedded in the problem definition and provides a way to approach the exact solution asymptotically. In this work, above mentioned electromechanical problem is modeled using VAM which breaks down the 3-D elasticity problem into two parts, namely a 2-D non-linear cross-sectional analysis and a 1-D non-linear analysis, along the reference curve. The recovery relations obtained as a by-product in the cross-sectional analysis earlier are used to obtain 3-D stresses, displacements and velocity contours. The piezo-composite laminate which is chosen for an initial phase of computational modeling is made up of commercially available Macro Fiber Composites (MFCs) stacked together in an arbitrary lay-up and applied with electrical voltages for actuation. The expressions of sectional forces and moments as obtained from cross-sectional analysis in closed-form show the electro-mechanical coupling and relative contribution of electric field in individual layers of the piezo-composite laminate. The spatial and temporal constitutive law as obtained from the cross-sectional analysis are substituted into 1-D fully intrinsic, geometrically exact equilibrium equations of motion and 1-D intrinsic kinematical equations to solve for all 1-D generalized variables as function of time and an along the reference curve co-ordinate, x(1).

Two player game variant of the Erdos-Szekeres problem

The classical Erdos-Szekeres theorem states that a convex k-gon exists in every sufficiently large point set. This problem has been well studied and finding tight asymptotic bounds is considered a challenging open problem. Several variants of the Erdos-Szekeres problem have been posed and studied in the last two decades. The well studied variants include the empty convex k-gon problem, convex k-gon with specified number of interior points and the chromatic variant. In this paper, we introduce the following two player game variant of the Erdos-Szekeres problem: Consider a two player game where each player playing in alternate turns, place points in the plane. The objective of the game is to avoid the formation of the convex k-gon among the placed points. The game ends when a convex k-gon is formed and the player who placed the last point loses the game. In our paper we show a winning strategy for the player who plays second in the convex 5-gon game and the empty convex 5-gon game by considering convex layer configurations at each step. We prove that the game always ends in the 9th step by showing that the game reaches a specific set of configurations.

Exploring the Consequences of a Representative ``Disallowed'' Conformation of Aib on a 3(10)-Helical Fold

The structural effects of a representative ``disallowed'' conformation of Aib on the 3(10)-helical fold of an octapeptidomimetic are explored. The 1D (H-1, C-13) & 2D NMR, FT-IR and CD data reveal that the octapeptide 1, adopts a 3(10)- helical conformation in solution, as it does in its crystal structure. The C-terminal methyl carboxylate (CO2Me) of 1 was modified into an 1,3-oxazine (Oxa) functional group in the peptidomimetic 2. This modification results in the stabilization of the backbone of the C-terminal Aib (Aib(star)-Oxa) of 2, in a conformation (phi, psi = 180, 0) that is natively disallowed to Aib. Consequent to the presence of this natively disallowed conformation, the 3(10)- helical fold is not disrupted in the body of the peptidomimetic 2. But the structural distortions that do occur in 2 are primarily in residues in the immediate vicinity of the natively disallowed conformation, rather than in the whole peptide body. Non-native electronic effects resulting from modifications in backbone functional groups can be at the origin of stabilizing residues in natively disallowed conformations. (C) 2014 Wiley Periodicals, Inc. Biopolymers

Exploring the Consequences of a Representative ``Disallowed'' Conformation of Aib on a 3(10)-Helical Fold

The structural effects of a representative ``disallowed'' conformation of Aib on the 3(10)-helical fold of an octapeptidomimetic are explored. The 1D (H-1, C-13) & 2D NMR, FT-IR and CD data reveal that the octapeptide 1, adopts a 3(10)- helical conformation in solution, as it does in its crystal structure. The C-terminal methyl carboxylate (CO2Me) of 1 was modified into an 1,3-oxazine (Oxa) functional group in the peptidomimetic 2. This modification results in the stabilization of the backbone of the C-terminal Aib (Aib(star)-Oxa) of 2, in a conformation (phi, psi = 180, 0) that is natively disallowed to Aib. Consequent to the presence of this natively disallowed conformation, the 3(10)- helical fold is not disrupted in the body of the peptidomimetic 2. But the structural distortions that do occur in 2 are primarily in residues in the immediate vicinity of the natively disallowed conformation, rather than in the whole peptide body. Non-native electronic effects resulting from modifications in backbone functional groups can be at the origin of stabilizing residues in natively disallowed conformations. (C) 2014 Wiley Periodicals, Inc. Biopolymers

Ultra-sensitive pressure dependence of bandgap of rutile-GeO2 revealed by many body perturbation theory

The reported values of bandgap of rutile GeO2 calculated by the standard density functional theory within local-density approximation (LDA)/generalized gradient approximation (GGA) show a wide variation (similar to 2 eV), whose origin remains unresolved. Here, we investigate the reasons for this variation by studying the electronic structure of rutile-GeO2 using many-body perturbation theory within the GW framework. The bandgap as well as valence bandwidth at Gamma-point of rutile phase shows a strong dependence on volume change, which is independent of bandgap underestimation problem of LDA/GGA. This strong dependence originates from a change in hybridization among O-p and Ge-(s and p) orbitals. Furthermore, the parabolic nature of first conduction band along X-Gamma-M direction changes towards a linear dispersion with volume expansion. (C) 2015 AIP Publishing LLC.

A Variant of the Hadwiger-Debrunner (p, q)-Problem in the Plane

Let X be a convex curve in the plane (say, the unit circle), and let be a family of planar convex bodies such that every two of them meet at a point of X. Then has a transversal of size at most . Suppose instead that only satisfies the following ``(p, 2)-condition'': Among every p elements of , there are two that meet at a common point of X. Then has a transversal of size . For comparison, the best known bound for the Hadwiger-Debrunner (p, q)-problem in the plane, with , is . Our result generalizes appropriately for if is, for example, the moment curve.

Parameterized Algorithms and Kernels for 3-Hitting Set with Parity Constraints

The 3-Hitting Set problem involves a family of subsets F of size at most three over an universe U. The goal is to find a subset of U of the smallest possible size that intersects every set in F. The version of the problem with parity constraints asks for a subset S of size at most k that, in addition to being a hitting set, also satisfies certain parity constraints on the sizes of the intersections of S with each set in the family F. In particular, an odd (even) set is a hitting set that hits every set at either one or three (two) elements, and a perfect code is a hitting set that intersects every set at exactly one element. These questions are of fundamental interest in many contexts for general set systems. Just as for Hitting Set, we find these questions to be interesting for the case of families consisting of sets of size at most three. In this work, we initiate an algorithmic study of these problems in this special case, focusing on a parameterized analysis. We show, for each problem, efficient fixed-parameter tractable algorithms using search trees that are tailor-made to the constraints in question, and also polynomial kernels using sunflower-like arguments in a manner that accounts for equivalence under the additional parity constraints.