Clan-based evolutionary approach for solving control problems

This paper presents a dan-based evolutionary approach for solving control problems. Three selected control problems, viz. linear-quadratic, harvest, and push-cart problems, are solved using the proposed approach. Results are compared with those of the evolutionary programming (EP) approach. In most of the cases, the proposed approach is successful in obtaining (near) optimal solutions for these selected problems.

Verifying scalings for bending rigidity of bilayer membranes using mesoscale models

The bending rigidity kappa of bilayer membranes was studied with coarse grained soft repulsive potentials using dissipative particle dynamics (DPD) simulations. Using a modified Andersen barostat to maintain the bilayers in a tensionless state, the bending rigidity was obtained from a Fourier analysis of the height fluctuations. From simulations carried out over a wide range of membrane thickness, the continuum scaling relation kappa proportional to d(2) was captured for both the L-alpha and L-beta phases. For membranes with 4 to 6 tail beads, the bending rigidity in the L-beta phase was found to be 10-15 times higher than that observed for the L-alpha phase. From the quadratic scalings obtained, a six fold increase in the area stretch modulus, k(A) was observed across the transition. The magnitude of increase in both kappa and k(A) from the L-alpha to the L-beta phase is consistent with current experimental observations in lipid bilayers and to our knowledge provides for the first time a direct evaluation of the mechanical properties in the L-beta phase.

Optimal blood glucose regulation of diabetic patients using single network adaptive critics

Diabetes is a long-term disease during which the body's production and use of insulin are impaired, causing glucose concentration level to increase in the bloodstream. Regulating blood glucose levels as close to normal as possible leads to a substantial decrease in long-term complications of diabetes. In this paper, an intelligent online feedback-treatment strategy is presented for the control of blood glucose levels in diabetic patients using single network adaptive critic (SNAC) neural networks (which is based on nonlinear optimal control theory). A recently developed mathematical model of the nonlinear dynamics of glucose and insulin interaction in the blood system has been revised and considered for synthesizing the neural network for feedback control. The idea is to replicate the function of pancreatic insulin, i.e. to have a fairly continuous measurement of blood glucose and a situation-dependent insulin injection to the body using an external device. Detailed studies are carried out to analyze the effectiveness of this adaptive critic-based feedback medication strategy. A comparison study with linear quadratic regulator (LQR) theory shows that the proposed nonlinear approach offers some important advantages such as quicker response, avoidance of hypoglycemia problems, etc. Robustness of the proposed approach is also demonstrated from a large number of simulations considering random initial conditions and parametric uncertainties. Copyright (C) 2009 John Wiley & Sons, Ltd.

Linear dynamics and active control of an elastically supported travelling string

In this paper, the linear dynamics and active control of a string travelling with uniform velocity is presented. Discrete elastic supports are introduced along the length of the string. Finite element formulation is adopted to obtain the governing equations of motion. The velocity of translation introduces gyroscopic terms in the system equations. The effect of translation and the discrete elastic supports on the free vibration solution is studied. The solution is utilized in actively controlling the string vibrations due to an initial disturbance. The control, affected in modal space, is optimal with respect to a quadratic performance index. Numerical results are presented to demonstrate the effectiveness of the control strategy in regulating the travelling string vibrations.

Brownian dynamics simulation of dense binary colloidal mixtures. II. Translational and bond-orientational order

We report the Brownian dynamics simulation results on the translational and bond-angle-orientational correlations for charged colloidal binary suspensions as the interparticle interactions are increased to form a crystalline (for a volume fraction phi = 0.2) or a glassy (phi = 0.3) state. The translational order is quantified in terms of the two- and four-point density autocorrelation functions whose comparisons show that there is no growing correlation length near the glass transition. The nearest-neighbor orientational order is determined in terms of the quadratic rotational invariant Q(l) and the bond-orientational correlation functions g(l)(t). The l dependence of Q(l) indicates that icosahedral (l = 6) order predominates at the cost of the cubic order (l = 4) near the glass as well as the crystal transition. The density and orientational correlation functions for a supercooled liquid freezing towards a glass fit well to the streched-exponential form exp[-(t/tau)(beta)]. The average relaxation times extracted from the fitted stretched-exponential functions as a function of effective temperatures T* obey the Arrhenius law for liquids freezing to a crystal whereas these obey the Vogel-Tamman-Fulcher law exp[AT(0)*/(T* - T-0*)] for supercooled Liquids tending towards a glassy state. The value of the parameter A suggests that the colloidal suspensions are ''fragile'' glass formers like the organic and molecular liquids.

