Search results

Classical continuum models, i.e. continuum models that do not incorporate an internal length scale, suffer from excessive mesh dependence when strain‐softening models are used in numerical analyses and cannot reproduce the size effect commonly observed in quasi‐brittle failure. In this contribution three different approaches will be scrutinized which may be used to remedy these two intimately related deficiencies of the classical theory, namely (i) the addition of higher‐order deformation gradients, (ii) the use of micropolar continuum models, and (iii) the addition of rate dependence. By means of a number of numerical simulations it will be investigated under which conditions these enriched continuum theories permit localization of deformation without losing ellipticity for static problems and hyperbolicity for dynamic problems. For the latter class of problems the crucial role of dispersion in wave propagation in strain‐softening media will also be highlighted.

Classical continuum models, i.e. continuum models that do not incorporate an internal length scale, suffer from pathological mesh‐dependence when strain‐softening models are employed in failure analyses. In this contribution the governing field equations are regularized by adding rotational degrees‐of‐freedom to the conventional translational degrees‐of‐freedom. This so‐called elasto‐plastic Cosserat continuum model, for which an efficient and accurate integration algorithm and a consistent tangent operator are also derived in this contribution, warrants convergence of the load—deflection curve to a unique solution upon mesh refinement and a finite width of the localization zone. This is demonstrated for an infinitely long shear layer and a biaxial test of a strain‐softening elasto‐plastic von Mises material.

We present the finite element simulations of reactive mineral‐carrying fluids mixing and mineralization in pore‐fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for lead, zinc and iron minerals in the regime of temperature‐gradient‐driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore‐fluid velocity and the equilibrium concentration of this particular kind of mineral. Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral‐carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

The dispersive behaviour of waves in softening problems is analysed. Attention is focused on the influence of the numerical scheme on the dispersion characteristics in the process of localization of deformation. Distinction has been made between softening models defined in a standard plasticity framework and in a gradient‐dependent plasticity theory. Waves in a standard softening plasticity continuum do not disperse but due to spatial discretization dispersion is introduced which results in a mesh size dependent length scale effect. On the other hand, wave propagation in a gradient‐dependent softening plasticity continuum is dispersive. By carrying out the dispersion analysis on the discretized system the influence of numerical dispersion on material dispersion can be quantified which enables us to determine the accuracy for the solution of the localization zone. For a modelling with and without the inclusion of strain gradients accuracy considerations with respect to mass discretization, finite element size, time integration scheme and time step have been carried out.

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to solve any new kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid‐saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore‐fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

This paper aims to develop a finite element analysis strategy, which is suitable for the analysis of progressive failure that occurs in pressure-dependent materials in practical engineering problems.

The numerical difficulties stemming from the strain-softening behaviour of the frictional material, which is represented by a non-associated Drucker–Prager material model, is tackled using the Cosserat continuum theory, while the mixed finite element formulation based on Hu–Washizu variational principle is adopted to allow the utilization of low-order finite elements.

The effectiveness and robustness of the low-order finite element are verified, and the simulation for a real-world landslide which occurred at the upstream side of Carsington embankment in Derbyshire reconfirms the advantages of the developed elastoplastic Cosserat continuum scheme in capturing the entire progressive failure process when the strain-softening and the non-associated plastic law are involved.

The permit of using low-order finite elements is of great importance to enhance computational efficiency for analysing large-scale engineering problems. The case study reconfirms the advantages of the developed elastoplastic Cosserat continuum scheme in capturing the entire progressive failure process when the strain-softening and the non-associated plastic law are involved.

In many scientific and engineering fields, large‐scale heat transfer problems with temperature‐dependent pore‐fluid densities are commonly encountered. For example, heat transfer from the mantle into the upper crust of the Earth is a typical problem of them. The main purpose of this paper is to develop and present a new combined methodology to solve large‐scale heat transfer problems with temperature‐dependent pore‐fluid densities in the lithosphere and crust scales.

The theoretical approach is used to determine the thickness and the related thermal boundary conditions of the continental crust on the lithospheric scale, so that some important information can be provided accurately for establishing a numerical model of the crustal scale. The numerical approach is then used to simulate the detailed structures and complicated geometries of the continental crust on the crustal scale. The main advantage in using the proposed combination method of the theoretical and numerical approaches is that if the thermal distribution in the crust is of the primary interest, the use of a reasonable numerical model on the crustal scale can result in a significant reduction in computer efforts.

From the ore body formation and mineralization points of view, the present analytical and numerical solutions have demonstrated that the conductive‐and‐advective lithosphere with variable pore‐fluid density is the most favorite lithosphere because it may result in the thinnest lithosphere so that the temperature at the near surface of the crust can be hot enough to generate the shallow ore deposits there. The upward throughflow (i.e. mantle mass flux) can have a significant effect on the thermal structure within the lithosphere. In addition, the emplacement of hot materials from the mantle may further reduce the thickness of the lithosphere.

The present analytical solutions can be used to: validate numerical methods for solving large‐scale heat transfer problems; provide correct thermal boundary conditions for numerically solving ore body formation and mineralization problems on the crustal scale; and investigate the fundamental issues related to thermal distributions within the lithosphere. The proposed finite element analysis can be effectively used to consider the geometrical and material complexities of large‐scale heat transfer problems with temperature‐dependent fluid densities.

Numerical methods are used to solve double diffusion driven reactive flow transport problems in deformable fluid‐saturated porous media. In particular, the temperature dependent reaction rate in the non‐equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAC and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore‐fluid flow, heat transfer and species transport/chemical reactions in deformable fluid‐saturated porous media has been developed. The coupled problem is divided into two sub‐problems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, it is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper. The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

We present a numerical methodology for the study of convective pore‐fluid, thermal and mass flow in fluid‐saturated porous rock basins. In particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore‐fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore‐fluid flow and mass transport in fluid‐saturated hydrothermal basins; (2) Convective pore‐fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence on the distribution patterns of convective pore‐fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore‐fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

The formulation of finite elements with incompatible discontinuous modes is examined rigorously. Both weak and strong discontinuities are considered. Starting from a careful elaboration of the kinematics for both types of discontinuities a comprehensive finite element formulation is derived based on a three‐field variational statement. Similarities and differences are highlighted between the various formulations which ensue.