• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.9
• /
• pp.789-795
• /
• 2014

In this study, the distribution characteristics of diffusers with various shapes that are installed in an open-type thermal storage system are numerically investigated. Four diffusers are designed to distribute a working fluid evenly through the holes on bifurcated pipes. Three-dimensional steady simulations of incompressible laminar flow are conducted using commercial software (ANSYS-FLUENT). The simulation results show that both the bidirectional header-type diffuser and the H-type diffuser distribute the working fluid evenly whereas both the unidirectional and the bidirectional diffusers distribute the working fluid unevenly. The results also show that the H-type diffuser requires a higher head of pump than the bidirectional header-type diffuser. Therefore, the bidirectional header-type diffuser is recommended for use because it enables even distribution of the working fluid and requires a low head of pump.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.28-40
• /
• 2002

Turbulent flows around two-dimensional wing sections in ground effect are analysed by incompressible RANS equations and a finite difference method. The Baldwin-Lomax algebraic turbulence model is used to simulate high Reynolds number flows. The main purpose of this study is to clarify the two-dimensional ground effect and its flow characteristics due to different ground boundary conditions, i.e., moving and fixed bottom boundary. As a first step, to validate the present numerical code, the computational result of Clark-Y(t/C 11.7%) is compared with published numerical results and experimental data. Then, NACA4412 section in ground effect is calculated for various ground clearances with two bottom boundary conditions. According to the computational results, the difference in the lift and moment simulated with the two bottom boundary conditions is negligible, but the drag force simulated by the fixed bottom is to some extent smaller than that by the moving bottom. Therefore, it can be concluded that the drag force measured in a wind tunnel with the fixed bottom could be smaller than that with the moving bottom.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.29 no.8
• /
• pp.419-428
• /
• 2017

In this study, energy loss due to ventilation load in the dock system was analyzed through simulation. Also, flow generated in the dock system of the warehouse was measured using manufactured measuring devices. Numerical simulation was conducted by simulating the most common picking tasks by examining the actual working environment. Incompressible and unsteady turbulent flows were assumed, and the turbulence model was the k-e standard model. Proper grid was selected through grid dependency test. Measurement was conducted using Honeywell and Vaisala sensors, and flow and temperature inside the warehouse were measured and compared with simulation results to validate simulation. When comparing amount of loss occurring in two hours and amount of loss occurring in 15 minutes, docking time of the former was eight times longer but energy loss was 3.8 times lower. Ventilation load occurring during the initial period after opening docking system accounted for a large proportion of total ventilation load. Also, comparing the load when the dock was closed and the load when the truck was parked, ventilation load was significantly higher than load due to heat conduction from the wall. Therefore, in improving the docking system, it is effective to reduce the gap by improving compatibility of the docking system and truck, rather than wall material.

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

• International Journal of Fluid Machinery and Systems
• /
• v.9 no.4
• /
• pp.370-381
• /
• 2016

As a single-channel pump is used for wastewater treatment, this particular pump type can prevent performance reduction or damage caused by foreign substances. However, the design methods for single-channel pumps are different and more difficult than those for general pumps. In this study, a design optimization method to improve the hydrodynamic performance of a single-channel pump impeller is implemented. Numerical analysis was carried out by solving three-dimensional steady-state incompressible Reynolds-averaged Navier-Stokes equations using the shear stress transport turbulence model. As a state-of-the-art impeller design method, two design variables related to controlling the internal cross-sectional flow area of a single-channel pump impeller were selected for optimization. Efficiency was used as the objective function and was numerically assessed at twelve design points selected by Latin hypercube sampling in the design space. An optimization process based on a radial basis neural network model was conducted systematically, and the performance of the optimum model was finally evaluated through an experimental test. Consequently, the optimum model showed improved performance compared with the base model, and the unstable flow components previously observed in the base model were suppressed remarkably well.

• Advances in aircraft and spacecraft science
• /
• v.3 no.3
• /
• pp.243-257
• /
• 2016

