SIMPLE REMESHING SCHEME FOR THE FINITE ELEMENT BASED SIMULATION OF DOPANT DIFFUSION IN SILICON.

Process simulation programs are valuable in generating accurate impurity profiles. Apart from accuracy the programs should also be efficient so as not to consume vast computer memory. This is especially true for devices and circuits of VLSI complexity. In this paper a remeshing scheme to make the finite element based solution of the non-linear diffusion equation more efficient is proposed. A remeshing scheme based on comparing the concentration values of adjacent node was then implemented and found to remove the problems of oscillation.

Centrifuge and discrete element modelling of tunnelling effects on pipelines

The paper presents centrifuge test data of the problem of tunnelling effects on buried pipelines and compares them to predictions made using DEM simulations. The paper focuses on the examination of pipeline bending moments, their distribution along the pipe, and their development with tunnel volume loss. Centrifuge results are obtained by PIV analysis and compared to results obtained using the DEM model. The DEM model was built to replicate the centrifuge model as closely as possible and included numerical features formulated specially for this task, such as structural elements to replicate the tunnel and pipeline. Results are extremely encouraging, with deviations between DEM and centrifuge test bending moment results being very small. © 2010 Taylor & Francis Group, London.

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.

Implementation of a PI phase angle controller for finite element analysis of the BDFM

Time-stepping finite element analysis of the BDFM for a specific load condition is shown to be a challenging problem because the excitation required cannot be predetermined and the BDFM is not open loops stable for all operating conditions. A simulation approach using feedback control to set the torque and stabilise the BDFM is presented together with implementation details. The performance of the simulation approach is demonstrated with an example and computed results are compared with measurements.

A subdivision-based implementation of the hierarchical b-spline finite element method

A novel technique is presented to facilitate the implementation of hierarchical b-splines and their interfacing with conventional finite element implementations. The discrete interpretation of the two-scale relation, as common in subdivision schemes, is used to establish algebraic relations between the basis functions and their coefficients on different levels of the hierarchical b-spline basis. The subdivision projection technique introduced allows us first to compute all element matrices and vectors using a fixed number of same-level basis functions. Their subsequent multiplication with subdivision matrices projects them, during the assembly stage, to the correct levels of the hierarchical b-spline basis. The proposed technique is applied to convergence studies of linear and geometrically nonlinear problems in one, two and three space dimensions. © 2012 Elsevier B.V.

3D modeling of high-T c superconductors by finite element software

A three-dimensional (3D) numerical model is proposed to solve the electromagnetic problems involving transport current and background field of a high-T c superconducting (HTS) system. The model is characterized by the E-J power law and H-formulation, and is successfully implemented using finite element software. We first discuss the model in detail, including the mesh methods, boundary conditions and computing time. To validate the 3D model, we calculate the ac loss and trapped field solution for a bulk material and compare the results with the previously verified 2D solutions and an analytical solution. We then apply our model to test some typical problems such as superconducting bulk array and twisted conductors, which cannot be tackled by the 2D models. The new 3D model could be a powerful tool for researchers and engineers to investigate problems with a greater level of complicity.

A finite element model of magnetization of superconducting bulks using a solid-state flux pump

Superconductors have a bright future; they are able to carry very high current densities, switch rapidly in electronic circuits, detect extremely small perturbations in magnetic fields, and sustain very high magnetic fields. Of most interest to large-scale electrical engineering applications are the ability to carry large currents and to provide large magnetic fields. There are many projects that use the first property, and these have concentrated on power generation, transmission, and utilization; however, there are relatively few, which are currently exploiting the ability to sustain high magnetic fields. The main reason for this is that high field wound magnets can and have been made from both BSCCO and YBCO, but currently, their cost is much higher than the alternative provided by low-Tc materials such as Nb3Sn and NbTi. An alternative form of the material is the bulk form, which can be magnetized to high fields. This paper explains the mechanism, which allows superconductors to be magnetized without the need for high field magnets to perform magnetization. A finite-element model is presented, which is based on the E-J current law. Results from this model show how magnetization of the superconductor builds up cycle upon cycle when a traveling magnetic wave is induced above the superconductor. © 2011 IEEE.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

Linear multiscale analysis and finite element validation of stretching and bending dominated lattice materials

The paper presents a multiscale procedure for the linear analysis of components made of lattice materials. The method allows the analysis of both pin-jointed and rigid-jointed microtruss materials with arbitrary topology of the unit cell. At the macroscopic level, the procedure enables to determine the lattice stiffness, while at the microscopic level the internal forces in the lattice elements are expressed in terms of the macroscopic strain applied to the lattice component. A numeric validation of the method is described. The procedure is completely automated and can be easily used within an optimization framework to find the optimal geometric parameters of a given lattice material. © 2011 Elsevier Ltd. All rights reserved.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.