• Proceedings of the Korean Geotechical Society Conference
• /
• 2008.03a
• /
• pp.676-679
• /
• 2008

Numerical simulations of fluid injection into particulate materials were conducted to observe cavity initiation and propagation using distinct element method. After generation of initial particles and wall elements, confining stress was applied by servo-control method. The fluid scheme solves the continuity and Navior-Stokes equations numerically, then derives pressure and velocity vectors for fixed grid by considering the existence of particles within the fluid cell. Fluid was injected as 7-step into the assembly in the x-direction from the inlet located at the center of the left boundary under confining stress condition, 0.1MPa and 0.5MPa, respectively. For each simulation, movement of particles, flow rate, fluid velocity, pressure history, wall stress including cavity initiation and propagation by interaction of flulid-paricles were analyzed.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.48 no.6
• /
• pp.489-493
• /
• 1999

In this paper, an analysis of microstrip-slotline transition is performed using a 3D vector Finite Element Method(FEM). Artificial anistropic absorber technique is employed to implement an matching boundary condition in FEM. On the base of the analysis, power divider/combiner is designed. The structure of the power combiner already developed are Branch-line coupler, Rat-race coupler, Wilkinson coupler, Lange coupler, etc. Which are all planar, If the frequency goes up, the coupling efficiency of these planar couplers is decreased on account of skin loss. Especially, in millimeter-wave band, the efficiency of more than two ways combiner is radically reduced, so that application in power amplifier circuit is almost impossible, Microstrip-slotline transition structure is a power combining technique integrated into wave-guide, so that the loss is small and the efficiency is high. Theoretically, we can mount several transistors into the power-combiner. This makes it possible to develop a high power amplifier. The numerically calculated performances of the device that is, we believe, the best are compared to the experimental results in Ka-Band(26.5GHz-40GHz).

• Smart Structures and Systems
• /
• v.18 no.4
• /
• pp.787-800
• /
• 2016

The aim of the paper is to analyze nonlinear transverse vibration of an embedded piezoelectric plate reinforced with single walled carbon nanotubes (SWCNTs). The system in rested in a Pasternak foundation. The micro-electro-mechanical model is employed to calculate mechanical and electrical properties of nanocomposite. Using nonlinear strain-displacement relations and considering charge equation for coupling between electrical and mechanical fields, the motion equations are derived based on energy method and Hamilton's principle. These equations can't be solved analytically due to their nonlinear terms. Hence, differential quadrature method (DQM) is employed to solve the governing differential equations for the case when all four ends are clamped supported and free electrical boundary condition. The influences of the elastic medium, volume fraction and orientation angle of the SWCNTs reinforcement and aspect ratio are shown on frequency of structure. The results indicate that with increasing volume fraction of SWCNTs, the frequency increases. This study might be useful for the design and smart control of nano/micro devices such as MEMS and NEMS.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.453-458
• /
• 2008

Under driving condition, A vehicle experiences various kinds of loads, which brings on the buckling and fracture of suspension systems. Lower control arm (LCA), which consists of 2 bush joints and 1 ball joint connection, is the one of the most important parts in the suspension system. The bush joints absorb the impact load and reduce the vibration from the road. When analyzing the LCA behavior, it is important to understand the material properties and boundary conditions of bushing systems correctly, because of the nonlinearity characteristics of the rubber. In this paper, in order to predict the large scale deformation of the LCA more precisely, three factors are newly suggested, that is, coupling of bush stiffness between translation and rotation, bush extraction force and maximum rotation angle of ball joint. LCA stiffness is estimated by CAE and component test. Analysis and test results are almost same and the validity of considering three factors in LCA analysis is verified.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.709-714
• /
• 2000

A splitting method for the direct numerical simulation of solid-liquid mixtures is presented, where a symmetric pressure equation is newly proposed. Through numerical experiment, it is found that the newly proposed splitting method works well with a matrix-free formulation fer some bench mark problems avoiding an erroneous pressure field which appears when using the conventional pressure equation of a splitting method. When deriving a typical pressure equation of a splitting method, the motion of a solid particle has to be approximated by the 'intermediate velocity' instead of treating it as unknowns since it is necessary as a boundary condition. Therefore, the motion of a solid particle is treated in such an explicit way that a particle moves by the known form drag (pressure drag) that is calculated from the pressure equation in the previous step. From the numerical experiment, it was shown that this method gives an erroneous pressure field even for the very small time step size as a particle velocity increases. In this paper, coupling the unknowns of particle velocities in the pressure equation is proposed, where the resulting matrix is reduced to the symmetric one by applying the projector of the combined formulation. It has been tested over some bench mark problems and gives reasonable pressure fields.

