• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.1
• /
• pp.49-57
• /
• 2008

One of the most important design concerns for undersea tunnels is to establish design water load and flow rate. These are greatly dependent on the hydraulic factors such as water head, cover depth, hydraulic boundary conditions. In this paper, the influence of the hydraulic design factors on the ground loading and the inflow rate was investigated using the coupled finite element method. A horse shoe-shaped tunnel constructed 30 m below sea bottom was adopted to evaluate the water head effect considering various water depth for varying hydraulic conditions and relative permeability between lining and ground. The effect of cover depth was analysed for varying cover depth with the water depth of 60 m. The results were considered in terms of pore water pressure, ground loading and flow rate. Ground loading increases with an increase in water head and cover depth without depending on hydraulic boundary conditions. This points out that in leaking tunnels an increase in water depth increases seepage force which consequently increases ground loading. Furthermore, it is identified that an increase in water head and cover depth increases the rate of inflow and a decrease in the permeability ratio reduces the rate of inflow considerably.

• Geomechanics and Engineering
• /
• v.23 no.3
• /
• pp.189-206
• /
• 2020

The accurate modeling of boundary conditions is important in simulations of the discrete element method (DEM) and can affect the numerical results significantly. In conventional triaxial compression (CTC) tests, the specimens are wrapped by flexible membranes allowing to deform freely. To accurately model the boundary conditions of CTC, new flexible boundary algorithms for 2D and 3D DEM simulations are proposed. The new algorithms are computationally efficient and easy to implement. Moreover, both horizontal and vertical component of confining pressure are considered in the 2D and 3D algorithms, which can ensure that the directions of confining pressure are always perpendicular to the specimen surfaces. Furthermore, the boundaries are continuous and closed in the new algorithms, which can prevent the escape of particles from the specimens. The effectiveness of the proposed algorithms is validated by biaxial and triaxial simulations of granular materials. The results show that the algorithms allow the boundaries to deform non-uniformly on the premise of maintaining high control accuracy of confining pressure. Meanwhile, the influences of different lateral boundary conditions on the numerical results are discussed. It is indicated that the flexible boundary is more appropriate for the models with large strain or significant localization than rigid boundary.

• Geomechanics and Engineering
• /
• v.21 no.5
• /
• pp.471-487
• /
• 2020

In this work, the buckling analysis of material sandwich plates based on a two-parameter elastic foundation under various boundary conditions is investigated on the basis of a new theory of refined trigonometric shear deformation. This theory includes indeterminate integral variables and contains only four unknowns in which any shear correction factor not used, with even less than the conventional theory of first shear strain (FSDT). Applying the principle of virtual displacements, the governing equations and boundary conditions are obtained. To solve the buckling problem for different boundary conditions, Galerkin's approach is utilized for symmetric EGM sandwich plates with six different boundary conditions. A detailed numerical study is carried out to examine the influence of plate aspect ratio, elastic foundation coefficients, ratio, side-to-thickness ratio and boundary conditions on the buckling response of FGM sandwich plates. A good agreement between the results obtained and the available solutions of existing shear deformation theories that have a greater number of unknowns proves to demonstrate the precision of the proposed theory.

• Nuclear Engineering and Technology
• /
• v.55 no.5
• /
• pp.1911-1923
• /
• 2023

The nuclear safety assessment involving large transient simulations is forcing the community to develop methods for coupling thermal-hydraulics and neutronic codes and three-dimensional (3D) Computational Fluid Dynamics (CFD) codes. In this paper a set of dynamic boundary conditions are implemented in OpenFOAM^Ⓡ in order to apply zero-dimensional (0D) approaches coupling with 3D thermal-hydraulic simulation in a single framework. This boundary conditions are applied to model pipelines, tanks, pumps, and heat exchangers. On a first stage, four tests are perform in order to assess the implementations. The results are compared with experimental data, full 3D CFD, and system code simulations, finding a general good agreement. The semi-implicit implementation nature of these boundary conditions has shown robustness and accuracy for large time steps. Finally, an application case, consisting of a simplified open pool with a cooling external circuit is solved to remark the capability of the tool to simulate thermal hydraulic systems commonly found in nuclear installations.

• 한국신재생에너지학회:학술대회논문집
• /
• 2005.06a
• /
• pp.701-704
• /
• 2005

Aquifer Thermal Energy Storage (ATES) can be a cost-effective and renewable geothermal energy source, depending on site-specific and thermohydraulic conditions. To design an effective ATES system having influenced by groundwater movement, understanding of thermo hydraulic processes is necessary. The heat transfer phenomena for an aquifer heat storage are simulated using FEFLOW with the scenario of heat pump operation with pumping and waste water reinjection in a two layered confined aquifer model. Temperature distribution of the aquifer model is generated, and hydraulic heads and temperature variations are monitored at the both wells during 365 days. The average groundwater velocities are determined with two hydraulic gradient sets according to boundary conditions, and the effect of groundwater flow are shown at the generated thermal distributions of three different depth slices. The generated temperature contour lines at the hydraulic gradient of 0.00 1 are shaped circular, and the center is moved less than 5m to the groundwater flow direction in 365 days simulation period. However at the hydraulic gradient of 0.01, the contour center of the temperature are moved to the end of boundary at each slice and the largest movement is at bottom slice. By the analysis of thermal interference data between two wells the efficiency of the heat pump system model is validated, and the variation of heads is monitored at injection, pumping and no operation mode.

