Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

Impact and energy absorption of empty and foam-filled conical tubes

Foam-filled conical tubes have recently emerged as efficient energy absorbing devices to mitigate the adverse effects of impacts. The primary aim of this thesis was to generate research and design information on the impact and energy absorption response of empty and foam-filled conical tubes, and to facilitate their application in energy absorbing systems under axial and oblique loading conditions representative of those typically encountered in crashworthiness and impact applications. Finite element techniques supported by experiments and existing results were used in the investigation. Major findings show that the energy absorption response can be effectively controlled by varying geometry and material parameters. A useful empirical formula was developed for providing engineering designers with an initial estimate of the load ratio and hence energy absorption performances of these devices. It was evident that foam-filled conical tubes enhance the energy absorption capacity and stabilise the crush response for both axial and oblique impact loading without a significant increase in the initial peak load. This is practically beneficial when higher kinetic energy needs to be absorbed, thus reducing the impact force transmitted to the protected structure and occupants. Such tubes also increase and maintain the energy absorption capacity under global bending as well as minimise the reduction of energy absorption capacity with increasing load angle. Furthermore, the results also highlight the feasibility of adding a foam-filled conical tube as a supplementary device in energy absorbing systems, since the overall energy absorption performance of such systems can be favourably enhanced by only including a relatively small energy absorbing device. Above all, the results demonstrate the superior performance of foam-filled conical tube for mitigating impact energy in impact and crashworthiness applications.

Structure and deformation correlation of closed-cell aluminium foam subject to uniaxial compression

We report the results of an experimental and numerical study conducted on a closed-cell aluminium foam that was subjected to uniaxial compression with lateral constraint. X-ray computed tomography was utilized to gain access into the three-dimensional (3-D) structure of the foam and some aspects of the deformation mechanisms. A series of advanced 3-D image analyses are conducted on the 3-D images aimed at characterizing the strain localization regions. We identify the morphological/geometrical features that are responsible for the collapse of the cells and the strain localization. A novel mathematical approach based on a Minkowski tensor analysis along with the mean intercept length technique were utilized to search for signatures of anisotropy across the foam sample and its evolution as a function of loading. Our results show that regions with higher degrees of anisotropy in the undeformed foam have a tendency to initiate the onset of cell collapse. Furthermore, we show that strain hardening occurs predominantly in regions with large cells and high anisotropy. We combine the finite element method with the tomographic images to simulate the mechanical response of the foam. We predict further deformation in regions where the foam is already deformed. Crown Copyright (C) 2012 Published by Elsevier Ltd. on behalf of Acta Materialia Inc. All rights reserved.

Impact of symmetrized and Burt-Foreman Hamiltonians on spurious solutions and energy levels of InAs/GaAs quantum dots

We present a systematic investigation of calculating quantum dots (QDs) energy levels using finite element method in the frame of eight-band k . p method. Numerical results including piezoelectricity, electron and hole levels, as yell as wave functions are achieved. In the calculation of energy levels, we do observe spurious solutions (SSs) no matter Burt-Foreman or symmetrized Hamiltonians are used. Different theories are used to analyse the SSs, we find that the ellipticity theory can give a better explanation for the origin of SSs and symmetrized Hamiltonian is easier to lead to SSs. The energy levels simulated with the two Hamiltonians are compared to each other after eliminating SSs, different Hamiltonians cause a larger difference on electron energy levels than that on hole energy levels and this difference decreases with the increase of QD size.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

Impact of finite element idealisation on the prediction of welded fuselage stiffened panel buckling

Lap joints are widely used in the manufacture of stiffened panels and influence local panel sub-component stability, defining buckling unit dimensions and boundary conditions. Using the Finite Element method it is possible to model joints in great detail and predict panel buckling behaviour with accuracy. However, when modelling large panel structures such detailed analysis becomes computationally expensive. Moreover, the impact of local behaviour on global panel performance may reduce as the scale of the modelled structure increases. Thus this study presents coupled computational and experimental analysis, aimed at developing relationships between modelling fidelity and the size of the modelled structure, when the global static load to cause initial buckling is the required analysis output. Small, medium and large specimens representing welded lap-joined fuselage panel structure are examined. Two element types, shell and solid-shell, are employed to model each specimen, highlighting the impact of idealisation on the prediction of welded stiffened panel initial skin buckling.

Finite element analysis of warpage in laminated aluminium alloy plates for machining of primary aeronautic parts

The aim of this paper consists in presenting a method of simulating the warpage in 7xxx series aluminium alloy plates. To perform this simulation finite element software MSC.Patran and MSC.Marc were used. Another result of this analysis will be the influence on material residual stresses induced on the raw material during the rolling process upon the warpage of primary aeronautic parts, fabricated through machining (milling) at Embraer. The method used to determinate the aluminium plate residual stress was Layer Removal Test. The numerical algorithm Modified Flavenot Method was used to convert layer removal and beam deflection in stress level. With such information about the level and profile of residual stresses become possible, during the step that anticipate the manufacturing to incorporate these values in the finite-element approach for modelling warpage parts. Based on that warpage parameter surely the products are manufactured with low relative vulnerability propitiating competitiveness and price. © 2007 American Institute of Physics.