First-order hyperpolarizabilities of octupolar aromatic molecules: symmetrically substituted triazines

The first hyperpolarizabilities of some symmetrically substituted triazines have been measured and compared with those of the corresponding symmetrically substituted benzenes. The octupolar triazines have higher quadratic polarizabilities than the corresponding octupolar benzenes. The triazine ring seems to be a better central acceptor than the benzene ring, but if it acts as a donor as in sym-triphenyl triazine, the nonlinearity improves further.

Fluctuation-dissipation relation for nonlinear Langevin equations

It is shown that the fluctuation-dissipation theorem is satisfied by the solutions of a general set of nonlinear Langevin equations with a quadratic free-energy functional (constant susceptibility) and field-dependent kinetic coefficients, provided the kinetic coefficients satisfy the Onsager reciprocal relations for the irreversible terms and the antisymmetry relations for the reversible terms. The analysis employs a perturbation expansion of the nonlinear terms, and a functional integral calculation of the correlation and response functions, and it is shown that the fluctuation-dissipation relation is satisfied at each order in the expansion.

Low-temperature magnetoresistance of ?-phase FexNi80-xCr20 alloys near the critical composition for ferromagnetism

We present a comprehensive study of magnetoresistance (MR) of the crystalline pseudobinary ?-phase Fe alloy series FexNi80-xCr20 (50?x?66). This alloy series shows exotic magnetic phases as the composition (x) is varied. It has a critical composition for ferromagnetism at x=xc?59�60. MR was measured in the temperature range 1.7�110 K and up to a field of 7 T. The observed MR was small and the change was ?1%. The temperature dependence of MR was found to contain a positive and a negative contribution. The positive term was found to be ?H2 and it dominates at high field and high temperatures. We explain this as a manifestation of Kohler�s rule. The negative MR was found to have a quadratic dependence on magnetization M. The magnitude of the negative MR reaches a maximum as x?xc.

A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

Magnetism in MnPSe3: a layered 3d(5) antiferromagnet with unusually large XY anisotropy

The anisotropic magnetic susceptibilities of single crystals of the isostructural layered antiferromagnets, MnPS3 (T-N = 78 K) and MnPSe3 (T-N = 74 K), have been measured as functions of temperature. In both compounds, divalent manganese is present in the high-spin S = 5/2 state. The anisotropies in the susceptibilities of the two are, however, very different; while the susceptibility of MnPS3 is isotropic, that of MnPSe3 shows a large XY anisotropy, unusual for a manganese compound. The anisotropic susceptibilities are described by the zero-field spin Hamiltonian: H = DSiz2 - Sigma J(ij).(S) over right arrow (S) over right arrow(j) with the quadratic single-ion anisotropy term introducing anisotropy in an otherwise isotropic situation. The exchange J and the single-ion zero-field-splitting (ZFS) parameter D were evaluated using the correlated effective-field theory of Lines. For MnPSe3, J/k = -5.29 K and D/k = 26.6 K, while for isotropic MnPS3, J/k = -8.1 K. It is suggested that the large value of the ZFS parameter for MnPSe3 as compared to MnPS3 could be due to the large ligand spin-orbit contribution of the heavier selenium.

Shape optimization of RC flexural members

A natural velocity field method for shape optimization of reinforced concrete (RC) flexural members has been demonstrated. The possibility of shape optimization by modifying the shape of an initially rectangular section, in addition to variation of breadth and depth along the length, has been explored. Necessary shape changes have been computed using the sequential quadratic programming (SQP) technique. Genetic algorithm (Goldberg and Samtani 1986) has been used to optimize the diameter and number of main reinforcement bars. A limit-state design approach has been adopted for the nonprismatic RC sections. Such relevant issues as formulation of optimization problem, finite-element modeling, and solution procedure have been described. Three design examples-a simply supported beam, a cantilever beam, and a two-span continuous beam, all under uniformly distributed loads-have been optimized. The results show a significant savings (40-56%) in material and cost and also result in aesthetically pleasing structures. This procedure will lead to considerable cost saving, particularly in cases of mass-produced precast members and a heavy cast-in-place member such as a bridge girder.

A fixed-grid finite element based enthalpy formulation for generalized phase change problems: role of superficial mushy region

