• Journal of the Korea Computer Graphics Society
• /
• v.10 no.4
• /
• pp.13-21
• /
• 2004

This paper presents a physically based technique for simulating complex multiphase fluids. This work is motivated by the "stable fluids" method developed by Stam to handle gaseous fluids. We extend this technique to water, which calls for the development of methods for modeling multiphase fluids and suppressing dissipation. We construct a multiphase fluid formulation by combining the Navier-Stokes equations with the level set method. By adopting constrained interpolation profile (CIP)-based advection, we reduce the numerical dissipation and diffusion significantly. We further reduce the dissipation by converting potential1y dissipative cel1s into droplets or bubbles that undergo Lagrangian motion. Due to the multiphase formulation, the proposed method properly simulates the interaction of water with surrounding air, instead of simulating water in a void space. Moreover, the introduction of the non-dissipative technique means that, in contrast to previous methods, the simulated water does not unnecessarily lose mass and its motion is not damped to an unphysical extent. Experiments showed that the proposed method is stable and runs fast. It is demonstrated that two-dimensional simulation runs in real-time.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.84-87
• /
• 2002

When computing the flow around complex three dimensional configurations, the generation of grid is the sunt time consuming part of any calculation. The object of this study is to develop the grid duster techniques capable of resolving complex flows with shock waves, expansion waves, shear layers, and cursive shapes, The Dot insert method of Non-Uniform Rational B-Splines is described as a id control method.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.152-156
• /
• 2001

Three dimensional orthogonal grid generation is able to control effectively the grid spacing near the boundaries, but there are some difficulty to meshing complex geometry. The mesh complex geometry by orthogonal grid generation method must divide block of geometry It is required a careful skill, and long time. Its also difficulty to make unstructured mesh on complex geometry. Particularly, three dimensional geometry must have more time and effort. Recently, there have been growing interests in mesh generation of complex grometry, aslike an automobile headlamp, the heart. The method of easily meshing complex geometry is resarched to solve them. We suggest octree grid into one among these methods. As octrce grid is automaticaly adapted at the boundaries by determine the level operations to control the grid spacing near the boundaries are unnecessary. In this paper we showed throe dimensional mesh generation, and heat-flow analysis on the octree mesh.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.19-22
• /
• 1998

Automatic mesh generation procedures applied to industrial now problems lead to complex mesh topologies where usually no special considerations to mesh resolution are taken. In the present study a fast and flexible solution algorithm in combination with generalized higher order discretization schemes is presented and its application to intake port calculation is demonstrated.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.41-42
• /
• 2006

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.493-495
• /
• 2010

A three-dimensional flow simulation is performed to investigate the flow field in the ski resort on complex terrain. The present paper aims to study the wind effects of mountainous terrain on the gondola safety. Strong wind happens in the ski resort on the mountain by complex terrain and it causes the dangerous accident of gondola. A digital map around the ski resort area is used to model the actual complex terrain for a 3-D analysis domain. Wind direction and speed to be used as a boundary condition are taken from local meteorological reports. The numerical results show details of the velocity distribution around a ski resort. From the results, we can suggest the modification of the installation of gondola for the safety due to strong wind.

• Korea-Australia Rheology Journal
• /
• v.17 no.4
• /
• pp.207-215
• /
• 2005

We present a finite difference solution for electrokinetic flow in rectangular microchannels encompassing Navier's fluid slip phenomena. The externally applied body force originated from between the nonlinear Poisson-Boltzmann field around the channel wall and the flow-induced electric field is employed in the equation of motion. The basic principle of net current conservation is applied in the ion transport. The effects of the slip length and the long-range repulsion upon the velocity profile are examined in conjunction with the friction factor. It is evident that the fluid slip counteracts the effect by the electric double layer and induces a larger flow rate. Particle streak imaging by fluorescent microscope and the data processing method developed ourselves are applied to straight channel designed to allow for flow visualization of dilute latex colloids underlying the condition of simple fluid. The reliability of the velocity profile determined by the flow imaging is justified by comparing with the finite difference solution. We recognized the behavior of fluid slip in velocity profiles at the hydrophobic surface of polydimethylsiloxane wall, from which the slip length was evaluated for different conditions.

