• Geomechanics and Engineering
• /
• v.9 no.6
• /
• pp.775-791
• /
• 2015

Due to the extension of communication ways (metro, highways, railways), hence, to improve traffic flow imposes often the difficult crossing that generally drive to the construction of underground works (tunnel, water conveyance tunnel...) plays a major role in the redevelopment of urban areas. This study is focused on the assessment of the interaction response of parallel tunnels, so this study uses the results from the simulation of two tunnels to illustrate a few observations that may aid in practical designs. In this article, simultaneous drilling of highway's twin tunnels is simulated by means of Finite Element Method (FEM) implemented in Plaxis program. So the treated subject appears in a setting of geotechnical where one can be to construct several tunnels sometimes in a ground of weak mechanical characteristics. The objective of this study is to simulate numerically the interaction effects caused by construction of two parallels tunnels. This is an important factor in the study of the total answer of the problem interaction between parallels underground works. The importance of the effects transmitted is function of several parameters as the type of the works, and the mechanical characteristics (tunnel size, depth, and the relative position between two tunnels, lining thickness...). This article describes numerical analyses of two parallels tunnels interaction. This study will be applied to a real case of a section tunnel T4 of the highway East-West (Algeria); the study presented below comprises a series of numerical simulations of two tunnels using the computer program Plaxis which is used in the analyses is based on Finite Element Method.

• Journal of Soil and Groundwater Environment
• /
• v.28 no.6
• /
• pp.9-15
• /
• 2023

In order to properly evaluate the spatio-temporal variations of groundwater flow, the data obtained in field experiments should be corroborated into numerical simulations. Particle tracking method is a simple simulation tool often employed in groundwater simulation to predict groundwater flow paths or solute transport paths. Particle tracking simulations visually show overall the particle flow path along the entire aquifer, but no previous simulation studies has yet described the parameter values at grid nodes around the particle path. Therefore, in this study, a new technical approach was proposed that enables acquisition of parameters associated with particle transport in grid nodes distributed in the center of the particle path in groundwater. Since the particle tracking path is commonly referred to as streamline, the algorithm and codes developed in this works designated streamline analysis method. The streamline analysis method can be applied in two-dimensional and three-dimensional finite element or finite difference grid networks, and can be utilized not only in the groundwater field but also in all fields that perform numerical modeling.

• Applied Chemistry for Engineering
• /
• v.30 no.4
• /
• pp.445-452
• /
• 2019

Shear thickening is an intriguing phenomenon in the fields of chemical engineering and rheology because it originates from complex situations with experimental and numerical measurements. This paper presents results from the numerical modeling of the particle-fluid dynamics of a two-dimensional mixture of colloidal particles immersed in a fluid. Our results reveal the characteristic particle behavior with an application of a shear force to the upper part of the fluid domain. By combining the lattice Boltzmann and discrete element methods with the calculation of the lubrication forces when particles approach or recede from each other, this study aims to reveal the behavior of the suspension, specifically shear thickening. The results show that the calculated suspension viscosity is in good agreement with the experimental results. Results describing the particle deviation, diffusivity, concentration, and contact numbers are also demonstrated.

• Korean Journal of Hydrosciences
• /
• v.6
• /
• pp.51-66
• /
• 1995

Various Eulerian-Lagerangian numerical models for the one-dimensional longtudinal dispersion equation are studied comparatively. In the models studied, the transport equation is decoupled into two component parts by the operator-splitting approach ; one part governing advection and the other dispersion. The advection equation has been solved using the method of characteristics following flud particles along the characteristic line and the result are interpolated onto an Eulerian grid on which the dispersion equation is solved by Crank-Nicholson type finite difference method. In solving the advection equation, various interpolation schemes are tested. Among those, Hermite interpo;ation po;ynomials are superor to Lagrange interpolation polynomials in reducing both dissipation and dispersion errors.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.1
• /
• pp.111-120
• /
• 1998

This paper is the first one of two-parted research efforts focusing on the modeling of rubber pad forming process. The rubber pad, driven by the pressurized fluid during the forming process, pushes the sheet metal to solid tool half and forms a part to final shape. In this part of the paper, a numerical procedure for the FE analysis of the rubber pad deformation is presented. The developed three-dimensional FE model is based on the total Lagrangian description of rubber maerial characterized by nearly incompressible hyper-elastic behavior under a large deformation assumption. Validity of the model as well as effects of different algorithms corresponding to incompresibility constraints and time integration methods on numerical solution responses are also demonstrated.

• Proceedings of the Korea Society for Simulation Conference
• /
• 2001.10a
• /
• pp.209-214
• /
• 2001

Acoustic field in steady-state is characterized by a Helmholtz equation. The transmission characteristics of acoustic wave devices is however influenced by the presence of the mean flow in the medium. The effect of the mean How introduces additional terms in the equation. In the present paper, two approaches are considered. One is that the equation is directly discretized by FEM for one-dimensional and the axisymmetric case. Another is that the equation is first transformed into the standard Helmholtz equation which is solved by BEM. The numerical demonstrations are made for the axisymmetric FEM and the three-dimensional BEM modeling. The numerical examination for a straight circular duct is first considered. The solutions are compared wish the analytical ones. The examination is then extended to the case when the mean How is locally present in a muser with expansion chamber.

