• International Journal of Naval Architecture and Ocean Engineering
• /
• v.4 no.3
• /
• pp.256-266
• /
• 2012

Cavitating flow simulation is of practical importance for many engineering systems, such as marine propellers, pump impellers, nozzles, torpedoes, etc. The present work has developed the base code to solve the cavitating flows past the axisymmetric bodies with several forebody shapes. The governing equation is the Navier-Stokes equation based on homogeneous mixture model. The momentum is in the mixture phase while the continuity equation is solved in liquid and vapor phase, separately. The solver employs an implicit preconditioning algorithm in curvilinear coordinates. The computations have been carried out for the cylinders with hemispherical, 1-caliber, and 0-caliber forebody and, then, compared with experiments and other numerical results. Fairly good agreements with experiments and numerical results have been achieved. It has been concluded that the present numerical code has successfully accounted for the cavitating flows past axisymmetric bodies. The present code has also shown the capability to simulate ventilated cavitation.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.11a
• /
• pp.187-193
• /
• 2000

In order to reduce the booming noise caused by first bending mode of drive shaft, this paper proposes a simulation program for prediction of the bending mode frequency of any tubular shaft. This program consists of a pre-processor for modeling of geometrical shape of drive shaft and applying the boundary conditions of various joints, a processor for constructing of global finite element matrices using beam elements and an eigen-solver based on MATLAB program. Using this simulation program, the effective and accurate FE model for a shaft attached in vehicle can be obtained by aid of database for stiffness of each joint. Thus the resonance frequencies and mode shapes of a shaft can be calculated accurately. Because the effect of the resonance on interior noise can be verified, more improved shaft can be proposed at the early stage of design.

• Ocean Systems Engineering
• /
• v.9 no.1
• /
• pp.43-77
• /
• 2019

A local-analytic-based Navier-Stokes solver has been employed in conjunction with a compound ocean structure motion analysis program for time-domain simulation of passing ship effects induced by multiple post-Panamax class ships in the exact condition of a real waterway. The exact seabed bathymetry was reproduced to the utmost precision attainable using the NOAA geophysical database for Virginia Beach, NOAA nautical charts for Hampton Roads and Norfolk harbor, and echo sounding data for the navigation channel and waterfront facilities. A parametric study consists of 112 simulation cases with various combinations of ship lanes, ship speeds, ship heading (inbound or outbound), channel depths, drift angles, and passing ship coupling (in head-on or overtaking encounters) were carried out for two waterfront facilities at NAVSTA Norfolk and Craney Island Fuel Terminal. The present paper provides detailed parametric study results at both locations to investigate the site-specific passing ship effects on the motion responses of ships moored at nearby piers.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.148-156
• /
• 2015

A continuous type side jet controller which has four nozzles with thrust control devices was considered. It is deployed to a missile for high maneuverability and fast controllability in the terminal guidance phase. However, it causes more complex aerodynamic jet interactions between the side jet and the supersonic free stream than does the conventional impulse type side jet with a small single thruster. In this paper, a numerical investigation of the jet interference effects for the missile equipped with a continuous type side jet thruster is presented. A three-dimensional flow field was simulated by using a commercial unstructured-based CFD solver. The numerical simulation method was validated through comparison with wind tunnel test results for the single jet. The method of defining jet direction for this type of side jet control to minimize simulation cases was also introduced. Flow fields investigation and jet interaction effects for various flow conditions, jet pressure ratios and defined jet direction conditions were performed. From the numerical simulation for the continuous type side jet, extensive aerodynamic interference data were obtained to construct an aerodynamic coefficients database for precise missile control.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.15 no.3
• /
• pp.114-121
• /
• 2007

In this study, a modified reverse flow type, one of the triple effect absorption cycles, is studied for performance improvement. The cycle simulation is carried out by using EES(Engineering Equation Solver) program for the working fluid of $H_2O/LiBr$ solution. The split-ratios of solution flow rate, UA of each component, pumping mass flow rate of solution are considered as key parameters. The results show that the optimal SRH (split ratio of high side) and SRL (split ratio of low side) values are 0.596 and 0.521, respectively. Under these conditions, the COP is maximized to 2.1. The optimal pumping mass flow rate is selected as 3 kg/s and the corresponding UAEV A is 121 kW/K in the present system. The present simulation results are compared to the other literature results from Kaita's (2002) and Cho's (1998) triple effect absorption systems. The present system has a lower solution temperature and a higher COP than the Kaita's modified reverse flow, and it also gives a higher COP than the Cho's parallel flow by adjusting split ratios.

• 한국연소학회:학술대회논문집
• /
• 2005.10a
• /
• pp.181-188
• /
• 2005