• Korea-Australia Rheology Journal
• /
• v.15 no.2
• /
• pp.83-90
• /
• 2003

In cases of the microfluidic channel, the electrokinetic influence on the transport behavior can be found. The externally applied body force originated from the electrostatic interaction between the nonlinear Poisson-Boltzmann field and the flow-induced electrical field is applied in the equation of motion. The electrostatic potential profile is computed a priori by applying the finite difference scheme, and an analytical solution to the Navier-Stokes equation of motion for slit-like microchannel is obtained via the Green's function. An explicit analytical expression for the induced electrokinetic potential is derived as functions of relevant physicochemical parameters. The effects of the electric double layer, the zeta potential of the solid surface, and the charge condition of the channel wall on the velocity profile as well as the electroviscous behavior are examined. With increases in either electric double layer or zeta potential, the average fluid velocity in the channel of same charge is entirely reduced, whereas the electroviscous effect becomes stronger. We observed an opposite behavior in the channel of opposite charge, where the attractive electrostatic interactions are presented.

• Korea-Australia Rheology Journal
• /
• v.14 no.4
• /
• pp.189-201
• /
• 2002

The interaction between hydrophobically modified hydroxyethyl cellulose (hmHEC), containing approximately 1 wt% side-alkyl chains of $C_{16}$, and an anionic sodium dodecyl sulphate (SDS) surfactant was investigated. For a semi-dilute solution of 0.5 wt% hmHEC, the previously observed behaviour of a maximum in solution viscosity at intermediate SDS concentrations, followed by a drop at higher SDS concentrations, until above the cmc of surfactant when the solution resembles that of the unsubstituted polymer, was confirmed. Additionally, a two-phase region containing a hydrogel phase and a water-like supernatant was found at low SDS concentrations up to 0.2 wt%, a concentration which is akin to the critical association concentration, cac, of SDS in the presence of hmHEC. Above this concentration, SDS molecules bind strongly to form mixed micellar aggregates with the polymer alkyl side-chains, thus strengthening the network junctions, resulting in the observed increase in viscosity and elastic modulus of the solution. The shear behaviour of this polymer-surfactant complex during steady and step stress experiments was examined In great detail. Between SDS concentrations of 0.2 and 0.25 wt%, the shear viscosity of the hmHEC-polymer complex network undergoes shear-induced thickening, followed by a two-stage shear-induced fracture or break-up of the network. The thickening is thought to be due to structural rearrangement, causing the network of flexible polymers to expand, enabling some polymer hydrophobic groups to be converted from intra- to inter-chain associations. At higher applied stress, a partial local break-up of the network occurs, while at even higher stress, above the critical or network yield stress, a complete fracture of the network into small microgel-like units, Is believed to occur. This second network rupture is progressive with time of shear and no steady state in viscosity was observed even after 300 s. The structure which was reformed after the cessation of shear is found to be significantly different from the original state.

• Analytical Science and Technology
• /
• v.7 no.3
• /
• pp.413-419
• /
• 1994

The determination of the neuromuscular blocking agents vecuronium bromide(VeBr) in biological fluids has been investigated. The method depends on the formation of insoluble red complex between vecuronium bromide and rose bengal in aqueous layer. The amount of vecuronium bromide was calculated from that of extracted rose bengal which was determined by spectrofluorimetry or HPLC/fluorescence detection method. It was possible to analyze VeBr in the range of $2{\sim}32{\mu}g/ml$(r=0.998 for water soln., 0.999 for urine, 0.996 for plasma). This method was applied to the analysis of VeBr in biological fluids, urine and plasma.