• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.60-68
• /
• 1991

The updated Lagrangian Finite Element Method is introduced to analyse rigid body-fluid impact problem which is characterized by incompressible Navier-Stokes equations and impact-contact conditions between free surface and rigid body. For the convenience of numerical computation, velocity fields are splinted into vicous and pressure parts, and then the governing equations and boundary conditions are decomposed in accordance with the decomposition. However, Viscous stresses acting an the solid boundaries are neglected on the assumption that very small velocity gradients may occur during extremely small time interval of the impact. Four coded quadrilateral elements are used to discretize the space domain and the fully explicit time-marching algorithm is employed with a reasonably small time step. At the beginning of each time step, contact velocity of the rigid body is computed from the momentum balance between the body and the fluid. The velocity field is then computed to satisfy the discretized equations of motions and incompressibility and contact constraints as well as an exact free surface boundary condition. At the end of each time step, the fluid domain is updated from the velocity field. In the present time stepping numerical analysis, behaviour of the free surface near the body can be observed without any difficulty which is very important in the water impact problem. The applicability of the algorithm is illustrated by a wedge type falling body problem. The numerical solutions for time-varying pressure distributions and impact loadings acting ion the surface are obtained.

• Advances in nano research
• /
• v.9 no.3
• /
• pp.133-145
• /
• 2020

In this study, nonlinear vibrations and dynamic instabilities of a smart embedded micro shell conveying varied fluid flow and subjected to the combined electro-thermo-mechanical loadings are investigated. With the aim of designing new hydraulic sensors and actuators, the piezoelectric materials are employed for the body and the effects of applying electric field on the stability of the system as well as the induced voltage due to the dynamic behavior of the system are studied. The nonlocal piezoelasticity theory and the nonlinear cylindrical shell model in conjunction with the energy approach are utilized to mathematically modeling of the structure. The fluid flow is assumed to be isentropic, incompressible and fully develop, and for more generality of the problem both steady and time dependent flow regimes are considered. The mathematical modeling of fluid flow is also carried out based on a scalar potential function, time mean Navier-Stokes equations and the theory of slip boundary condition. Employing the modified Lagrange equations for open systems, the nonlinear coupled governing equations of motion are achieved and solved via the state space problem; forth order numerical integration and Bolotin's method. In the numerical results, a comprehensive discussion is made on the dynamical instabilities of the system (such as divergence, flutter and parametric resonance). We found that applying positive electric potential field will improve the stability of the system as an actuator or vibration amplitude controller in the micro electro mechanical systems.

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.7
• /
• pp.56-66
• /
• 1999

Thermoacoustic oscillation is a significant problem in cylindrical-type combustors such as common internal combustion engines, industrial furnaces, gas turbine, etc. This kind of low frequency oscillation can lead to serious consequences such as destruction of the combustor and production of strong noise. The accurate numerical simulation of thermoacoustic phenomena is a complex and challenging problem, especially when considering the chemical reaction of mixtures. As with other simulations of aerodynamics and aeroacoustics, the direct computation of thermoacoustic phenomena requires that Navier-Stokes equations be solved using accurate numerical differentiation and time-marching schemes, with non-reflecting boundary conditions. The numerical approach used here aims at qualitative analysis and efficient prediction of those problems, not at the development of an accurate scheme. The numerical prediction developed in this work is shown to be reasonably matched with experimental result.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.8
• /
• pp.661-669
• /
• 2007

Study on the feasibility assessment for applying the overlap grid method to numerical calculations on the ground effect has been performed. The objective of the present study is to settle the problem in the grid generation process. A low Mach number preconditioned turbulent flow solver using the overlap grid and the multi-block grid methods has been developed and applied to the ground effect simulation around the RAE 101 airfoil. It has been verified that the overlap grid method not only can provide sufficiently accurate solutions but also work out the grid generation problem in the ground effect simulations.

• Wind and Structures
• /
• v.15 no.3
• /
• pp.189-208
• /
• 2012

