Convection-diffusion in MHD: Comparison of primitive variable and vector potential solutions of the magnetic induction equation

Abstract not available

Parallelisation of a 3-D electromagnetic scattering code

A high-frequency time domain finite element scattering code using a combination of edge and piecewise constant elements on unstructured tetrahedral meshes is described. A comparison of computation with theory is given for scattering from a sphere. A parallel implementation making use of the bulk synchronous parallel (BSP) programming model is described in detail; a BSP performance model of the parallelized field calculation is derived and compared to timing measurements on up to 128 processors on a Cray T3D.

Optimal static and sequential design : a critical review

The aim of this thesis is to review and augment the theory and methods of optimal experimental design. In Chapter I the scene is set by considering the possible aims of an experimenter prior to an experiment, the statistical methods one might use to achieve those aims and how experimental design might aid this procedure. It is indicated that, given a criterion for design, a priori optimal design will only be possible in certain instances and, otherwise, some form of sequential procedure would seem to be indicated. In Chapter 2 an exact experimental design problem is formulated mathematically and is compared with its continuous analogue. Motivation is provided for the solution of this continuous problem, and the remainder of the chapter concerns this problem. A necessary and sufficient condition for optimality of a design measure is given. Problems which might arise in testing this condition are discussed, in particular with respect to possible non-differentiability of the criterion function at the design being tested. Several examples are given of optimal designs which may be found analytically and which illustrate the points discussed earlier in the chapter. In Chapter 3 numerical methods of solution of the continuous optimal design problem are reviewed. A new algorithm is presented with illustrations of how it should be used in practice. It is shown that, for reasonably large sample size, continuously optimal designs may be approximated to well by an exact design. In situations where this is not satisfactory algorithms for improvement of this design are reviewed. Chapter 4 consists of a discussion of sequentially designed experiments, with regard to both the philosophies underlying, and the application of the methods of, statistical inference. In Chapter 5 we criticise constructively previous suggestions for fully sequential design procedures. Alternative suggestions are made along with conjectures as to how these might improve performance. Chapter 6 presents a simulation study, the aim of which is to investigate the conjectures of Chapter 5. The results of this study provide empirical support for these conjectures. In Chapter 7 examples are analysed. These suggest aids to sequential experimentation by means of reduction of the dimension of the design space and the possibility of experimenting semi-sequentially. Further examples are considered which stress the importance of the use of prior information in situations of this type. Finally we consider the design of experiments when semi-sequential experimentation is mandatory because of the necessity of taking batches of observations at the same time. In Chapter 8 we look at some of the assumptions which have been made and indicate what may go wrong where these assumptions no longer hold.

Fractional differential equations with dependence on the Caputo-Katugampola derivative

In this paper we present a new type of fractional operator, the Caputo–Katugampola derivative. The Caputo and the Caputo–Hadamard fractional derivatives are special cases of this new operator. An existence and uniqueness theorem for a fractional Cauchy type problem, with dependence on the Caputo–Katugampola derivative, is proven. A decomposition formula for the Caputo–Katugampola derivative is obtained. This formula allows us to provide a simple numerical procedure to solve the fractional differential equation.

Dimensionamento de uma estrutura de contenção de terras flexível com vários níveis de apoio

