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.

Tailoring the second mode of Euler-Bernoulli beams: an analytical approach

In this paper, we study the inverse mode shape problem for an Euler-Bernoulli beam, using an analytical approach. The mass and stiffness variations are determined for a beam, having various boundary conditions, which has a prescribed polynomial second mode shape with an internal node. It is found that physically feasible rectangular cross-section beams which satisfy the inverse problem exist for a variety of boundary conditions. The effect of the location of the internal node on the mass and stiffness variations and on the deflection of the beam is studied. The derived functions are used to verify the p-version finite element code, for the cantilever boundary condition. The paper also presents the bounds on the location of the internal node, for a valid mass and stiffness variation, for any given boundary condition. The derived property variations, corresponding to a given mode shape and boundary condition, also provides a simple closed-form solution for a class of non-uniform Euler-Bernoulli beams. These closed-form solutions can also be used to check optimization algorithms proposed for modal tailoring.

Spiral-wave dynamics in ionically realistic mathematical models for human ventricular tissue: the effects of periodic deformation

We carry out an extensive numerical study of the dynamics of spiral waves of electrical activation, in the presence of periodic deformation (PD) in two-dimensional simulation domains, in the biophysically realistic mathematical models of human ventricular tissue due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, and Panfilov (the TNNPO4 model). We first consider simulations in cable-type domains, in which we calculate the conduction velocity theta and the wavelength lambda of a plane wave; we show that PD leads to a periodic, spatial modulation of theta and a temporally periodic modulation of lambda; both these modulations depend on the amplitude and frequency of the PD. We then examine three types of initial conditions for both TP06 and TNNPO4 models and show that the imposition of PD leads to a rich variety of spatiotemporal patterns in the transmembrane potential including states with a single rotating spiral (RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state with multiple broken spirals, and a state SA in which all spirals are absorbed at the boundaries of our simulation domain. We find, for both TP06 and TNNPO4 models, that spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and the initial condition. We examine how these different types of spiral-wave states can be eliminated in the presence of PD by the application of low-amplitude pulses by square- and rectangular-mesh suppression techniques. We suggest specific experiments that can test the results of our simulations.

Controlling the cross-section of waveguides inscribed on bulk chalcogenide glasses using fast laser

Waveguides have been fabricated on melt-quenched, bulk chalcogenide glasses using the femto-second laser inscription technique at low repetition rates in the single scan regime. The inscribed waveguides have been characterized by butt-coupling method and the diameter of the waveguide calculated using the mode-field image of the waveguide. The waveguide cross-section symmetry is analyzed using the heat diffusion model by relating the energy and translation speed of the laser. The net-fluence and symmetry of the waveguides are correlated based on the theoretical values and experimental results of guiding cross-section.

Study of structures, energies and vibrational frequencies of (O-2)(n)(+) (n=2-5) clusters by GGA and meta-GGA density functional methods

Using Generalized Gradient Approximation (GGA) and meta-GGA density functional methods, structures, binding energies and harmonic vibrational frequencies for the clusters O-4(+), O-6(+), O-8(+) and O-10(+) have been calculated. The stable structures of O-4(+), O-6(+), O-8(+) and O-10(+) have point groups D-2h, D-3h, D-4h, and D-5h optimized on the quartet, sextet, octet and dectet potential energy surfaces, respectively. Rectangular (D-2h) O-4(+) has been found to be more stable compared to trans-planar (C-2h) on the quartet potential energy surface. Cyclic structure (D-3h) of CA cluster ion has been calculated to be more stable than other structures. Binding energy (B.E.) of the cyclic O-6(+) is in good agreement with experimental measurement. The zero-point corrected B.E. of O-8(+) with D4h symmetry on the octet potential energy surface and zero-point corrected B.E. of O-10(+) with D-5h symmetry on the dectet potential energy surface are also in good agreement with experimental values. The B.E. value for O-4(+) is close to the experimental value when single point energy is calculated by Brueckner coupled-cluster method, BD(T). (C) 2014 Elsevier B.V. All rights reserved.

Non-uniform beams and stiff strings isospectral to axially loaded uniform beams and piano strings

In this paper, we derive analytical expressions for mass and stiffness functions of transversely vibrating clamped-clamped non-uniform beams under no axial loads, which are isospectral to a given uniform axially loaded beam. Examples of such axially loaded beams are beam columns (compressive axial load) and piano strings (tensile axial load). The Barcilon-Gottlieb transformation is invoked to transform the non-uniform beam equation into the axially loaded uniform beam equation. The coupled ODEs involved in this transformation are solved for two specific cases (pq (z) = k (0) and q = q (0)), and analytical solutions for mass and stiffness are obtained. Examples of beams having a rectangular cross section are shown as a practical application of the analysis. Some non-uniform beams are found whose frequencies are known exactly since uniform axially loaded beams with clamped ends have closed-form solutions. In addition, we show that the tension required in a stiff piano string with hinged ends can be adjusted by changing the mass and stiffness functions of a stiff string, retaining its natural frequencies.

A matrix model for QCD

Gribov's observation that global gauge fixing is impossible has led to suggestions that there may be a deep connection between gauge fixing and confinement. We find an unexpected relation between the topological nontriviality of the gauge bundle and colored states in SU(N) Yang-Mills theory, and show that such states are necessarily impure. We approximate QCD by a rectangular matrix model that captures the essential topological features of the gauge bundle, and demonstrate the impure nature of colored states explicitly. Our matrix model also allows the inclusion of the QCD theta-term, as well as to perform explicit computations of low-lying glueball masses. This mass spectrum is gapped. Since an impure state cannot evolve to a pure one by a unitary transformation, our result shows that the solution to the confinement problem in pure QCD is fundamentally quantum information-theoretic.

CROSS-SECTIONAL ANALYSIS OF PRE-TWISTED THICK BEAMS USING VARIATIONAL ASYMPTOTIC METHOD

An asymptotically-exact methodology is presented for obtaining the cross-sectional stiffness matrix of a pre-twisted moderately-thick beam having rectangular cross sections and made of transversely isotropic materials. The anisotropic beam is modeled from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy of the beam is computed making use of the constitutive law and the kinematical relations derived with the inclusion of geometrical nonlinearities and initial twist. Large displacements and rotations are allowed, but small strain is assumed. The Variational Asymptotic Method is used to minimize the energy functional, thereby reducing the cross section to a point on the reference line with appropriate properties, yielding a 1-D constitutive law. In this method as applied herein, the 2-D cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged as orders of the small parameters. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that renders the 1-D strain measures well-defined. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

Boxicity and cubicity of product graphs

The boxicity (cubicity) of a graph G is the minimum natural number k such that G can be represented as an intersection graph of axis-parallel rectangular boxes (axis-parallel unit cubes) in R-k. In this article, we give estimates on the boxicity and the cubicity of Cartesian, strong and direct products of graphs in terms of invariants of the component graphs. In particular, we study the growth, as a function of d, of the boxicity and the cubicity of the dth power of a graph with respect to the three products. Among others, we show a surprising result that the boxicity and the cubicity of the dth Cartesian power of any given finite graph is, respectively, in O(log d/ log log d) and circle dot(d/ log d). On the other hand, we show that there cannot exist any sublinear bound on the growth of the boxicity of powers of a general graph with respect to strong and direct products. (C) 2015 Elsevier Ltd. All rights reserved.

Electromechanical Analysis of Tapered Piezoelectric Bimorph at High Electric Field

Piezoelectric bimorph laminar actuator of tapered width exhibits better performance for out-of-plane deflection compared to the rectangular surface area, while consuming equal surface area. This paper contains electromechanical analysis and modeling of a tapered width piezoelectric bimorph laminar actuator at high electric field in static state. The analysis is based on the second order constitutive equations of piezoelectric material, assuming small strain and large electric field to capture its behavior at high electric field. Analytical expressions are developed for block force, output strain energy, output energy density, input electrical energy, capacitance and energy efficiency at high electric field. The analytical expressions show that for fixed length, thickness, and surface area of the actuator, how the block force and output strain energy gets improved in a tapered surface actuator compared to a rectangular surface. Constant thickness, constant length and constant surface area of the actuator ensure constant mass, and constant electrical capacitance. We consider high electric field in both series and parallel electrical connection for the analysis. Part of the analytical results is validated with the experimental results, which are reported in earlier literature.

Analysis of Transverse Shear Strains in Pre-Twisted Thick Beams using Variational Asymptotic Method

The cross-sectional stiffness matrix is derived for a pre-twisted, moderately thick beam made of transversely isotropic materials and having rectangular cross sections. An asymptotically-exact methodology is used to model the anisotropic beam from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy is computed making use of the beam constitutive law and kinematical relations derived with the inclusion of geometrical nonlinearities and an initial twist. The energy functional is minimized making use of the Variational Asymptotic Method (VAM), thereby reducing the cross section to a point on the beam reference line with appropriate properties, forming a 1-D constitutive law. VAM is a mathematical technique employed in the current problem to rigorously split the 3-D analysis of beams into two: a 2-D analysis over the beam cross-sectional domain, which provides a compact semi-analytical form of the properties of the cross sections, and a nonlinear 1-D analysis of the beam reference curve. In this method, as applied herein, the cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged in orders of the small parameters. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that render the 1-D strain measures well-defined. The zeroth-order 3-D warping field thus yielded is then used to integrate the 3-D strain energy density over the cross section, resulting in the 1-D strain energy density, which in turn helps identify the corresponding cross-sectional stiffness matrix. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

A Statistical Model of Current Loops and Magnetic Monopoles

We formulate a natural model of loops and isolated vertices for arbitrary planar graphs, which we call the monopole-dimer model. We show that the partition function of this model can be expressed as a determinant. We then extend the method of Kasteleyn and Temperley-Fisher to calculate the partition function exactly in the case of rectangular grids. This partition function turns out to be a square of a polynomial with positive integer coefficients when the grid lengths are even. Finally, we analyse this formula in the infinite volume limit and show that the local monopole density, free energy and entropy can be expressed in terms of well-known elliptic functions. Our technique is a novel determinantal formula for the partition function of a model of isolated vertices and loops for arbitrary graphs.

Development of cup shaped microneedle array for transdermal drug delivery

Microneedle technology is one of the attractive methods in transdermal drug delivery. However, the clinical applications of this method are limited owing to: complexity in the preparation of multiple coating solutions, drug leakage while inserting the microneedles into the skin and the outer walls of the solid microneedle can hold limited quantity of drug. Here, the authors present the fabrication of an array of rectangular cup shaped silicon microneedles, which provide for reduced drug leakage resulting in improvement of efficiency of drug delivery and possibility of introducing multiple drugs. The fabricated solid microneedles with rectangular cup shaped tip have a total height of 200 mu m. These cup shaped tips have dimensions: 60 x 60 mu m (length x breadth) with a depth of 60 mu m. The cups are filled with drug using a novel in-house built drop coating system. Successful drug dissolution was observed when the coated microneedle was used on mice. Also, using the above method, it is possible to fill the cups selectively with different drugs, which enables simultaneous multiple drug delivery. (C) 2015 American Vacuum Society.

Interplay of Substrate Conductivity, Cellular Microenvironment, and Pulsatile Electrical Stimulation toward Osteogenesis of Human Mesenchymal Stem Cells in Vitro

The influences of physical stimuli such as surface elasticity, topography, and chemistry over mesenchymal stem cell proliferation and differentiation are well investigated. In this context, a fundamentally different approach was adopted, and we have demonstrated the interplay of inherent substrate conductivity, defined chemical composition of cellular microenvironment, and intermittent delivery of electric pulses to drive mesenchymal stem cell differentiation toward osteogenesis. For this, conducting polyaniline (PANI) substrates were coated with collagen type 1 (Coll) alone or in association with sulfated hyaluronan (sHya) to form artificial extracellular matrix (aECM), which mimics the native microenvironment of bone tissue. Further, bone marrow derived human mesenchymal stem cells (hMSCs) were cultured on these moderately conductive (10(-4)10(-3) S/cm) aECM coated PANI substrates and exposed intermittently to pulsed electric field (PEF) generated through transformer-like coupling (TLC) approach over 28 days. On the basis of critical analysis over an array of end points, it was inferred that Coll/sHya coated PANI (PANI/Coll/sHya) substrates had enhanced proliferative capacity of hMSCs up to 28 days in culture, even in the absence of PEF stimulation. On the contrary, the adopted PEF stimulation protocol (7 ms rectangular pulses, 3.6 mV/cm, 10 Hz) is shown to enhance osteogenic differentiation potential of hMSCs. Additionally, PEF stimulated hMSCs had also displayed different morphological characteristics as their nonstimulated counterparts. Concomitantly, earlier onset of ALP activity was also observed on PANI/Coll/sHya substrates and resulted in more calcium deposition. Moreover, real-time polymerase chain reaction results indicated higher mRNA levels of alkaline phosphatase and osteocalcin, whereas the expression of other osteogenic markers such as Runt-related transcription factor 2, Col1A, and osteopontin exhibited a dynamic pattern similar to control cells that are cultured in osteogenic medium. Taken together, our experimental results illustrate the interplay of multiple parameters such as substrate conductivity, electric field stimulation, and aECM coating on the modulation of hMSC proliferation and differentiation in vitro.

Thermal Modeling of Organic Light-Emitting Diode Display Panels

Power densities required to operate active-matrix organic-light-emitting diode (AMOLED) based displays for high luminance applications, lead to temperature rise due to self heating. Temperature rise leads to significant degradation and consequent reduction in life time. In this work numerical techniques based computational fluid dynamics (CFD) is used to determine the temperature rise and its distribution for an AMOLED based display for a given power density and size. Passive cooling option in form of protruded rectangular fins is implemented to reduce the display temperature.