• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.272-275
• /
• 2003

A mixed analytical/numerical method is developed here to solve the low Reynolds number kepsilon turbulence model. In this method the advection-diffusion part is solved numerically, while the source terms are split into two parts: one part is solved analytically and the next is solved numerically.

• Structural Engineering and Mechanics
• /
• 제3권2호
• /
• pp.135-144
• /
• 1995

Since structural systems may fail in any one of several failure modes, computation of system reliability is always difficult. A method using numerical quadrature for computing structural system reliability with either one or more than one failure mode is presented in this paper. Statistically correlated safety margin equations are transformed into a group of uncorrelated variables and the joint density function of these uncorrelated variables can be generated by using the Maximum Entropy Method. Structural system reliability is then obtained by integrating the joint density function with the transformed safety domain enclosed within a set of linear equations. The Gaussian numerical integration method is introduced in order to improve computational accuracy. This method can be used to evaluate structural system reliability for Gaussian or non-Gaussian variables with either linear or nonlinear safety boundaries. It is also valid for implicit safety margins such as computer programs. Both the theory and the examples show that this method is simple in concept and easy to implement.

• Structural Engineering and Mechanics
• /
• 제81권4호
• /
• pp.523-527
• /
• 2022

The objective of this study is to simulate the flow in pipes with various boundary conditions. Free-pressure fluid model, is used in the pipe based on Navier-Stokes equation. The models are solved by using the numerical method. A problem called "stability of pipes" is used in order to compare frequency and critical fluid velocity. When the initial conditions of problem satisfied the instability conditions, the free-pressure model could accurately predict discontinuities in the solution field. Employing nonlinear strains-displacements, stress-strain energy method the governing equations were derived using Hamilton's principal. Differential quadrature method (DQM) is used for obtaining the frequency and critical fluid velocity. The results of this paper are analyzed by hyperbolic numerical method. Results show that the level of numerical diffusion in the solution field and the range of well-posedness are two important criteria for selecting the two-fluid models. The solutions for predicting the flow variables is approximately equal to the two-pressure model 2. Therefore, the predicted pressure changes profile in the two-pressure model is more consistent with actual physics. Therefore, in numerical modeling of gas-liquid two-phase flows in the vertical pipe, the present model can be applied.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2001년도 봄 학술발표회 논문집
• /
• pp.105-112
• /
• 2001

The purpose of this study is to assess the stability of tunnel for a high speed railway crossing the fault zone. The area where the tunnel crossed the fault zone can be unstable during construction and operation. Geotechnical investigations have been conducted to determine an optimum excavation method by obtaining the material properties around the fault zone and to check the stability of the tunnel. For the numerical analysis, the FLAC, numerical analysis code based on finite difference method, was utilized to analyze the behavior of the fault at three points having typical ground conditions. Based on the results of numerical analysis, the combinations of compaction grouting and LW grouting were determined as suitable methods for pre-excavation Improvement of the ground surrounding the tunnel opening. In conclusion, the stability of the tunnel construction for the high speed railway within the fault zone may be obtained by adopting the optimum excavation method and the reinforcement method. The numerical analysis based on FLAC program contains errors caused by assumptions used in numerical analysis, therefore constant monitoring with respect to the change of ground condition and groundwater is highly recommended to minimize the numerical error and the possibility of damage to tunnel.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2000년도 봄 학술발표회 논문집
• /
• pp.467-474
• /
• 2000

The application of Terzaghi's theory of consolidation for analysing the settlement of multi-layered soils is not strictly valid because the theory involves an assumption that the soil is homogeneous. The settlement of stratified soils with confined aquifer can be analysed using numerical techniques whereby the governing differential equation is replaced by 2-dimensional finite difference approximations. The problems of discontinuous layer interface are very important in the algorithm and programming for the analysis of multi-layered consolidation using a numerical analysis, finite difference method(F.D.M.). Better results can be obtained by the process for discontinuous layer interface, since it can help consolidation analysis to model the actual ground The purpose of this paper provides an efficient computer algorithm based on numerical analysis using finite difference method(F.D.M) which account for multi-layered soils with confined aquifer to determine the degree of consolidation and excess pore pressures relative to time and positions more realistically.