O presente trabalho, no âmbito de projeto final de curso de metrado em Engenharia da Construção, teve como objetivo o estudo do comportamento de estruturas de suporte de terras flexíveis multi-apoiadas (com diferentes tipos de apoio) para dois tipos solos homogéneos. Recorreu-se às teorias clássicas, como a de Rankine, desenvolvidas para estruturas de suporte de terras rígidas. Às teorias semi-empíricas de Terzaghi & Peck que culminaram nos diagramas de Terzaghi & Peck. Apesar de os digramas de Terzaghi & Peck serem diagramas de pressões de terras a usar em estruturas de suporte de terras flexíveis, apresentam algumas limitações importantes, como a sua aplicação apenas em solos heterogéneos, com presença ou não de níveis freáticos, e sem fornecer distribuição das pressões de terras na zona passiva (zona enterrada). Como na atualidade os modelos de elementos finitos permitem simular de modo muito mais rigoroso os problemas da engenharia. O presente trabalho esteve focado em analisar um caso prático em diferentes solos e com diferentes tipos de apoios. Será estudado mediante os métodos analíticos usando as teorias clássicas e posteriormente métodos numéricos (com diferentes programas de cálculo). Finalmente serão comparados os resultados obtidos mediante os diferentes métodos usados. As estruturas foram inicialmente pré-dimensionadas usando os métodos clássicos. Assim foram usados os diagramas de pressões de terras de Terzaghi & Peck para a zona ativa (zona em escavação) e a teoria de Rankine para conhecer as pressões de terras na zona enterrada da cortina (parede moldada) e recorrendo ao software Ftool para a obtenção dos parâmetros de dimensionamento de estruturas de suporte de terras objeto de estudo. Posteriormente utilizaram-se os programas de cálculo automático CYPE 2015 k, e o programa de cálculo de elementos finitos PLAXIS Introductory 2010. Estes programas permitem simular o faseamento construtivo do muro. Para estudar a influência de algúns parâmetros no comportamento da Resumo IV cortina o estudo foi realizado com dois solos distintos, um solo argiloso mole e um solo arenoso denso. Assim como para dois tipos de apoios distintos, ancoragens ativas e escoras passivas. Foram analisados diferentes parâmetros na estrutura de suporte; pressões horizontais das terras, deslocamentos horizontais, esforço axial, transverso e momento fletor.

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.

Mathematical analysis of maximum power generated by photovoltaic systems and fitting curves for standard test conditions

The rural electrification is characterized by geographical dispersion of the population, low consumption, high investment by consumers and high cost. Moreover, solar radiation constitutes an inexhaustible source of energy and in its conversion into electricity photovoltaic panels are used. In this study, equations were adjusted to field conditions presented by the manufacturer for current and power of small photovoltaic systems. The mathematical analysis was performed on the photovoltaic rural system I- 100 from ISOFOTON, with power 300 Wp, located at the Experimental Farm Lageado of FCA/UNESP. For the development of such equations, the circuitry of photovoltaic cells has been studied to apply iterative numerical methods for the determination of electrical parameters and possible errors in the appropriate equations in the literature to reality. Therefore, a simulation of a photovoltaic panel was proposed through mathematical equations that were adjusted according to the data of local radiation. The results have presented equations that provide real answers to the user and may assist in the design of these systems, once calculated that the maximum power limit ensures a supply of energy generated. This real sizing helps establishing the possible applications of solar energy to the rural producer and informing the real possibilities of generating electricity from the sun.

Directional moisture transfer through a wild silkworm cocoon wall

A silkworm cocoon is a porous biological structure with multiple protective functions. In the current work, the authors have used both experimental and numerical methods to reveal the unique moisture transfer characteristics through a wild Antheraea pernyi silkworm cocoon wall, in comparison with the long-domesticated Bombyx mori silkworm cocoon walls. The water vapor transmission and water vapor permeability (WVP) properties show that the A. pernyi cocoons exhibit directional moisture transfer behavior, with easier moisture transfer from inside out than outside in [e.g., the average WVP is 0.057 g/(h m bar) from inside out and is 0.034 g/(h m bar) from outside in]. Numerical analysis shows that the cubic mineral crystals in the outer section of the A. pernyi cocoon wall create a rough surface that facilitates air turbulence and promotes disturbance amplitude of the flow field, leading to lengthened water vapor transfer path and increased tortuosity of the moist air. It also indicates the vortex of water vapor can be generated in the outer section of cocoon wall, which increases the diffusion distance of water vapor and enhances the turbulence kinetic energy and turbulence eddy dissipation, signifying higher moisture resistance in the outer section. The difference in moisture resistance of the multiple A. pernyi cocoon layers is largely responsible for the unique directional moisture transfer behavior of this wild silkworm cocoon. These findings may inspire a biomimicry approach to develop novel lightweight moisture management materials and structures.

