• Geomechanics and Engineering
• /
• v.23 no.5
• /
• pp.455-466
• /
• 2020

This paper presents analysis of the interaction between tunnel and Qanat with a particular interest for the optimization of Qanat shape using the discrete element code, PFC2D, and the results will be compared with the FEM results of PLAXIS2D. For these concerns, using software PFC2D based on Discrete Element Method (DEM), a model with dimension of 100m ^⁎ 100 m was prepared. A circular tunnel with dimension of 9 m was situated 20 m below the ground surface. Also one Qanat was situated perpendicularly above the tunnel roof. Distance between Qanat center and ground surface was 8 m. Five different shapes for Qanat were selected i.e., square, semi-circular, vertical ellipse, circular and horizontal ellipse. Confining pressure of 5 MPa was applied to the model. The vertical displacement of balls situated in ground surface was picked up to measure the ground subsidence. Also two measuring circles were situated at the tunnel roof and at the Qanat roof to check the vertical displacements. The properties of the alluvial soil of Tehran city are: γ_dry=19 (KN/㎥), E= 750 (kg/㎠), ν=0.35, c=0.3(kg/㎠), φ=34°. In order to validate the DEM results, a comparison between the numerical results (obtained in this study) and analytical and field monitoring have been done. The PFC2D results are compared with the FEM results. The results shows that when Qanat has rectangular shape, the tensile stress concentration at the Qanat corners has maximum value while it has minimum value for vertical ellipse shape. The ground subsidence for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. The vertical displacements at the tunnel roof for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. Historical shape of Qante approved the finding of this research.

• Proceeding of Spring/Autumn Annual Conference of KHA
• /
• 2006.11a
• /
• pp.316-321
• /
• 2006

Multiplex use facility emphasized the importance of indoor air environment because many people using that place. People using parking lot of multiplex use facility in little time, but a pollutant of in indoor air menacing user health. This study measure and analyze the CO, $CO_2$ for better air quality of an indoor parking lot in the multiplex user facility. As a result, two-way opened above ground parking lot is higher numerical value of CO, $CO_2$ than an underground parking lot. And a parking lot of square form is higher numerical value of $CO_2$ than a parking lot of rectangle form.

• Publications of The Korean Astronomical Society
• /
• v.15 no.spc2
• /
• pp.65-71
• /
• 2000

We have investigated the numerical methods to calculate model atmosphere for the analysis of spectral lines emitted from the sun and stars. Basic equations used in our calculations are radiative transfer, statistical equilibrium and charge-particle conservations. Transfer equation has been solved to get emitting spectral line profile as an initial value problem using Adams-Bashforth-Moulton method with accuracy as high as 12th order. And we have calculated above non linear differential equations simultaneously as a boundary value problem by finite difference method of 3 points approximation through Feautrier elimination scheme. It is found that all computing programs coded by above numerical methods work successfully for our model atmosphere.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.10a
• /
• pp.243-249
• /
• 2003

The head-on collision of two solitary waves are examined using a boundary element method. Attachment, detachment times and alplitudes and maximum run-up times and amplitudes are computed. Consolidation times show local minimum value if two waves are of equal amplitudes are colliding. Attachment times show local maximum value if the amplitudes of two waves are the same. The detachment time show local maximum if two wves are the same. The detachment amplitude show local minimum values if the amplitude e(=a/h) is greater than 0.3.

• Journal of the Korean Graphic Arts Communication Society
• /
• v.28 no.1
• /
• pp.43-50
• /
• 2010

The 200mesh screen was tensioned and fixed on a frame. When applied tension to 5N and 10N per unit area, a side length of opening of the screen was 163.223${\mu}m$ and 168.224${\mu}m$, respective. But side length not tensioned was 158.879${\mu}m$. We knew that a side length of opening of the screen rarely changes with tension applied to the screen. The appearances that a side length of opening of the screen expand little are due to a decreasing diameter of thread by means of tension. In a thickness measurement of screen, While the high density mesh screen that had a lot of knots that crossed a line of latitude and longitude per unit area appeared a higher numerical value, the low density mesh screen that had a few knots appeared a low numerical value.

