• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.7
• /
• pp.455-462
• /
• 2017

The magnetic interaction between elliptic Janus magnetic particles are investigated using a direct simulation method. Each particle is a one-to-one mixture of paramagnetic and nonmagnetic materials. The fluid is assumed to be incompressible Newtonian and nonmagnetic. A uniform magnetic field is applied externally in a horizontal direction. A finite-element-based fictitious domain method is employed to solve the magnetic particulate flow in the creeping flow regime. In the magnetic problem, the magnetic field in the entire domain, including the particles and the fluid, is obtained by solving the governing equation for the magnetic potential. Then, the magnetic forces acting on the particles are calculated via a Maxwell stress tensor formulation. In a single particle problem, it is found that the orientation angle at equilibrium is affected by the aspect ratio of the particle. As for the two-particle interaction, the dynamics and the final conformation of the particles are significantly influenced by the aspect ratio, the orientation, and the spatial positions of the particles. For the given positions of the particles, the fluid flow is also influenced by the orientation of each particle. The self-assembly structure of the particles is not a fixed one, but it varies with the above-mentioned factors.

• The Journal of the Acoustical Society of Korea
• /
• v.37 no.4
• /
• pp.163-173
• /
• 2018

This paper performed the prediction of the acoustic loads applied to the surface of the flight vehicle during flight. Acoustic loads during flight arise from the pressure fluctuations on the surface of body. The conventional method of predicting the acoustic loads in flight uses semi-empirical method derived from theoretical and experimental results. However, there is a limit in obtaining the flow characteristics and the boundary layer parameters of the flight vehicle which are used as the input values of the empirical equation through experiments. Therefore, in this paper, we use the hybrid method which combines the results of CFD (Computational Fluid Dynamics) with semi-empirical methods to predict the acoustic loads acting on flight vehicle during flight. For the flight vehicle with cone-cylinder-flare shape, acoustic loads were estimated for the subsonic, transonic, supersonic, and Max-q (Maximum dynamic pressure) condition flight. For the hybrid method, two kind of boundary layer edge estimation methods based on CFD results are compared and the acoustic loads prediction results were compared according to empirical equations presented by various researchers.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.5
• /
• pp.425-433
• /
• 2015

In this paper, we have numerically investigated two-dimensional dielectrophoretic (DEP) motions of a single particle suspended freely in a viscous fluid, interacting with a nearby nonconducting planar wall, under an externally applied uniform direct-current electric field. Particularly, we solve the Maxwell equation with a large sharp jump in the electric conductivity at the particle-fluid interface and then integrate the Maxwell stress tensor to compute the DEP force on the particle. Results show that, under an electric field parallel to the wall, one particle is always repelled to move far away from the wall and the motion depends strongly on the particle-wall spacing and the particle conductivity. The motion strength vanishes when the particle is as conductive as the fluid and increases as the conductivity deviates further from that of the fluid.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.3
• /
• pp.9-18
• /
• 1985

The numerical calculation for solving boundary-value problem related to potential flows with a free surface is carried out by application of the localized finite element method. Only forced motion of 2-D body in infinitely deep fluid is considered, although this schemes is equally applicable to any first order time-harmonic problems of similar nature. The infinite domain of the fluid is separated into the inner flow field and the outer flow field with common inter-surface boundary. The finite element method is applied to obtain the solution in the inner flow field and the Green functions are utilized to represent the solution in the outer flow field. At the inter-surface boundary, the continuity of the value of potential and the normal derivative of the potential(i.e. matching condition) is conserved. The present method has better computational efficiency than the previous LFEM and the integral equation method of Frank. This enhanced computational efficiency is presumably due to the fact that the present method gives a symmetric coefficient matrix and requires less computational time in calculating the influence coefficient matrix of Green function than the integral equation method. And the irregular frequency desen't exist because the uniqueness of the solution is assured by the such that the exact free surface condition is satisfied on the boundary of the localized finite element region(i.e. inner region). As an example of the above method, the hydrodynamic forces for the circular cylinder and the rectangular cylinders are calculated. In the computed results, the small number of singularity distribution segments($3{\sim}6$) give good result relative to Ursell's and Vugts'.

• Applied Chemistry for Engineering
• /
• v.10 no.2
• /
• pp.324-329
• /
• 1999

In this work, high pressure binary phase equilibria data of carbon dioxide and isopropyl alcohol were obtained by experiment. A static type experimental apparatus was made to measure temperature, pressure and phase equilibria composition. The experimental apparatus was tested by comparing the measured phase equilibria data of the carbon dioxide-isopropyl alcohol system at $80^{\circ}C$ with those of Rodosz. The binary phase behavior data of carbon dioxide-isopropyl alcohol system were measured in range of 41 to 133 bar and at temperatures of 40, 60, 80, 100 and $120^{\circ}C$. The solubility of isopropyl alcohol increases as the temperatures increases at constant pressure. Also, these carbon dioxide-alcohol solute system have critical-mixture curves that exhibit maxima in pressure at temperatures between the critical temperatures of carbon dioxide and isopropyl alcohol. The experimental data obtained in this study were modeled using the statistical associating fluid theory(SAFT) equation of state. A good fit of the data was obtained with SAFT using two adjustable parameters for the carbon dioxide-isopropyl alcohol system.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1487-1503
• /
• 2019