The enthalpy method is primarily developed for studying phase change in a multicomponent material, characterized by a continuous liquid volume fraction (phi(1)) vs temperature (T) relationship. Using the Galerkin finite element method we obtain solutions to the enthalpy formulation for phase change in 1D slabs of pure material, by assuming a superficial phase change region (linear (phi(1) vs T) around the discontinuity at the melting point. Errors between the computed and analytical solutions are evaluated for the fluxes at, and positions of, the freezing front, for different widths of the superficial phase change region and spatial discretizations with linear and quadratic basis functions. For Stefan number (St) varying between 0.1 and 10 the method is relatively insensitive to spatial discretization and widths of the superficial phase change region. Greater sensitivity is observed at St = 0.01, where the variation in the enthalpy is large. In general the width of the superficial phase change region should span at least 2-3 Gauss quadrature points for the enthalpy to be computed accurately. The method is applied to study conventional melting of slabs of frozen brine and ice. Regardless of the forms for the phi(1) vs T relationships, the thawing times were found to scale as the square of the slab thickness. The ability of the method to efficiently capture multiple thawing fronts which may originate at any spatial location within the sample, is illustrated with the microwave thawing of slabs and 2D cylinders. (C) 2002 Elsevier Science Ltd. All rights reserved.

Band Structures of Bilayer Graphene Superlattices

We formulate a low energy effective Hamiltonian to study superlattices in bilayer graphene (BLG) using a minimal model which supports quadratic band touching points. We show that a one dimensional (1D) periodic modulation of the chemical potential or the electric field perpendicular to the layers leads to the generation of zero-energy anisotropic massless Dirac fermions and finite energy Dirac points with tunable velocities. The electric field superlattice maps onto a coupled chain model comprised of ``topological'' edge modes. 2D superlattice modulations are shown to lead to gaps on the mini-Brillouin zone boundary but do not, for certain symmetries, gap out the quadratic band touching point. Such potential variations, induced by impurities and rippling in biased BLG, could lead to subgap modes which are argued to be relevant to understanding transport measurements.

Viscoelastic phenomenology based structure assignment for closed-loop vibration control of a beam with sensors and actuators

In this paper we incorporate a novel approach to synthesize a class of closed-loop feedback control, based on the variational structure assignment. Properties of a viscoelastic system are used to design an active feedback controller for an undamped structural system with distributed sensor, actuator and controller. Wave dispersion properties of onedimensional beam system have been studied. Efficiency of the chosen viscoelastic model in enhancing damping and stability properties of one-dimensional viscoelastic bar have been analyzed. The variational structure is projected on a solution space of a closed-loop system involving a weakly damped structure with distributed sensor and actuator with controller. These assign the phenomenology based internal strain rate damping parameter of a viscoelastic system to the usual elastic structure but with active control. In the formulation a model of cantilever beam with non-collocated actuator and sensor has been considered. The formulation leads to the matrix identification problem of two dynamic stiffness matrices. The method has been simplified to obtain control system gains for the free vibration control of a cantilever beam system with collocated actuator-sensor, using quadratic optimal control and pole-placement methods.

Fast edge splitting and Edmonds' arborescence construction for unweighted graphs

Given an unweighted undirected or directed graph with n vertices, m edges and edge connectivity c, we present a new deterministic algorithm for edge splitting. Our algorithm splits-off any specified subset S of vertices satisfying standard conditions (even degree for the undirected case and in-degree ≥ out-degree for the directed case) while maintaining connectivity c for vertices outside S in Õ(m+nc2) time for an undirected graph and Õ(mc) time for a directed graph. This improves the current best deterministic time bounds due to Gabow [8], who splits-off a single vertex in Õ(nc2+m) time for an undirected graph and Õ(mc) time for a directed graph. Further, for appropriate ranges of n, c, |S| it improves the current best randomized bounds due to Benczúr and Karger [2], who split-off a single vertex in an undirected graph in Õ(n2) Monte Carlo time. We give two applications of our edge splitting algorithms. Our first application is a sub-quadratic (in n) algorithm to construct Edmonds' arborescences. A classical result of Edmonds [5] shows that an unweighted directed graph with c edge-disjoint paths from any particular vertex r to every other vertex has exactly c edge-disjoint arborescences rooted at r. For a c edge connected unweighted undirected graph, the same theorem holds on the digraph obtained by replacing each undirected edge by two directed edges, one in each direction. The current fastest construction of these arborescences by Gabow [7] takes Õ(n2c2) time. Our algorithm takes Õ(nc3+m) time for the undirected case and Õ(nc4+mc) time for the directed case. The second application of our splitting algorithm is a new Steiner edge connectivity algorithm for undirected graphs which matches the best known bound of Õ(nc2 + m) time due to Bhalgat et al [3]. Finally, our algorithm can also be viewed as an alternative proof for existential edge splitting theorems due to Lovász [9] and Mader [11].