• Journal of applied mathematics & informatics
• /
• v.26 no.1_2
• /
• pp.45-60
• /
• 2008

Here we consider the evaluation of the the dynamic component of the second order force due to wave diffraction by a circular cylinder analytically and numerically. The cylinder is fixed, vertical, surface piercing in water of finite uniform depth. The formulation of the wave-structure interaction is based on the assumption of a homogeneous, ideal, incompressible, and inviscid fluid. The nonlinearity in the wave-structure interaction problem arises from the free surface boundary conditions, namely, dynamic and kinematic free surface boundary conditions. We expand the velocity potential and free surface elevation functions in terms of a small parameter and then consider the second order diffraction problem. After deriving the pressure using Bernoulli's equation, we obtain the analytical expression for the dynamic component of the second order force on the cylinder by integrating the pressure over the wetted surface. The computation of the dynamic force component requires only the first order velocity potential. Numerical results for the dynamic force component are presented.

• Fibers and Polymers
• /
• v.7 no.4
• /
• pp.424-431
• /
• 2006

Roll drafting, a mechanical operation attenuating fiber bundles to an appropriate thickness, is an important operation unit for manufacturing staple yams. It influences not only the linear density regularity of the slivers or staple yams that are produced, but also the quality of the textile product and the efficiency of the thereafter processes. In this research, the dynamic states of the fiber bundle in the roll drafting zone were analyzed by simulation, based on the mathematical model that describes the dynamic behavior of the flowing bundle. The state variables are the linear density and velocity of the fiber bundles and we simulated the dynamics states of the bundle flow, e.g., the profiles of the linear density and velocity in the draft zone for various values of the model parameters and boundary conditions, including the initial conditions to obtain their influence on the dynamic state. Results showed that the mean velocity profile of the fiber bundle was strongly influenced by draft ratio and process speed, while the input sliver linear density has hardly affected the process dynamics. Velocity variance of individual fibers that could be supposed to be a disturbing factor in drafting was also influenced by the process speed. But the major disturbance occurred due to the velocity slope discontinuity at the front roll, which was strongly influenced by the process speed. Thickness of input sliver didn't play any important role in the process dynamics.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• 한국항공운항학회:학술대회논문집
• /
• 2004.11a
• /
• pp.223-226
• /
• 2004

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.6
• /
• pp.2267-2275
• /
• 2013

Existing conventional model for analysis of shallow water flow just assumed the internal boundary condition as free-slip, which resulted in the wrong prediction about the velocity, vorticity, water level, shear stress distribution, and time variation of drag and lift force around a structure. In this study, a finite element model that can predict flow characteristics around the structure accurately was developed and internal boundary conditions were generalized as partial slip condition using slip length concept. Laminar flow characteristics behind circular cylinder were analyzed by varying the internal boundary conditions. The simulation results of (1) time variations of longitudinal and transverse velocities, and vorticity; (2) wake length; (3) vortex shedding phenomena by slip length; (4) and mass conservation showed that the vortex shedding had never observed and laminar flow like creeping motion was occurred under free-slip condition. Assignment of partial slip condition changed the velocity distribution on the cylinder surface and influenced the magnitude of the shear stress and the occurrence of vorticity so that the period of vortex shedding was reduced compared with the case of no slip condition. The maximum mass conservation error occurred in the case of no slip condition, which had the value of 0.73%, and there was 0.21 % reduction in the maximum mass conservation error by changing the internal boundary condition from no slip to partial slip condition.

• 한국연소학회:학술대회논문집
• /
• 2006.04a
• /
• pp.32-38
• /
• 2006

A premixed flame propagating in a tube suffers strong variation in its shape and structure depending on boundary conditions. The effects of thermal boundary conditions and flow fields on flame propagation are numerically investigated. Navier-Stokes equations and species equations are solved with a one-step irreversible global reaction model of methane-air mixture. Finite volume method using an adaptive grid method is applied to investigate the flame structure. In the case of an adiabatic wall, friction force on the wall significantly affected the flame structure while in the case of an isothermal wall, local quenching near the wall dominated flame shapes and propagation. In both cases, variations of flow fields occurred not only in the near field of the flame but also within the flame itself, which affected propagation velocities. This study provides an overview of the characteristics of flames in small tubes at a steady state.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.701-708
• /
• 1997

