• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.1
• /
• pp.58-66
• /
• 2014

Numerical simulations for the flows containing free surface remain difficult problems because the drastic differences of physical properties of water and air, The difference of densities makes the solution instable in particular. For the stabilities of the solutions, the most typical methods to simulate free surface flows, such as Volume Of Fluid(VOF) and Level-Set(LS) methods, impose transient zones where the physical prosperities are continuously distributed. The thickness of the transient zone is the source of the numerical errors. The other side, marker-density method does not use such a transient zone. In the traditional marker-density method, however, the air velocities of free surface cells are extrapolated from the water velocity, and the pressures on the free surface are extrapolated from the air pressures for the stability of the solution. In this study, the marker-density method is modified for the decrease of such numerical errors. That is, the pressure on the free surface is determined to coincide with the pressure gradient terms of the governing equations, and the velocity of free surface cells are calculated with the governing equations. Two-dimensional steady spilling breakers behind of a submersed hydrofoil and three-dimensional spilling breaker near a wedge shaped ship model are simulated using INHAWAVE-II including the modified marker-density(MMD) method. The results are compared with the results of Fluent V6.3 including VOF method and several published research results.

• Journal of computational fluids engineering
• /
• v.18 no.4
• /
• pp.25-34
• /
• 2013

A numerical study of a laminar natural convection of the CuO-water nanofluid in a square cavity using the Buongiorno's nonhomogeneous model is presented. All the governing equations including the volume fraction equation are discretized on a cell-centered, non-uniform grid employing the finite-volume method with a primitive variable formulation. Calculations are performed over a range of Rayleigh numbers and volume fractions of the nanopartile. From the computed results, it is shown that both the homogeneous and nonhomogeneous models predict the deterioration of the natural convection heat transfer well with an increase of the volume fraction of nanoparticle at the same Rayleigh number, which was observed in the previous experimental studies. It is also shown that the differences in the computed results of the average Nusselt number at the wall between the homogeneous and nonhomogeneous models are very small, and this indicates that the slip mechanism of the Brown diffusion and thermophoresis effects are negligible in the laminar natural convection of the nanofluid. The degradation of the heat transfer with an increase of the volume fraction of the nanoparticle in the natural convection of nanofluid is due to the increase of the viscosity and the decrease of the thermal expansion coefficient and the specific heat. It is clarified in the present study that the previous controversies between the numerical and experimental studies are owing to the different definitions of the Nusselt number.

• Korean Journal of Radiology
• /
• v.24 no.6
• /
• pp.564-573
• /
• 2023

Objective: To investigate the feasibility of assessing the viscoelastic properties of the brain using magnetic resonance elastography (MRE) and a novel MRE transducer to determine the relationship between the viscoelastic properties and glymphatic function in neurologically normal individuals. Materials and Methods: This prospective study included 47 neurologically normal individuals aged 23-74 years (male-to-female ratio, 21:26). The MRE was acquired using a gravitational transducer based on a rotational eccentric mass as the driving system. The magnitude of the complex shear modulus |G^*| and the phase angle 𝛗 were measured in the centrum semiovale area. To evaluate glymphatic function, the Diffusion Tensor Image Analysis Along the Perivascular Space (DTI-ALPS) method was utilized and the ALPS index was calculated. Univariable and multivariable (variables with P < 0.2 from the univariable analysis) linear regression analyses were performed for |G^*| and 𝛗 and included sex, age, normalized white matter hyperintensity (WMH) volume, brain parenchymal volume, and ALPS index as covariates. Results: In the univariable analysis for |G^*|, age (P = 0.005), brain parenchymal volume (P = 0.152), normalized WMH volume (P = 0.011), and ALPS index (P = 0.005) were identified as candidates with P < 0.2. In the multivariable analysis, only the ALPS index was independently associated with |G^*|, showing a positive relationship (β = 0.300, P = 0.029). For 𝛗, normalized WMH volume (P = 0.128) and ALPS index (P = 0.015) were identified as candidates for multivariable analysis, and only the ALPS index was independently associated with 𝛗 (β = 0.057, P = 0.039). Conclusion: Brain MRE using a gravitational transducer is feasible in neurologically normal individuals over a wide age range. The significant correlation between the viscoelastic properties of the brain and glymphatic function suggests that a more organized or preserved microenvironment of the brain parenchyma is associated with a more unimpeded glymphatic fluid flow.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.6
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.11
• /
• pp.689-695
• /
• 2016

The performance of proton exchange membrane fuel cells (PEMFC) is strongly related to the water flow and accumulation in the gas diffusion layer (GDL) and catalyst layer. Understanding the behavior of fluid from the characteristics of the media is crucial for the improvement of the performance and design of the GDL. In this paper, a numerical method is proposed to calculate the design parameters of the GDL, i.e., permeability, tortuosity, and effective diffusivity. The fluid flow in a channel filled with randomly packed hard spheres is simulated to validate the method. The flow simulation was performed by lattice Boltzmann method with bounce back condition for the solid volume fraction in the porous media, with different values of porosities. Permeability, which affects the flow, was calculated from the average pressure drop and the velocity in the porous media. Tortuosity, calculated by the ratio the average path length of the randomly injected massless particles to the thickness of the porous media, and the resultant effective diffusivity were in good agreement with the theoretical model. The suggested method can be used to calculate the parameters of real GDL accurately without any modification.

• Journal of Advanced Marine Engineering and Technology
• /
• v.38 no.1
• /
• pp.1-7
• /
• 2014

