• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.8
• /
• pp.734-741
• /
• 2013

In contrast of external excitations, parametric excitations can produce a large response when the excitation frequency is away from the linear natural frequencies. The Mathieu equation is the simplest differential equation with periodic coefficients, which lead to the parametric excitation. The Mathieu equation may have the unbounded solutions. This work conducted the stability analysis for the Mathieu equation, using Floquet theory and numerical method. Using Lindstedt's perturbation method, harmonic solutions of the Mathieu equation and transition curves separating stable from unstable motions were obtained. Using Floquet theory with numerical method, stable and unstable regions were calculated. The numerical method had the same transition curves as the perturbation method. Increased stable regions due to the inclusion of damping were calculated.

• The Journal of the Acoustical Society of Korea
• /
• v.17 no.3E
• /
• pp.35-42
• /
• 1998

Transient acoustic and elastic wave propagation in inhomogeneous media are studied by using the Fourier method. It is known that the fourier method has advantages in memory requirements and computing speed over conventional methods such as FDM and FEM, because the Fourier method needs less grid points for achieving the same accuracy. To verify the proposed numerical scheme, several examples having analytic solutions are considered, where two different semi-infinite media are in contact along a plane boundary. The comparisons of numerical results by the Fourier method and analytic solutions show good agreements. In addition, the fourier method is applied to a layered half-plane, in which an elastic semi-infinite medium is covered by an elastic layer of finite thickness. It is showed how to derive the analytic solutions by using the Cagniard-de Hoop method. The numerical solutions are in excellent agreements with analytic results.

• Journal of Advanced Marine Engineering and Technology
• /
• v.33 no.6
• /
• pp.852-859
• /
• 2009

For the analysis of solids containing a number of microinclusions or microvoids, in which the mechanical effect of each inclusion or void, a numerical approach is need to be developed to understand the mechanical behavior of damaged solids containing these defects. In this study, the simulation method using the natural element method is proposed for the analysis of effective elastic moduli. The mechanical effect of each inclusion or void is considered by controlling the material constants for Gaussian points. The relationship between area fraction of microinclusions or microvoids and effective elastic moduli is studied to verify the validity of the proposed method. The obtained results are in good agreement with the theoretical results such as differential method, self-consistent method, Mori-Tanaka method, as well as the numerical results by rigid body spring model.

• International journal of steel structures
• /
• v.18 no.5
• /
• pp.1699-1709
• /
• 2018

Suspen-dome structures are extensively used due to their superiority over traditional structures. The friction between cable and joints may severely influence the distribution of cable force, especially during the pre-stressing construction period. An accurate and efficient numerical method has not yet been developed that can be used for estimating the influence of friction on cable force distribution. Thus, this study proposes an efficient friction element to simulate friction between cable and joint. A flowchart for estimating the value of friction force is introduced. These novel numerical methods were adopted to estimate the influence of friction on cable force distribution. The accuracy and efficiency of these numerical methods were validated through numerical tests.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.808-813
• /
• 2007

We present the numerical characteristics of a new true stress-strain curve acquisition method over a large range of strains by the tensile test and a finite element method through comparing the results obtained by various finite element mesh systems. The method is introduced in detail. The effects of the finite element mesh systems on the results are investigated to show its numerical characteristics of the new method. It is shown that the method is quite robust, implying that it can be used as a special function of the tensile test machines.

• Journal of applied mathematics & informatics
• /
• v.20 no.1_2
• /
• pp.611-621
• /
• 2006

This paper is to discuss the numerical conformal mapping from the unit disk onto Jordan region, which can be solved by Theodorsen equation. Wegmann's method has been known as the most efficient one for the Theodorsen equation. However, we found divergence through numerical experiments by the iterative method of Wegmann. The divergence occurs especially when some degree of difficulty is high. We analyze the cause of divergence and propose an improved method by applying a low frequency pass filter to Wegmann's method. By this proposed method we can get a stable convergence for all the problems which was unstable with the Wegmann's method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.6
• /
• pp.443-453
• /
• 2009

The Rayleigh-Taylor instability, the bubble rising in both partially and fully filled containers and the droplet splash are simulated by an in-house solution code(PowerCFD), which are typical benchmark problems among multiphase flows with material interface due to density difference. The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present method simulates efficiently and accurately complex free surface flows such as multiphase flows with material interface due to both density difference and instability.

• Transactions of the Korean Society of Pressure Vessels and Piping
• /
• v.11 no.2
• /
• pp.31-39
• /
• 2015

Pipe wall-thinning by flow-accelerated corrosion (FAC) and various types of erosion is a significant and costly damage phenomenon in secondary piping systems of nuclear power plants (NPPs). Most NPPs have management programs to ensure pipe integrity due to wall-thinning that includes periodic measurements for pipe wall thicknesses using ultrasonic tests (UTs). Nevertheless, thinning evaluations are not easy because the amount of thickness reduction being measured is often quite small compared to the accuracy of the inspection technique. U.S. Electric Power Research Institute (EPRI) had proposed Total Point Method (TPM) as a thinning occurrence evaluation method, which is a very useful method for detecting locally thinned pipes or fittings. However, evaluation engineers have to discern manually the measurement data because there are no numerical algorithm for TPM. In this study, numerical algorithms were developed based on non-parametric and parametric statistical method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1335-1342
• /
• 2000

Three-dimensional flow analysis and numerical optimization methods are presented for the design of an axial-flow fan. Steady, incompressible, three-dimensional Reynolds-averaged Navier-Stokes equations are used as governing equations, and standard k- ${\varepsilon}$ turbulence model is chosen as a turbulence model. Governing equations are discretized using finite volume method. Steepest descent method, conjugate gradient method and BFGS method are compared to determine the searching directions. Golden section method and quadratic fit-sectioning method are tested for one dimensional search. Objective function is defined as a ratio of generation rate of the turbulent kinetic energy to pressure head. Two variables concerning sweep angle distribution are selected as the design variables. Performance of the final fan designed by the optimization was tested experimentally.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.