• 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.

• Advances in biomechanics and applications
• /
• v.1 no.2
• /
• pp.127-141
• /
• 2014

A preliminary study of a new pulsatile pump that will work to a frequency greater than 1 Hz, is presented. The fluid-structure interaction between a Newtonian blood flow and a piston drive that moves with periodic speed is simulated. The mechanism is of double effect and has four valves, two at the input flow and two at the output flow; the valves are simulated with specified velocity of closing and reopening. The simulation is made with finite elements software named COMSOL Multiphysics 3.3 to resolve the flow in a preliminary planar configuration. The geometry is 2D to determine areas of high speeds and high shear stresses that can cause hemolysis and platelet aggregation. The opening and closing valves are modelled by solid structure interacting with flow, the rhythmic opening and closing are synchronized with the piston harmonic movement. The boundary conditions at the input and output areas are only normal traction with reference pressure. On the other hand, the fluid structure interactions are manifested due to the non-slip boundary conditions over the piston moving surfaces, moving valve contours and fix pump walls. The non-physiologic frequency pulsatile pump, from the viewpoint of fluid flow analysis, is predicted feasible and with characteristic of low hemolysis and low thrombogenesis, because the stress tension and resident time are smaller than the limit and the vortices are destroyed for the periodic flow.

• The Korean Journal of Rheology
• /
• v.11 no.1
• /
• pp.57-65
• /
• 1999

In order to measure yield stress of PIM feedstocks simply and effectively, a yield stress measuring technique was developed by a vane method. The vane method had an advantage that there was no wall-slip, while it had a drawback that it could not measure viscosity change at various shear rates. A Newtonian fluid was tested for the appropriateness of the measuring technique. The end effect of a vane was checked to produce an acceptable error. The torque peak has been considered to be developed at yielding of non-Newtonian fluids with yield stress. However, it was influenced very much by control system of the instrument so that the torque value at the stable region was taken to calculate yield stress. Torque at zero rotational speed was obtained by extrapolating the torque values at various speeds to remove the effect of the rotational drag. As general verification, yield stress of feedstocks made of Tungsten carbide powder with wax-based binder was measured at different temperatures and various powder concentrations.

• Textile Coloration and Finishing
• /
• v.22 no.3
• /
• pp.187-193
• /
• 2010

