Optimal design of Functionally Graded Materials using a procedural model and Particle Swarm Optimization

A new method for the optimal design of Functionally Graded Materials (FGM) is proposed in this paper. Instead of using the widely used explicit functional models, a feature tree based procedural model is proposed to represent generic material heterogeneities. A procedural model of this sort allows more than one explicit function to be incorporated to describe versatile material gradations and the material composition at a given location is no longer computed by simple evaluation of an analytic function, but obtained by execution of customizable procedures. This enables generic and diverse types of material variations to be represented, and most importantly, by a reasonably small number of design variables. The descriptive flexibility in the material heterogeneity formulation as well as the low dimensionality of the design vectors help facilitate the optimal design of functionally graded materials. Using the nature-inspired Particle Swarm Optimization (PSO) method, functionally graded materials with generic distributions can be efficiently optimized. We demonstrate, for the first time, that a PSO based optimizer outperforms classical mathematical programming based methods, such as active set and trust region algorithms, in the optimal design of functionally graded materials. The underlying reason for this performance boost is also elucidated with the help of benchmarked examples. © 2011 Elsevier Ltd. All rights reserved.

基于GIS的区域土壤侵蚀模型设计

区域土壤侵蚀模型是国家和区域土壤侵蚀调查、水土保持宏观规划的支持工具。借鉴国外区域土壤侵蚀模型,以DEM栅格为空间单元,对区域土壤侵蚀的单元模型进行定量表达,包括降雨、植被截留、入渗,微地形填洼等,并利用GIS功能完成径流传递和汇集部分的计算,建立了区域土壤侵蚀模型。所建立模型在延河流域的试运算,结果接近现实,能反映土壤侵蚀时空分布趋势。模型的建立可为水土保持的宏观决策提供支持。

A four-phase confocal elliptical cylinder model for predicting the effective thermal conductivity of coated fibre composites

A four-phase confocal elliptical cylinder model is proposed from which a generalised self-consistent method is developed for predicting the thermal conductivity of coated fibre reinforced composites. The method can account for the influence of the fibre section shape ratio on conductivity, and the physical reasonableness of the model is demonstrated by using the fibre distribution function. An exact solution is obtained for thermal conductivity by applying conformal mapping and Laurent series expansion techniques of the analytic function. The solution to the three-phase confocal elliptical model, which simulates composites with idealised fibre-matrix interfaces, is arrived at as the degenerated case. A comparison with other available micromechanics methods, Hashin and Shtrikman's bounds and experimental data shows that the present method provides convergent and reasonable results for a full range of variations in fibre section shapes and for a complete spectrum of the fibre volume fraction. Numerical results show the dependence of the effective conductivities of composites on the aspect ratio of coated fibres and demonstrate that a coating is effective in enhancing the thermal transport property of a composite. The present solutions are helpful to analysis and design of composites.

Boundary triplets and M-functions for non-selfadjoint operators, with applications to elliptic PDEs and block operator matrices

B.M. Brown, M. Marletta, S. Naboko, I. Wood: Boundary triplets and M-functions for non-selfadjoint operators, with applications to elliptic PDEs and block operator matrices, J. London Math. Soc., June 2008; 77: 700-718. The full text of this article will be made available in this repository in June 2009 Sponsorship: EPSRC,INTAS

A Hamiltonian approach to simulation of acoustic systems involving nonlinear interactions

Nonlinear interactions take place in most systems that arise in music acoustics, usually as a result of player-instrument coupling. Several time-stepping methods exist for the numerical simulation of such systems. These methods generally involve the discretization of the Newtonian description of the system. However, it is not always possible to prove the stability of the resulting algorithms, especially when dealing with systems where the underlying force is a non-analytic function of the phase space variables. On the other hand, if the discretization is carried out on the Hamiltonian description of the system, it is possible to prove the stability of the derived numerical schemes. This Hamiltonian approach is applied to a series of test models of single or multiple nonlinear collisions and the energetic properties of the derived schemes are discussed. After establishing that the schemes respect the principle of conservation of energy, a nonlinear single-reed model is formulated and coupled to a digital bore, in order to synthesize clarinet-like sounds.