Towards simulation of NASA35 axial compressor

Purpose – This paper aims to validate and analyse the NASA35 axial compressor performance based on a numerical approach. Design/methodology/approach – Knowledge about flow property change during compressor operation at high and relatively low speed is still limited. This work provides a numerical approach to address these problems. Validation of numerical methods is proposed to generate confidence the numerical approach adopted, and after that, analysis of compressor performance at different operation conditions is carried out. Findings – The numerical methods proposed are proved capable in predicting compressor performance. Changes of flow property during compressor operation are discussed and explained. Research limitations/implications – The current numerical work is carried out based on the first stage of the NASA35 axial compressor, where the interactive effects from adjacent stage are not counted in. Furthermore, the steady-state simulation enforces an averaging of flow at rotor-stator interface, where the transient rotor-stator interaction is removed. Practical implications – This work validates the numerical methods used in the prediction of NASA35 axial compressor performance, and a similar numerical approach can be used for other turbomachinery simulation cases. Originality/value – This work reinforces the understanding of axial compressor operation and provides reliable results for further investigation of a similar type of compressor. In addition, details of flow field within the NASA35 compressor during operation are given and explained which experiments still have difficult to achieve.

The Influence of Healed Intrablock Rockmass Structure on the Behaviour of Deep Excavations in Complex Rockmasses

Conventional rockmass characterization and analysis methods for geotechnical assessment in mining, civil tunnelling, and other excavations consider only the intact rock properties and the discrete fractures that are present and form blocks within rockmasses. Field logging and classification protocols are based on historically useful but highly simplified design techniques, including direct empirical design and empirical strength assessment for simplified ground reaction and support analysis. As modern underground excavations go deeper and enter into more high stress environments with complex excavation geometries and associated stress paths, healed structures within initially intact rock blocks such as sedimentary nodule boundaries and hydrothermal veins, veinlets and stockwork (termed intrablock structure) are having an increasing influence on rockmass behaviour and should be included in modern geotechnical design. Due to the reliance on geotechnical classification methods which predate computer aided analysis, these complexities are ignored in conventional design. Given the comparatively complex, sophisticated and powerful numerical simulation and analysis techniques now practically available to the geotechnical engineer, this research is driven by the need for enhanced characterization of intrablock structure for application to numerical methods. Intrablock structure governs stress-driven behaviour at depth, gravity driven disintegration for large shallow spans, and controls ultimate fragmentation. This research addresses the characterization of intrablock structure and the understanding of its behaviour at laboratory testing and excavation scales, and presents new methodologies and tools to incorporate intrablock structure into geotechnical design practice. A new field characterization tool, the Composite Geological Strength Index, is used for outcrop or excavation face evaluation and provides direct input to continuum numerical models with implicit rockmass structure. A brittle overbreak estimation tool for complex rockmasses is developed using field observations. New methods to evaluate geometrical and mechanical properties of intrablock structure are developed. Finally, laboratory direct shear testing protocols for interblock structure are critically evaluated and extended to intrablock structure for the purpose of determining input parameters for numerical models with explicit structure.

A Higher order and stable method for the numerical integration of Random Differential Equations

Trabalho apresentado no Congresso Nacional de Matemática Aplicada à Indústria, 18 a 21 de novembro de 2014, Caldas Novas - Goiás

Selecting the individual numerical scale and prioritization method in the analytic hierarchy process : a 2-tuple fuzzy linguistic approach

