• Journal of Ocean Engineering and Technology
• /
• v.25 no.4
• /
• pp.36-41
• /
• 2011

In the present paper, a direct forcing/fictitious domain (DF/FD) level set method is proposed to simulate the FSI (fluid-solid interaction) in two-phase flow. The main idea is to combine the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. The DF/FD method is a non-Lagrange-multiplier version of a distributed Lagrange multiplier/fictitious domain (DLM/FD) method. This method does not sacrifice the accuracy and robustness by employing a discrete ${\delta}$ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, easy implementation, and utilization of the original governing equation without modification. Simulations of various water-entry problems have been conducted to validate the capability and accuracy of the present method in solving the FSI in two-phase flow. Consequently, the present results are found to be in good agreement with those of previous studies.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.

• Tribology and Lubricants
• /
• v.32 no.1
• /
• pp.1-8
• /
• 2016

This paper focuses on the thermal deformation induced preload change in the tilting pad journal bearing, using a three-dimensional (3D) thermo-hydro-dynamic (THD) approach. Preload is considered as a critical factor in designing the tilting pad journal bearing. The initial preload measured under nil external load and nil thermal gradient is influenced by two factors, namely, the thermal deformation and elastic deformation. Thermal deformation is due to a temperature distribution in the bearing pads, whereas the elastic deformation is due to fluid forces acting on the pads. This study focuses on the changes induced in preload and film clearance due to thermal deformation. The generalized Reynolds equation is used to evaluate the force of the fluid and the 3D energy equation is used to calculate the temperature of the lubricant. The abovementioned equations are combined by establishing a relationship between viscosity and temperature. The heat transfer within the bearing pads, the lubricant, and the spinning journal is calculated using the heat flux boundary condition. The 3D Finite Element Method (FEM) is used in modeling the (1) heat conduction in the spinning journal and bearing pads, (2) thermal gradient induced thermal distortion of the spinning journal and pads, and (3) viscous shearing, and heat conduction and convection in a thin film. This evaluation method has an increased fidelity, and it can prove to be a cost-effective tool that can be used by designers to predict the dynamic behavior of a bearing.

• Tribology and Lubricants
• /
• v.33 no.5
• /
• pp.207-213
• /
• 2017

In sliding bearings, viscous friction due to high shear acting on the bearing surface raises the oil temperature. One of the mechanisms responsible for generating the load-carrying capacity in parallel surfaces is known as the viscosity wedge effect. In this paper, we investigate the effect of film-temperature boundary conditions on the thermohydrodynamic (THD) lubrication of parallel slider bearings. For this purpose, the continuity equation, Navier-Stokes equation, and the energy equation with temperature-viscosity-density relations are numerically analyzed using the commercial computational fluid dynamics (CFD) code FLUENT. Two different film-temperature boundary conditions are adopted to investigate the pressure generation mechanism. The temperature and viscosity distributions in the film thickness and flow directions were obtained, and the factors related to the pressure generation in the equation of motion were examined in detail. It was confirmed that the temperature gradients in the film and flow directions contribute heavily to the thermal wedge effect, due to which parallel slider bearing can not only support a considerable load but also reduce the frictional force, and its effect is significantly changed with the film-temperature boundary conditions. The present results can be used as basic data for THD analysis of surface-textured sliding bearings; however, further studies on various film-temperature boundary conditions are required.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.4
• /
• pp.6-11
• /
• 2009

In this study, the internal waves in two-density layered fluids were analyzed using the Numerical Wave Tank (NWT) technique in the frequency domain. The NWT is based on a two-domain Boundary Element Method with the potential fluids using the whole-domain matrix scheme. From the mathematical solution of the two-domain boundary integral equation, two different wave modes could be classified: a surface wave mode and an internal wave mode, and each mode were shown to have a wave number determined by a respective dispersion relation. The magnitudes of the internal waves against surface waves were investigated for various fluid densities and water depths. The calculated results are compared with available theoretical data.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.13 no.2
• /
• pp.83-90
• /
• 2010

