Selectivity, pulse fishing and endogenous lifespan in Beverton-Holt models

Optimal management in a multi-cohort Beverton-Holt model with any number of age classes and imperfect selectivity is equivalent to finding the optimal fish lifespan by chosen fallow cycles. Optimal policy differs in two main ways from the optimal lifespan rule with perfect selectivity. First, weight gain is valued in terms of the whole population structure. Second, the cost of waiting is the interest rate adjusted for the increase in the pulse length. This point is especially relevant for assessing the role of selectivity. Imperfect selectivity reduces the optimal lifespan and the optimal pulse length. We illustrate our theoretical findings with a numerical example. Results obtained using global numerical methods select the optimal pulse length predicted by the optimal lifespan rule.

A Mechanical Model For The Quantification Of The Effect Of Laser Quenching On Ctod In Steels

The technology of laser quenching is widely used to improve the surface properties of steels in surface engineering. Generally, laser quenching of steels can lead to two important results. One is the generation of residual stress in the surface layer. In general, the residual stress varies from the surface to the interior along the quenched track depth direction, and the residual stress variation is termed as residual stress gradient effect in this work. The other is the change of mechanical properties of the surface layer, such as the increases of the micro-hardness, resulting from the changes of the microstructure of the surface layer. In this work, a mechanical model of a laser-quenched specimen with a crack in the middle of the quenched layer is developed to quantify the effect of residual stress gradient and the average micro-hardness over the crack length on crack tip opening displacement (CTOD). It is assumed that the crack in the middle of the quenched layer is created after laser quenching, and the crack can be a pre-crack or a defect due to some reasons, such as a void, cavity or a micro-crack. Based on the elastic-plastic fracture mechanics theory and using the relationship between the micro-hardness and yield strength, a concise analytical solution, which can be used to quantify the effect of residual stress gradient and the average micro-hardness over the crack length resulting from laser quenching on CTOD, is obtained. The concise analytical solution obtained in this work, cannot only be used as a means to predict the crack driving force in terms of the CTOD, but also serve as a baseline for further experimental investigation of the effect after laser-quenching treatment on fracture toughness in terms of the critical CTOD of a specimen, accounting for the laser-quenching effect. A numerical example presented in this work shows that the CTOD of the quenched can be significantly decreased in comparison with that of the unquenched. (C) 2008 Elsevier B.V. All rights reserved.

Optimal uncertainty quantification via convex optimization and relaxation

Many engineering applications face the problem of bounding the expected value of a quantity of interest (performance, risk, cost, etc.) that depends on stochastic uncertainties whose probability distribution is not known exactly. Optimal uncertainty quantification (OUQ) is a framework that aims at obtaining the best bound in these situations by explicitly incorporating available information about the distribution. Unfortunately, this often leads to non-convex optimization problems that are numerically expensive to solve.

This thesis emphasizes on efficient numerical algorithms for OUQ problems. It begins by investigating several classes of OUQ problems that can be reformulated as convex optimization problems. Conditions on the objective function and information constraints under which a convex formulation exists are presented. Since the size of the optimization problem can become quite large, solutions for scaling up are also discussed. Finally, the capability of analyzing a practical system through such convex formulations is demonstrated by a numerical example of energy storage placement in power grids.

When an equivalent convex formulation is unavailable, it is possible to find a convex problem that provides a meaningful bound for the original problem, also known as a convex relaxation. As an example, the thesis investigates the setting used in Hoeffding's inequality. The naive formulation requires solving a collection of non-convex polynomial optimization problems whose number grows doubly exponentially. After structures such as symmetry are exploited, it is shown that both the number and the size of the polynomial optimization problems can be reduced significantly. Each polynomial optimization problem is then bounded by its convex relaxation using sums-of-squares. These bounds are found to be tight in all the numerical examples tested in the thesis and are significantly better than Hoeffding's bounds.

Environmental taxes and the Double dividend hypothesis

In this End of Degree Dissertation I built a basic economic model, where two perfectly competitive markets interact, as a tool to understand, illustrate and evaluate environmentalfiscal reforms that have already been implemented or have been announced in many developed countries. Thus, in this dissertation I try to explain as simply as possible the theorical aspects of the “double dividend” hyphotesis and provide a numerical example with illustrative purpose.

