• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.155-158
• /
• 2009

We developed a numerical method to analyze the contact-line problems, incorporating a dynamic contact angle model. We used level set method to capture free surface. The method is applied to the analysis of dynamic behavior of a droplet in DC electrowetting. The result is compared with an experimental data and result of perturbation method.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.9-10
• /
• 2006

A numerical method for simulating tree surface flows including the surface tension is presented. Numerical scheme is based an a fractional-step method with a finite volume formulation and the interface between liquid and gas is tracked by Volume of Fluid (VOF) method. Piecewise Linear Interface Calculation (PLIC) method is used to reconstruct the interface and the surface tension is considered using a Continuum Surface Force (CSF) model. Several free surface flow phenomena were simulated to show its effectiveness to find such phenomena.

• Journal of the Korean Mathematical Society
• /
• v.33 no.4
• /
• pp.865-878
• /
• 1996

We consider numerical methods for finding a solution to a nonlinear system of algebraic equations $$ (1) f(x) = 0, $$ where the function $f : R^n \to R^n$ is ain $x \in R^n$. In [10], a quasi-Gauss-Newton method is proposed and shown the computational efficiency over SQRT algorithm by numerical experiments. The convergence rate of the method has not been proved theoretically. In this paper, we show theoretically that the iterate $x_k$ obtained from the quasi-Gauss-Newton method for the problem (1) actually converges to a root by Kantorovich-type convergence analysis. We also show the rate of convergence of the method is superlinear.

• Geomechanics and Engineering
• /
• v.16 no.4
• /
• pp.399-413
• /
• 2018

This paper addresses the issue of field measurement of excavation damage zone (EDZ) and its numerical simulation method considering both excavation unloading and blasting load effects. Firstly, a 2000 m-deep rock cavern in China is focused. A detailed analysis is conducted on the field measurement data regarding the mechanical response of rock masses subjected to excavation and blasting operation. The extent of EDZ is revealed 3.6 m-4.0 m, accounting for 28.6% of the cavern span, so it is significantly larger than rock caverns at conventional overburden depth. The rock mass mechanical response subjected to excavation and blasting is time-independent. Afterwards, based on findings of the field measurement data, a numerical evaluation method for EDZ determination considering both excavation unloading and blasting load effects is presented. The basic idea and general procedures are illustrated. It features a calibration operation of damage constant, which is defined in an elasto-plastic damage constitutive model, and a regression process of blasting load using field blasting vibration monitoring data. The numerical simulation results are basically consistent with the field measurement results. Further, some issues regarding the blasting loads, applicability of proposed numerical method, and some other factors are discussed. In conclusion, the field measurement data collected from the 2000 m-deep rock cavern and the corresponding findings will broaden the understanding of tunnel behavior subjected to excavation and blasting at great depth. Meanwhile, the presented numerical simulation method for EDZ determination considering both excavation unloading and blasting load effects can be used to evaluate rock caverns with similar characteristics.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.3
• /
• pp.248-257
• /
• 2011

In the present study, we proposed a new numerical wave tank model to analyze the vertical tension-leg circular floating bodies, using a 2-D Navier-Stokes solver. An IBM(Immersed Boundary Method) capable of handling interactions between waves and moving structures with complex geometry on a standard regular Cartesian grid system is coupled to the VOF(Volume of Fluid) method for tracking the free surface. Present numerical results for the motions of the floating body were compared with existing experimental data as well as numerical results based on FAVOR(Fractional Area Volume Obstacle Representation) algorithm. For detailed examinations of the present model, the additional hydraulic experiments for floating motions and free surface transformations were conducted. Further, the versatility of the proposed numerical model was verified via the numerical and physical experiments for the general rectangular floating bodies. Numerical results were compared with experiments and good agreement was archived.

• Journal of the Korean Wood Science and Technology
• /
• v.36 no.3
• /
• pp.24-29
• /
• 2008

Numerical analysis model was used to analyse the behaviors of multiple bolted connections. Axial-bending element was supposed as basic model, and the effects of frame members and steel fasteners were classified for the behavior analysis. In the condition only two bolts were used, the traditional analytical methods, which show somewhat accuracy, have some advantages more than numerical analysis that need many time consuming, However, more many bolts were used in practical field condition. Also, it is impossible to analyse the behaviors of various bolts layouts and arrangements conditions by traditional analytical methods. Therefore, there is only numerical analysis method for the accurate behavioranalysis on the practical bolted connection condition. Therefore, numerical analysis method was applied on the various multiple bolted connections. On the result exactness and the reflection of connection condition, numerical analysis method showed the superiority more than widely used traditional empirical analysis methods as yield model.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.35 no.2
• /
• pp.41-48
• /
• 2023

In this study, simplified linear numerical method that can simulate wave generation and transformation by a moving bottom is introduced. Numerical analysis is conducted in wave number domain after continuity equation, linear dynamic and kinematic free surface boundary conditions and linear kinematic bottom boundary condition are Fourier transformed, and the results are expressed in space domain by an inverse Fourier transform. In the wavenumber domain, the dynamic free water surface boundary condition and the kinematic free water surface boundary condition are numerically calculated, and the velocity potential in the mean water level (z = 0) satisfies the continuity equation and the kinematic bottom boundary condition. Wave generation and transformation are investigated when the triangular and rectangular shape of bottoms move periodically. The results of the simplified numerical method are compared with the results of previous analytical solutions and agree well with them. Stability of numerical results according to the calculation time interval (Δt) and the calculation wave number interval (Δk) was also investigated. It was found that the numerical results were appropriate when Δt ≤ T(period)/1000 and Δk ≤ π/100.

• 한국연소학회:학술대회논문집
• /
• 2006.10a
• /
• pp.202-208
• /
• 2006

The purpose of this research is to model numerically the turbulent gas-particle flows in a rectangular chamber using Eulerian-Eulerian Method. A computer code using the ${\kappa}-{\varepsilon}-Ap$ two-phase turbulence model is developed for the numerical study. This code and the Eulerian multiphase model in FLUENT were used for the numerical simulations of the two-phase flow in a rectangular chamber. The numerical results calculated by the two different turbulent gas-particle codes have shown that the ${\kappa}-{\varepsilon}-Ap$ model results in a stronger diffusion of the flow momentum in the gas-particle turbulence interaction than the Eulerian multiphase model in FLUENT.

• Earthquakes and Structures
• /
• v.7 no.2
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2005.04a
• /
• pp.89-96
• /
• 2005

Numerical procedure of design optimality evaluation is studied for caisson structural systems. Two kinds of evaluation methods can be considered; mathematical optimality criteria method (MOCM) and numerical optimization method (NOM). The choice of the method depends on the available information of the system MOCM can be used only when the information of all function values, gradients and Lagrange multipliers is available, which may not be realistic in practice. Therefore, in this study, NOMs are applied for the structural optimality evaluation, where only design variables are necessary. To this end, Metropolis genetic algorithm (MGA) is advantageously used and applied for a standard optimization model of caisson composite breakwater. In the numerical example, cost and constraint functions are assumed to be changed from the orignal design situation and their effects are evaluated for optimality. From the theoretical consideration and numerical experience, it is found that the proposed optimality evaluation procedure with MGA-based NOM is efficient and practically applicable.