Efficient parametric uncertainty analysis within the hybrid Finite Element/Statistical Energy Analysis method

This paper is concerned with the development of efficient algorithms for propagating parametric uncertainty within the context of the hybrid Finite Element/Statistical Energy Analysis (FE/SEA) approach to the analysis of complex vibro-acoustic systems. This approach models the system as a combination of SEA subsystems and FE components; it is assumed that the FE components have fully deterministic properties, while the SEA subsystems have a high degree of randomness. The method has been recently generalised by allowing the FE components to possess parametric uncertainty, leading to two ensembles of uncertainty: a non-parametric one (SEA subsystems) and a parametric one (FE components). The SEA subsystems ensemble is dealt with analytically, while the effect of the additional FE components ensemble can be dealt with by Monte Carlo Simulations. However, this approach can be computationally intensive when applied to complex engineering systems having many uncertain parameters. Two different strategies are proposed: (i) the combination of the hybrid FE/SEA method with the First Order Reliability Method which allows the probability of the non-parametric ensemble average of a response variable exceeding a barrier to be calculated and (ii) the combination of the hybrid FE/SEA method with Laplace's method which allows the evaluation of the probability of a response variable exceeding a limit value. The proposed approaches are illustrated using two built-up plate systems with uncertain properties and the results are validated against direct integration, Monte Carlo simulations of the FE and of the hybrid FE/SEA models. © 2013 Elsevier Ltd.

Finite element beam propagation method for analysis of plasmonic waveguide

The basic idea of the finite element beam propagation method (FE-BPM) is described. It is applied to calculate the fundamental mode of a channel plasmonic polariton (CPP) waveguide to confirm its validity. Both the field distribution and the effective index of the, fundamental mode are given by the method. The convergence speed shows the advantage and stability of this method. Then a plasmonic waveguide with a dielectric strip deposited on a metal substrate is investigated, and the group velocity is negative for the fundamental mode of this kind of waveguide. The numerical result shows that the power flow direction is reverse to that of phase velocity.

Interaction of linear waves with infinitely long horizontal cylinders studied by boundary element method

The two-dimensional problems concerning the interaction of linear water waves with cylinders of arbitrary shape in two-layer deep water are investigated by use of the Boundary Integral Equation method (BIEM). Simpler new expressions for the Green functions are derived, and verified by comparison of results obtained by BIEM with these by an analytical method. Examined are the radiation and scattering of linear waves by two typical configurations of cylinders in two-layer deep water. Hydrodynamic behaviors including hydrodynamic coefficients, wave forces, reflection and transmission coefficients and energies are analyzed in detail, and some interesting physical phenomena are observed.

Analysis of Si/GaAs Bonding Stresses with the Finite Element Method

In conjunction with ANSYS, we use the finite element method to analyze the bonding stresses of Si/GaAs. We also apply a numerical model to investigate a contour map and the distribution of normal stress,shearing stress,and peeling stress,taking into full consideration the thermal expansion coefficient as a function of temperature. Novel bonding structures are proposed for reducing the effect of thermal stress as compared with conventional structures. Calculations show the validity of this new structure.

Elastic Strain Field Distribution of a Self-Organized Growth Lensed-Shaped Quantum Dot by Finite Element Method

The stress and strain fields in self-organized growth coherent quantum dots (QD) structures are investigated in detail by two-dimension and three-dimension finite element analyses for lensed-shaped QDs. The nonobjective isolate quantum dot system is used. The calculated results can be directly used to evaluate the conductive band and valence band confinement potential and strain introduced by the effective mass of the charge carriers in strain QD.

Stress analysis of silica-based arrayed waveguide grating by a finite element method

The stress distribution in silica optical waveguides on silicon is calculated by using finite element method (FEM). The waveguides are mainly subjected to compressive stress along the x direction and the z direction, and it is accumulated near the interfaces between the core and cladding layers. The shift of central wavelength of silica arrayed waveguide grating (AWG) on silicon-substrate with the designed wavelength and the polarization dependence are caused by the stress in the silica waveguides.

U-transformation-finite element method in 3-dimensional analysis of orthotropic laminates

For an orthotropic laminate, an equivalent system with doubly cyclic periodicity is introduced. Then a 3-dimensional finite element model for the equivalent system is transformed into the unitary space, where the large finite element matrix equation is decoupled into some small matrix equations. Such a decoupling very efficiently reduces the computational effort. For an orthotropic laminate with four clamped edges, no exact elasticity solution is available, and the deflection values predicted by different methods have a considerable difference each other for a small length-to-thickness ratio. The present predictions are the largest because the present method is a full 3-dimensional finite element analysis without superfluous constraints. Illustrative numerical examples are presented to observe the distributions of stresses through the thickness of the laminates. (C) 2010 Elsevier Ltd. All rights reserved.

3D quench simulation using finite element method for CR superconducting magnet's strand

The CR superconducting magnet is a dipole of the FAIR project of GSI in Germany. The quench of the strand is simulated using FEM software ANSYS. From the simulation, the quench propagation can be visualized. Programming with APDL, the value of propagation velocity of normal zone is calculated. Also the voltage increasing over time of the strand is computed and pictured. Furthermore, the Minimum Propagation Zone (MPZ) is studied. At last, the relation between the current and the propagation velocity of normal zone, and the influence of initial temperature on quench propagation are studied.

