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.

Comparison of finite element and finite volume methods application in geometrically nonlinear stress analysis

A novel three-dimensional finite volume (FV) procedure is described in detail for the analysis of geometrically nonlinear problems. The FV procedure is compared with the conventional finite element (FE) Galerkin approach. FV can be considered to be a particular case of the weighted residual method with a unit weighting function, where in the FE Galerkin method we use the shape function as weighting function. A Fortran code has been developed based on the finite volume cell vertex formulation. The formulation is tested on a number of geometrically nonlinear problems. In comparison with FE, the results reveal that FV can reach the FE results in a higher mesh density.

Numerical modelling and validation of Marangoni and surface tension phenomena using the finite volume method

Surface tension induced flow is implemented into a numerical modelling framework and validated for a number of test cases. Finite volume unstructured mesh techniques are used to discretize the mass, momentum and energy conservation equations in three dimensions. An explicit approach is used to include the effect of surface tension forces on the flow profile and final shape of a liquid domain. Validation of this approach is made against both analytical and experimental data. Finally, the method is used to model the wetting balance test for solder alloy material, where model predictions are used to gain a greater insight into this process. Copyright © 2000 John Wiley & Sons, Ltd.

Viscoelastic flow through a planar contraction using a semi-Lagrangian finite volume method

A semi-Lagrangian finite volume scheme for solving viscoelastic flow problems is presented. A staggered grid arrangement is used in which the dependent variables are located at different mesh points in the computational domain. The convection terms in the momentum and constitutive equations are treated using a semi-Lagrangian approach in which particles on a regular grid are traced backwards over a single time-step. The method is applied to the 4 : 1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions. The development of vortex behaviour with increasing values of We is analyzed.

A semi-Lagrangian finite volume method for Newtonian contraction flows

A new finite volume method for solving the incompressible Navier--Stokes equations is presented. The main features of this method are the location of the velocity components and pressure on different staggered grids and a semi-Lagrangian method for the treatment of convection. An interpolation procedure based on area-weighting is used for the convection part of the computation. The method is applied to flow through a constricted channel, and results are obtained for Reynolds numbers, based on half the flow rate, up to 1000. The behavior of the vortex in the salient corner is investigated qualitatively and quantitatively, and excellent agreement is found with the numerical results of Dennis and Smith [Proc. Roy. Soc. London A, 372 (1980), pp. 393-414] and the asymptotic theory of Smith [J. Fluid Mech., 90 (1979), pp. 725-754].

Combined vertex-based - cell-centred finite volume method for flows in complex geometries

CFD modelling of 'real-life' processes often requires solutions in complex three dimensional geometries, which can often result in meshes where aspects of it are badly distorted. Cell-centred finite volume methods, typical of most commercial CFD tools, are computationally efficient, but can lead to convergence problems on meshes which feature cells with high non-orthogonal shapes. The vertex-based finite volume method handles distorted meshes with relative ease, but is computationally expensive. A combined vertex-based - cell-centred (VB-CC) technique, detailed in this paper, allows solutions on distorted meshes that defeat purely cell-centred physical models to be employed in the solution of other transported quantities. The VB-CC method is validated with benchmark solutions for thermally driven flow and turbulent flow. An early application of this hybrid technique is to three-dimensional flow over an aircraft wing, although it is planned to use it in a wide variety of processing applications in the future.

A cell-centred finite volume method for modelling viscoelastic flow

An unstructured cell-centred finite volume method for modelling viscoelastic flow is presented. The method is applied to the flow through a planar channel and the 4:1 planar contraction for creeping flow of an Oldroyd-B fluid. The results are presented for a range of Weissenberg numbers. In the case of the planar channel results are compared with analytical solutions. For the 4:1 planar contraction benchmark problem the convection terms in the constitutive equations are approximated using both first and second order differencing schemes to compare the techniques and the effect of mesh refinement on the solution is investigated. This is the first time that a fully unstructured, cell-centredfinitevolume technique has been used to model the Oldroyd-B fluid for the test cases presented in this paper.

Numerical simulation of a patent technology for sealing of deep-sea oil wells using nonlinear finite element method

For over 50 years bridge plugs and cement have been used for well abandonment and work over and are still the material of choice. However the failures of cement abandonments using bridge plugs has been reported on many occasions, some of which have resulted in fatal consequences. A new patented product is designed to address the shortcomings associated with using bridge plugs and cement. The new developed tools use an alloy based on bismuth that is melted in situ using Thermite reaction. The tool uses the expansion properties of bismuth to seal the well. Testing the new technology in real field under more than 2 km deep sea water can be expensive. Virtual simulation of the new device under simulated thermal and mechanical environment can be achieved using nonlinear finite element method to validate the product and reduce cost. Experimental testing in the lab is performed to measure heat generated due to thermite reaction. Then, a sequential thermal mechanical explicit/implicit finite element solver is used to simulate the device under both testing lab and deep water conditions.