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

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.1
• /
• pp.23-30
• /
• 2004

In the Lagrangian vortex particle method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations, a numerical scheme for calculating pressure fields is presented. Implementation of the numerical method is directly connected with the well-established surface panel methods, just by dealing with the dynamic coupling among vorticity field. Assuming the vorticity and the velocity fields are to be calculated in time domain analysis, the pressure calculation for a complete set of solution at present time step is performed in a similar way to the one used in the Eulerian description. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsive started circular cylinder for Reynolds number 550. The comparative study with the Eulerian finite Volume method provides an extensive understanding and application of the mesh-free Lagrangian vortex methods for numerical simulation of viscous flows around arbitrary bodies of general shape.

• Atmosphere
• /
• v.30 no.4
• /
• pp.347-360
• /
• 2020

The automatic algorithm optimized for the Korean Peninsula was developed to detect and track the center of typhoon based on a geometrical method using high-resolution retrieved WISSDOM (WInd Syntheses System using DOppler Measurements) wind and reflectivity data. This algorithm analyzes the center of typhoon by detecting the geometric circular structure of the typhoon's eye in radar reflectivity and vorticity 2D field data. For optimizing the algorithm, the main factors of the algorithm were selected and the optimal thresholds were determined through sensitivity experiments for each factor. The center of typhoon was detected for 5 typhoon cases that approached or landed on Korean Peninsula. The performance was verified by comparing and analyzing from the best track of Korea Meteorological Administration (KMA). The detection rate for vorticity use was 15% higher on average than that for reflectivity use. The detection rate for vorticity use was up to 90% for DIANMU case in 2010. The difference between the detected locations and best tracks of KMA was 0.2° on average when using reflectivity and vorticity. After the optimization, the detection rate was improved overall, especially the detection rate more increased when using reflectivity than using vorticity. And the difference of location was reduced to 0.18° on average, increasing the accuracy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1262-1272
• /
• 2003

A two-dimensional direct numerical simulation is performed to investigate the dynamic behaviors of a single vortex in counter reacting and non-reacting flow field. A predictor-corrector-type numerical scheme with a low Mach number approximation is used in this simulation. A 16-step augmented reduced mechanism is adopted to treat the chemical reaction. The budget of the vorticity transport equation is examined to reveal a mechanism leading to the formation, destruction and transport of a single vortex according to the direction of vortex generation in reacting and non-reacting flows. The results show that air-side vortex has more larger strength than that of fuel-side vortex in both non-reacting and reacting flows. In reacting flow, the vortex is more dissipated than that in non-reacting flow as the vortex approach the flame. The total circulation in reacting flow, however, is larger than that in non-reacting flow because the convection transport of vorticity becomes much large by the increased velocity near the flame region. It is also found that the stretching and the convection terms mainly generate vorticity in non-reacting and reacting flows. The baroclinic torque term generates vorticity, while the viscous and the volumetric expansion terms attenuate vorticity in reacting flow. Furthermore, the contribution of volumetric expansion term on total circulation for air-side vortex is much larger than that of fuel-side vortex. It is also estimated that the difference of total circulation near stagnation plane according to the direction of vortex generation mainly attributes to the convection term.

• Journal of Astronomy and Space Sciences
• /
• v.25 no.4
• /
• pp.405-414
• /
• 2008

To better understand the physical processes that control the high-latitude lower thermospheric dynamics, we analyze the divergence and vorticity of the high-latitude neutral wind field in the lower thermosphere during the southern summertime for different IMF conditions. For this study the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEG CM) is used. The analysis of the large-scale vorticity and divergence provides basic understanding flow configurations to help elucidate the momentum sources that ulti-mately determine the total wind field in the lower polar thermosphere and provides insight into the relative strengths of the different sources of momentum responsible for driving winds. The mean neutral wind pattern in the high-latitude lower thermosphere is dominated by rotational flow, imparted primarily through the ion drag force, rather than by divergent flow, imparted primarily through Joule and solar heating. The difference vorticity, obtained by subtracting values with zero IMF from those with non-zero IMF, in the high-latitude lower thermosphere is much larger than the difference divergence for all IMF conditions, indicating that a larger response of the thermospheric wind system to enhancement in the momentum input generating the rotational motion with elevated IMF than the corresponding energy input generating the divergent motion. the difference vorticity in the high-latitude lower thermosphere depends on the direction of the IMF. The difference vorticity for negative and positive $B_y$ shows positive and negative, respectively, at higher magnetic latitudes than $-70^{\circ}$. For negative $B_z$, the difference vorticities have positive in the dusk sector and negative in the dawn sector. The difference vorticities for positive $B_z$ have opposite sign. Negative IMF $B_z$ has a stronger effect on the vorticity than does positive $B_z$.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.36 no.6
• /
• pp.1037-1048
• /
• 2016

