• Coupled systems mechanics
• /
• v.9 no.3
• /
• pp.289-309
• /
• 2020

The paper studies the fluid flow profile contained between the orthotropic plate and rigid wall under the action of the moving load on the plate and main attention is focused on the fluid velocity profile in the load moving direction. It is assumed that the plate material is orthotropic one and the fluid is viscous and barotropic compressible. The plane-strain state in the plate and the plane flow of the fluid is considered. The motion of the plate is described by utilizing the exact equations of elastodynamics for anisotropic bodies, however, the flow of the fluid by utilizing the linearized Navier-Stokes equations. For the solution of the corresponding boundary value problem, the moving coordinate system associated with the moving load is introduced, after which the exponential Fourier transformation is employed with respect to the coordinate which indicates the distance of the material points from the moving load. The exact analytical expressions for the Fourier transforms of the sought values are obtained, the originals of which are determined numerically. Presented numerical results and their analyses are focused on the question of how the moving load acting on the face plane of the plate which is not in the contact with the fluid can cause the fluid flow and what type profile has this flow along the thickness direction of the strip filled by the fluid and, finally, how this profile changes ahead and behind with the distance of the moving load.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.10
• /
• pp.754-760
• /
• 2008

Computation of moving interface by the level set method typically requires the reinitialization of level set function. An inaccurate estimation of level set function $\phi$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, an accurate and robust reinitialization process is essential to the simulation of free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates level set function directly using a normal vector on the interface without solving there-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1 splitting/SUPG (Streamline Upwind Petrov-Galerkin) FEM are adopted to discretize advection equation of the level set function and the incompressible Navier-Stokes equation, respectively. Advection equation and re-initialization process of free surface capturing are validated with benchmark problems, i.e., a broken dam flow and timereversed single vortex flow. The simulation results are in good agreement with the existing results.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.15-25
• /
• 2004

This paper outlines a new convergence acceleration de-signed to solve scramjet flowfields with zones of re-circulation. Named the “marching-window”, the algorithm consists of performing pseudo-time iterations on a minimal width subdomain composed of a sequence of cross-stream planes of nodes. The upstream boundary of the subdomain is positioned such that all nodes upstream exhibit a residual smaller than the user-specified convergence threshold. The advancement of the downstream boundary follows the advancement of the upstream boundary, except in zones of significant streamwise ellipticity where a streamwise ellipticity sensor ensures its continuous progress. Compared to the standard pseudo-time marching approach, the march-ing-window is here seen to decrease the work required for convergence by up to 24 times for supersonic flows with little streamwise ellipticity and by up to 8 times for supersonic flows with large streamwise separated regions. The memory requirements are observed to be reduced sixfold by not allocating memory to the nodes not included in the computational subdomain. The marching-window satisfies the same convergence criterion as the standard pseudo-time stepping methods, hence resulting in the same converged solution within the tolerance of the user-specified convergence threshold. The extension of the marching-window to the weakly-ionized Navier-Stokes equations is also discussed.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.95-104
• /
• 2002

Numerical solution of the full Navier-Stokes equation as well as the energy equation has been obtained for the unsteady natural convection in a rectangular enclosure. One side wall was maintained at very high temperature simulating fires. Especially the effect of surface radiation was taken into account. While the enclosed air was assumed to be transparent, the internal walls directly interacted one another through the surface radiation. Due to a significant temperature difference in the flow field, the equation of state was used instead of the Boussinesq approximation. It was found that the rapid heating of the adiabatic ceiling and floor by the incoming radiation from the hot wall made the evolution at thermo-fluid field highly unstable in the initial period. Therefore, the secondary cells brought about at the floor region greatly affected the heat transfer mechanism inside the enclosure. The heat transfer rate was augmented by the radiation, resulting in requiring less time for the flow to reach the steady state. At the steady state neglecting radiation two internal hydraulic jumps were clearly observed in upper/left as well as in lower/right comer. However, the hydraulic jump in the lower/right comer could not be observed for the case including radiation due to its high momentum flow over the bottom wall. Radiation resulted in a faster establishment of the steady state phenomena.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.1 no.1
• /
• pp.1-12
• /
• 2009

Wing-in-Ground vehicles and aerodynamically assisted boats take advantage of increased lift and reduced drag of wing sections in the ground proximity. At relatively low speeds or heavy payloads of these craft, a flap at the wing trailing-edge can be applied to boost the aerodynamic lift. The influence of a flap on the two-dimensional NACA 4412 airfoil in viscous ground-effect flow is numerically investigated in this study. The computational method consists of a steady-state, incompressible, finite volume method utilizing the Spalart-Allmaras turbulence model. Grid generation and solution of the Navier-Stokes equations are completed using computer program Fluent. The code is validated against published experimental and numerical results of unbounded flow with a flap, as well as ground-effect motion without a flap. Aerodynamic forces are calculated, and the effects of angle of attack, Reynolds number, ground height, and flap deflection are presented for a split and plain flap. Changes in the flow introduced with the flap addition are also discussed. Overall, the use of a flap on wings with small attack angles is found to be beneficial for small flap deflections up to 5% of the chord, where the contribution of lift augmentation exceeds the drag increase, yielding an augmented lift-to-drag ratio.

• Ocean Systems Engineering
• /
• v.9 no.2
• /
• pp.111-133
• /
• 2019

