• 한국전산구조공학회논문집
• /
• 제23권6호
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• 한국지반환경공학회 논문집
• /
• 제12권7호
• /
• pp.57-65
• /
• 2011

본 연구에서는 3차원 개별요소해석 코드인 Particle Flow Code, $PFC^{3D}$(Itasca)를 이용, 조립재료의 실내 삼축압축시험에 대한 개별요소 수치 모델링을 수행하였으며, 해석 모델과 개별요소를 대상으로 다양한 상사 조건에 대한 개별요소 수치 모델링을 수행, 그 결과를 통해 각각의 스케일 조건이 최종 수치 모델링 결과에 미치는 영향을 분석하였다. 3차원 개별요소 수치 모델링은 기존 2차원 모델링 대비 별도의 간극률 환산 없이 정확한 초기 조건 구현이 가능했으며, 응력-변형 및 체적변화 거동, 강도정수등에 있어 실내시험 결과와 유사한 수치 해석적 예측이 가능하였다. 해석 모델과 개별요소에 대한 다양한 상사비 조건별 수치 모델링 결과, 3차원 해석 시의 안정적 예측결과 도출 및 수치 시험실 활용에 대한 적정성을 확보하되, 해석시간 단축 및 해석 효율성 확보를 위해서는 해석 모델과 개별요소에 대한 적정 상사비 결정이 필요함을 알 수 있었다. 해석 모델의 크기와 개별요소의 입경크기를 변화시켜 개별요소 수치모델링을 수행한 결과, 대부분의 경우 전체적인 응력-변형 거동에 차이가 발생하였지만, 점착력과 내부 마찰각의 강도정수는 $D_{mod}/D_{gmax}$ < 10 조건에 유사한 결과를 보였으며, 개별요소 방법이 수치 시험실 기법을 이용한 강도정수 산정에 효과적으로 적용될 수 있음을 확인할 수 있었다.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2005년도 춘계 학술발표회 논문집
• /
• pp.264-271
• /
• 2005

Sand pile method, such as sand drain method and sand compaction pile method, has been popularly used as an improvement method for soft clay grounds. The effect of accelerating consolidation of soft clay grounds has been evaluated with Barron's solution. By the way, the consolidation behavior of soft clay ground with sand piles strongly depends on both the nonlinear mechanical interaction between sand piles and surrounding clays and the degradation permeability of clays. In this paper, the method for determining consolidation parameters of soft clay ground with sand drains by using Barron's solution was proposed, through a series of numerical simulations. Through the method, the change in both volume compressibility and permeability during consolidation was reasonably evaluated.

• 대한산업공학회지
• /
• 제14권2호
• /
• pp.115-129
• /
• 1988

This paper discusses an analytical method to derive the efficiency and the average storage level of an unreliable two-machine, two-stage automatic transfer line with interstage buffer storage. Extensive numerical experiments and simulations are also conducted to see the behavior of the system under the wide range of parameter values as well as various distributions.

• East Asian mathematical journal
• /
• 제35권5호
• /
• pp.597-612
• /
• 2019

In this study, we propose a new logarithmic barrier approach to solve linear programming problem using the projective method of Karmarkar. We are interested in computation of the direction by Newton's method and of the step-size using minorant functions instead of line search methods in order to reduce the computation cost. Our new approach is even more beneficial than classical line search methods. We reinforce our purpose by many interesting numerical simulations proved the effectiveness of the algorithm developed in this work.

• Geomechanics and Engineering
• /
• 제30권5호
• /
• pp.423-435
• /
• 2022

