Finite element solution of Navier Stokes equations adapted to a priori error estimates

The paper is based on qualitative properties of the solution of the Navier-Stokes equations for incompressible fluid, and on properties of their finite element solution. In problems with corner-like singularities (e.g. on the well-known L-shaped domain) usually some adaptive strategy is used. In this paper we present an alternative approach. For flow problems on domains with corner singularities we use the a priori error estimates and asymptotic expansion of the solution to derive an algorithm for refining the mesh near the corner. It gives very precise solution in a cheap way. We present some numerical results.

Moderate energy impact analysis combining phenomenological contact law with localised damage and integral equation method

A computational impact analysis methodology has been developed, based on modal analysis and a local contact force-deflection model. The contact law is based on Hertz contact theory while contact stresses are elastic, defines a modified contact theory to take account of local permanent indentation, and considers elastic recovery during unloading. The model was validated experimentally through impact testing of glass-carbon hybrid braided composite panels. Specimens were mounted in a support frame and the contact force was inferred from the deceleration of the impactor, measured by high-speed photography. A Finite Element analysis of the panel and support frame assembly was performed to compute the modal responses. The new contact model performed well in predicting the peak forces and impact durations for moderate energy impacts (15 J), where contact stresses locally exceed the linear elastic limit and damage may be deemed to have occurred. C-scan measurements revealed substantial damage for impact energies in the range of 30-50 J. For this regime the new model predictions might be improved by characterisation of the contact law hysteresis during the unloading phase, and a modification of the elastic vibration response in line with damage levels acquired during the impact. © 2011 Elsevier Ltd. All rights reserved.

Branching toughens fibrous networks.

Fibrous collagenous networks are not only stiff but also tough, due to their complex microstructures. This stiff yet tough behavior is desirable for both medical and military applications but it is difficult to reproduce in engineering materials. While the nonlinear hyperelastic behavior of fibrous networks has been extensively studied, the understanding of toughness is still incomplete. Here, we identify a microstructure mimicking the branched bundles of a natural type I collagen network, in which partially cross-linked long fibers give rise to novel combinations of stiffness and toughness. Finite element analysis shows that the stiffness of fully cross-linked fibrous networks is amplified by increasing the fibril length and cross-link density. However, a trade-off of such stiff networks is reduced toughness. By having partially cross-linked networks with long fibrils, the networks have comparable stiffness and improved toughness as compared to the fully cross-linked networks. Further, the partially cross-linked networks avoid the formation of kinks, which cause fibril rupture during deformation. As a result, the branching allows the networks to have stiff yet tough behavior.

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.

Fracture toughness of fibrous membranes

Random fibrous networks exist in both natural biological and engineering materials. While the nonlinear deformation of fibrous networks has been extensively studied, the understanding of their fracture behaviour is still incomplete. To study the fracture toughness of fibrous materials, the near-tip region is crucial because failure mechanisms such as fibril rupture occur in this region. The consideration of this region in fracture studies is, however, a difficult task because it involves microscopic mechanical responses at a small length scale. This paper extends our previous finite element analysis by incorporating the microscopic responses into a macroscopic domain by using a submodeling technique. The detailed study of microstructures at crack tips show a stochastic toughness of membranes due to the random nature of fibrous networks. Further, the sizes of crack tip region, which are sufficient to provide a reasonable prediction of fracture behaviour in a specific type of fibrous network, were presented. Future work includes improving the current linear assumption in the macroscopic models to become nonlinear.

Finite element modelling of the sismic behaviour of water front structures

Water front structures have suffered significant damage in many of the recent earthquakes. These include gravity type quay walls, vertically composite walls, cantilever retaining walls, anchored bulkheads and similar structures. One of the primary causes for the poor performance of these classes of structures is the liquefaction of the foundation soil and in some instances liquefaction of the backfill soil. The liquefaction of the soil in-front of the quay wall tends to cause large lateral displacements and rotation of the wall. Often such gravity walls are placed on rubble mound deposited onto the sea bed.This paper presents finite element analyses of such a problem in which strength degradation of the foundation soil and the backfill material will be modelled using PZ mark III constitutive model. The performance of the wall in terms of its lateral displacement, vertical settlement and/or the rotation suffered by the wall will be presented. In addition, the contours of the horizontal and vertical effective stresses and the excess pore pressure ratio will be presented at different time instants together with hyrdraulic gradients. Immediately after the earthquake, the hydraulic gradients indicate migration of pore water into the region below the wall, suggesting further softening of the foundation soil below the wall.