In this study, nonlinear wave interaction in the presence of a uniform current is studied using numerical model, named CADMAS-SURF which is based on the Navier-Stokes equations coupled with Volume of Fluid for tracking free surface deformation. The original CADMAS-SURF developed for interaction of wave with structure is modified/extended to simulate nonlinear fluid dynamic motions within wave-current coexisting field. The capability of Numerical Wave-Current Tank (NWCT) in this study is validated by comparing with available existing laboratory experiments for both wave-following and wave-opposing current. The numerical results for interaction between wave and current are shown to be in good agreement with experimental data. Then, this study focused on the dynamic motions of the water velocity, surface elevation and vorticity within combined wave-current field in demonstrating complex nonlinear physical phenomena due to interaction between wave and current. In addition, NWCT is applied to simulate a more complex wave-current-structure field for wave propagating over a submerged breakwater associated with current. Detailed discussion including characteristics of velocity and vorticity fields and the relation between free surface and vorticity are given.

• Journal of the Korean Society of Visualization
• /
• v.3 no.2
• /
• pp.45-50
• /
• 2005

Electro-osmotic flow (EOF) instability in a microchannel has been experimentally investigated using a micro-PIV system. The micro-PIV system consisting of a two-head Nd:Yag laser and cooled CCD camera was used to measure instantaneous velocity fields and vorticity contours of the EOF instability in a T-shape glass microchannel. The electrokinetic flow instability occurs in the presence of electric conductivity gradients. Charge accumulation at the interface of conductivity gradients leads to electric body forces, driving the coupled flow and electric field into an unstable dynamics. The threshold electric field above which the flow becomes unstable and rapid mixing occurs is about 1000V/cm. As the electric field increases, the flow pattern becomes unstable and vortical motion is enhanced. This kind of instability is a key factor limiting the robust performance of complex electrokinetic bio-analytical devices, but can also be used for rapid mixing and effective flow control fer micro-scale bio-chips.

• Journal of KSNVE
• /
• v.7 no.2
• /
• pp.215-222
• /
• 1997

The flow and acoustic fields due to a vortex ring interaction with a rigid sphere are simulated numerically. The flow field is regarded as three-dimensional inviscid and incompressible. The vorticity is assumed to be concentrated inside the finite core of vortex filament. The vortex filament curve, described by parabolic blending curve function, is used to effectively solve the modified Biot-Savart equation. The interaction between a vortex ring and a rigid sphere using the parabolic blending curve is calculated. The trajectory of the vortex ring is obtained with several different initial positions between the ring and the sphere. The force variations acting on the sphere are calculated by using the boundary integral method. Finally, we can also obtain the acoustic signals at the far field observation positions from the force variations acting on the rigid surface. We can find that the dipole axis of the directivity patterns are rotated during the interacting phenomena.

• Journal of the Chungcheong Mathematical Society
• /
• v.27 no.4
• /
• pp.635-641
• /
• 2014

The existence and the regularity of the particle trajectory flow X(x, t) along a velocity field u on $\mathbb{R}^n$ are discussed under the BMO-blow-up condition: $${\int}_{0}^{T}{\parallel}{\omega}({\tau}){\parallel}_{BMO}d{\tau}<{\infty}$$ of the vorticity ${\omega}{\equiv}{\nabla}{\times}u$. A comment on our result related with the mystery of turbulence is presented.

• 한국연소학회:학술대회논문집
• /
• 2003.12a
• /
• pp.185-193
• /
• 2003

Basic flow fields of two dimensional V flames were examined as a preliminary work to study the instability of premixed flame with vorticity generation. Laminar premixed propane and methane flame were anchored by electrically heated wire to make two dimensional V flames. Flow fields were measured mainly by PIV(Particle Image Velocimetray) and the results were compared with those obtained by LDV(Laser Doppler Velocimetry) to confirm their reliability. Because the curvatures of V flames are so small, V flames were locally assumed to be inclined planar flames in gravitational field. The measured flow fields were locally compared with those of analytical solutions, which showed the qualitatively similar results. In downstream region, the vorticity fields were nearly constant except region near the center line, which support the assumption of locally one dimensional flame. Besides it was tried to find experimentally the similarity of flow fields in downstream region. Finally, stability diagram of propane and methane flames were drawn for the equivalence ratio less than one and the wide range of mean velocity.