• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.803-806
• /
• 2002

The dispersion of Lagrangian fluid particles in a turbulent channel flow is studied by a direct numerical simulation. Four points Hermite interpolation in the homogeneous direction and Chebyshev polynomials in the inhomogeneous direction is adopted by assesing the acceleration of fluid particles. In order to characterize the inhomogeneous Lagrangian statistics, accurate single particle Lagrangian statistics are obtained along the wall normal direction. Integral time scales of Lagrangian velocity can be normalized by Eulerian mean shear stresses.

• Journal of Ocean Engineering and Technology
• /
• v.7 no.1
• /
• pp.33-55
• /
• 1993

선수충격파의 문제를 푸는데 있어서 Boundary Integral Method(BIM)의 여러가지 수치 해석방법이 검토되었으며, 특히 여러가지 Time stepping scheme, Green function, far-field 조건등에 따른 수치해석안정성과 정확성의 상관관계가 연구되었다. von Neumann 안정성해석과 matrix 안정성해석 등을 이용한 선형 안정성해석을 기초로하여, 수치해석방법의 안정성 여부를 체계적으로 조사할 수 있는 parameter(Free Surface Stability number)를 설정하고, 이 parameter의 변화에 따른 비선형 운동해석을 연구하였다. 그 결과 비선형성이 심하지 않은 기진파의 경우에서는 비선형 운동해석의 수치해석 안정성의 선형 수치해석 안정성과 큰 차이가 없음을 알 수 있게 된다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• Journal of Korean Society for Atmospheric Environment
• /
• v.9 no.E
• /
• pp.347-357
• /
• 1993

The Lagrangian grid-free numerical method, the discrete vortex method, was applied to solve the Navier-Stokes euqations. This method avoids the introduction of numerical viscosity swamping the real physical viscosity at high Reynolds number, unlike Eulerian method, e.g. finite difference and element methods. The boundary integral equation method for the potential flow solution was included to make the discrete vortex method more feasible for complex geometries. The fast adaptive multipole expansion method was incorporated to reduce the computational time from $O(N^2)$ to O(N) for the computations of vortex-vortex interactions. The test problems were air flow around one circular cylinder and two circular cylinders in tandem with various gaps. The numerical results were in excellent gareement with the experimental and other computational results. The applicabilty of the method was discussed with the indoor and the outdoor air pollution problems, especially the contaminant transport in the recirculation regions.

• The Bulletin of The Korean Astronomical Society
• /
• v.41 no.1
• /
• pp.50.1-50.1
• /
• 2016

The incorporation of radiation from massive stars is essential for modeling the dynamics and chemistry of star-forming clouds, yet it is a computationally demanding task for three-dimensional problems. We describe the implementation and tests of radiative transfer module due to point sources on a three-dimensional Cartesian grid in the Eulerian MHD code Athena. To solve the integral form of the radiation transfer equation, we adopt a widely-used long characteristics method with spatially adaptive ray tracing in which rays are split when sampling of cells becomes coarse. We use a completely asynchronous communication pattern between processors to accelerate transport of rays through a computational domain, a major source of performance bottleneck. The results of strong and weak scaling tests show that our code performs well with a large number of processors. We apply our radiation hydrodynamics code to some test problems involving dynamical expansion of HII regions.

• Advances in Computational Design
• /
• v.5 no.4
• /
• pp.417-443
• /
• 2020

Tunnels have been an integral part of human civilization. Due to complexity in its design and structure, the stability of underground structures under extreme loading conditions has utmost importance. Increased terrorism and geo-political conflicts have forced the engineers and researchers to study the response of underground structures, especially tunnels under blast loading. The present study has been carried out to seek the response of tunnel structures under blast load using the finite element technique. The tunnel has been considered in quartzite rock of northern India. The Mohr-Coulomb constitutive model has been adopted for the elastoplastic behaviour of rock. The rock model surrounding the tunnel has dimensions of 30 m x 30 m x 35 m. Both unlined and lined (concrete) tunnel has been studied. Concrete Damage Plasticity model has been considered for the concrete lining. Four different parameters (i.e., tunnel diameter, liners thickness, overburden depth and mass of explosive) have been varied to observe the behaviour under different condition. To carry out blast analysis, Coupled-Eulerian-Lagrangian (CEL) modelling has been adopted for modelling of TNT (Trinitrotoluene) and enclosed air. JWL (Jones-Wilkins-Lee) model has been considered for TNT explosive modelling. The paper concludes that deformations in lined tunnels follow a logarithmic pattern while in unlined tunnels an exponential pattern has been observed. The stability of the tunnel has increased with an increase in overburden depth in both lined and unlined tunnels. Furthermore, the tunnel lining thickness also has a significant effect on the stability of the tunnel, but in smaller diameter tunnel, the increase in tunnel lining thickness has not much significance. The deformations in the rock tunnel have been decreased with an increase in the diameter of the tunnel.

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.355-372
• /
• 2011

Based on the fully nonlinear velocity potential theory, the liquid sloshing in a three dimensional tank under random excitation is studied. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing scheme, B-spline curve, is applied to both the longitudinal and transverse directions of the tank to eliminate the possible saw-tooth instabilities. When the tank is undergoing one dimensional regular motion of small amplitude, the calculated results are found to be in very good agreement with linear analytical solution. In the simulation, the normal standing waves, travelling waves and bores are observed. The extensive calculation has been made for the tank undergoing specified random oscillation. The nonlinear effect of random sloshing wave is studied and the effect of peak frequency used for the generation of random oscillation is investigated. It is found that, even as the peak value of spectrum for oscillation becomes smaller, the maximum wave elevation on the side wall becomes bigger when the peak frequency is closer to the natural frequency.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.