• Tribology and Lubricants
• /
• v.5 no.2
• /
• pp.42-47
• /
• 1989

The thermal behaviour of layerd solids, typified in practice by surface coated materials, is evaluated for the specific case of a fast moving heat source. This is intended to represent the particular instance of solids in sliding contact and the consequences of friction. The finite difference method has been utilised to establish the temperature distributions at the surface and also the sub-surface region for coating materials which are either less conductive or more conductive than the substrate to which they are attached. The effects of variation in layer thickness, and also the load, speed and friction coefficient, are evaluated.

• Ocean Systems Engineering
• /
• v.8 no.2
• /
• pp.167-182
• /
• 2018

RTI (Rayleigh-Taylor instability) is investigated by a multi-liquid MPS (Moving Particle Semi-implicit) method for both viscous and inviscid flows for various density differences, initial-disturbance amplitudes, viscosities, and surface tensions. The MPS simulation can be continued up to the late stage of high nonlinearity with complicated patterns and its initial developments agree well with the linear theoretical results. According to the relevant linear theory, the difference between inviscid and viscous fluids is the rising velocity at which upward-mushroom-like RTI flow with vortex formation is generated. However, with the developed MPS program, significant differences in both growing patters and developing speeds are observed. Also, more dispersion can be observed in the inviscid case. With larger Atwood (AT) number, stronger RTI flows are developed earlier, as expected, with higher potential-energy differences. With larger initial disturbances, quite different patterns of RTI-development are observed compared to the small-initial-disturbance case. If AT number is small, the surface tension tends to delay and suppress the RTI development when it is sufficiently large. Interestingly, at high AT number, the RTI-suppressions by increased surface tension become less effective.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.382-385
• /
• 2006

A high-order spectral/boundary-element method is newly adapted as an efficient numerical tool. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved fly using the high-order spectral method and body potential is solved fly using the high-order boundary element method. Through the combination of these two methods, the wave-making problems fly a submerged sphere moving with the large amplitude oscillation are solved in time-domain. With the example calculations, nonlinear effects on free-surface profiles and hydrodynamic forces are shown and discussed.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.21-32
• /
• 1992

Numerical solutions are presented for solving the free surface flow created by a three-dimensional body moving beneath the free surface with constant velocity at an angle of attack. The solution is obtained using a panel method based on the perturbation potential, which employs Havelock sources and normal dipoles distributed on the body surface and Havelock normal dipoles in the wake downstream of the trailing edge. A pressure Kutta condition with an iterative solution procedure is implemented to satisfy equal pressure condition on the upper and lower surfaces at the trailing edge. Numerical calculation examples in the present paper include an ellipsoid at zero angle of attack, a rectangular planform wing at a small angle of attack in the limit of zero Froude number and then free surface flows and hydrodynamic forces acting on the submerged spheroid and parabolic strut are calculated. Discussions are made about the validity of the present method.

• Journal of the Korean Society for Precision Engineering
• /
• v.22 no.5 s.170
• /
• pp.205-213
• /
• 2005

This paper addresses a method for numerical simulation in the pressing process of hot glass. Updated Lagrangian finite element formulations are employed for the flow and energy equations to accommodate moving meshes. The model is assumed axi-symmetric and creep flow is assumed due to the high viscosity. Commercial software ANSYS is used to solve the coupled flow and energy equations. Moving contact points as well as free surface during the pressing are effectively calculated and updated by utilizing API functions of CAD software Unigraphics. The mesh distortion problem near the wall is overcome by automatic remeshing, and the temperatures of the new mesh are conveniently interpolated by using a unique function of ANSYS. The developed model is applied to the pressing process of TV glasses. In conclusion, the presented method shows that the pressing process accompanying moving boundary can be simulated by effectively combining general purpose software without resorting to special dedicated codes.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.6
• /
• pp.388-396
• /
• 2017

A Lagrangian approach based computational fluid dynamics (CFD) was used to simulate large and/or sharp deformations and fragmentations of interfaces, including free surfaces, through tracing each particle with physical quantities. According to the concept of the particle-based CFD method, it is possible to apply it to both fluid particles and solid particles such as sand, gravel, and rock. However, the presence of more than two different phases in the same domain can make it complicated to calculate the interaction between different phases. In order to solve multiphase problems, particle interaction models for multiphase problems, including surface tension, buoyancy-correction, and interface boundary condition models, were newly adopted into the moving particle semi-implicit (MPS) method. The newly developed MPS method was used to simulate a typical validation problem involving dam breaking. Because the soil and other particles, excluding the water, may have different viscosities, various viscosity coefficients were applied in the simulations for validation. The newly developed and validated MPS method was used to simulate the mobile beds induced by broken dam flows. The effects of the viscosity on soil particles were also investigated.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.720-731
• /
• 2002

Flow with moving free surfaces is analyzed with an the Eulerian coordinate system. This study proposes a semi-implicit filling algorithm using VOF in which the PLIC (Piecewise Linear Interface Calculation) -type interface reconstruction method and the donor-acceptor-type front advancing scheme are adopted. Also, a new scheme using extrapolation of the stream function is proposed to find the velocity of the node that newly enters the computational domain. The effect of wall boundary conditions on the flow field and temperature field is examined by numerically solving a two-dimensional casting process.

• Journal of Ocean Engineering and Technology
• /
• v.28 no.2
• /
• pp.93-101
• /
• 2014

In general, particle simulation methods such as the MPS(Moving Particle Simulation) or SPH(Smoothed Particle Hydrodynamics) methods have some serious drawbacks for pressure solutions. The pressure field shows spurious high fluctuations both temporally and spatially. It is well known that pressure fluctuation primarily occurs because of the numerical approximation of the partial differential operators. The MPS and SPH methods employ a pre-defined kernel function in the approximation of the gradient and Laplacian operators. Because this kernel function is constructed artificially, an accurate solution cannot be guaranteed, especially when the distribution of particles is irregular. In this paper, we propose a particle simulation method based on the moving least-square technique for solving the partial differential operators using a Taylor-series expansion. The developed method was applied to the hydro-static pressure and dam-broken problems to validate it.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.

• Journal of Mechanical Science and Technology
• /
• v.20 no.7
• /
• pp.1043-1050
• /
• 2006

The aerodynamic interaction between a wing and a rail is investigated using a boundary-element method. The source and doublet singularities are distributed on the wing and its guide-way rail surface. The unknown strengths of the singularities are determined by inverting the aerodynamic influence coefficient matrices. Present method is validated by comparing computed results with the other numerical data. Rail width and rail height affect the aerodynamic characteristics of the wing only if the rail is narrower than the wing span. Although the present results are limited to the inviscid, irrotational flows, it is believed that the present method can be applied to the conceptual design of the high speed ground transporters moving over the rail.