Boundary Element Method (BEM) Analysis for Galvanic Corrosion of Hot Dip Galvanized Steel Immersed in Seawater

A numerical analysis of galvanic corrosion of hot-dip galvanized steel immersed in seawater was presented. The analysis was based on the boundary element methods (BEMs) coupled with Newton-Raphson iterative technique to treat the nonlinear boundary conditions, which were determined by the experimental polarization curves. Results showed that galvanic current density concentrates on the boundary of steel substrate and zinc coating, and the sacrificial protection of zinc coating to steel substrate results in overprotection of steel cathode. Not only oxygen reduction but also hydrogen reduction could occur as cathode reactions, which probably led up to the adsorption and absorption of hydrogen atoms. Flat galvanized steel tensile sample shows a brittle behavior similar to hydrogen embrittlement according to the SSRT (show strain rate test) in seawater.

An extension of Purcell's vector method with applications to panel element equations

The classical Purcell's vector method, for the construction of solutions to dense systems of linear equations is extended to a flexible orthogonalisation procedure. Some properties are revealed of the orthogonalisation procedure in relation to the classical Gauss-Jordan elimination with or without pivoting. Additional properties that are not shared by the classical Gauss-Jordan elimination are exploited. Further properties related to distributed computing are discussed with applications to panel element equations in subsonic compressible aerodynamics. Using an orthogonalisation procedure within panel methods enables a functional decomposition of the sequential panel methods and leads to a two-level parallelism.

Shape-optimized mesh partitioning and load balancing for parallel adaptive FEM

We present a dynamic distributed load balancing algorithm for parallel, adaptive Finite Element simulations in which we use preconditioned Conjugate Gradient solvers based on domain-decomposition. The load balancing is designed to maintain good partition aspect ratio and we show that cut size is not always the appropriate measure in load balancing. Furthermore, we attempt to answer the question why the aspect ratio of partitions plays an important role for certain solvers. We define and rate different kinds of aspect ratio and present a new center-based partitioning method of calculating the initial distribution which implicitly optimizes this measure. During the adaptive simulation, the load balancer calculates a balancing flow using different versions of the diffusion algorithm and a variant of breadth first search. Elements to be migrated are chosen according to a cost function aiming at the optimization of subdomain shapes. Experimental results for Bramble's preconditioner and comparisons to state-of-the-art load balancers show the benefits of the construction.

Dynamic load-balancing of finite element applications with the DRAMA library

The DRAMA library, developed within the European Commission funded (ESPRIT) project DRAMA, supports dynamic load-balancing for parallel (message-passing) mesh-based applications. The target applications are those with dynamic and solution-adaptive features. The focus within the DRAMA project was on finite element simulation codes for structural mechanics. An introduction to the DRAMA library will illustrate that the very general cost model and the interface designed specifically for application requirements provide simplified and effective access to a range of parallel partitioners. The main body of the paper will demonstrate the ability to provide dynamic load-balancing for parallel FEM problems that include: adaptive meshing, re-meshing, the need for multi-phase partitioning.

Micro-mechanical parameterisations for continuum modelling of granular material using the discrete element method

The present work uses the discrete element method (DEM) to describe assemblies of particulate bulk materials. Working numerical descriptions of entire processes using this scheme are infeasible because of the very large number of elements (1012 or more in a moderately sized industrial silo). However it is possible to capture much of the essential bulk mechanics through selective DEM on important regions of an assembly, thereafter using the information in continuum numerical descriptions of particulate processes. The continuum numerical model uses population balances of the various components in bulk solid mixtures. It depends on constitutive relationships for the internal transfer, creation and/or destruction of components within the mixture. In this paper we show the means of generating such relationships for two important flow phenomena – segregation whereby particles differing in some important property (often size) separate into discrete phases, and degradation, whereby particles break into sub-elements, through impact on each other or shearing. We perform DEM simulations under a range of representative conditions, extracting the important parameters for the relevant transfer, creation and/or destruction of particles in certain classes within the assembly over time. Continuum predictions of segregation and degradation using this scheme are currently being successfully validated against bulk experimental data and are beginning to be used in schemes to improve the design and operation of bulk solids process plant.

A hybrid Laplace transform/finite difference boundary element method for diffusion problems

The solution process for diffusion problems usually involves the time development separately from the space solution. A finite difference algorithm in time requires a sequential time development in which all previous values must be determined prior to the current value. The Stehfest Laplace transform algorithm, however, allows time solutions without the knowledge of prior values. It is of interest to be able to develop a time-domain decomposition suitable for implementation in a parallel environment. One such possibility is to use the Laplace transform to develop coarse-grained solutions which act as the initial values for a set of fine-grained solutions. The independence of the Laplace transform solutions means that we do indeed have a time-domain decomposition process. Any suitable time solver can be used for the fine-grained solution. To illustrate the technique we shall use an Euler solver in time together with the dual reciprocity boundary element method for the space solution