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.

The objective of this paper is to provide a better understanding of the effect of irregular columnar jointed structure on the permeability and flow characteristics of rock masses.

An efficient numerical procedure is proposed to investigate the permeability and fluid flow in columnar jointed rock masses (CJRMs), of which the columnar jointed networks are generated by a modified constrained centroid Voronoi algorithm according to the field statistical results. The fractures are represented explicitly by using the lower-dimensional zero thickness elements. And the modeling scheme is validated by a benchmark test for flow in fractured porous media. The effective permeability and representative elementary volume (REV) size of CJRMs are estimated using finite element method (FEM). The influences of joint density and variation coefficient of columnar joint structure on the permeability of the rock mass are discussed.

The simulation results indicate that the permeability is scale-dependent and tends to be stable with increase of model size. The hydraulic REV size is determined as 3.5 m for CJRMs in the present study. Moreover, the joint density is a dominant factor affecting the permeability of CJRMs. The average permeability of columnar jointed structures increases linearly with the joint density under the same REV size, while the influence from the coefficient of variation can be neglected.

The present paper investigates the REV size of the CJRMs and the effect of joint parameters on the permeability. The proposed method and the results obtained are useful on understanding the hydraulic characteristic of the irregular CJRMs in engineering projects.

In this paper, three research topics are presented referring to different aspects of multifield problems in civil engineering. The first example deals with long term behaviour of wood under multiaxial states of stress and the effect of moisture changes on the deformation behaviour of wood. The second example refers to the application of a three‐phase model for soils to the numerical simulation of dewatering of soils by means of compressed air. The soil is modelled as a three phase‐material, consisting of the deformable soil skeleton and the fluid phases – water and compressed air. The third example is concerned with computational durability mechanics of concrete structures. As a particular example of chemically corrosive mechanisms, the material degradation due to the dissolution of calcium and external loading is addressed.

The stress–strain behaviors of rockfill materials in dams are significantly affected by the anisotropy and grain crushing. However, these factors are rarely considered in numerical simulations of high rockfill dams. This study intends to develop a reasonable and practical constitutive model for rockfill materials to overcome the above problems.

The effects of anisotropy and grain crushing are comprehensively considered by the spatial position of the reference state line. After the improved generalized plasticity model for rockfill materials (referred to as the PZR model) is developed and verified by laboratory tests, it is used with the finite element method to simulate the stress–strain behaviors of the Nuozhadu high core rockfill dam.

The simulated results agree well with the laboratory tests data and the situ monitoring data, verifying the reliability and practicability of the developed PZR model.

A new anisotropic state parameter is proposed to reflect the nonmonotonic variation in the strength as the major principal stress direction angle varies. This advantage is verified by the simulation of a set of conventional triaxial tests with different inclination angles of the compaction plane. 2) This is the first time that the elastoplastic model is verified by the situ monitoring data of high core rockfill dams. The numerical simulation results show that the PZR model can well reflect the stress–strain characteristics of rockfill materials in high core rockfill dams and is better than the traditional EB model.