The increase of air traffic volume has brought an increasing amount of issues related to carbon and NOx emissions and noise pollution. Aircraft manufacturers are concentrating their efforts to develop technologies to increase aircraft efficiency and consequently to reduce pollutant discharge and noise emission. Ultra High By-Pass Ratio engine concepts provide reduction of fuel consumption and noise emission thanks to a decrease of the jet velocity exhausting from the engine nozzles. In order to keep same thrust, mass flow and therefore section of fan/nacelle diameter should be increased to compensate velocity reduction. Such feature will lead to close-coupled architectures for engine installation under the wing. A strong jet-wing interaction resulting in a change of turbulent mixing in the aeroacoustic field as well as noise enhancement due to reflection phenomena are therefore expected. On the other hand, pressure fluctuations on the wing as well as on the fuselage represent the forcing loads, which stress panels causing vibrations. Some of these vibrations are re-emitted in the aeroacoustic field as vibration noise, some of them are transmitted in the cockpit as interior noise. In the present work, the interaction between a jet and wing or fuselage is reproduced by a flat surface tangential to an incompressible jet at different radial distances from the nozzle axis. The change in the aerodynamic field due to the presence of the rigid plate was studied by hot wire anemometric measurements, which provided a characterization of mean and fluctuating velocity fields in the jet plume. Pressure fluctuations acting on the flat plate were studied by cavity-mounted microphones which provided point-wise measurements in stream-wise and spanwise directions. Statistical description of velocity and wall pressure fields are determined in terms of Fourier-domain quantities. Scaling laws for pressure auto-spectra and coherence functions are also presented.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.80-86
• /
• 1994

A numerical method for simulations of inviscid incompressible flow fields around a ship advancing on the free surface is developed. A body fitted coordinate system, generated by numerically solving elliptic type partial differential equations is used to conform the ship and free surface configurations. Three dimensional Euler equations transformed to the non-staggered body fitted coordinate system are discretised by finite difference method. Time and spatial derivatives are discretised by forward and centered differencings, respectively, and artificial dissipations are added to discretised convection terms for improvements of numerical stability. At each time steps, free surface elevations are recomputed to satisfy nonlinear free surface conditions. Poisson equations for pressure field are solved iteratively and the velocity field for next time step is extrapolated. To verify the developed numerical method, flow fields around a Wigley model are simulated(Fn=0.250-0.408) and compared with experimental data to show good agreements.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.1
• /
• pp.1-9
• /
• 2021

In the present study, the model-scale Propeller Open Water (POW) tests for the propeller of 176K bulk carrier and 8600TEU container ship were conducted through Computational Fluid Dynamics (CFD) simulation. In order to solve the incompressible viscous flow field, the Reynolds-averaged Navier-Stokes (RaNS) equations were employed as the governing equations. The γ-Reθ(gamma-Re-theta) transition model combined with the SST k-ωturbulence model was introduced to describe the laminar-turbulence transition considering the low Reynolds number of model-scale. Firstly, the flow simulation developing over a flat plate was performed to verify the transition modeling, in which the wall shear stresses were compared with experiments and other numerical results. Then, to investigate the effect of the model, the CFD simulation for the POW test was performed and the simulated propeller performance was validated through comparison with the experiment conducted at Korea Research Institute of Ships & Ocean Engineering (KRISO).

• Wind and Structures
• /
• v.34 no.5
• /
• pp.407-419
• /
• 2022

Many factors should be considered by architects and designers for designing a tall building. Wind load is one of these important factors that govern the design of tall building structures and can become a serious challenge when buildings tend to be built very tall and slender. On the other hand, through the initial stages of a design process, choosing the design geometry greatly affects the wind-induced forces on a tall building. With this respect, geometric shapes with 3-fold rotational symmetry are one of the applied plan shapes in tall buildings. This study, therefore, aims to investigate the aerodynamic characteristics of 8 different geometrical shapes using Computational Fluid Dynamics (CFD) by measuring the drag and lift forces. A case study approach was conducted in which different building shape models have the same total gross area and the same height of 300 meters. The simulation was an incompressible transient flow that ran 1700 timesteps (85 seconds on the real-time scale). The results show a great difference between wind-induced force performance of buildings with different plan shapes. Generally, it is stated that the shapes with the same area, but with smaller perimeters, are better choices for reducing the drag force on buildings. Applying the lift force, the results show that the buildings with plan shapes that have rounded corners act better in crosswind flow while, those with sharp corners induce larger forces in the same direction. This study delivers more analytical understanding of building shapes and their behavior against the wind force through the parametric modelling.