Design of anisotropic reflector with birefringent thin films

A novel design for dielectric anisotropic mirrors with birefringent thin films for normal incidence is presented. This mirror consists of a stack of quarter-wave biaxial layers. The biaxial anisotropic layers can be fabricated by oblique deposition. The reflectance is different for two linear polarizations of light incidence on the mirrors. As a numerical example, the design is carried out on glass with TiO2 and ZrO2. These thin films could be applied to anisotropic reflective devices for lasers.

Propagation properties of beams generated by Gaussian mirror resonator in uniaxial crystals

Based on the paraxial vectorial theory of beams propagating in uniaxially anisotropic media, we have derived the analytical propagation equations of beams generated by Gaussian mirror resonator (GMR) in uniaxial crystals, and given the typical numerical example to illustrate our analytical results. Due to the anisotropy crystals, the ordinary and extraordinary beams originated by incident beams generated by GMR propagate with different diffraction lengths, thus the linear polarization state and axial symmetry of the incident beams generated by GMR do not remain during propagating in crystals. (c) 2006 Elsevier GmbH. All rights reserved.

Propagation and polarization properties of hollow Gaussian beams in uniaxial crystals

Based on the paraxial vectorial theory of beams propagating in uniaxially anisotropic media, we have derived the analytical propagation equations of hollow Gaussian beams (HGBs) in uniaxial crystals, and given the typical numerical example to illustrate our analytical results. Due to the anisotropy crystals, the ordinary and extraordinary beams originated by incident HGBs propagate with different diffraction lengths, thus the linear polarization state and axial symmetry of incident HGBs do not remain during propagating in crystals. (c) 2007 Published by Elsevier B.V.

Towards identification of a general model of damping

Characterization of damping forces in a vibrating structure has long been an active area of research in structural dynamics. In spite of a large amount of research, understanding of damping mechanisms is not well developed. A major reason for this is that unlike inertia and stiffness forces it is not in general clear what are the state variables that govern the damping forces. The most common approach is to use `viscous damping' where the instantaneous generalized velocities are the only relevant state variables. However, viscous damping by no means the only damping model within the scope of linear analysis. Any model which makes the energy dissipation functional non-negative is a possible candidate for a valid damping model. This paper is devoted to develop methodologies for identification of such general damping models responsible for energy dissipation in a vibrating structure. The method uses experimentally identified complex modes and complex natural frequencies and does not a-priori assume any fixed damping model (eg., viscous damping) but seeks to determine parameters of a general damping model described by the so called `relaxation function'. The proposed method and several related issues are discussed by considering a numerical example of a linear array of damped spring-mass oscillators.

An analysis of the cooperative mechano-sensitive feedback between intracellular signaling, focal adhesion development, and stress fiber contractility

Cells communicate with their external environment via focal adhesions and generate activation signals that in turn trigger the activity of the intracellular contractile machinery. These signals can be triggered by mechanical loading that gives rise to a cooperative feedback loop among signaling, focal adhesion formation, and cytoskeletal contractility, which in turn equilibrates with the applied mechanical loads. We devise a signaling model that couples stress fiber contractility and mechano-sensitive focal adhesion models to complete this above mentioned feedback loop. The signaling model is based on a biochemical pathway where IP3 molecules are generated when focal adhesions grow. These IP3 molecules diffuse through the cytosol leading to the opening of ion channels that disgorge Ca2+ from the endoplasmic reticulum leading to the activation of the actin/myosin contractile machinery. A simple numerical example is presented where a one-dimensional cell adhered to a rigid substrate is pulled at one end, and the evolution of the stress fiber activation signal, stress fiber concentrations, and focal adhesion distributions are investigated. We demonstrate that while it is sufficient to approximate the activation signal as spatially uniform due to the rapid diffusion of the IP3 through the cytosol, the level of the activation signal is sensitive to the rate of application of the mechanical loads. This suggests that ad hoc signaling models may not be able to capture the mechanical response of cells to a wide range of mechanical loading events. © 2011 American Society of Mechanical Engineers.