A direct boundary element method(DBEM) is developed for thin bodies whose surfaces are rigid or compliant. Th eHelmholtz integral equation and its normal derivative integral equation are adopted simultaneously to calculate the pressure on both sides of the thin body, instead of the jump values across it, to account for the different surface conditions of each side. Unlike the usual assumption, the normal velocity is assumed to be discontinuous across the thin body. In this approach, only the neutral surface of the thin body has to be discontinuous across the thin body. In this approach, only the neural surface of the thin body has to be discretized. The method is validated by comparison with analytic and/or numerical results for acoustic scattering and radiation from several surface conditions of the thin body; the surfaces are rigid when stationary or vibrating, and part of the interior surface is lined with a sound-absorbing material.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.31 no.5
• /
• pp.479-486
• /
• 2011

In image-based CFD modeling of carotid bifurcation hemodynamics, it is often not possible (or at least not convenient) to impose measured velocity profiles at the common carotid artery inlet. Instead, fully-developed velocity profiles are usually imposed based on measured flow rates. However, some studies reported a pronounced influence of inflow boundary conditions that were based on actual velocity profiles measured by magnetic resonance imaging which showing the unusual presence of a high velocity band in the middle of the vessel during early diastole inconsistent with a Dean-type velocity profile. We demonstrated that those velocity profiles were induced by the presence of modest secondary curvature of the inlet and set about to test whether such more "realistic" velocity profiles might indeed have a more pronounced influence on the carotid bifurcation hemodynamics. We found that inlet boundary condition with axisymmetric fully-developed velocity profile(Womersley flow) is reasonable as long as sufficient CCA inlet length of realistic geometry is applied.

• Journal of the Korean Society of Mechanical Technology
• /
• v.13 no.2
• /
• pp.55-61
• /
• 2011

This study analyzed on the characteristics of temperature distribution in an active regeneration DPF using computer simulation. In order to verify the boundary condition of analysis, results of temperature distribution in DPF are compared between experimental and computer simulation. Using this boundary condition, temperature distribution and filter's durability in DPF analyzed according to various operating conditions. The results of computational analysis are agreed well with experimental ones from the tendency of temperature distribution of axis and radius direction. The temperature increases and the axial temperature gradients in DPF according to velocity of exhaust gas are lowered as the high velocity of exhaust gas. But the temperature gradients of radius direction at exit side in DPF are grown as the high velocity of exhaust gas. The results according to inlet temperature of exhaust gas show that the increase ratios of temperature in DPF are grown as the high temperature of exhaust gas.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.10a
• /
• pp.81-88
• /
• 2000

The elastic wave equation is solved using the finite-difference method in 3D space to simulate the seismic wave propagation. It is based on the velocity-stress formulation of the equation of motion on a staggered grid. The nonreflecting boundary conditions are used to attenuate the wave field close to the numerical boundary. To satisfy the stress-free conditions at the free-surface boundary, a new formulation combining the zero-stress formalism with the vacuum one is applied. The effective media parameters are employed to satisfy the traction continuity condition across the media interface. With use of the moment-tensor components, the wide range of source mechanism parameters can be specified. The numerical experiments are carried out in order to test the applicability and accuracy of this scheme and to understand the fundamental features of the wave propagation under the generalized elastic media structure. Computational results show that the scheme is sufficiently accurate for modeling wave propagation in 3D elastic media and generates all the possible phases appropriately in under the given heterogeneous velocity structure. Also the characteristics of the ground motion in an sedimentary basin such as the amplification, trapping, and focusing of the elastic wave energy are well represented. These results demonstrate the use of this simulation method will be helpful for modeling the ground motion of seismological and engineering purpose like earthquake hazard assessment, seismic design, city planning, and etc..