Active vibration suppression of beams with discrete actuators using optimal dynamic inversion

Most of the structural elements like beams, cables etc. are flexible and should be modeled as distributed parameter systems (DPS) to represent the reality better. For large structures, the usual approach of 'modal representation' is not an accurate representation. Moreover, for excessive vibrations (possibly due to strong wind, earthquake etc.), external power source (controller) is needed to suppress it, as the natural damping of these structures is usually small. In this paper, we propose to use a recently developed optinial dynamic inversion technique to design a set of discrete controllers for this purpose. We assume that the control force to the structure is applied through finite number of actuators, which are located at predefined locations in the spatial domain. The method used in this paper determines control forces directly from the partial differential equation (PDE) model of the system. The formulation has better practical significance, both because it leads to a closed form solution of the controller (hence avoids computational issues) as well as because a set of discrete actuators along the spatial domain can be implemented with relative ease (as compared to a continuous actuator).

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Intermolecular interactions in pi -conjugated systems: Application to polyenes

A computational scheme has been developed for strongly interacting systems wherein the intermolecular interaction is introduced as a charge-induced-dipole term. Within this approximation, the model Hamiltonian is exactly solved using a valence-bond basis. The validity of the scheme has been checked by use of exact calculations on small model systems. The method has been applied to finite polyenes to study the shifts in the ground-state energies and dipole-allowed excited-state energies in the presence of neighbors. Our calculations show a red shift in the optical gap of the infinite polyene by 0.124 eV, which is rather small compared to the experimental red shift. This is traced to the larger inaccuracy in the calculated shift in the excited state. The calculated shift in the ground-state energies are more accurate and hence the method is better suited for studying the effect of intermolecular interactions on the properties of the ground state.

Photophysics of conjugated polymers: interplay between Forster energy migration and defect concentration in shaping a photochemical funnel in PPV