Discrete commutative difference operator theory

The object of this thesis is to formulate a basic commutative difference operator theory for functions defined on a basic sequence, and a bibasic commutative difference operator theory for functions defined on a bibasic sequence of points, which can be applied to the solution of basic and bibasic difference equations. in this thesis a brief survey of the work done in this field in the classical case, as well as a review of the development of q~difference equations, q—analytic function theory, bibasic analytic function theory, bianalytic function theory, discrete pseudoanalytic function theory and finally a summary of results of this thesis

Studies in the geometry of the discrete plane

This thesis is an attempt to initiate the development of a discrete geometry of the discrete plane H = {(qmxo,qnyo); m,n e Z - the set of integers}, where q s (0,1) is fixed and (xO,yO) is a fixed point in the first quadrant of the complex plane, xo,y0 ¢ 0. The discrete plane was first considered by Harman in 1972, to evolve a discrete analytic function theory for geometric difference functions. We shall mention briefly, through various sections, the principle of discretization, an outline of discrete a alytic function theory, the concept of geometry of space and also summary of work done in this thesis

Numerical conformal mapping and mesh generation for polygonal and multiply-connected regions

Details are given of a boundary-fitted mesh generation method for use in modelling free surface flow and water quality. A numerical method has been developed for generating conformal meshes for curvilinear polygonal and multiply-connected regions. The method is based on the Cauchy-Riemann conditions for the analytic function and is able to map a curvilinear polygonal region directly onto a regular polygonal region, with horizontal and vertical sides. A set of equations have been derived for determining the lengths of these sides and the least-squares method has been used in solving the equations. Several numerical examples are presented to illustrate the method.

Invariants of binary differential equations