Nanofluids, suspending nano-particles of various materials, have a good heat transfer characteristics compared with pure base fluids. For the reason, nanofluids have been considered as one of the measure to improve energy efficiency, and studied to apply on a working fluid of a high performance heat exchanger. This study tested thermal conductivities of nanofluids applying diamond nano-powder on DI water, ethylene glycol, and ethyl alcohol. Nanofluids are fabricated by matrix synthetic method, and the volume percent of diamond nano-powder contained in the base fluid are 0.1, 0.3, 0.5, and 1vol%. As a result, thermal conductivities are enhanced with applying diamond nano-power. Especially, the conductivity is highly increased up to 23% at 1vol% nanofluid applying diamond nano-powder on DI water.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1570-1579
• /
• 2000

Finite element analysis of fluid flow with moving free surface has been carried out in two and tree dimensions. The new VOF-based numerical algorithm that has been proposed by the present authors was applied to several 2-D and 3-D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools that have been newly introduced by the present authots; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2-D and 3-D cavity filling and sloshing problems, which demonstrated versatility and effectiveness of the new free surface tracking scheme as well as the overall solution procedure. The proposed numerical algorithm resolved successfully the interacting free surface with each other. The simulated results demonstrated the applicability of proposed numerical algorithm to the practical problems of large free surface motion. Also, it has been demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non-uniform grid systems and can be extended to the 3-D free surface flow problem without additional efforts.

• International Journal of Fluid Machinery and Systems
• /
• v.5 no.3
• /
• pp.109-116
• /
• 2012

Water hammer pumps can effectively use the water hammer phenomenon for water pumping. They are capable of providing an effective fluid transport method in regions without a well-developed social infrastructure. The results of experiments examining the effect of the geometric form of water hammer pumps by considering their major dimensions have been reported. However, these conventional studies have not fully evaluated pump performance in terms of pump head and flow rate, common measures of pump performance. The authors have focused on the effects on the pump performance of various geometric form factors in water hammer pumps. The previous study examined how the hydrodynamic characteristics was affected by the inner diameter ratio of the drive and lift pipes and the angle of the drive pipe, basic form factors of water hammer pumps. The previous papers also showed that the behavior of water hammer pump operation could be divided into four characteristic phases. The behavior of temporal changes in valve chamber and air chamber pressures according to the air volume in the air chamber located downstream of the lift valve was also clarified in connection with changes in water hammer pump performance. In addition, the effects on water hammer pump performance of the length of the spring attached to the drain valve and the drain pipe angle, form factors around the drain valve, were examined experimentally. This study focuses on the form of the lift valve, a major component of water hammer pumps, and examines the effects of the size of the lift valve opening area on water hammer pump performance. It also clarifies the behavior of flow in the valve chamber during water hammer pump operation.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.435-449
• /
• 2014

The effects of the gravity field and the free surface on the cavity shape and the drag are investigated through a numerical analysis for the steady supercavitating flow past a simple two-dimensional body underneath the free surface. The continuity and the RANS equations are numerically solved for an incompressible fluid using a $k-{\epsilon}$ turbulence model and a mixture fluid model has been applied for calculating the multiphase flow of air, water and vapor using the method of volume of fluid and the Schnerr-Sauer cavitation model. Numerical solutions have been obtained for the supercavitating flow about a two-dimensional $30^{\circ}$ wedge in wide range of depths of submergence and inflow velocities. The results are presented for the cavity shape, especially the length and the width, and the drag of the wedge in comparison with those of the case for the infinite fluid flow neglecting the gravity and the free surface. The influences of the gravity field and the free surface on the aforementioned quantities are discussed. The length and the width of the supercavity are reduced and the centerline of the cavity rises toward the free surface due to the effects of the gravity field and the free surface. The drag coefficient of the wedge, however, is about the same except for shallow depths of submergence. As the supercavitating wedge is approaching very close to the free surface, it is found the length and the width of a cavity are shorten even though the cavitation number is reduced. Also the present result suggests that, under the influence of the gravity field and the free surface, the length of the supercavity for a certain cavitation number varies and moreover is proportional to the inverse of the submergence depth Froude number.

• International Journal of Fluid Machinery and Systems
• /
• v.9 no.1
• /
• pp.66-74
• /
• 2016

Concerning the numerical simulation of high-speed water jet with intensive cavitation this paper presents a practical compressible mixture flow method by coupling a simplified estimation of bubble cavitation and a compressible mixture flow computation. The mean flow of two-phase mixture is calculated by URANS for compressible fluid. The intensity of cavitation in a local field is evaluated by the volume fraction of gas phase varying with the mean flow, and the effect of cavitation on the flow turbulence is considered by applying a density correction to the evaluation of eddy viscosity. High-speed submerged water jets issuing from a sheathed sharp-edge orifice nozzle are treated when the cavitation number, ${\sigma}=0.1$, and the computation result is compared with experimental data The result reveals that cavitation occurs initially at the entrance of orifice and bubble cloud develops gradually while flowing downstream along the shear layer. Developed bubble cloud breaks up and then sheds downstream periodically near the sheath exit. The pattern of cavitation cloud shedding evaluated by simulation agrees experimental one, and the possibility to capture the unsteadily shedding of cavitation clouds is demonstrated. The decay of core velocity in cavitating jet is delayed greatly compared to that in no-activation jet, and the effect of the nozzle sheath is demonstrated.