The effects of molecular weight (MW) and concentration on the rheological properties of poly(vinyl alcohol) (PVA) solutions in dimethyl sulfoxide (DMSO) were investigated at $30^{\circ}C$. Ubbelohde viscometer and rotational rheometer were employed for dilute and concentrated regime, respectively. In the dilute regime, the Mark-Houwink exponent ($\alpha$) of the solutions determined from three different MWs proved 0.73. The critical concentration (C*), in which the entanglement and overlap of polymer molecules began to take place, decreased with increasing the MW of PVA. Huggins constant ($K_H$) values ranged from 0.33 to 0.45 over the MW examined. In the log-log plot of $\eta_{sp}$ versus [$\eta$]C, the PVA with higher degree of polymerization (DP) gave a greater slope exhibiting the inflection point in the vicinity of C*. In the dynamic viscosity ($\eta'$) curve, the PVA solutions of DP 1700 presented Newtonian fluid behavior over most of the frequency range examined. However, the lower Newtonian flow region reduced with increasing the DP. As the PVA concentration increased, $\eta'$ was increased and the onset shear rate for pseudoplasticity was decreased. In the Cole-Cole plot, PVA solutions showed almost a single master curve in a slope of ca. 1.65 regardless of the DP. However, the increase of the concentration from 8 to 12 wt% for PVA solutions of DP 5000 decreased the slope from 1.73 to 1.57. In the tan $\delta$ curve, the onset frequency for sol-gel transition was shifted from 154 to 92 rad/s with increasing the DP from 3300 to 5000 and from 192 to 46 rad/s with increasing the concentration from 8 to 12 wt%. In addition, longer relaxation time ($\lambda$) was observed with increasing the DP and concentration.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1465-1477
• /
• 1993

Powder Injection Molding(PM) is an advanced and complicated technology for manufacturing ceramic or metal products making use of a conventional injection molding process, which is generally used for plastic products. Among many technologies involved in the successful PIM, injection molding process is one of the key steps to form a desired shape out of powder/binder mixtures. Thus, it is of great importance to have a numerical tool to predict the powder injection molding filling process. In this regard, a finite element analysis system has been developed for numerical simulations of filling process of powder injection molding. Powder/polymer mixtures during the filling pro cess of injection molding can be rheologically characterized as Non-Newtonian fluids with a so called yield phenomena and have a peculiar feature of apparent slip phenomena on the wall boundaries surrounding mold cavity. Therefore, in the present study, a physical modeling of the filling process of powder/polymer mixtures was developed to take into account both the yield stress and slip phenomena and a finite element formulation was developed accordingly. The numerical analysis scheme for filling simulation is accomplished by combining a finite element method with control volume technique to simulate the movement of flow front and a finite difference method to calculate the temperature distribution. The present study presents the modeling, numerical scheme and some numerical analysis results showing the effect of the yield stress and slip phenomena.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.2 no.1
• /
• pp.37-48
• /
• 1990

The drag and heat transfer reduction phenomena and degradation effects of drag reducing polymer solutions which are known as the viscoelastic fluids are investigated experimentally for the turbulent circular tube flows. Two stainless steel tubes are used for the experimental flow loops. Aqueous solutions of Polyacrylamide Separan AP-273 with concentrations from 300 to 1000 wppm are used as working fluids. Flow loops are set up to measure the friction factors and heat transfer coefficients of test tubes in the once-through system and the recirculating flow system. Test tubes are heated by power supply directly to apply constant heat flux boundary conditions on the wall. Capillary tube viscometer and falling ball viscometer are used to measure the viscous characteristics of fluids and the characteristic relaxation time of a fluid is determined by the Powell-Eyring model. The order of magnidude of the thermal entrance length of a drag reducing polymer solution is close to the order of magnitude of the laminar entrance length of Newtonian fluids. Dimensionless heat transfer coefficients of the viscoelastic non-Newtonian fluids may be represented as a function of flow behavior index n and newly defined viscoelastic Graetz number. As degradation continues viscosity and the characteristic relaxation time of the testing fluids decrease and heat transfer coefficients increase. The characteristic relaxation time is used to define the Weissenberg number and variations of friction factors and heat transfer coefficients due to degradation are presented in terms of the Weissenberg number.

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

In this study, the rheological properties of nitromethane gel propellants on nano/micron sized gelling agent are investigated. Silicon dioxide is used as the gellant with 5 wt%, 6.5 wt% and 8 wt% concentration, respectively, where the measurements are conducted under steady-state shear flow conditions using a rotational rheometer. The nitromethane/silicon dioxide gel showed non-Newtonian flow behavior for the entire experimental shear rate ranges. The gel fuels with nano-sized gellant had a slightly higher viscosity than the gel fuels with micron-sized one for low shear rate range. Additionally, it was found that Herschel-Bulkley model can hardly describe the rheological behavior of nitromethane gel propellant, but the NM model(by Teipel and Forter-Barth) is better suited to explain the rheological behavior of nitromethane gel propellant.

• Journal of the Korean earth science society
• /
• v.45 no.4
• /
• pp.292-303
• /
• 2024

We study quasi-spherical, supersonic accretion flows around black holes using high-accuracy numerical simulations. We describe a code, the Lagrangian Total Variation Diminishing (TVD), and a remap routine to address a specific issue in the Advection Dominated Accretion Flow (ADAF) that is, appropriately handling the angular momentum even near the inner boundary. The Lagrangian TVD code is based on an explicit finite difference scheme on mass-volume grids to track fluid particles with time. The consequences are remapped on fixed grids using the explicit Eulerian finite-difference algorithm with a third-order accuracy. Test results show that one can successfully handle flows and resolve shocks within two to three computational cells. Especially, the calculation of a hydrodynamical accretion disk without viscosity around a black hole shows that one can conserve nearly 100% of specific a ngular momentum in one-and two-dimensional cylindrical coordinates. Thus, we apply this code to obtain a numerically similar ADAF solution. We perform simulations, including viscosity terms in one-dimensional spherical geometry on the non-uniform grids, to obtain greater quantitative consequences and to save computational time. The error of specific angular momentum in Newtonian potential is less than 1% between r~10r_s and r~10⁴ r_s, where r_s is sink size. As Narayan et al. (1997) suggested, the ADAFs in pseudo-Newtonian potential become supersonic flows near the black hole, and the sonic point is r_sonic~5.3r_g for flow with α =0.3 and γ=1 .5. Such simulations indicate that even the ADAF with γ=5/3 is differentially rotating, as Ogilvie (1999) indicated. Hence, we conclude that the Lagrangian TVD and remap code treat the role of viscosity more precisely than the other scheme, even near the inner boundary in a rotating accretion flow around a nonrotating black hole.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2951-2954
• /
• 2007

Tumor hemodynamics in vascular state is numerically simulated using pressure node solution. The tumor angiogenesis pattern in our previous study is used for the geometry of vessel networks. For tumor angiogenesis, the equation that governed angiogenesis comprises a tumor angiogenesis factor (TAF) conservation equation in time and space, which is solved numerically using the Galerkin finite element method. A stochastic process model is used to simulate vessel formation and vessel. In this study, we use a two-dimensional model with planar vessel structure. Hemodynamics in vessel is assumed as incompressible steady flow with Newtonian fluid properties. In parent vessel, arterial pressure is assigned as a boundary condition whereas a constant terminal pressure is specified in tumor inside. Kirchhoff's law is applied to each pressure node to simulate the pressure distribution in vessel networks. Transient pressure distribution along with angiogenesis pattern is presented to investigate the effect of tumor growth in tumor hemodynamics.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2002.10b
• /
• pp.91-92
• /
• 2002

Most studies of elsatohydrodynamic lubrication are oriented only to the measurement of film thickness itself with optical interferometer. In order to exactly investigate the characteristics of a certain lubricant under the condition of additives. especially for traction performance. it is also important to get the information of traction force as well. In this work. we developed the device for measuring friction force of EHL contact condition, which can trace the film thickness over the contact area with optical interferometer. To verify the validity of the measuring system, the friction force and film thickness under EHL condition are measured with the variation of additive ratios of viscosity Index improvers.