• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.112-119
• /
• 2016

NEXTSat-1 is the next-generation small-size artificial satellite system planed by the Satellite Technology Research Center(SatTReC) in Korea Advanced Institute of Science and Technology(KAIST). For the control of attitude and transition of the orbit, the system has adopted a RHM(Resisto-jet Head Module), which has a very simple geometry with a reasonable efficiency. An axisymmetric model is devised with two coil-resistance heaters using xenon(Xe) gas, and the minimum required specific impulse is 60 seconds under the thrust more than 30 milli-Newton. To design the module, seven basic parameters should be decided: the nozzle shape, the power distribution of heater, the pressure drop of filter, the diameter of nozzle throat, the slant length and the angle of nozzle, and the size of reservoir, etc. After quasi one-dimensional analysis, a theoretical value of specific impulse is calculated, and the optima of parameters are found out from the baseline with a series of multi-physical numerical simulations based on the compressible Navier-Stokes equations for gas and the heat conduction energy equation for solid. A commercial code, COMSOL Multiphysics is used for the computation with a FEM (finite element method) based numerical scheme. The final values of design parameters indicate 5.8% better performance than those of baseline design after the verification with all the tuned parameters. The present method should be effective to reduce the time cost of trial and error in the development of RHM, the thruster of NEXTSat-1.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.5
• /
• pp.16-21
• /
• 2010

During the loading/offloading operation of a liquefied natural gas carrier (LNGC) that is moored in a side-by-side configuration with an offshore plant, sloshing that occurs due to the partially filled LNG tank and the interactive effect between the two floating bodies are important factors that affect safety and operability. Therefore, a time-domain software program, called CHARM3D, was developed to consider the interactions between sloshing and the motion of a floating body, as well as the interactions between multiple bodies using the potential-viscous hybrid method. For the simulation of a floating body in the time domain, hydrodynamic coefficients and wave forces were calculated in the frequency domain using the 3D radiation/diffraction panel program based on potential theory. The calculated values were used for the simulation of a floating body in the time domain by convolution integrals. The liquid sloshing in the inner tanks is solved by the 3D-FDM Navier-Stokes solver that includes the consideration of free-surface non-linearity through the SURF scheme. The computed sloshing forces and moments were fed into the time integration of the ship's motion, and the updated motion was, in turn, used as the excitation force for liquid sloshing, which is repeated for the ensuing time steps. For comparison, a sloshing motion coupled analysis program based on linear potential theory in the frequency domain was developed. The computer programs that were developed were applied to the side-by-side offloading operation between the offshore plant and the LNGC. The frequency-domain results reproduced the coupling effects qualitatively, but, in general, the peaks were over-predicted compared to experimental and time-domain results. The interactive effects between the sloshing liquid and the motion of the vessel can be intensified further in the case of multiple floating bodies.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.6
• /
• pp.15-20
• /
• 2005

Two-dimensional incompressible Navier-Stokes flow solver is developed using SIMPLER method to study the unsteady viscous flow physics over two-dimensional ellipses. Unsteady viscous flows over various thickness-to-chord ratios of 0.6, 0.8, 1.0, and 1.2 elliptic cylinders are simulated at different Reynolds numbers of 200, 400, and 1,000. This study is focused on the understanding the effects of Reynolds number and elliptic cylinder thickness on the drag and lift forces. The present numerical solutions are compared with available experimental and numerical results and show a good agreement. Through this study, it is observed that the Reynolds number and the cylinder thickness affect not only the frequency of the force oscillations but also the mean values and the amplitudes of the total drag and lift forces significantly.

• Wind and Structures
• /
• v.15 no.6
• /
• pp.481-494
• /
• 2012

The greenhouse type metal structures are increasingly used in modern construction of livestock farms because they are less laborious to construct and they provide a more favorable microclimate for the growth of animals compared to conventional livestock structures. A key stress factor for metal structures is the wind. The external pressure coefficient ($c_{pe}$) is used for the calculation of the wind effect on the structures. A high pressure coefficient value leads to an increase of the construction weight and subsequently to an increase in the construction cost. The EC1 in conjunction with EN 13031-1:2001, which is specialized for greenhouses, gives values for this coefficient. This value must satisfy two requirements: the safety of the structure and a reduced construction cost. In this paper, the Navier - Stokes and continuity equations are solved numerically with the finite element method (Galerkin Method) in order to simulate the two dimensional, incompressible, viscous air flow over the vaulted roofs of single span and twin-span with eaves livestock greenhouses' structures, with a height of 4.5 meters and with length of span of 9.6 and 14 m. The simulation was carried out in a wind tunnel. The numerical results of pressure coefficients, as well as, the distribution of them are presented and compared with data from Eurocodes for wind actions (EC1, EN 13031-1:2001). The results of the numerical experiment were close to the values given by the Eurocodes mainly on the leeward area of the roof while on the windward area a further segmentation is suggested.