Ground vibration generated repeatedly in blasting tunnel excavation sites is known to be one of the major hazards induced by blasting operations. Various studies have been conducted to minimize these hazards, both theoretical and empirical methods using electronic detonator, the deck charge method, the center-cut method among others Among these various existing methods for controlling the ground vibration, in this study, we investigated the cut method. In particular, we analyzed and compared the V-cut method, which is commonly used in tunnel blasting, to the double-drilled parallel method, which has recently been introduced in tunnel excavation site. To understand the rock fragmentation efficiency as well as the ground vibration controllability of the two methods, we performed in-situ field blasting tests with both cut methods at a tunnel excavation site. Additionally, numerical analysis by FLAC3D has been executed for a better understanding of fracture propagation pattern and ground vibration generation by each cut method. Ground vibration levels, by PPVs measured in field blasting tests and PPVs estimated in numerical simulations, showed a lower value in the double-drilled parallel compared with the V-cut method, although the exact values are quite different in field measurement and numerical estimation.

• Ocean Systems Engineering
• /
• 제7권4호
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Interaction and multiscale mechanics
• /
• 제5권1호
• /
• pp.27-40
• /
• 2012

The two-dimensional incompressible flow of a linear viscoelastic fluid we considered in this research has rapidly oscillating initial conditions which contain both the large scale and small scale information. In order to grasp this double-scale phenomenon of the complex flow, a multiscale analysis method is developed based on the mathematical homogenization theory. For the incompressible flow of a linear viscoelastic Maxwell fluid, a well-posed multiscale system, including averaged equations and cell problems, is derived by employing the appropriate multiple scale asymptotic expansions to approximate the velocity, pressure and stress fields. And then, this multiscale system is solved numerically using the pseudospectral algorithm based on a time-splitting semi-implicit influence matrix method. The comparisons between the multiscale solutions and the direct numerical simulations demonstrate that the multiscale model not only captures large scale features accurately, but also reflects kinetic interactions between the large and small scale of the incompressible flow of a linear viscoelastic fluid.

• 한국전산유체공학회지
• /
• 제1권1호
• /
• pp.19-25
• /
• 1996

A multiblock grid generation has been developed to be reliably used for a Navier-Stokes simulation of the turbomachinery flow-fields A multiblock structure simplifies the creation of structured H-grids about complex turbomachinery geometries and facilitate the creation of a grid in the tip flow region. The numerical algorithm adopts the combination of the algebraic and elliptic method to create the internal grids efficiently and quickly. The grid refinement process is enhanced by developing strategies to utilized Bezier curves and splines along with weighted transfinite interpolation technique and by formulating the grid-imbedding method for the viscous boundary-layer meshes. For purposes of illustration, the grid generator is applied to the high turning turbine rotor blades. Two different types of computational grids are provided to be compared with respect to the grid adaptation to the flow simulations. Extension to three-dimensions was done to show the possibility of its application to the tip-flow simulations. The grid quality of the multiblock structure is good in the passages, with gloval orthogonality and adequate smoothness.

• Bulletin of the Korean Chemical Society
• /
• 제33권3호
• /
• pp.779-789
• /
• 2012

In the kinetic theory of dense fluids the many-particle collision bracket integral is given in terms of a classical collision operator defined in the phase space. To find an algorithm to compute the collision bracket integrals, we revisit the eigenvalue problem of the Liouville operator and re-examine the method previously reported [Chem. Phys. 1977, 20, 93]. Then we apply the notion and concept of the eigenfunctions of the Liouville operator and knowledge acquired in the study of the eigenfunctions to cast collision bracket integrals into more convenient and suitable forms for numerical simulations. One of the alternative forms is given in the form of time correlation function. This form, on a further manipulation, assumes a form reminiscent of the Chapman- Enskog collision bracket integrals, but for dense gases and liquids as well as solids. In the dilute gas limit it would give rise precisely to the Chapman-Enskog collision bracket integrals for two-particle collision. The alternative forms obtained are more readily amenable to numerical simulation methods than the collision bracket integrals expressed in terms of a classical collision operator, which requires solution of classical Lippmann-Schwinger integral equations. This way, the aforementioned kinetic theory of dense fluids is made fully accessible by numerical computation/simulation methods, and the transport coefficients thereof are made computationally as accessible as those in the linear response theory.