• 대한기계학회논문집A
• /
• 제21권10호
• /
• pp.1702-1711
• /
• 1997

A numerical study for the number of terms of a power series stress function and the effect of hybrid element size on stress analysis around a hole in loaded orthotropic composites is presented. The hybrid method coupling experimental and/or theoretical inputs and complex variable formulations involving conformal mappings and analytical continuity is used to calculate tangential stress on the boundary of the hole in uniaxially loaded, finite width glass epoxy tensile plate. The tests are done by rarying the number of terms, element size and nodal locations on the external boundary of the hybrid region. The numerical results indicate that the hybrid method is accurate and powerful in both experimental and numerical stress analysis.

• Nuclear Engineering and Technology
• /
• 제44권8호
• /
• pp.847-854
• /
• 2012

Despite extensive research efforts, the mechanism of the nucleate boiling phenomena is still not clear. A direct numerical simulation of the boiling phenomena is one of the promising approaches in order to clarify its heat transfer characteristics and discuss their mechanism. Therefore, many DNS procedures have been developed based on recent highly advancing computer technologies. This brief review focuses on the state of the art in direct numerical simulation of the pool boiling phenomena over the past two decades. In this review, the fundamentals of the boiling phenomena and the bubble departure and micro-layer models are briefly introduced, and then the numerical procedures for tracking or capturing interface/surface shape such as the front tracking method, level set method, volume of fluid treatments, and other methods (Lattice Boltzmann method, phase-field method and so on) are briefly reviewed.

• 한국가시화정보학회:학술대회논문집
• /
• 한국가시화정보학회 2006년도 추계학술대회 논문집
• /
• pp.86-91
• /
• 2006

In this parer, we present the bubble forming and motion in the micro channel by using the two-dimensional numerical computation and experiment. In the numerical computation, The Lattice Boltzmann method(LBM) and free-energy model is used to treat the interfacial force and deformation of binary fluid system, drawn in to a micro channel and a numerical simulation is carried out by using the parallel computation method. The urn in this investigation is to examine the applicability of LBM to numerical analysis and experimental method of binary fluid separation and motion in the micro channel.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제27권4호
• /
• pp.272-280
• /
• 2023

We study the numerical analysis for the Cahn-Hilliard (CH) equation using the decoupled projection (DP) method. The CH equation is a fourth order nonlinear partial differential equation that is hard to solve. Therefore, various of numerical schemes have been proposed to solve the CH equation. To verify the relation of each existing scheme for the CH equation, we consider the DP method for linear convex splitting schemes. We present the numerical experiments to demonstrate our analysis. Throughout this study, it is expected to construct a novel numerical scheme using the relation with existing numerical schemes.

• Journal of information and communication convergence engineering
• /
• 제12권1호
• /
• pp.39-45
• /
• 2014

A numerical method for solving Hamiltonian equations is said to be symplectic if it preserves the symplectic structure associated with the equations. Various symplectic methods are widely used in many fields of science and technology. A symplectic method preserves an approximate Hamiltonian perturbed from the original Hamiltonian. It theoretically supports the effectiveness of symplectic methods for long-term integration. Although it is also related to long-term integration, numerical stability of symplectic methods have received little attention. In this paper, we consider explicit symplectic methods defined for Hamiltonian equations with Hamiltonians of the special form, and study their numerical stability using the harmonic oscillator as a test equation. We propose a new stability criterion and clarify the stability of some existing methods that are visually based on the criterion. We also derive a new method that is better than the existing methods with respect to a Courant-Friedrichs-Lewy condition for hyperbolic equations; this new method is tested through a numerical experiment with a nonlinear wave equation.