• Ocean Systems Engineering
• /
• v.3 no.3
• /
• pp.219-236
• /
• 2013

The accurate prediction of motion in waves of a marine vehicle is essential to assess the maximum sea state vs. operational requirements. This is particularly true for small crafts, such as Autonomous Surface Vessels (ASV). Two different numerical methods to predict motions of a SWATH-ASV are considered: an inviscid strip theory initially developed at MIT for catamarans and then adapted for SWATHs and new a hybrid strip theory, based on the numerical solution of the radiation forces by an unsteady viscous, non-linear free surface flow solver. Motion predictions obtained by the viscous flow method are critically discussed against those obtained by potential flow strip theory. Effects of viscosity are analyzed by comparison of sectional added mass and damping calculated at different frequencies and for different sections, RAOs and motions response in irregular waves at zero speed. Some relevant conclusions can be drawn from this study: influence of viscosity is definitely non negligible for SWATH vessels like the one presented: amplitude of the pitch and heave motions predicted at the resonance frequency differ of 20% respectively and 50%; in this respect, the hybrid method with fully non-linear, viscous free surface calculation of the radiation forces turns out to be a very valuable tool to improve the accuracy of traditional strip theories, without the burden of long computational times requested by fully viscous time domain three dimensional simulations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.9
• /
• pp.2944-2952
• /
• 1996

In this paper the CSCM type upwind flux difference splitting Navier-Stokes method has been applied to study the ARL-SL19 transoni $c^ersonic compressor cascade flow. First, the general characteristics of baseline cascade flow were analyzed. At freestream Mach n.1.612 and exit/inlet pressure ratio 2.15, the results from current laminar flow were compared well in suction surface with the experiment; however, not well in pressure surface. Second, numerical study of the transoni $c^ersonic compressor cascade flow demonstrated the effectiveness of a passive control by the various size cavities. A cavity under the shock foot point at the suction surface of the blades was used as a passive control. The passive control of shock-boundary layer interaction by a cavity reduced total pressure losses. The effect of cavity length and depth was studied. The total pressure loss was reduced by about 10% and the isentropic efficiency was improved slightly. The effect of cavity depth in current study(d/l = 0.05, 0.02) was not found strong. Further adequate turbulence modeling and TVD schemes would help to capture the shock more accurately and increase the effectiveness of the current shock-boundary layer interaction study using upwind flux difference splitting computational methods.thods.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.4
• /
• pp.33-38
• /
• 2003

The Physics of the Plume-induced shock and separation Particularly at a high Plume to exit pressure ratio and supersonic speeds up to Mach 3.0 with and without a passive control method, porous extension, were studied using computational techniques. Mass-averaged Navier-Stokes equations with the RNG $\kappa$-$\varepsilon$ turbulence model were solved using a fully implicit finite volume scheme and a 4-stage Runge-Kutta method. The control methodology for plume-afterbody interactions is to use a perforated wall attached at either the nozzle exit or the edge of the missile base. The Effect of porous wall length on plume interference is also investigated The computational results show the main effect of the porous extension on plume-afterbody interactions is to restrain the plume from strongly underexpanding during a change in flight conditions. With control, a change in porous extension length has no significant effect rut plume interference.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.4
• /
• pp.90-97
• /
• 2013

Oil booms are one of the most widely used types of equipment for the protection of coastal areas against oil spills. In some situations, however, there are several types of oil leaks from the oil boom. Important factors regarding these phenomena include the surrounding ocean environment, such as waves, the density and viscosity of oil, the length of the oil boom skirt, etc. To estimate the performance of the oil boom, it is necessary to predict the behavior of the spilled oil and oil boom. In the present study, the prediction of oil boom performance in waves was carried out using the Pusan-National-University-modified Moving Particle Semi-implicit (PNU-MPS) method, which is an improved version of the original MPS proposed by Koshizuka and Oka (1996). The governing equations, which consist of continuity and Navier-Stokes equations, are solved by Lagrangian moving particles, and all terms expressed by differential operators in the governing equations are replaced by the particle interaction models based on a kernel function. The simulation results were validated through a comparison with the results of Violeau et al. (2007)..

• Wind and Structures
• /
• v.17 no.6
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.127-138
• /
• 1999

In this study steady state solutions of cavity flows driven by two moving walls are studied. The north and east walls of the cavity are movable where as the remaining two walls are fixed in space. Numerical experiments for three different driving schemes for moving walls are done at two different Reynolds numbers of Re=40 and 400. The first scheme is to accelerate north and east walls simultaneously. In the second one, the north wall is started first and the east wall is accelerated later. In the third one the east wall starts first. It is usually expected that all these three cases yield the same steady state solution after sufficiently long time. However, present numerical experiments show that such a usual belief is valid only when the Reynolds number is low enough (Re=40). At higher Reynolds number (Re=400), the flow develops to three different steady states depending on the history of the boundary condition change.