Recent single molecule experiments have suggested the existence of a photochemical funnel in the photophysics of conjugated polymers, like poly[2-methoxy-5-(2'-ethylhexyl)oxy-1,4-phenylenevinylene] (MEH-PPV). The funnel is believed to be a consequence of the presence of conformational or chemical defects along the polymer chain and efficient non-radiative energy transfer among different chromophore segments. Here we address the effect of the excitation energy dynamics on the photophysics of PPV. The PPV chain is modeled as a polymer with the length distribution of chromophores given either by a Gaussian or by a Poisson distribution. We observe that the Poisson distribution of the segment lengths explains the photophysics of PPV better than the Gaussian distribution. A recently proposed version of an extended particle-in-a-box' model is used to calculate the exciton energies and the transition dipole moments of the chromophores, and a master equation to describe the excitation energy transfer among different chromophores. The rate of energy transfer is assumed to be given here, as a first approximation, by the well-known Forster expression. The observed excitation population dynamics confirms the photochemical funneling of excitation energy from shorter to longer chromophores of the polymer chain. The time scale of spectral shift and energy transfer for our model polymer, with realistic values of optical parameters, is in the range of 200-300 ps. We find that the excitation energy may not always migrate towards the longest chromophore segments in the polymer chain as the efficiency of energy transfer between chromophores depends on the separation distance between the two and their relative orientation.

Dynamical cluster approximation: Nonlocal dynamics of correlated electron systems

We recently introduced the dynamical cluster approximation (DCA), a technique that includes short-ranged dynamical correlations in addition to the local dynamics of the dynamical mean-field approximation while preserving causality. The technique is based on an iterative self-consistency scheme on a finite-size periodic cluster. The dynamical mean-field approximation (exact result) is obtained by taking the cluster to a single site (the thermodynamic limit). Here, we provide details of our method, explicitly show that it is causal, systematic, Phi derivable, and that it becomes conserving as the cluster size increases. We demonstrate the DCA by applying it to a quantum Monte Carlo and exact enumeration study of the two-dimensional Falicov-Kimball model. The resulting spectral functions preserve causality, and the spectra and the charge-density-wave transition temperature converge quickly and systematically to the thermodynamic limit as the cluster size increases.

Polarization relaxation, dielectric dispersion, and solvation dynamics in dense dipolar liquid

A unified treatment of polarization relaxation, dielectric dispersion and solvation dynamics in a dense, dipolar liquid is presented. It is shown that the information of solvent polarization relaxation that is obtained by macroscopic dielectric dispersion experiments is not sufficient to understand dynamics of solvation of a newly created ion or dipole. In solvation, a significant contribution comes from intermediate wave vector processes which depend critically on the short range (nearest‐neighbor) spatial and orientational order that are present in a dense, dipolar liquid. An analytic expression is obtained for the time dependent solvation energy that depends, in addition to the translational and rotational diffusion coefficients of the liquid, on the ratio of solute–solvent molecular sizes and on the microscopic structure of the polar liquid. Mean spherical approximation (MSA) theory is used to obtain numerical results for polarization relaxation, for wave vector and frequency dependent dielectric function and for time dependent solvation energy. We find that in the absence of translational contribution, the solvation of an ion is, in general, nonexponential. In this case, the short time decay is dominated by the longitudinal relaxation time but the long time decay is dominated by much slower large wave vector processes involving nearest‐neighbor molecules. The presence of a significant translational contribution drastically alters the decay behavior. Now, the long‐time behavior is given by the longitudinal relaxation time constant and the short time dynamics is controlled by the large wave vector processes. Thus, although the continuum model itself is conceptually wrong, a continuum model like result is recovered in the presence of a sizeable translational contribution. The continuum model result is also recovered in the limit of large solute to solvent size ratio. In the opposite limit of small solute size, the decay is markedly nonexponential (if the translational contribution is not very large) and a complete breakdown of the continuum model takes place. The significance of these results is discussed.

Solvation of an ion and of a dipole in a dipolar liquid: How different are the dynamics

Theoretical expressions for the time-dependent solvation energy of an ion and of a dipole in a dense dipolar liquid are derived from microscopic considerations. We show that in contradiction to the prediction of the continuum models, the dynamics of these two species are significantly different from each other. Especially, the zero wavevector contribution, which is significant for ions, is totally absent for dipoles. Dipolar solvation may be profoundly influenced by the translational modes of the host solvent.

Optimal infinite-horizon feedback laws for a general class of constrained discrete-time systems: stability and moving-horizon approximations

Stability results are given for a class of feedback systems arising from the regulation of time-varying discrete-time systems using optimal infinite-horizon and moving-horizon feedback laws. The class is characterized by joint constraints on the state and the control, a general nonlinear cost function and nonlinear equations of motion possessing two special properties. It is shown that weak conditions on the cost function and the constraints are sufficient to guarantee uniform asymptotic stability of both the optimal infinite-horizon and movinghorizon feedback systems. The infinite-horizon cost associated with the moving-horizon feedback law approaches the optimal infinite-horizon cost as the moving horizon is extended.

Approximation of Internal Refractive Index Variation Improves Image Guided Diffuse Optical Tomography of Breast

Effective usage of image guidance by incorporating the refractive index (RI) variation in computational modeling of light propagation in tissue is investigated to assess its impact on optical-property estimation. With the aid of realistic patient breast three-dimensional models, the variation in RI for different regions of tissue under investigation is shown to influence the estimation of optical properties in image-guided diffuse optical tomography (IG-DOT) using numerical simulations. It is also shown that by assuming identical RI for all regions of tissue would lead to erroneous estimation of optical properties. The a priori knowledge of the RI for the segmented regions of tissue in IG-DOT, which is difficult to obtain for the in vivo cases, leads to more accurate estimates of optical properties. Even inclusion of approximated RI values, obtained from the literature, for the regions of tissue resulted in better estimates of optical properties, with values comparable to that of having the correct knowledge of RI for different regions of tissue.

Unified Galerkin- and DAE-Based Approximation of Fractional Order Systems

We consider numerical solutions of nonlinear multiterm fractional integrodifferential equations, where the order of the highest derivative is fractional and positive but is otherwise arbitrary. Here, we extend and unify our previous work, where a Galerkin method was developed for efficiently approximating fractional order operators and where elements of the present differential algebraic equation (DAE) formulation were introduced. The DAE system developed here for arbitrary orders of the fractional derivative includes an added block of equations for each fractional order operator, as well as forcing terms arising from nonzero initial conditions. We motivate and explain the structure of the DAE in detail. We explain how nonzero initial conditions should be incorporated within the approximation. We point out that our approach approximates the system and not a specific solution. Consequently, some questions not easily accessible to solvers of initial value problems, such as stability analyses, can be tackled using our approach. Numerical examples show excellent accuracy. DOI: 10.1115/1.4002516]

Influence of an additional discrete elastic support on the stability of a Leipholz column

The stability characteristics of a conservatively loaded structure are expected to improve if additional supports are provided to the structure. The same, however, may not be said of a non-conservatively loaded structure; several factors, such as the location and stiffness of supports, type of structure and loading, have a significant influence on the stability characteristics. The influence of an arbitrarily located elastic support on the stability characteristics of a Leipholz column is examined.