The hydrodynamic interaction of two floating bodies in waves freely floating or connected by a rigid link is studied by using a boundary element method in the frequency-domain. The exact two-body hydrodynamic coefficients of added mass, wave damping and exciting force are calculated from the radiation-diffraction potential solution of the improved Green integral equation associated with the free surface Green function. The irregular frequencies in the conventional Green integral equation make it difficult to predict the physical resonance of the fluid in the gap between two bodies floating side by side. However, the improved Green integral equation employed in this study is free of irregular frequencies and always yields the exact solution of the multi-body radiation-diffraction potential boundary value problem. The 6 degree-of-freedom motions of two bodies freely floating side by side or connected parallel by a rigid link have been calculated for the incident wave frequencies ranging from 0.1 to 5 radians per second in head, left and right bow quartering seas. The 6-component load of the rigid link have also been presented.

• Structural Engineering and Mechanics
• /
• v.68 no.4
• /
• pp.443-457
• /
• 2018

In this article, a comparative study has been made to investigate the scattering behaviour of three layered structure model on torsional surface wave. For such model intermediate layer is taken as fiber reinforced composite, resting over a dry sandy Gibson substratum and underlying by different anelastic media. We consider two distinct mediums for topmost layer. In the first case, topmost layer has been taken as fluid saturated homogeneous porous layer, while in the second case the fluid saturated porous layer has been replaced by a transversely isotropic layer. Simple form expression for the secular equation of torsional surface wave has been worked out in both the cases by executing specific boundary conditions, which comprises Whittaker's function and its derivative, for imminent result that have been elaborated asymptotically. Some special cases have been constituted which are in excellent compliance with recorded literatures. For the sake of comparative study, numerical estimation and graphical illustration have been accomplished to identify the effects of the width ratio of the layers, Biot's gravity parameter, sandy parameter, porosity parameter and other heterogeneity parameters corresponding to the layers and half spaces, horizontal compressive and tensile initial stress on the phase velocity of torsional surface wave.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.103-108
• /
• 2015

Fluid-body interaction analysis of floating body with six degree-of-freedom motion is presented. In this study, three-dimensional incompressible Navier-Stokes equations are employed as a governing equation. The numerical method is based on a finite-volume approach on a cartesian grid together with a fractional-step method. To represent the body motion, the immersed boundary method for direct forcing is employed. In order to simulate the coupled six degree-of-freedom motion, Euler's equations based on rigid body dynamics are utilized. To represent the complex body shape, level-set based algorithm is utilized. In order to describe the free surface motion, the volume of fluid method utilizing the tangent of hyperbola for interface capturing scheme is employed. This study showed three different continuums(air, water and body) are simultaneously simulated by newly developed code. To demonstrate the applicability of the current approach, two different problems(dam-breaking with stationary obstacle and water entry) are simulated and all results are validated.

• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.327-337
• /
• 2020

The stability and dynamic performance of a flapper-nozzle servo valve depend on several factors, such as the motion of the armature component and the deformation of the spring tube. As the only connection between the armature component and the fixed end, the spring tube plays a decisive role in the dynamic response of the entire system. Aiming at predicting the vibration characteristics of the servo valves to combine them with the control algorithm, an innovative dynamic stiffness based on a distributed parameter model (DPM) is proposed that can reflect the dynamic deformation of the spring tube and a suitable discrete method is applied according to the working condition of the spring tube. With the motion equation derived by DPM, which includes the impact of inertia, damping, and stiffness force, the mathematical model of the spring tube dynamic stiffness is established. Subsequently, a suitable program for this model is confirmed that guarantees the simulation accuracy while controlling the time consumption. Ultimately, the transient response of the spring tube is also evaluated by a finite element method (FEM). The agreement between the simulation results of the two methods shows that dynamic stiffness based on DPM is suitable for predicting the transient response of the spring tube.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.14 no.3
• /
• pp.244-252
• /
• 2004

The theoretical method is developed to investigate the vibration characteristics for the partially fluid-filled continuous cylindrical shells with the intermediate supports. The intermediate supports are simulated by two types of artificial springs : the translational spring for the translation for each direction and the rotational spring for a rotation. The springs are continuously distributed along the circumferential direction. By allowing the spring stiffness to become very high compared to the stiffness of the structure, the rigid intermediate supports are approximated. In the theoretical procedure, the Love's thin shell theory is adopted to formulate the theoretical model. The frequency equation of the continuous cylindrical shell is derived by the Rayleigh-Ritz approach based on the energy method. Comparison and convergence studies are carried out to verify and establish the appropriate number of series term and the artificial spring stiffness to produce results with an acceptable order of accuracy. The effect of intermediate supports, their positions and fluid level on the natural frequencies and mode shapes are studied.