Lock-in and drag amplification phenomena are studied for a flexible cantilever using a simplified fluid-structure interaction approach. Instead of solving the 3D domain, a simplified setup is devised, in which 2D flow problems are solved on a number of planes parallel to the wind direction and transversal to the structure. On such planes, the incompressible Navier-Stokes equations are solved to estimate the fluid action at different positions of the line-like structure. The fluid flow on each plane is coupled with the structural deformation at the corresponding position, affecting the dynamic behaviour of the system. An Arbitrary Lagrangian-Eulerian (ALE) approach is used to take in account the deformation of the domain, and a fractional-step scheme is used to solve the fluid field. The stabilization of incompressibility and convection is achieved through orthogonal quasi-static subscales, an approach that is believed to provide a first step towards turbulence modelling. In order to model the structural problem, a special one-dimensional element for thin walled cross-section beam is implemented. The standard second-order Bossak method is used for the time integration of the structural problem.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.2
• /
• pp.267-278
• /
• 1996

Steady, laminar, forced convection flow and heat transfer in the entrance region of an internally finned circular duct with a finite thermal conductivity has been analyzed numerically. The problem under investigation is a three-dimensional boundary layer problem, and is solved by employing a marching-type procedure which involves solution of a series of 2-dimensional elliptic problems in the cross-stream plane. Two types of inlet hydrodynamic conditions are considered : (a) uniform velocity flow and (b) fully developed flow. From the above inlet conditions, the effects of the fin height(h), fin number(N) and conductivity ratio($k_r$) on the flow and thermal characteristics are investigated. The numerical results show that the height and number of fins, and ratio of the solid to fluid thermal conductivity have pronounced effect on the solution. Considering pressure drop, optimized dimensionless fin height is 0.4.

• Korea-Australia Rheology Journal
• /
• v.13 no.3
• /
• pp.141-148
• /
• 2001

In this research, experimental studies leave been performed on the hydrodynamic interaction between a spherical particle and a plane wall by measuring the force between the particle and wall. To approach the system as a resistance problem, a servo-driving system was set-up by assembling a microstepping motor, a ball screw and a linear motion guide for the particle motion. Glycerin and dilute solution of polyacrylamide in glycerin were used as Newtonian and non-Newtonian fluids, respectively. The polymer solution behaves like a Boger fluid when the concentration is 1,000 ppm or less. The experimental results were compared with the asymptotic solution of Stokes equation. The result shows that fluid inertia plays all important role in the particle-wall interaction in Newtonian fluid. This implies that the motion of two particles in suspension is not reversible even in Newtonian fluid. In non-Newtonian fluid, normal stress difference and viscoelasticity play important roles as expected. In the dilute solution weak shear thinning and the migration of polymer molecules in the inhomogeneous flow field also affect the physic of the problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.7
• /
• pp.763-770
• /
• 2004

The problem of edge overcoating developed near the edge of the steel strip is studied quantitatively in the gas wiping process of continuous hot-dip galvanizing. It has been assumed that the edge overcoating occurs due to the reduced impact pressure of wiping gas on the strip edge and it is one of detrimental problems to the quality of coating products. In order to analyse the edge overcoating problem numerically, three-dimensional unsteady flows due to the gas wiping are calculated by using a commercial code, STAR-CD. Standard $\kappa$-$\varepsilon$ model is used as a turbulence model. The 1D code for calculation of coating thickness is constructed by using continuity and Navier-Stokes equations. The calculation results have shown good agreement with measurements of edge overcoating thickness, taken from galvanizing line trials. Therefore it is conformed that the major cause of edge overcoating is the reduced impact pressure of wiping gas on the strip surface.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.1
• /
• pp.51-57
• /
• 2009

In this study, a new topology optimization method is developed to design heat-dissipating structure with forced convection. To cool down electrical devices or mechanical machines, two types of convection models have been widely used: the natural convection model with a large Archimedes number and the forced convection with a small Archimedes number. In these days, lots of engineering application areas such as electrochemical conversion devices (Fuel cell) or rocket propulsion engines adopt the forced convection to dissipate the generated heat. Therefore, to our knowledge, it becomes an important issue to design flow channels inside which the generated heat dissipate. Thus, this paper studies optimal topological designs considering fluid-heat interactions. To consider the effect of the advection in the heat transfer problem, the incompressible Navier-stokes equation is solved. This paper numerically studies the coupling phenomena and presents optimal channel design considering forced convection.