• Korea-Australia Rheology Journal
• /
• v.20 no.1
• /
• pp.15-26
• /
• 2008

It has been known that there is a viscoelastic instability in the channel flow past a cylinder at high Deborah (De) number. Some of our numerical simulations and a boundary layer analysis indicated that this instability is related to the shear flow in the gap between the cylinder and the channel walls in our previous work. The critical condition for instability initiation may be related to an inflection velocity profile generated by the normal stress near the cylinder surface. At high De, the elastic normal stress coupling with the streamline curvature is responsible for the shear instability, which has been recognized by the community. In this study, an instability criterion for the flow problem is proposed based on the analysis on the pressure gradient and some supporting numerical simulations. The critical De number for various model fluids is given. It increases with the geometrical aspect ratio h/R (half channel width/cylinder radius) and depends on a viscosity ratio ${\beta}$(polymer viscosity/total viscosity) of the model. A shear thinning first normal stress coefficient will delay the instability. An excellent agreement between the predicted critical Deborah number and reported experiments is obtained.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.388-391
• /
• 2017

In this study, heat transfer analysis of variable thrust control system have been conducted by using commercial CFD code and FEM code. We Carried out computational fluid dynamics analysis to obtain the temperature and convective heat transfer coefficient of hot gas of variable thrust control system. Data are used as boundary condition for heat transfer analysis using mapping method. Temperature of O-ring for sealing was predicted

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.667-679
• /
• 2020

As the LNG carrier operates in ice covered waters, it is key to ensure the overall safety, which is related to the coupling effect of ice-breaking process and internal liquid sloshing. This paper focuses on the sloshing simulation of the ice-breaking LNG carrier, and the numerical method is proposed using Circumferential Crack Method (CCM) and Volume of Vluid (VOF) with two main key factors (velocity νx and force Fx). The ship motion analysis is carried out by CCM when the ship navigates in the ice-covered waters with a constant propulsion power. The velocity νx is gained, which is the initial excitation condition for the calculation of internal sloshing force Fx. Then, the ship motion is modified based on iterative computations under the union action of ice-breaking force and liquid sloshing load. The sloshing simulation under the LNG tank is studied with the modified ship motion. Moreover, an ice-breaking LNG ship with three-leaf type tank is used for case study. The internal LNG sloshing is simulated with three different liquid heights, including free surface shape and sloshing pressure distribution at a given moment, pressure curves at monitoring points on the bulkhead. This present method is effective to solve the sloshing simulation during ice-breaking process, which could be a good reference for the design of the polar ice-breaking LNG carrier.

• Journal of Environmental Science International
• /
• v.19 no.1
• /
• pp.27-37
• /
• 2010

A system coupled the prognostic WRF mesoscale model and CALMET diagnostic model has been employed for predicting high-resolution wind field over complex coastal area. WRF has three nested grids down to from during two days from 24 August 2007 to 26 August 2007. CALMET simulation is performed using both initial meteorological field from WRF coarsest results and surface boundary condition that is Shuttle Radar Topography Mission (SRTM) 90m topography and Environmental Geographic Information System (EGIS) 30m landuse during same periods above. Four Automatic Weather System (AWS) and a Sonic Detection And Ranging (SODAR) are used to verify modeled wind fields. Horizontal wind fields in CM_100m is not only more complex but better simulated than WRF_1km results at Backwoon and Geumho in which there are shown stagnation, blocking effects and orographically driven winds. Being increased in horizontal grid spacing, CM_100m is well matched with vertically wind profile compared SODAR. This also mentions the importance of high-resolution surface boundary conditions when horizontal grid spacing is increased to produce detailed wind fields over complex terrain features.

• Journal of the Korea Computer Graphics Society
• /
• v.20 no.4
• /
• pp.9-16
• /
• 2014

This paper proposes a novel method that couples fluid and deformation/fracture. Our method considers two interaction types: fluid-object interaction and fluid-fluid interaction. In fluid-fluid interaction, we simulate water and smoke separately and blend their velocities in the intersecting region depend on their densities. Our method separates projection process into two steps for each of water and smoke. This reduces the number of grid cells required for projection in order to optimize the number of iterations for convergence and improve stability of the simulation. In water projection step, smoke region regarded as the cells with Dirichlet boundary condition. The smoke projection step solves water region with Neumann boundary condition. To take care of fluid-object interaction, we make use of the fluid pressure to update velocities of the each of the mass points so that the object can deform or fracture. Although our method doesn't provide physically accurate results, the various examples show that our method generate appealing visuals with good performance.