The methods and performance of a pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT employing the drift-flux model are presented. This code aims at providing an accurate yet fast core thermal-hydraulics solution capability to high-fidelity multiphysics core analysis systems targeting massively parallel computing platforms. The four equation drift-flux model is adopted for two-phase calculations, and numerical solutions are obtained by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models involving turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analyses. ESCOT is parallelized by a domain decomposition scheme that involves both radial and axial decomposition to enable highly parallelized execution. The ESCOT solutions are validated through the applications to various experiments which include CNEN $4{\times}4$, Weiss et al. two assemblies, PNNL $2{\times}6$, RPI $2{\times}2$ air-water, and PSBT covering single/two-phase and unheated/heated conditions. The parameters of interest for validation include various flow characteristics such as turbulent mixing, spacer grid pressure drop, cross-flow, reverse flow, buoyancy effect, void drift, and bubble generation. For all the validation tests, ESCOT shows good agreements with measured data in the extent comparable to those of other subchannel-scale codes: COBRA-TF, MATRA and/or CUPID. The execution performance is examined with a mini-sized whole core consisting of 89 fuel assemblies and for an OPR1000 core. It turns out that it is about 1.5 times faster than a subchannel code based on the two-fluid three field model and the axial domain decomposition scheme works as well as the radial one yielding a steady-state solution for the OPR1000 core within 30 s with 104 processors.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.377-387
• /
• 2010

In this study, the turbulent free surface around KVLCC1 employed in the circular motion test simulation is numerically calculated using a commercial CFD(Computational Fluid Dynamics) code, FLUENT. Also, hydrodynamic forces and yaw moments around a ship model are calculated during the steady turning. Numerical simulations of the turbulent flows with free surface around KVLCC1 have been carried out by use of RANS equation based on calculation of hydrodynamic forces and yaw moments exerted upon the ship hull. Wave elevation is simulated by using the VOF method. VOF method is known as one of the most effective numerical techniques handling two-fluid domains of different density simultaneously. Boundary layer thickness and wake field are changed various yaw velocities of ship model during the steady turning. The calculated hydrodynamic forces are compared with those obtained by model tests.

• The KSFM Journal of Fluid Machinery
• /
• v.4 no.1 s.10
• /
• pp.22-29
• /
• 2001

Centrifugal fans are widely used in industrial practices but the noise generated by these machines causes one of the most serious problems. In general, the centrifugal fan noise is often dominated by tones at BPF(blade passage frequency) and its higher harmonics. This is a consequence of the strong interaction between the flow discharged from the impeller and the cutoff in the easing. However, only a few researches have been carried out on predicting the noise because of the difficulty in obtaining detailed information about the flow field and casing effects on noise radiation. The objective of this study is to develop a prediction method for the unsteady flow field and the acoustic pressure field of a centrifugal fan, and to calculate the effects of small vanes that are attached in original impeller - Splitter impeller. We assume that the impeller rotates with a constant angular velocity and the flow field around the impeller is incompressible and inviscid. So, a discrete vortex method (DVM) is used to model the centrifugal fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. The splitter impeller changes the acoustic characteristics as well as performance. Two-splitter type impeller and splitter impeller which splitter locates in jet region are good for acoustic characteristics.

• International Journal of Fluid Machinery and Systems
• /
• v.8 no.3
• /
• pp.169-182
• /
• 2015

When simulating the dynamic behaviour of a hydro power plant, it is essential to have a good representation of the turbine behaviour. The pressure transients in the system occurs because the flow changes, which the turbine defines. The flow through the turbine is a function of the pressure, the speed of rotation and the wicket gate opening and is, most often described in a performance diagram or Hill diagram. In the Hill diagram, the efficiency is drawn like contour lines, hence the name. A turbines Hill diagram is obtained by performance tests on scaled model in a laboratory. However, system dynamic simulations have to be performed in the early stage of a project, before the turbine manufacturer has been chosen and the Hill diagram is known. Therefore one have to rely on diagrams for a turbine with similar speed number. The Hill diagram is drawn through measured points, so for using the diagram in a simulation program, one have to iterate in the diagram based on curve fitting of the measured points. This paper describes an alternative method. By means of the Euler turbine equation, it is possible to set up two differential equations which represents the turbine performance with good enough accuracy for the dynamic simulations. The only input is the turbine's main geometry, the runner blade in- and outlet angle and the guide vane angle at best efficiency point of operation (BEP). In the paper, simulated turbine characteristics for a high head Francis turbine, and for a reversible pump turbine are compared with laboratory measured characteristics.

• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1583-1602
• /
• 2011

A special system of partial differential equations (PDEs) occur in a natural way while studying a class of irrotational inviscid fluid flow problems involving infinite channels. Certain aspects of solutions of such PDEs are analyzed in the context of flow problems involving multiple layers of fluids of different constant densities in a channel associated with arbitrary bottom topography. The whole analysis is divided into two parts-part A and part B. In part A the linearized theory is employed along with the standard Fourier analysis to understand such flow problems and physical quantities of interest are derived analytically. In part B, the same set of problems handled in part A are examined in the light of a weakly non-linear theory involving perturbation in terms of a small parameter and it is shown that the original problems can be cast into KdV type of nonlinear PDEs involving the bottom topography occurring in one of the coefficients of these equations. Special cases of bottom topography are worked out in detail and expressions for quantities of physical importance are derived.