Flow through a scaled horizontal axis marine current turbine was numerically simulated after validation and the turbine design was optimized. The computational fluid dynamics (CFD) code Ansys-CFX 16.1 for numerical modeling, an in-house blade element momentum (BEM) code for analytical modeling and an in-house surrogate-based optimization (SBO) code were used to find an optimal turbine design. The blade-pitch angle (${\theta}$) and the number of rotor blades (NR) were taken as design variables. A single objective optimization approach was utilized in the present work. The defined objective function was the turbine's power coefficient ($C_P$). A $3{\times}3$ full-factorial sampling technique was used to define the sample space. This sampling technique gave different turbine designs, which were further evaluated for the objective function by solving the Reynolds-Averaged Navier-Stokes equations (RANS). Finally, the SBO technique with search algorithm produced an optimal design. It is found that the optimal design has improved the objective function by 26.5%. This article presents the solution approach, analysis of the turbine flow field and the predictability of various surrogate based techniques.

• Coupled systems mechanics
• /
• v.8 no.3
• /
• pp.199-218
• /
• 2019

This paper studies the forced vibration of the hydro-elastic system consisting of the anisotropic (orthotropic) plate, compressible viscous fluid and rigid wall within the scope of the exact equations and relations of elastodynamics for anisotropic bodies for describing of the plate motion, and with utilizing the linearized exact Navier-Stokes equations for describing of the fluid flow. For solution of the corresponding boundary value problem it is employed time-harmonic presentation of the sought values with respect to time and the Fourier transform with respect to the space coordinate on the coordinate axis directed along the plate length. Numerical results on the pressure acting on the interface plane between the plate and fluid are presented and discussed. The main aim in this discussion is focused on the study of the influence of the plate material anisotropy on the frequency response of the mentioned pressure. In particular, it is established that under fixed values of the shear modulus of the plate material a decrease in the values of the modulus of elasticity of the plate material in the direction of plate length causes to increase of the absolute values of the interface pressure. The numerical results are presented not only for the viscous fluid case but also for the inviscid fluid case.

• Ocean Systems Engineering
• /
• v.12 no.1
• /
• pp.23-38
• /
• 2022

In this paper, a pitch-type wave energy converter (WEC-rotor) is investigated in irregular wave conditions for the real sea testing at the west coast of Jeju Island, South Korea. The present research builds on and extends our previous work on regular waves to irregular waves. The hydrodynamic characteristics of the WEC-rotor are assessed by establishing a quasi-two-dimensional numerical wave tank using computational fluid dynamics by solving the Reynolds-averaged Navier-Stokes equation. The numerical solution is validated with physical experiments, and the comparison shows good agreement. Furthermore, the hydrodynamic performance of the WEC-rotor is explored by investigating the effect of the power take-off (PTO) loading torque by one-way and two-way systems, the wave height, the wave period, operational and high sea wave conditions. Irrespective of the sea wave conditions, the absorbed power is quadratic in nature with the one-way and two-way PTO loading systems. The power absorption increases with the wave height, and the increment is rapid and mild in the two-way and one-way PTO loading torques, respectively. The pitch response amplitude operator increases as the wave period increases until the maximum value and then decreases. For a fixed PTO loading, the power and efficiency are higher in the two-way PTO loading system than in the one-way PTO loading system at different wave periods.

• Ocean Systems Engineering
• /
• v.4 no.3
• /
• pp.185-200
• /
• 2014

The objective of the present simulations is to evaluate the applicability of the standard $k-{\varepsilon}$ turbulence model in engineering practice in the subcritical to supercritical flow regimes. Two-dimensional numerical simulations of flow around a circular cylinder at $Re=1{\times}10^5$, $5{\times}10^5$ and $1{\times}10^6$, had been performed using Unsteady Reynolds-Averaged Navier Stokes (URANS) equations with the standard $k-{\varepsilon}$ turbulence model. Solution verification had been studied by evaluating grid and time step size convergence. For each Reynolds number, several meshes with different grid and time step size resolutions were chosen to calculate the hydrodynamic quantities such as the time-averaged drag coefficient, root-mean square value of lift coefficient, Strouhal number, the coefficient of pressure on the downstream point of the cylinder, the separation angle. By comparing the values of these quantities of adjacent grid or time step size resolutions, convergence study has been performed. Solution validation is obtained by comparing the converged results with published numerical and experimental data. The deviations of the values of present simulated quantities from those corresponding experimental data become smaller as Reynolds numbers increases from $1{\times}10^5$ to $1{\times}10^6$. This may show that the standard $k-{\varepsilon}$ model with enhanced wall treatment appears to be applicable for higher Reynolds number turbulence flow.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.3B
• /
• pp.293-303
• /
• 2011

An efficient diagonalized approximate factorization algorithm (DAF) is developed for the solution of three-dimensional incompressible viscous flows. The pressure-based, artificial compressibility (AC) method is used for calculating steady incompressible Navier-Stokes equations. The AC form of the governing equations is discretized in space using a second-order-accurate finite volume method. The present DAF method is applied to derive a second-order accurate splitting of the discrete system of equations. The primary objective of this study is to investigate the computational efficiency of the present DAF method. The solutions of the DAF method are evaluated relative to those of well-known four-stage Runge-Kutta (RK4) method for fully developed and developing laminar flows in curved square ducts and a laminar flow in a cavity. While converged solutions obtained by DAF and RK4 methods on the same computational meshes are essentially identical because of employing the same discrete schemes in space, both algorithms shows significant discrepancy in the computing efficiency. The results reveal that the DAF method requires substantially at least two times less computational time than RK4 to solve all applied flow fields. The increase in computational efficiency of the DAF methods is achieved with no increase in computational resources and coding complexity.