• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• Journal of Aerospace System Engineering
• /
• v.16 no.4
• /
• pp.35-44
• /
• 2022

The study on the wake vortex behavior during the aircraft's take-off and landing flight phase is critical to the flight safety of the aircraft, following close behind and the economy of the airport. The study on the wake vortex behavior should include the understanding of the ground effect on the behavior of the multiple wake vortices, generated from aircraft during the take-off and landing flight phase. In thia study, numerical schemes that can consider the ground effect were devised, by applying a vorticity boundary condition and an image method into the existing two-dimensional Fourier-spectral method. The present method was validated by comparing the present results, with the computed and measured data in the published literature. It was shown that the present method can predict the generation and behavior of the secondary vortex near the ground with reasonable accuracy. In future, the effect of the atmospheric conditions such as the stratification and the wind shear on the behavior of the vortex pair will be studied.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.2B
• /
• pp.179-189
• /
• 2010

This paper presents a numerical investigation of vorticity generation in fully vegetated open-channel flows. The Reynolds stress model is used for the turbulence closure. Open-channel flows with rough bed-smooth sidewalls and smooth bed-rough sidewalls are simulated. The computed vectors show that in channel flows with rough bed and rough sidewalls, the free-surface secondary currents become relatively smaller and larger, respectively, compared with that of plain channel flows. Also, open-channel flows over vegetation are simulated. The computed bottom vortex occupies the entire water depth, while the free-surface vortex is reduced. The contours of turbulent anisotropy and Reynolds stress are presented with different density of vegetation. The budget analysis of vorticity equation is carried out to investigate the generation mechanism of secondary currents. The results of the budget analysis show that in plain open-channel flow, the production by anisotropy is important in the vicinity of the wall and free-surface boundaries, and the production by Reynolds stress is important in the region away from the boundaries. However, this rule is not effective in vegetated channel flows. Also, in plain channel flows, the vorticity is generated mainly in the vicinity of the free-surface and the bottom, while in vegetated channel flows, the regions of the bottom and vegetation height are important to generate the vorticity.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.6
• /
• pp.2267-2275
• /
• 2013

Existing conventional model for analysis of shallow water flow just assumed the internal boundary condition as free-slip, which resulted in the wrong prediction about the velocity, vorticity, water level, shear stress distribution, and time variation of drag and lift force around a structure. In this study, a finite element model that can predict flow characteristics around the structure accurately was developed and internal boundary conditions were generalized as partial slip condition using slip length concept. Laminar flow characteristics behind circular cylinder were analyzed by varying the internal boundary conditions. The simulation results of (1) time variations of longitudinal and transverse velocities, and vorticity; (2) wake length; (3) vortex shedding phenomena by slip length; (4) and mass conservation showed that the vortex shedding had never observed and laminar flow like creeping motion was occurred under free-slip condition. Assignment of partial slip condition changed the velocity distribution on the cylinder surface and influenced the magnitude of the shear stress and the occurrence of vorticity so that the period of vortex shedding was reduced compared with the case of no slip condition. The maximum mass conservation error occurred in the case of no slip condition, which had the value of 0.73%, and there was 0.21 % reduction in the maximum mass conservation error by changing the internal boundary condition from no slip to partial slip condition.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.6
• /
• pp.1163-1170
• /
• 1992

A computer program was developed to analyze the two-dimensional unsteady incompressible viscous flow behind a rectangular bluff body between two parallel plates. The Peaceman-Rachford alternating direction implicit numerical method and Wachspress parameter were adopted to solve the governing equations in vorticity-transport and stream function formulation. The steady state flow and the vortex flow behind a rectangular bluff body in a chemical were investigated for Reynolds numbers of 200 and 500. The vortex shedding was generated by a physical pertubation numerically imposed at the center of the flow field for a short time. It was observed that the perturbed flow became periodic after a transient period.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.2
• /
• pp.67-76
• /
• 1988

A fully non-linear two-dimensional model for studying basic tidal dynamics within the Arabian Gulf has been formulated. The model extends to the east as far as continental shelf edge and has a mesh resolution of 1/12 degree latitude by 1/10 degree longitude. The developed model was used to investigate the distribution of dominant $M_2$, $S_2$, $K_1$, $O_1$ and $N_2$ tides adopting the open boundary tidal information from Schwiderski's global ocean tidal model. Computed co amplitude and cophase charts were presented independently and comparison between these charts and existing Admiralty tidal charts were made. Maximum bottom stress vector during the semidiurnal period due to $M_2$ and $M_4$ tides are also presented to suggest the possible direction of sediment movement in the region.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.17 no.1
• /
• pp.29-37
• /
• 2011

PIV measurements of the velocity field, pressure field, vorticity, and turbulent intensity in the rear of curved section of an oil fence with current speed showed that the flow directions in the rear of flow boundary area were similar to those in the front of it. As the current speed increased, the patterns of pressure distribution were changed, and the turbulent flow became more irregular. CFD simulations under the same conditions as the PIV tests showed that the flow patterns of the wake were similar to those by PIV tests in speed of 0.3 m/s and less, but were distinctively deviated from those in 0.4 m/s due to the flexibility of the oil fence, which was not properly taken care of in CFD modeling.

• Journal of Environmental Science International
• /
• v.6 no.4
• /
• pp.323-333
• /
• 1997

An objective method for the generation of velocity streamfunction is presented for dealing with discretely sampled oceauc data. The method treats a Poisson equation (forced by vorticity) derived from Helmholtz theorem In which streamfunction is obtained by isolating the non-divergent part of the two-dimensional flow field. With a mixed boundary condition and vorticity field estimated from observed field, the method Is Implemented over the Texas-Louisiana show based on the current meter data of the Texas-Louisiana Shelf Circulation and Transport Processes Study (LATEX) measured at 31 moorings for 32 months (April 1992 - November 1994). The resulting streamfunction pattern is quote consistent with observations. The streamfunction field by this method presents an opportunity to initiauze and to verier computer models for local forecasts of enoronmental flow conditions for ell spill, nutrient and plankton transports as well as opportuuty to understand shelf-wlde low-frequency currents.