• Journal of Ocean Engineering and Technology
• /
• v.23 no.1
• /
• pp.38-42
• /
• 2009

Recently, various particle based simulation techniques, which solve the Navier Stokes and continuity equations, have been developed and applied to complicated engineering problems. However, although progress is being made on their visualization or rendering techniques, these are still insufficient. In this study, to render a smooth configuration for a free surface, a rendering technique was developed that included the generation of density fields from the location information for simulated particles and the creation model for a polygonal surface. The developed rendering technique was applied to the visualization of a dynamic free surface flow interacting with a structure using a particle based simulation technique.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.395-406
• /
• 2011

In this research a numerical simulation method is developed for moving body in free surface flows using fixed staggered rectangular grid system. The non-linear free surface near the body is defined by marker-density method. The body boundary is defined by line segment connecting the points where the body surface and grid line meet. Continuity equation and Navier-Stokes equations are used as governing equations and the equations are coupled with two-step projection method. The velocities and pressures of body boundary and free surface cells are calculated with simultaneous iterative method. To treat a body movement in a fixed grid system, the volume displaced by moving body is added to the divergence of the body boundary cell. For the verification of the present numerical method. vortex shedding period of advancing cylinder is calculated and the period is compared with existing experiment results. Moreover, added mass and damping coefficients of a vertically excited box are calculated and the computed results are compared with published experiment results. Impulsive pressure and water level variation due to sloshing phenomenon are simulated and the results are compared with published experiment results. Varying the plunger shape, the waves generated by plunging type wave maker are compared with the 2nd order Stokes wave theory The plunger shape generating the wave that shows the best agreement with the theory is represented.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.05a
• /
• pp.392-399
• /
• 1997

An approximation technique was developed for the simulation of free surface flows in non-orthogonal coordinates. The main idea of this approach is to approximate VOF by the second order linear equation in the transformed domain on the assumption that the continuity of free surface would be maintained. The method was justified through a set of numerical test to examine if its original shape could be maintained when the circles are convected in uniform velocity in horizontal direction in curvilinear coordinates. Finally a simple problem was solved by applying the method to CFX4.1 general purpose CFDS code.

• Journal of the Korean Society of Visualization
• /
• v.16 no.1
• /
• pp.54-60
• /
• 2018

It is important to develop more efficient and productive casting processes for an automated high precision molten-metal casting system. Detailed analysis of molten-metal flow in the casting process by the numerical approach will help to optimize the control of a ladle. In this study, the smoothed particle hydrodynamics method was applied to analyze casting flow characteristics with different tilting angular speed and initial molten-metal level. The smoothed particle hydrodynamics technique has advantages to easily handle non-linear free surface behavior with the absence of a computational mesh. We found that tilting angular speed has relatively greater effect on the casting flowrate and that the effect of the initial molten-metal level is only minor. Further extensive study will be necessary to find an optimal condition for high efficient casting system.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.05a
• /
• pp.43-46
• /
• 2005

Micro molding technology is a promising mass production technology for polymer based microstructures. Mass production technologies such as the micro injection/compression molding, hot embossing, and micro reaction molding are already in use. In the present study, we have developed a numerical analysis system to simulate three-dimensional non-isothermal cavity filling for hot embossing, with a special emphasis on the free surface capturing. Precise free surface capturing has been successfully accomplished with the level set method, which is solved by means of the Runge-Kutta discontinuous Galerkin (RKDG) method. The RKDG method turns out to be excellent from the viewpoint of both numerical stability and accuracy of volume conservation. The Stokes equations are solved by the stabilized finite element method using the equal order tri-linear interpolation function. To prevent possible numerical oscillation in temperature Held we employ the streamline upwind Petrov-Galerkin (SUPG) method. With the developed code we investigated the detailed change of free surface shape in time during the mold filling. In the filling simulation of a simple rectangular cavity with repeating protruded parts, we find out that filling patterns are significantly influenced by the geometric characteristics such as the thickness of base plate and the aspect ratio and pitch of repeating microstructures. The numerical analysis system enables us to understand the basic flow and material deformation taking place during the cavity filling stage in microstructure fabrications.

• Earthquakes and Structures
• /
• v.4 no.5
• /
• pp.557-585
• /
• 2013

The recent huge earthquake ground motion records in Japan result in the reconsideration of seismic design forces for nuclear power stations from the view point of seismological research. In addition, the seismic design force should be defined also from the view point of structural engineering. In this paper it is shown that one of the occurrence mechanisms of such large acceleration in recent seismic records (recorded in or near massive structures and not free-field ground motions) is due to the interaction between a massive building and its surrounding soil which induces amplification of local mode in the surface soil. Furthermore on-site investigation after earthquakes in the nuclear power stations reveals some damages of soil around the building (cracks, settlement and sand boiling). The influence of plastic behavior of soil is investigated in the context of interaction between the structure and the surrounding soil. Moreover the amplification property of the surface soil is investigated from the seismic records of the Suruga-gulf earthquake in 2009 and the 2011 off the Pacific coast of Tohoku earthquake in 2011. Two methods are introduced for the analysis of the non-stationary process of ground motions. It is shown that the non-stationary Fourier spectra can detect the temporal change of frequency contents of ground motions and the displacement profile integrated from its acceleration profile is useful to evaluate the seismic behavior of the building and the surrounding soil.

