This study demonstrates a three-dimensional numerical simulation of earth pressure balance (EPB) shield tunnelling using a coupled discrete element method (DEM) and a finite difference method (FDM). The analysis adopted the actual size of a spoke-type EPB shield tunnel boring machine (TBM) consisting of a cutter head with cutting tools, working chamber, screw conveyor, and shield. For the coupled model to reproduce the in situ ground condition, the ground formation was generated partially using the DEM (for the limited domain influenced by excavation), with the rest of the domain being composed of FDM grids. In the DEM domain, contact parameters of particles were calibrated via a series of large-scale triaxial test analyses. The model simulated tunnelling as the TBM operational conditions were controlled. The penetration rate and the rotational speed of the screw conveyor were automatically adjusted as the TBM advanced to prevent the generation of excessive or insufficient torque, thrust force, or chamber pressure. Accordingly, these parameters were maintained consistently around their set operational ranges during excavation. The simulation results show that the proposed numerical model based on DEM–FDM coupling could reasonably simulate EPB driving while considering the TBM operational conditions.

The presence of retained colloidal particles causes the retardation of contaminant transport when the contaminant is favorably adsorbed to colloidal particles. Although the particle size distribution affects the retention behavior of colloidal particles, the impact of particle size distribution on contaminant transport has not been reported to date. This study investigates the impact of the particle size distribution of the colloidal particles on contaminant transport through numerical simulation by representing the particle size distribution as a lognormal distribution function. In addition, the bed efficiency and contaminant saturation of simulated breakthrough curves were calculated, and a contaminant transport model with the Langmuir isotherm for the reaction between the contaminant–sand and contaminant–colloidal particle was introduced and validated with experimental data. The simulated breakthrough curves, bed efficiency, and contaminant saturation indicated that an increase in the mean and standard deviation of the particle size distribution causes the retardation of contaminant transport.

In this paper, a two-dimensional (2-D) Finite Element (FE) analysis of a geothermal well was performed with respect to five different cross-sections corresponding to the design specifications for the geothermal well that is currently constructed in Pohang, South Korea. Among the essential components (such as ground formation, casing, and cementing) of a geothermal well, the thermal and mechanical stability of the cementing component was discussed based on a series of parametric studies with consideration of the thermal conductivity and Young’s modulus of the cementing component. With increasing number of casing layers, the cementing component experiences less stress concentration. In addition, the lower thermal conductivity of the cementing material is advantageous for effectively controlling radial displacement. Consequently, it should be noted in geothermal well cementing construction that long-term strength degradation of the cementing might cause the severe structural instability of an entire geothermal well.

As weathered granite soils are extensively distributed in the Korean Peninsula, determining the optimal foam conditioning parameters for weathered granite soils can provide a high workability, low permeability, and high compressibility of weathered granite soils for successful operation of an earth pressure balance (EPB) shield tunnel boring machine (TBM). In this study, a laboratory-scale excavation apparatus was developed to conduct various laboratory excavation simulations for measuring the excavation-induced torque and the abrasivity of cutter bits. In addition, slump values, compressibility, and permeability were measured to evaluate the workability of foam-conditioned weathered granite soils. From the experimental results, the optimal foaming conditions for the weathered granite soil was obtained at foam injection ratios (F/Rs) of 60–70%, foam expansion ratios (FERs) of 15–18, and foaming agent concentrations (C_f) of 3–4%. The experimental results presented in this study can provide the rheological, mechanical, and hydraulic behavior of foam-conditioned weathered granite soils for successful EPB TBM operation.

The artificial ground freezing (AGF) method has been used in many geotechnical engineering applications such as temporary excavation support, underpinning and groundwater cutoff. The AGF method utilizes a refrigerant such as liquid nitrogen or brine, circulating through embedded freezing pipes in order to freeze the ground. In this paper, two in-situ cryogenic freezing experiments (i.e., single freezing-pipe test and frozen-wall formation test) were performed using liquid nitrogen to simulate the AGF in a Korean marine clay deposit, in which the freezing rate was evaluated. The thermal conductivity of frozen and unfrozen marine clay was evaluated by performing typical laboratory experiments. In addition, the strength and stiffness of frozen-thawed deposits were comparatively measured by sounding tests (i.e., piezocone penetration test and lateral loading test). The freezing rate of the frozen-wall formation test in the Korean marine clay deposit was approximately twice as high as that of the single freezing pipe test. Compared to the original marine clay deposit, the frozen-thawed marine clay showed a significant reduction in strength and stiffness.