• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.59-72
• /
• 1997

A numerical model is constructed to simulate the interactions of oblique shock wave / turbulent boundary layer on a strongly heated wall. The heated wall temperature is two times higher than the adiabatic wall temperature and the shock wave is strong enough to induce boundary layer separation. The numerical diffusion in the finite volume method is reduced by the use of a higher order convection scheme(UMIST scheme) which is a TVD version of QUICK scheme. The turbulence model is Chen-Kim two time scale model. The comparison of the wall pressure distribution with the experimental data ensures the validity of this numerical model. The effect of strong wall heating enlarges the separation region upstream and downstream. In order to eliminate the separation, wall suction is applied at the shock foot position. The bleeding slot width is about same as the upstream boundary layer thickness and suction mass flow is 10% of the flow rate in the upstream boundary layer. The final configuration of the shock reflection pattern and the wall pressure distribution approach to the non-viscous value when wall suction is applied.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.41-49
• /
• 1997

Computational results from Navier-Stokes equations are presented for the Stokes-wave/flat-plate boundary-layer and wake for small wave steepness(Ak=0.01), including exact and approximate treatments of the viscous free-surface boundary conditions. The macro-scale flow indicate that the variations of the external-flow pressure gradients cause acceleration or deceleration of the streamwise velocity component and alternating direction of the cross flow. Remarkably, the wake displays a greater response, i.e., a bias with regard to favorable as compared to adverse pressure gradients. The micro-scale flow indicates that the free-surface boundary conditions have a profound influence over the boundary layer and near/intermediate wake. Order-of-magnitude estimates are conformed to the computational results. And appreciable errors are introduced through approximations to the free-surface boundary conditions.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.2
• /
• pp.100-106
• /
• 1997

Viscous flow around an actual ship is calculated by an use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stress is modelled by using k-$\varepsilon$ turbulence model and the law of wall is applied near the body. Body fitted coordinates are introduced for the treatment of the complex boundary of the ship hull form. The transformed equations in the computational domain are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implcit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the disssssssscretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). The subject ship model of actual calculation is 4,410 TEU class container carrier. For 4 geosim models the calculated viscous resistancce values are compared with the model test results and analyzed on their componentss. The resistance performance of an actual ship is predicted very resonably, so this mothod may be utilized as a design tool of hull form.

• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.515-526
• /
• 2017

Dynamic response of functionally graded Carbon nanotubes (FG-CNT) reinforced pipes conveying viscous fluid under accelerated moving load is presented. The mixture rule is used for obtaining the material properties of nano-composite pipe. The radial force induced by viscous fluid is calculated by Navier-Stokes equation. The material properties of pipe are considered temperature-dependent. The structure is simulated by Reddy higher-order shear deformation shell theory and the corresponding motion equations are derived by Hamilton's principal. Differential quadrature (DQ) method and the Integral Quadrature (IQ) are applied for analogizing the motion equations and then the Newmark time integration scheme is used for obtaining the dynamic response of structure. The effects of different parameters such as boundary conditions, geometrical parameters, velocity and acceleration of moving load, CNT volume percent and distribution type are shown on the dynamic response of pipe. Results indicate that increasing CNTs leads to decrease in transient deflection of structure. In accelerated motion of the moving load, the maximum displacement is occurred later with respect to decelerated motion of moving load.

• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.66-69
• /
• 2006

In the present study, we propose a virtual boundary method for simulation of massive inextensible flexible strings immersed in viscous fluid flow. The fluid motion is governed by the Navier-Stokes equations and a momentum forcing is added in order to bring the fluid to move at the same velocity with the immersed surface. A massive inextensible flexible string model is described by another set of equations with an additional momentum forcing which is a result of the fluid viscosity and the pressure difference across the string. The momentum forcing is calculated by a feedback loop. Simulations of several numerical examples are carried out, inlcuding a hanging string which starts moving under gravity without ambient fluid, a string swimming within a uniform flow and a uniform flow over two side-by side strings. The numerical results agree well with the theoretical analysis and previous experimental observations. Preliminary results of a swimming elongated fishlike body will also be presented.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.8 s.239
• /
• pp.940-947
• /
• 2005

An immersed boundary method for simulation of density-stratified flows has been developed and applied to computation of viscous flows past three different types of obstacle under table density stratification, namely laminar flows past a vertical barrier, a cosine hill, and a sphere, respectively. Density forcing is introduced on the body surface or inside the body. Significant changes in flow characteristics are observed depending on Fr. The numerical results are in good agreement with other authors' experimental and numerical results currently available, and shed light on computation of density-stratified flows in complex geometries.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1991.04a
• /
• pp.8-14
• /
• 1991

The finite element technique incorporating infinite elements is applied to analyzing the general three dimensional wave-structure interaction problems within the limits of linear wave theory. The hydrodynamic farces are assumed to be inertially dominated, and viscous effects are neglected. In order to analyze the corresponding boundary value problems efficiently, two types of elements are developed. One is the infinite element for modeling the radiation condition at infinity, and the other is the fictitious bottom boundary element for the case of deep water. To validate those elements, numerical analyses are performed for several floating structures. Comparisons with the results from culler available solution methods show that the present method incorporating tile infinite and the fictitious bottom boundary elements gives good results.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.1
• /
• pp.39-61
• /
• 2017

In this present article, we analyzed the heat and mass transfer characteristics of the nonlinear unsteady radiative MHD flow of a viscous, incompressible and electrically conducting fluid past an infinite vertical porous plate under the influence of Soret and chemical reaction effects. The effect of physical parameters are accounted for two distinct types of thermal boundary conditions namely prescribed uniform wall temperature thermal boundary condition and prescribed heat flux thermal boundary condition. Based on the flow nature, the dimensionless flow governing equations are resolved to harmonic and non harmonic parts. In particular skin friction coefficient, Nusselt number and Sherwood number are found to evolve into their steady state case in the large time limit. Parametric study of the solutions are conducted and discussed.

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

The complete, fully-elliptic Reynolds-averaged Navier-Stokes equations have been solved using a two-layer model, in the $\kappa-\varepsilon$ turbulence model, for the axisymmetric body. Numerically generated boundary-fitted coordinate system and the finite analytic methods are used to solve the governing equations. Calculations are started after the middle body with given inlet conditions. The velocities and the turbulent quantities at the inlet section are specified by solving the boundary layer equations or by standard flat-plate boundary profiles. The effects of the inlet conditions on the solution are investigated.