• Journal of applied mathematics & informatics
• /
• v.13 no.1_2
• /
• pp.53-65
• /
• 2003

Here we present a study of small-amplitude, shallow water waves on sloping beds. The beds considered in this analysis are linear and quadratic in nature. First we start with stating the relevant governing equations and boundary conditions for the theory of water waves. Once the complete prescription of the water-wave problem is available based on some assumptions (like inviscid, irrotational flow), we normalize it by introducing a suitable set of non-dimensional variables and then we scale the variables with respect to the amplitude parameter. This helps us to characterize the various types of approximation. In the process, a summary of equations that represent different approximations of the water-wave problem is stated. All the relevant equations are presented in rectangular Cartesian coordinates. Then we derive the equations and boundary conditions for small-amplitude and shallow water waves. Two specific types of bed are considered for our calculations. One is a bed with constant slope and the other bed has a quadratic form of surface. These are solved by using separation of variables method.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2012.02a
• /
• pp.128-128
• /
• 2012

We investigated the atomic and electronic properties of graphene grown by Pd silicidation and intercalation using LEED, STM, and ARPES. Pd was deposited on the 6H-SiC(0001) surface at RT. The formation of Pd silicide gives rise to breaking of Si-C bonds of the SiC crystal, which enables to release C atoms at low temperature. The C atoms are transformed into graphene from $860^{\circ}C$ according to the LEED patterns as a function of annealing temperature. Even though the graphene spots were observed in the LEED pattern and the Fourier transformed STM images after annealing at $870^{\circ}C$, the topography images showed various superstructures so that graphene is covered with Pd silicide residue. After annealing at $950^{\circ}C$, monolayer graphene was revealed at the surface. The growth of graphene is not limited by surface obstacles such as steps and defects. In addition, we observed that six protrusions consisting of the honeycomb network of graphene has same intensity meaning non-broken AB-symmetry of graphene. The ARPES results in the vicinity of K point showed the non-doped linear ${\pi}$ band structure indicating monolayer graphene decoupled from the SiC substrate electronically. Note that the charge neutrality of graphene grown by Pd silicidation and intercalation was sustained regardless of annealing temperature in contrast with quasi-free- standing graphene induced by H and Au intercalation. Further annealing above $1,000^{\circ}C$ accelerates sublimation of the Pd silicide layer underneath graphene. This results in appearance of the $(6r3x6r3)R30^{\circ}$ structure and dissolution of the ${\pi}$ bands for quasi-free-standing graphene.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.219-235
• /
• 2014

With the development of the technology of underwater moving bodies, the need for developing the knowledge of surface effect interaction of free surface and underwater moving bodies is increased. Hence, the two-phase flow is a subject which is interesting for many researchers all around the world. In this paper, the non-linear free surface deformations which occur during the water-exit of a circular cylinder due to its buoyancy are solved using finite volume discretization based code, and using Volume of Fluid (VOF) scheme for solving two phase flow. Dynamic mesh model is used to simulate dynamic motion of the cylinder. In addition, the effect of cylinder mass in presence of an external force is studied. Moreover, the oblique exit and entry of a circular cylinder with two exit angles is simulated. At last, water-exit of a circular cylinder in six degrees of freedom is simulated in 3D using parallel processing. The simulation errors of present work (using VOF method) for maximum velocity and height of a circular cylinder are less than the corresponding errors of level set method reported by previous researchers. Oblique exit shows interesting results; formation of waves caused by exit of the cylinder, wave motion in horizontal direction and the air trapped between the waves are observable. In 3D simulation the visualization of water motion on the top surface of the cylinder and the free surface breaking on the front and back faces of the 3D cylinder at the exit phase are observed which cannot be seen in 2D simulation. Comparing the results, 3D simulation shows better agreement with experimental data, specially in the maximum height position of the cylinder.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.6
• /
• pp.23-28
• /
• 2011

In this paper, a new generation method for ship induced waves is proposed using the fully non-linear 3-D numerical model with non-reflected wave generation system (LES-WASS-3D). A ship induced wave generated by the newly proposed method is examined in comparison with that obtained by an empirical formula. It is then shown that there is a good agreement in free surface the elevation between them. As a result, it is revealed that a ship induced wave in a 3-D numerical wave field can be simulated well using LES-WASS-3D.