• Structural Engineering and Mechanics
• /
• v.80 no.5
• /
• pp.553-562
• /
• 2021

In this paper, a hybrid scaled boundary finite element and finite volume method (SBFEM-FVM) is proposed for simulating hydraulic-fracture propagation in brittle concrete materials. As a semi-analytical method, the scaled boundary finite element method is introduced for modelling concrete crack propagation under both an external force and water pressure. The finite volume method is employed to model the water within the crack and consider the relationship between the water pressure and the crack opening distance. The cohesive crack model is used to analyse the non-linear fracture process zone. The numerical results are compared with experimental data, indicating that the F-CMOD curves and water pressure changes under different loading conditions are approximately the same. Different types of water pressure distributions are also studied with the proposed coupled model, and the results show that the internal water pressure distribution has an important influence on crack propagation.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.1
• /
• pp.37-47
• /
• 2008

Design and construction of tunnels require understanding the influence of groundwater. Particularly, it is essential to know how the drainage conditions at the tunnel boundary affect flow behavior of ground adjacent to the tunnels. In this study flow behavior of a leaking tunnel was investigated using physical model tests for tunnel depths and various hydraulic boundary conditions. Particular concerns were given to flow lines toward tunnels. Test results showed that the boundary conditions hardly influence on flow patterns and time required to reach steady state conditions. It is revealed that with an increase in water depth, flow lines concentrated to the drain holes. The physical tests were numerically simulated. Numerical results showed that the flow behavior was represented appropriately by considering filter-drain hole drainage rather than boundary drainage all over the lining.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.10
• /
• pp.105-116
• /
• 2011

No analytical solution exists for evaluating in-situ hydraulic conductivity of vertical cutoff walls by analyzing slug test results. Recently, an analytical solution to interpret slug tests has been proposed for a partially penetrated well in an aquifer. However, this analytical solution cannot be directly applied to the cutoff wall because the solution has been developed exclusively for an infinite aquifer instead of a narrow cutoff wall. To consider the cutoff wall boundary conditions, the analytical solution has been modified in this study to take into account the narrow boundaries by introducing the imaginary well theory. Two boundary conditions are considered according to the existence of filter cakes: constant head boundary and no flux boundary. Generalized steady-state shape factors are presented for each geometric condition, which can be used for evaluating the in-situ hydraulic conductivity of cutoff walls. The constant head boundary condition provides higher shape factors and no flux boundary condition provides lower shape factors than the infinite aquifer, which enables to adjust the in-situ hydraulic conductivity of the cutoff wall. The hydraulic conductivities calculated from the analytical solution in this paper give about 1.2~1.7 times higher than those from the Bouwer and Rice method, one of the semi-empirical formulas. Considering the compressibility of the backfill material, the analytical solution developed in this study was proved to correspond to the case of incompressible backfill materials.

• Journal of the Korean Geotechnical Society
• /
• v.37 no.12
• /
• pp.25-31
• /
• 2021

Boundary conditions for thermal-hydraulic problems of soils play an essential role in the numerical accuracy. This study presents a boundary condition considering the thermo-hydraulic interaction between the ground and the atmosphere. Ground surface energy balance consists of solar radiation, ground radiation, wind convection, latent heat from water evaporation, and heat conduction to the ground. Equations for each heat flux are presented, and numerical analyses are performed in conjunction with the FEM program for the thermal-hydraulic phenomenon of unsaturated soils. Numerical results using the weather data at the Ulsan Meteorological Observatory are similar to the measured surface temperature. Latent heat caused by water evaporation during the daytime lowers the surface temperature of the bare soil, and a thermal equilibrium is reached at nighttime when the effect of the ground condition is significantly reduced. The temperature change of the surface ground is diminished at the deeper ground due to its thermal diffusion. Numerical analysis where the surface ground temperature is the primary concern requires considering the thermo-hydraulic interaction between the ground and the atmosphere.

• Journal of Soil and Groundwater Environment
• /
• v.25 no.4
• /
• pp.87-97
• /
• 2020

In this study, the applicability of the geostatistical evolution strategy as an inverse analysis method of estimating hydraulic properties of small-scale basin was tested. The geostatistical evolution strategy is a type of data assimilation method that can effectively estimate aquifer hydraulic conductivity by combining a global optimization model of the evolution strategy and a local optimization model of the ensemble Kalman filtering. In the applicability test, the geometry, hydraulic boundary conditions, and the distribution of groundwater monitoring wells of Hanlim-Eup were employed. On the other hand, a synthetic hydraulic conductivity distribution was generated and used as the reference property for ease of estimation quality assessment. In the estimations, two different cases were tested where, in Case I, both groundwater levels and hydraulic conductivity measurements were assumed to be available, and only the groundwater levels were available, in Case II. In both cases, the reference and estimated hydraulic conductivity fields were found to show reasonable similarity, even though the prior information for estimation was not accurate. The ability to estimate hydraulic conductivity without accurate prior information suggests that this method can be used effectively to estimate mathematical properties in real-world cases, many of which little prior information is available for the aquifer conditions.