Large eddy simulation(LES) methodology used to model the isothermal swirling flows in a dump combustor and the turbulent premixed flame in a model gas turbine combustor. The LES solver was implemented on parallel computer consisting 16 processors. In isothermal flow simulation, the results was compared with that of ${\kappa}-{\varepsilon}$ model as well as experimental data, in order to verify the capability of LES code. To model the turbulent premixed flame in a gas turbine, the G-equation flamelet model was used. The results showd that LES and RANS well predicted the mean velocity field of a non-swirling flow. However, in swirling flow, LES showed a better performance in predicting the mean axial and azimuthal velocities, and the central recirculation zone than those of RANS. In a model gas turbine combustor, the operation condition of high pressure and temperature induced the different phenomena, such as flame length and flow-field information, comparing with the condition of ambient pressure and temperature. Finally, it was identified that the flame and heat release oscillations are related to the vortex shedding generated by swirl flow and pressure wave propagation.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.35D no.6
• /
• pp.37-45
• /
• 1998

In this paper, we report the first three-dimensional simulation result of the transient enhanced diffusion(TED) of dopants in the ion-implanted silicon by employing our 3D semiconductor process simulator, INPROS system. In order to simulate three-dimensional TED redistribution of dopants in silicon, the dopant distributions after the ion implantation was calculated by Monte Carlo(MC) method, followed by finite element(FE) numerical solver for thermal annealing. Excellent agreement between the simulated 3D profile and the SIMS data has been obtained for ion-implanted arsenic and phosphorus after annealing the boron marker layer at 75$0^{\circ}C$ for 2 hours. Our three-dimensional TED simulation could successfully explain the reverse short channel effect(RSCE) by taking the 3D point defect distribution into account. A coupled TED simulation and device simulation allows reverse short channel effect on threshold to be accurately predicted.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.6
• /
• pp.487-497
• /
• 2011

Numerical analysis has been performed to investigate the rotor speed and pitch range variations when the airspeed is increased in autorotation. Transient Simulation Method(TSM) was used to obtain the steady states of autorotation. The rotor blade was analyzed by the two-dimensional compressible Navier-Stokes solver in order to adapt to the airspeed increase and the results were used in the transient simulation method. Meanwhile, the Pitt/Peters inflow theory was used to supply the induced velocity fields. For the prescribed torque equilibrium state, the combinations of velocity, shaft angle, and pitch angle were produced to investigate the rotor speeds and variable ranges. The rotor tip Mach number and rotor speed were correlated and the trim range of pitch angle was observed with respect to the shaft angle decrease.

• Nuclear Engineering and Technology
• /
• v.51 no.8
• /
• pp.1916-1938
• /
• 2019

The objective of this paper is to discuss the modeling principles of phenomena governing core degradation/melting and in-vessel melt relocation during severe accidents in light water reactors. The proposed modeling approach has been applied in the development of a new accident simulation package, COMPASS (COre Meltdown Progression Accident Simulation Software). COMPASS can be used either as a stand-alone tool to simulate in-vessel meltdown progression up to and including RPV failure, or as a component of an integrated simulation package being developed in Korea for the APR1400 reactor. Interestingly, since the emphasis in the development of COMPASS modeling framework has been on capturing generic mechanistic aspects of accident progression in light water reactors, several parts of the overall model should be useful for future accident studies of other reactor designs, both PWRs and BWRs. The issues discussed in the paper include the overall structure of the model, the rationale behind the formulation of the governing equations and the associated simplifying assumptions, as well as the methodology used to verify both the physical and numerical consistencies of the overall solver. Furthermore, the results of COMPASS validation against two experimental data sets (CORA and PHEBUS) are shown, as well as of the predicted accident progression at TMI-2 reactor.

• Wind and Structures
• /
• v.20 no.1
• /
• pp.75-93
• /
• 2015

This work is focused on a parametric numerical study of the barrier's bar inclination shelter effect in crosswind scenario. The parametric study combines mesh morphing and design of experiments in automated manner. Radial Basis Functions (RBF) method is used for mesh morphing and Ansys Workbench is used as an automation platform. Wind barrier consists of five bars where each bar angle is parameterized. Design points are defined using the design of experiments (DOE) technique to accurately represent the entire design space. Three-dimensional RANS numerical simulation was utilized with commercial software Ansys Fluent 14.5. In addition to the numerical study, experimental measurement of the aerodynamic forces acting on a vehicle is performed in order to define the critical wind disturbance scenario. The wind barrier optimization method combines morphing, an advanced CFD solver, high performance computing, and process automaters. The goal is to present a parametric aerodynamic simulation methodology for the wind barrier shelter that integrates accuracy and an extended design space in an automated manner. In addition, goal driven optimization is conducted for the most influential parameters for the wind barrier shelter.