In this paper, we study binary differential equations a(x, y)dy (2) + 2b(x, y) dx dy + c(x, y)dx (2) = 0, where a, b, and c are real analytic functions. Following the geometric approach of Bruce and Tari in their work on multiplicity of implicit differential equations, we introduce a definition of the index for this class of equations that coincides with the classical Hopf`s definition for positive binary differential equations. Our results also apply to implicit differential equations F(x, y, p) = 0, where F is an analytic function, p = dy/dx, F (p) = 0, and F (pp) not equal aEuro parts per thousand 0 at the singular point. For these equations, we relate the index of the equation at the singular point with the index of the gradient of F and index of the 1-form omega = dy -aEuro parts per thousand pdx defined on the singular surface F = 0.

A Time Efficient Optical Model for GATE Simulation of a LYSO Scintillation Matrix Used in PET Applications

A time efficient optical model is proposed for GATE simulation of a LYSO scintillation matrix coupled to a photomultiplier. The purpose is to avoid the excessively long computation time when activating the optical processes in GATE. The usefulness of the model is demonstrated by comparing the simulated and experimental energy spectra obtained with the dual planar head equipment for dosimetry with a positron emission tomograph ( DoPET). The procedure to apply the model is divided in two steps. Firstly, a simplified simulation of a single crystal element of DoPET is used to fit an analytic function that models the optical attenuation inside the crystal. In a second step, the model is employed to calculate the influence of this attenuation in the energy registered by the tomograph. The use of the proposed optical model is around three orders of magnitude faster than a GATE simulation with optical processes enabled. A good agreement was found between the experimental and simulated data using the optical model. The results indicate that optical interactions inside the crystal elements play an important role on the energy resolution and induce a considerable degradation of the spectra information acquired by DoPET. Finally, the same approach employed by the proposed optical model could be useful to simulate a scintillation matrix coupled to a photomultiplier using single or dual readout scheme.

Caracterização analítica e geométrica da metodologia geral de determinação de distribuições de Weibull para o regime eólico e suas aplicações

O regime eólico de uma região pode ser descrito por distribuição de frequências que fornecem informações e características extremamente necessárias para uma possível implantação de sistemas eólicos de captação de energia na região e consequentes aplicações no meio rural em regiões afastadas. Estas características, tais como a velocidade média anual, a variância das velocidades registradas e a densidade da potência eólica média horária, podem ser obtidas pela frequência de ocorrências de determinada velocidade, que por sua vez deve ser estudada através de expressões analíticas. A função analítica mais adequada para distribuições eólicas é a função de densidade de Weibull, que pode ser determinada por métodos numéricos e regressões lineares. O objetivo deste trabalho é caracterizar analítica e geometricamente todos os procedimentos metodológicos necessários para a realização de uma caracterização completa do regime eólico de uma região e suas aplicações na região de Botucatu - SP, visando a determinar o potencial energético para implementação de turbinas eólicas. Assim, foi possível estabelecer teoremas relacionados com a forma de caracterização do regime eólico, estabelecendo a metodologia concisa analiticamente para a definição dos parâmetros eólicos de qualquer região a ser estudada. Para o desenvolvimento desta pesquisa, utilizou-se um anemômetro da CAMPBELL.

Homoclinic bifurcations in reversible Hamiltonian systems

We study the existence of homoclic solutions for reversible Hamiltonian systems taking the family of differential equations u(iv) + au - u +f(u, b) = 0 as a model, where fis an analytic function and a, b real parameters. These equations are important in several physical situations such as solitons and in the existence of finite energy stationary states of partial differential equations, but no assumptions of any kind of discrete symmetry is made and the analysis here developed can be extended to others Hamiltonian systems and successfully employed in situations where standard methods fail. We reduce the problem of computing these orbits to that of finding the intersection of the unstable manifold with a suitable set and then apply it to concrete situations. We also plot the homoclinic values configuration in parameters space, giving a picture of the structural distribution and a geometrical view of homoclinic bifurcations. (c) 2005 Published by Elsevier B.V.

Modelagem dinâmica de transformadores de corrente para uso em estudos de qualidade da energia e proteção do sistema elétrico

Pós-graduação em Engenharia Elétrica - FEB

Frações contínuas que correspondem a séries de potências em dois pontos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

On the Luminous and Dark Matter Distribution in Early-Type Galaxies

The way mass is distributed in galaxies plays a major role in shaping their evolution across cosmic time. The galaxy's total mass is usually determined by tracing the motion of stars in its potential, which can be probed observationally by measuring stellar spectra at different distances from the galactic centre, whose kinematics is used to constrain dynamical models. A class of such models, commonly used to accurately determine the distribution of luminous and dark matter in galaxies, is that of equilibrium models. In this Thesis, a novel approach to the design of equilibrium dynamical models, in which the distribution function is an analytic function of the action integrals, is presented. Axisymmetric and rotating models are used to explain observations of a sample of nearby early-type galaxies in the Calar Alto Legacy Integral Field Area survey. Photometric and spectroscopic data for round and flattened galaxies are well fitted by the models, which are then used to get the galaxies' total mass distribution and orbital anisotropy. The time evolution of massive early-type galaxies is also investigated with numerical models. Their structural properties (mass, size, velocity dispersion) are observed to evolve, on average, with redshift. In particular, they appear to be significantly more compact at higher redshift, at fixed stellar mass, so it is interesting to investigate what drives such evolution. This Thesis focuses on the role played by dark-matter haloes: their mass-size and mass-velocity dispersion correlations evolve similarly to the analogous correlations of ellipticals; at fixed halo mass, the haloes are more compact at higher redshift, similarly to massive galaxies; a simple model, in which all the galaxy's size and velocity-dispersion evolution is due to the cosmological evolution of the underlying halo population, reproduces the observed size and velocity-dispersion of massive compact early-type galaxies up to redshift of about 2.