Search results

The objective of this paper is to establish a solution strategy for obtaining dual solutions, namely trivial (conventional) and nontrivial (unconventional) solutions, of coupled pore-fluid flow and chemical dissolution problems in heterogeneous porous media.

Through applying a perturbation to the pore-fluid velocity, original governing partial differential equations of a coupled pore-fluid flow and chemical dissolution problem in heterogeneous porous media are transformed into perturbed ones, which are then solved by using the semi-analytical finite element method. Through switching off and on the applied perturbation terms in the resulting perturbed governing partial differential equations, both the trivial and nontrivial solutions can be obtained for the original governing partial differential equations of the coupled pore-fluid flow and chemical dissolution problem in fluid-saturated heterogeneous porous media.

When a coupled pore-fluid flow and chemical dissolution system is in a stable state, the trivial and nontrivial solutions of the system are identical. However, if a coupled pore-fluid flow and chemical dissolution system is in an unstable state, then the trivial and nontrivial solutions of the system are totally different. This recognition can be equally used to judge whether a coupled pore-fluid flow and chemical dissolution system involving heterogeneous porous media is in a stable state or in an unstable state. The proposed solution strategy can produce dual solutions for simulating coupled pore-fluid flow and chemical dissolution problems in fluid-saturated heterogeneous porous media.

A solution strategy is proposed to obtain the nontrivial solution, which is often overlooked in the computational simulation of coupled pore-fluid flow and chemical dissolution problems in fluid-saturated heterogeneous porous media. The proposed solution strategy provides a useful way for understanding the underlying dynamic mechanisms of the chemical damage effect associated with the stability of structures that are built on soil foundations.

The objective of this paper is to develop a semi-analytical finite element method for solving chemical dissolution-front instability problems in fluid-saturated porous media.

The porosity, horizontal and vertical components of the pore-fluid velocity and solute concentration are selected as four fundamental unknown variables for describing chemical dissolution-front instability problems in fluid-saturated porous media. To avoid the use of numerical integration, analytical solutions for the property matrices of a rectangular element are precisely derived in a purely mathematical manner. This means that the proposed finite element method is a kind of semi-analytical method. The column pivot element solver is used to solve the resulting finite element equations of the chemical dissolution-front instability problem.

The direct use of horizontal and vertical components of the pore-fluid velocity as fundamental unknown variables can improve the accuracy of the related numerical solution. The column pivot element solver is useful for solving the finite element equations of a chemical dissolution-front instability problem. The proposed semi-analytical finite element method can produce highly accurate numerical solutions for simulating chemical dissolution-front instability problems in fluid-saturated porous media.

Analytical solutions for the property matrices of a rectangular element are precisely derived for solving chemical dissolution-front instability problems in fluid-saturated porous media. The proposed semi-analytical finite element method provides a useful way for understanding the underlying dynamic mechanisms of the washing land method involved in the contaminated land remediation.

WE write on the eve of an Annual Meeting of the Library Association. We expect many interesting things from it, for although it is not the first meeting under the new constitution, it is the first in which all the sections will be actively engaged. From a membership of eight hundred in 1927 we are, in 1930, within measurable distance of a membership of three thousand; and, although we have not reached that figure by a few hundreds—and those few will be the most difficult to obtain quickly—this is a really memorable achievement. There are certain necessary results of the Association's expansion. In the former days it was possible for every member, if he desired, to attend all the meetings; today parallel meetings are necessary in order to represent all interests, and members must make a selection amongst the good things offered. Large meetings are not entirely desirable; discussion of any effective sort is impossible in them; and the speakers are usually those who always speak, and who possess more nerve than the rest of us. This does not mean that they are not worth a hearing. Nevertheless, seeing that at least 1,000 will be at Cambridge, small sectional meetings in which no one who has anything to say need be afraid of saying it, are an ideal to which we are forced by the growth of our numbers.

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 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.

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.

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 service-profit chain model (Heskett, Jones, Loverman, Sasser, & Schlesinger, 1994) highlights the well-documented relationship between employee and customer attitudes suggesting that employees who are satisfied and engaged with their work provide better customer service resulting in higher levels of customer satisfaction and, ultimately, driving firm revenue. The authors propose an expansion of the service-profit margin identifying the leadership behaviors that create positive employee attitudes and engagement. Specifically, the authors suggest that leaders who focus on recognition, involvement, growth and development, health and safety, and teamwork (Kelloway, Nielsen, & Dimoff, 2017) create a psychologically healthy workplace for customer service providers and, ultimately, an enhanced customer experience.

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.

The purpose of this paper is to present an upscale theory of the thermal-mechanical coupling particle simulation for non-isothermal problems in two-dimensional quasi-static system, under which a small length-scale particle model can exactly reproduce the same mechanical and thermal results with that of a large length-scale one.

The objective is achieved by extending the upscale theory of particle simulation for two-dimensional quasi-static problems from an isothermal system to a non-isothermal one.

Five similarity criteria, namely geometric, material (mechanical and thermal) properties, gravity acceleration, (mechanical and thermal) time steps, thermal initial and boundary conditions (Dirichlet/Neumann boundary conditions), under which a small-length-scale particle model can exactly reproduce both the mechanical and thermal behavior with that of a large length-scale model for non-isothermal problems in a two-dimensional quasi-static system are proposed. Furthermore, to test the proposed upscale theory, two typical examples subjected to different thermal boundary conditions are simulated using two particle models of different length scale.

The paper provides some important theoretical guidances to modeling thermal-mechanical coupled problems at both the engineering length scale (i.e. the meter scale) and the geological length scale (i.e. the kilometer scale) using the particle simulation method directly. The related simulation results from two typical examples of significantly different length scales (i.e. a meter scale and a kilometer scale) have demonstrated the usefulness and correctness of the proposed upscale theory for simulating non-isothermal problems in two-dimensional quasi-static system.