Finite-element discretization of 3D energy-transport equations for semiconductors

In this thesis a mathematical model was derived that describes the charge and energy transport in semiconductor devices like transistors. Moreover, numerical simulations of these physical processes are performed. In order to accomplish this, methods of theoretical physics, functional analysis, numerical mathematics and computer programming are applied. After an introduction to the status quo of semiconductor device simulation methods and a brief review of historical facts up to now, the attention is shifted to the construction of a model, which serves as the basis of the subsequent derivations in the thesis. Thereby the starting point is an important equation of the theory of dilute gases. From this equation the model equations are derived and specified by means of a series expansion method. This is done in a multi-stage derivation process, which is mainly taken from a scientific paper and which does not constitute the focus of this thesis. In the following phase we specify the mathematical setting and make precise the model assumptions. Thereby we make use of methods of functional analysis. Since the equations we deal with are coupled, we are concerned with a nonstandard problem. In contrary, the theory of scalar elliptic equations is established meanwhile. Subsequently, we are preoccupied with the numerical discretization of the equations. A special finite-element method is used for the discretization. This special approach has to be done in order to make the numerical results appropriate for practical application. By a series of transformations from the discrete model we derive a system of algebraic equations that are eligible for numerical evaluation. Using self-made computer programs we solve the equations to get approximate solutions. These programs are based on new and specialized iteration procedures that are developed and thoroughly tested within the frame of this research work. Due to their importance and their novel status, they are explained and demonstrated in detail. We compare these new iterations with a standard method that is complemented by a feature to fit in the current context. A further innovation is the computation of solutions in three-dimensional domains, which are still rare. Special attention is paid to applicability of the 3D simulation tools. The programs are designed to have justifiable working complexity. The simulation results of some models of contemporary semiconductor devices are shown and detailed comments on the results are given. Eventually, we make a prospect on future development and enhancements of the models and of the algorithms that we used.

Numerical analysis of roadside design (NARD). Volume II: programmers manual. Final report.

Federal Highway Administration, Office of Safety and Traffic Operations Research and Development, McLean, Va.

Numerical analysis of roadside design (NARD). Volume I: users manual. Final report.

Federal Highway Administration, Office of Safety and Traffic Operations Research and Development, McLean, Va.

Characterisation and investigation of local variations in mechanical behaviour in cast aluminium using gradient solidification, Digital Image Correlation and finite element simulation

Due to design and process-related factors, there are local variations in the microstructure and mechanical behaviour of cast components. This work establishes a Digital Image Correlation (DIC) based method for characterisation and investigation of the effects of such local variations on the behaviour of a high pressure, die cast (HPDC) aluminium alloy. Plastic behaviour is studied using gradient solidified samples and characterisation models for the parameters of the Hollomon equation are developed, based on microstructural refinement. Samples with controlled microstructural variations are produced and the observed DIC strain field is compared with Finite Element Method (FEM) simulation results. The results show that the DIC based method can be applied to characterise local mechanical behaviour with high accuracy. The microstructural variations are observed to cause a redistribution of strain during tensile loading. This redistribution of strain can be predicted in the FEM simulation by incorporating local mechanical behaviour using the developed characterization model. A homogeneous FEM simulation is unable to predict the observed behaviour. The results motivate the application of a previously proposed simulation strategy, which is able to predict and incorporate local variations in mechanical behaviour into FEM simulations already in the design process for cast components.

Dynamic energy absorption characteristics of foam-filled conical tubes under oblique impact loading

This paper treats the crush behaviour and energy absorption response of foam-filled conical tubes subjected to oblique impact loading. Dynamic computer simulation techniques validated by experimental testing are used to carry out a parametric study of such devices. The study aims at quantifying the energy absorption of empty and foam-filled conical tubes under oblique impact loading, for variations in the load angle and geometry parameters of the tube. It is evident that foam-filled conical tubes are preferable as impact energy absorbers due to their ability to withstand oblique impact loads as effectively as axial impact loads. Furthermore, it is found that the energy absorption capacity of filled tubes is better maintained compared to that of empty tubes as the load orientation increases. The primary outcome of this study is design information for the use of foam-filled conical tubes as energy absorbers where oblique impact loading is expected.

Coupling of energy from quantum emitters to the plasmonic mode of V groove waveguides: A numerical study

This work is a theoretical investigation into the coupling of a single excited quantum emitter to the plasmon mode of a V groove waveguide. The V groove waveguide consists of a triangular channel milled in gold and the emitter is modeled as a dipole emitter, and could represent a quantum dot, nitrogen vacancy in diamond, or similar. In this work the dependence of coupling efficiency of emitter to plasmon mode is determined for various geometrical parameters of the emitter-waveguide system. Using the finite element method, the effect on coupling efficiency of the emitter position and orientation, groove angle, groove depth, and tip radius, is studied in detail. We demonstrate that all parameters, with the exception of groove depth, have a significant impact on the attainable coupling efficiency. Understanding the effect of various geometrical parameters on the coupling between emitters and the plasmonic mode of the waveguide is essential for the design and optimization of quantum dot–V groove devices.