• Journal of applied mathematics & informatics
• /
• v.28 no.1_2
• /
• pp.411-423
• /
• 2010

In this paper, a numerical method for singular perturbation problems arising in chemical reactor theory for general singularly perturbed two point boundary value problems with boundary layer at one end(left or right) of the underlying interval is presented. The original second order differential equation is replaced by an approximate first order differential equation with a small deviating argument. By using the trapezoidal formula we obtain a three term recurrence relation, which is solved using Thomas Algorithm. To demonstrate the applicability of the method, we have solved four linear (two left and two right end boundary layer) and one nonlinear problems. From the results, it is observed that the present method approximates the exact or the asymptotic expansion solution very well.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.19 no.3
• /
• pp.4462-4471
• /
• 1977

The governing differential equations for the tapered circular arch with fixed ends have been derived, and a numerical procedure for the solution of these equations have been developed. The governing differential equations were solved numerically by an initial value integration procedure and Shooting Methods for boundary value problems. The Rungekutta fourth order integration technique was used. The methods was programmed for a Cyber 73-18 computer System, and all esults were obtained on this computer. A detailed study has been made for a fixed arch with an angle of opening equal to 0.7 radian, and the results are presented in detail in tables and curves. It is hoped that the results presented herein is applied to the deformations of gives point from the tri-axial direction of tapered circular arch with fixed ends, bending moment, and torsional moment, and that at the same time results to be used for archwise structures in steel structure.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.1
• /
• pp.29-41
• /
• 1996

The two-stage numerical model was used to study the relation between three-dimensional local wind model, advection/diffusion model of random walk method and second moment method on Pusan coastal area. The first stage is three dimensional time-dependent local wind model which gives the wind field and vertical dirrusion coefficient. The second stage is advection/diffusion model which uses the results of the first stage as input data. First, wind fields on Pusan coastal area for none synoptic scale wind showed typical land and sea breeze circulation, and convergence zone occured at 1200LST in northern of domain, in succession, moved northward of domain. Emissions from Sinpyeong industrial district were trasnported toward the inland by sea breeze during daytime, and reached the end part of domain about 1800LST. During nighttime, emissions return to sea by land breeze and vertical diffusion also contributes to upward transport. In order to use this model for forecast of air pollution concentration on the Pusan coastal area, it is necessary that computed value must be compared with measured value and wind fields model must also be dealt in detail.

• Tribology and Lubricants
• /
• v.12 no.3
• /
• pp.63-71
• /
• 1996

The connecting rod bearing, which is subjected to periodical dynamic loading, is an impoRant component of the reciprocating engine. In the operation of this bearing, significant parameters are the oil film thickness and the film pressure. Peak film pressures of 20-30 MPa are not uncommon. So the elastic deformation of the bearing housing can have a significant effect on the bearing performance. In this study, a numerical analysis of connecting rod bearing is investigated. Elastic deformation of the bearing housing is considered in the analysis. Separate hydrodynamic and structural analysis are coupled through a direct iterative process. It is shown that as the result of the elastic deformation of the bearing housing, the eccentricity ratio is increased, and the minimum value of the minimum film thickness and the maximum value of the maximum film pressure are decreased. The variations of rotational speed and cylinder pressure affect the minimum film thickness and the maximum film pressure variations of the connecting rod bearing.

• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1273-1287
• /
• 2008

In this paper, the numerical integration method for general singularly perturbed two point boundary value problems with mixed boundary conditions of both left and right end boundary layer is presented. The original second order differential equation is replaced by an approximate first order differential equation with a small deviating argument. By using the trapezoidal formula we obtain a three term recurrence relation, which is solved using Thomas Algorithm. To demonstrate the applicability of the method, we have solved four linear (two left and two right end boundary layer) and one nonlinear problems. From the results, it is observed that the present method approximates the exact or the asymptotic expansion solution very well.