A finite element-spherical harmonics model for radiative transfer in inhomogeneous clouds. Part II. some applications

EVENT has been used to examine the effects of 3D cloud structure, distribution, and inhomogeneity on the scattering of visible solar radiation and the resulting 3D radiation field. Large eddy simulation and aircraft measurements are used to create realistic cloud fields which are continuous or broken with smooth or uneven tops. The values, patterns and variance in the resulting downwelling and upwelling radiation from incident visible solar radiation at different angles are then examined and compared to measurements. The results from EVENT confirm that 3D cloud structure is important in determining the visible radiation field, and that these results are strongly influenced by the solar zenith angle. The results match those from other models using visible solar radiation, and are supported by aircraft measurements of visible radiation, providing confidence in the new model.

Noether-type discrete conserved quantities arising from a finite element approximation of a variational problem

In this work, we prove a weak Noether-type Theorem for a class of variational problems that admit broken extremals. We use this result to prove discrete Noether-type conservation laws for a conforming finite element discretisation of a model elliptic problem. In addition, we study how well the finite element scheme satisfies the continuous conservation laws arising from the application of Noether’s first theorem (1918). We summarise extensive numerical tests, illustrating the conservation of the discrete Noether law using the p-Laplacian as an example and derive a geometric-based adaptive algorithm where an appropriate Noether quantity is the goal functional.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

A finite element simulation of micro-mechanical frictional behaviour in metal forming

Friction is a critical factor for sheet metal forming (SMF). The Coulomb friction model is usually used in most finite element (FE) simulation for SMF. However, friction is a function of the local contact deformation conditions, such as local pressure, roughness and relative velocity. Frictional behaviour between contact surfaces can be based on three cases: boundary, hydrodynamic and mixed lubrication. In our microscopic friction model based on the finite element method (FEM), the case of dry contact between sheet and tool has been considered. In the view of microscopic geometry, roughness depends upon amplitude and wavelength of surface asperities of sheet and tool. The mean pressure applied on the surface differs from the pressure over the actual contact area. The effect of roughness (microscopic geometric condition) and relative speed of contact surfaces on friction coefficient was examined in the FE model for the microscopic friction behaviour. The analysis was performed using an explicit FE formulation. In this study, it was found that the roughness of deformable sheet decreases during sliding and the coefficient of friction increases with increasing roughness of contact surfaces. Also, the coefficient of friction increases with the increase of relative velocity and adhesive friction coefficient between contact surfaces.

Finite element modelling of metallic tubular crash structures with an explicit code

Numerous experimental studies have been carried out to investigate the collapse of tubular metallic crash structures under axial compression. Some simple theoretical models have been developed but these often assume one type of progressive collapse, which is not always representative of the real situation. Finite Element (FE) models, when further refined, have the potential to predict the actual collapse mode and how it influences the load-displacement and energy absorption characteristics. This paper describes an FE modelling investigation with the explicit code LS−DYNA. An automatic mesh generation programme written by the authors is used to set up shell and solid element tube models. Mesh specification issues and features relating to the contact and friction models are discussed in detail. The crush modes, load-deflection characteristics and energy absorption values found in the simulations are compared with a reasonable degree of correlation to those observed in a physical testing programme; however, improvements are still required.

Object oriented finite element analysis: a distributed approach to mesh generation

The object-oriented finite element method (OOFEM) has attracted the attention of many researchers. Compared with the traditional finite element method, OOFEM software has the advantages of maintenance and reuse. Moreover, it is easier to expand the architecture to a distributed one. In this paper, we introduce a distributed architecture of a object-oriented finite element preprocessor. A comparison between the distributed system and the centralised system shows that the former, presented in the paper, greatly improves the performance of mesh generation. Other finite element analysis modules could be expanded according to this architecture.

The design and development of an integrated system for object-oriented finite element computing

Finite Element Method (FEM) is widely used in Science and Engineering since 1960’s. The vast majority of FEM software is procedure-oriented. However, this conventional style of designing FEM software encounters problems in maintenance, reuse, and expansion of the software. Recently the object-oriented finite element method attracts the attention of lots of researchers, and now there is a growing interest in this method. In this paper, the object-oriented finite element (OOFE) is briefly introduced. Then the design and development of an integrated OOFE system is described. A comparison of the integrated OOFE system and a procedure-oriented system shows that our OOFE system has many advantages.

Finite element simulation of deformation processes of micro-components

The aim of this paper is to improve the understanding of deformation of micro medical needle and thread during assembly and then to develop an economical and flexible deformation method. Therefore, the swaging process is computationally simulated with the finite element method in this paper. A commercially available explicit nonlinear finite element analysis code, LS-Dyna, is used to model the 3-D deformation and contact problem. As the firmness of the assembly on the needle depends on the contact force and friction, the contact and the slide between the needle and thread are taken into account in the simulation. The general surface-to-surface contact algorithm (STS) is used to simulate the contact. The paper provides an insight into the deformation of the micro products.

Three-dimensional Finite Element modelling of laser shock peening process

Laser shock peening (LSP) is an innovative surface treatment technique for metal alloys, with the great improvement of their fatigue, corrosion and wear resistance performance. Finite element method has been widely applied to simulate the LSP to provide the theoretically predictive assessment and optimally parametric design. In the current work, 3-D numerical modelling approaches, combining the explicit dynamic analysis, static equilibrium analysis algorithms and different plasticity models for the high strain rate exceeding 106s-1, are further developed. To verify the proposed methods, 3-D static and dynamic FEA of AA7075-T7351 rods subject to two-sided laser shock peening are performed using the FEA package–ABAQUS. The dynamic and residual stress fields, shock wave propagation and surface deformation of the treated metal from different material modelling approaches have a good agreement.