• Structural Engineering and Mechanics
• /
• v.83 no.4
• /
• pp.435-449
• /
• 2022

Constitutive modeling that could reasonably predict and effectively evaluate the post-peak structural behavior while eliminating the mesh-dependency in numerical simulation remains to be developed for general engineering applications. Based on the previous work, a simple one-dimensional modeling procedure is proposed to predict and evaluate the post-peak response, as characterized by the evolution of localized strain field, of a steel member to monotonically uniaxial tension. The proposed model extends the classic one-dimensional softening with localization model as introduced by (Schreyer and Chen 1986) to account for the localization length, and bifurcation and rupture points. The new findings of this research are as follows. Two types of strain-softening functions (bilinear and nonlinear) are proposed for comparison. The new failure criterion corresponding to the constitutive modeling is formulated based on the engineering strain inside the localization zone at rupture. Furthermore, a new mathematical expression is developed, based on the strain rate inside and outside the localization zone, to describe the displacement field at which bifurcation occurs. The model solutions are compared with the experimental data on four low-carbon cylindrical steel bars of different lengths. For engineering applications, the model solutions are also compared to the experimental data of a cylindrical steel bar system (three steel bars arranged in series). It is shown that the bilinear and nonlinear softening models can predict the energy dissipation in the post-peak regime with an average difference of only 4%.

• Computers and Concrete
• /
• v.16 no.5
• /
• pp.759-774
• /
• 2015

This study simulates the flexural behavior of ultra-high-performance fiber-reinforced concrete (UHPFRC) beams reinforced with steel and glass fiber-reinforced polymer (GFRP) rebars. For this, micromechanics-based modeling was first carried out on the basis of single fiber pullout models considering inclination angle. Two different tension-softening curves (TSCs) with the assumptions of 2-dimensional (2-D) and 3-dimensional (3-D) random fiber orientations were obtained from the micromechanics-based modeling, and linear elastic compressive and tensile models before the occurrence of cracks were obtained from the mechanical tests and rule of mixture. Finite element analysis incorporating smeared crack model was used due to the multiple cracking behaviors of structural UHPFRC beams, and the characteristic length of two times the element width (or two times the average crack spacing at the peak load) was suggested as a result of parametric study. Analytical results showed that the assumption of 2-D random fiber orientation is appropriate to a non-reinforced UHPFRC beam, whereas the assumption of 3-D random fiber orientation is suitable for UHPFRC beams reinforced with steel and GFRP rebars due to disorder of fiber alignment from the internal reinforcements. The micromechanics-based finite element analysis also well predicted the serviceability deflections of UHPFRC beams with GFRP rebars and hybrid reinforcements.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2006.04a
• /
• pp.293-297
• /
• 2006

Submarine groundwater discharge (SGD) and the interface between seawater and freshwater in an unconfined coastal aquifer was evaluated by numerical modeling. A two-dimensional vertical cross section of the aquifer was constructed. Coupled flow and salinity transport modeling were peformed by using a numerical code FEFLOW In this study, we investigated the changes in groundwater flow and salinity transport in coastal aquifer with hydraulic condition such as the magnitude of recharge flux, hydraulic conductivity. Especially, transient simulation considering tidal effect and seasonal change of recharge rate was simulated to compare the difference between quasi-steady state and transient state. Results show that SGD flux is in proportion to the recharge rate and hydraulic conductivity, and the interface between the seawater and the freshwater shows somewhat retreat toward the seaside as recharge flux increases. Considered tidal effect, SGD flux and flow directions are affected by continuous change of the sea level and the interface shows more dispersed pattern affected by velocity variation. The cases which represent variable daily recharge rate instead of annual average value also shows remarkably different result from the quasi-steady case, implying the importance of transient state simulation.

• ETRI Journal
• /
• v.19 no.2
• /
• pp.59-70
• /
• 1997

A study on the operation of a heat pipe with two heat sources has been performed to optimize the heat distribution of satellite equipment. A numerical modeling is used to predict the temperature profile for the heat pipe assuming cylindrical two-dimensional laminar flow for the vapor, and the conduction heat transfer for the wall and wick. An experimental study using the copper-water heat pipe with the length of 0.45 m has been performed to evaluate the numerical model and to compare the temperature distribution at the outer wall for the non-uniform heat distribution. The results on temperature profiles for the heat input range from 29 W to 47 W on each heater are presented. Also the correlation between the heat input and the temperature increase is presented for the optimum distribution on two heaters. The result shows that the outer wall temperature can be controlled by redistribution of heat sources. It is also concluded that the heat source closer to the condenser can carry more heat while maintaining lower temperatures at the outer wall.