Application of the fixed point method for solution in time stepping finite element analysis using the inverse vector Jiles-Atherton model

An implementation of the inverse vector Jiles-Atherton model for the solution of non-linear hysteretic finite element problems is presented. The implementation applies the fixed point method with differential reluctivity values obtained from the Jiles-Atherton model. Differential reluctivities are usually computed using numerical differentiation, which is ill-posed and amplifies small perturbations causing large sudden increases or decreases of differential reluctivity values, which may cause numerical problems. A rule based algorithm for conditioning differential reluctivity values is presented. Unwanted perturbations on the computed differential reluctivity values are eliminated or reduced with the aim to guarantee convergence. Details of the algorithm are presented together with an evaluation of the algorithm by a numerical example. The algorithm is shown to guarantee convergence, although the rate of convergence depends on the choice of algorithm parameters. © 2011 IEEE.

Eccentricity effects on acoustic radiation from a spherical source suspended within a thermoviscous fluid sphere.

Acoustic radiation from a spherical source undergoing angularly periodic axisymmetric harmonic surface vibrations while eccentrically suspended within a thermoviscous fluid sphere, which is immersed in a viscous thermally conducting unbounded fluid medium, is analyzed in an exact fashion. The formulation uses the appropriate wave-harmonic field expansions along with the translational addition theorem for spherical wave functions and the relevant boundary conditions to develop a closed-form solution in form of infinite series. The analytical results are illustrated with a numerical example in which the vibrating source is eccentrically positioned within a chemical fluid sphere submerged in water. The modal acoustic radiation impedance load on the source and the radiated far-field pressure are evaluated and discussed for representative values of the parameters characterizing the system. The proposed model can lead to a better understanding of dynamic response of an underwater acoustic lens. It is equally applicable in miniature transducer analysis and design with applications in medical ultrasonics.

Air traffic control with an expected value criterion

in this contribution we discuss a stochastic framework for air traffic conflict resolution. The conflict resolution task is posed as the problem of optimizing an expected value criterion. Optimization is carried out by Monte Carlo Markov Chain (MCMC) simulation. A numerical example illustrates the proposed strategy. Copyright © 2005 IFAC.

Mathematical relationships between representations of structure in linear interconnected dynamical systems

A dynamical system can exhibit structure on multiple levels. Different system representations can capture different elements of a dynamical system's structure. We consider LTI input-output dynamical systems and present four representations of structure: complete computational structure, subsystem structure, signal structure, and input output sparsity structure. We then explore some of the mathematical relationships that relate these different representations of structure. In particular, we show that signal and subsystem structure are fundamentally different ways of representing system structure. A signal structure does not always specify a unique subsystem structure nor does subsystem structure always specify a unique signal structure. We illustrate these concepts with a numerical example. © 2011 AACC American Automatic Control Council.

Designing output-feedback predictive controllers by reverse-engineering existing LTI controllers

An approach to designing a constrained output-feedback predictive controller that has the same small-signal properties as a pre-existing output-feedback linear time invariant controller is proposed. Systematic guidelines are proposed to select an appropriate (non-unique) realization of the resulting state observer. A method is proposed to transform a class of offset-free reference tracking controllers into the combination of an observer, steady-state target calculator and predictive controller. The procedure is demonstrated with a numerical example. © 2013 IEEE.

Reconfigurable predictive control for redundantly actuated systems with parameterised input constraints

A method is proposed for on-line reconfiguration of the terminal constraint used to provide theoretical nominal stability guarantees in linear model predictive control (MPC). By parameterising the terminal constraint, its complete reconstruction is avoided when input constraints are modified to accommodate faults. To enlarge the region of feasibility of the terminal control law for a certain class of input faults with redundantly actuated plants, the linear terminal controller is defined in terms of virtual commands. A suitable terminal cost weighting for the reconfigurable MPC is obtained by means of an upper bound on the cost for all feasible realisations of the virtual commands from the terminal controller. Conditions are proposed that guarantee feasibility recovery for a defined subset of faults. The proposed method is demonstrated by means of a numerical example. © 2013 Elsevier B.V. All rights reserved.