The validity of the priority vector used in the analytic hierarchy process (AHP) relies on two factors: the selection of a numerical scale and the selection of a prioritization method. The traditional AHP selects only one numerical scale (e.g., the Saaty scale) and one prioritization method (e.g., the eigenvector method) for each particular problem. For this traditional selection approach, there is disagreement on which numerical scale and prioritization method is better in deriving a priority vector. In fact, the best numerical scale and the best prioritization method both rely on the content of the pairwise comparison data provided by the AHP decision makers. By defining a set of concepts regarding the scale function and the linguistic pairwise comparison matrices (LPCMs) of the priority vector and by using LPCMs to unify the format of the input and output of AHP, this paper extends the AHP prioritization process under the 2-tuple fuzzy linguistic model. Based on the extended AHP prioritization process, we present two performance measure criteria to evaluate the effect of the numerical scales and prioritization methods. We also use the performance measure criteria to develop a 2-tuple fuzzy linguistic multicriteria approach to select the best numerical scales and the best prioritization methods for different LPCMs. In this paper, we call this type of selection the individual selection of the numerical scale and prioritization method. We also compare this individual selection with traditional selection by using both random and real data and show better results with individual selection.

A Numerical Solution Using an Adaptively Preconditioned Lanczos Method for a Class of Linear Systems Related with the Fractional Poisson Equation

This study considers the solution of a class of linear systems related with the fractional Poisson equation (FPE) (−∇2)α/2φ=g(x,y) with nonhomogeneous boundary conditions on a bounded domain. A numerical approximation to FPE is derived using a matrix representation of the Laplacian to generate a linear system of equations with its matrix A raised to the fractional power α/2. The solution of the linear system then requires the action of the matrix function f(A)=A−α/2 on a vector b. For large, sparse, and symmetric positive definite matrices, the Lanczos approximation generates f(A)b≈β0Vmf(Tm)e1. This method works well when both the analytic grade of A with respect to b and the residual for the linear system are sufficiently small. Memory constraints often require restarting the Lanczos decomposition; however this is not straightforward in the context of matrix function approximation. In this paper, we use the idea of thick-restart and adaptive preconditioning for solving linear systems to improve convergence of the Lanczos approximation. We give an error bound for the new method and illustrate its role in solving FPE. Numerical results are provided to gauge the performance of the proposed method relative to exact analytic solutions.

A numerical method to quantify differential axial shortening in concrete buildings

Differential distortion comprising axial shortening and consequent rotation in concrete buildings is caused by the time dependent effects of “shrinkage”, “creep” and “elastic” deformation. Reinforcement content, variable concrete modulus, volume to surface area ratio of elements and environmental conditions influence these distortions and their detrimental effects escalate with increasing height and geometric complexity of structure and non vertical load paths. Differential distortion has a significant impact on building envelopes, building services, secondary systems and the life time serviceability and performance of a building. Existing methods for quantifying these effects are unable to capture the complexity of such time dependent effects. This paper develops a numerical procedure that can accurately quantify the differential axial shortening that contributes significantly to total distortion in concrete buildings by taking into consideration (i) construction sequence and (ii) time varying values of Young’s Modulus of reinforced concrete and creep and shrinkage. Finite element techniques are used with time history analysis to simulate the response to staged construction. This procedure is discussed herein and illustrated through an example.

Look before you leap : a confidence-based method for selecting species criticality while avoiding negative populations in τ-leaping

The stochastic simulation algorithm was introduced by Gillespie and in a different form by Kurtz. There have been many attempts at accelerating the algorithm without deviating from the behavior of the simulated system. The crux of the explicit τ-leaping procedure is the use of Poisson random variables to approximate the number of occurrences of each type of reaction event during a carefully selected time period, τ. This method is acceptable providing the leap condition, that no propensity function changes “significantly” during any time-step, is met. Using this method there is a possibility that species numbers can, artificially, become negative. Several recent papers have demonstrated methods that avoid this situation. One such method classifies, as critical, those reactions in danger of sending species populations negative. At most, one of these critical reactions is allowed to occur in the next time-step. We argue that the criticality of a reactant species and its dependent reaction channels should be related to the probability of the species number becoming negative. This way only reactions that, if fired, produce a high probability of driving a reactant population negative are labeled critical. The number of firings of more reaction channels can be approximated using Poisson random variables thus speeding up the simulation while maintaining the accuracy. In implementing this revised method of criticality selection we make use of the probability distribution from which the random variable describing the change in species number is drawn. We give several numerical examples to demonstrate the effectiveness of our new method.