• 한국항공우주학회지
• /
• 제30권8호
• /
• pp.112-119
• /
• 2002

주기적으로 배열된 액체 로켓엔진의 연료 분사기 배열에 대한 3차원 열전달 수치해석을 수행하였다. 해석에는 SIMPLE 알고리즘을 사용한 유한체적법이 적용되었으며 2차원 후향 계단에 대해서 기존의 문헌과 일치하는 열전달 결과를 주었다. 분사기 요소 사이의 거리가 가까울수록 누셀 수와 압력강하가 증가하였다. 레이놀즈 수가 증가할 경우, 누셀 수는 증가 하였으며 무차원 압력강하는 약간 감소하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.10-13
• /
• 2010

The flow field in road tunnel is influenced by some facts such as piston effect of vehicle's move, operation of ventilation facilities, natural wind and buoyancy effect of fire plume. Among those, jet fan is one of main ventilation facilities especially in longitudinal ventilation system of tunnel. In this study to analyze tunnel flow induced by operation of jet fan, numerical simulation has been carried out. The velocity distributions and streamlines in tunnel are examined to consider the three-dimensional characteristics of tunnel flow caused by jet fan.

• 한국자동차공학회논문집
• /
• 제10권4호
• /
• pp.223-233
• /
• 2002

The air distribution duct with multiple outlets is an essential part of automotive air-conditioning system In a bus. The estimation of airflow rate in an automotive air-conditioning duct is typically very complicate due to large variations in cross-sectional area and abrupt changes in flow direction, as well as unbalanced distribution of the flow. In this paper, the flow characteristic in a duct with multiple outlets is investigated through experiment, CFD simulation and a one-dimensional simulation. Numerical simulations have been performed for two simplified air conditioning ducts with multiple outlets used in a medium bus. The three dimensional Navier-Stokes code was used to evaluate the overall pressure, velocity Held, and distribution rate at each diffuser according to the change of various design parameters such as ratio of cross-sectional area and radius of bifurcated region. In addition, a one-dimensional program based on Bernoulli equation was developed to obtain optimized diffuser area required to equalize discharge flow rate at each outlet. As a result of this study, optimized diffuser area of design variable by one-dimensional program was very reasonable as compared to the trend deduced from CFD Simulation. Therefore, the simple and convenient one-dimensional analysis developed in this study can be applied in practical design procedure for air-conditioning duct.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.684-689
• /
• 2001

In this study, three-dimensional viscous flow analysis and optimization are presented for the design of a mixed-flow fan. Steady, imcompressible, three-dimensional Reynolds averaged Navier-Stokes equations are used as governing equations, and standard $k-{\varepsilon}$ turbulence model is chosen as a turbulence model. Governimg equations are discretized using finite volume method. Upwind difference scheme is used for the discretization of the convective term and SIMPLEC algorithm is used as a velocity-pressure correction procedure. The computational results are compared with the results obtained by TASCflow. For the numerical optimization of the design, objective function is defined as a ratio of generation of the turbulent energy to pressure head. Sweep angles are used as design variables.

• 대한기계학회논문집A
• /
• 제22권1호
• /
• pp.121-131
• /
• 1998

This paper applies the FE analysis procedure, developed in the Part I of the companion article, to the three-dimensional rubber pad deformation during rubber-pad forming process. Effects of different algorithms corresponding to incompressibility constraint and time integration methods on numerical solution responses are investigated. Laboratory scale experiments support the validity of the developed FE procedure an demonstrate the accuracy of the numerical models. Full scale model responses are also predicted using the reasonable method and parameters obtained in laboratory modeling.

• 대한기계학회논문집A
• /
• 제22권1호
• /
• pp.111-120
• /
• 1998

This paper is the first one of two-parted research efforts focusing on the modeling of rubber pad forming process. The rubber pad, driven by the pressurized fluid during the forming process, pushes the sheet metal to solid tool half and forms a part to final shape. In this part of the paper, a numerical procedure for the FE analysis of the rubber pad deformation is presented. The developed three-dimensional FE model is based on the total Lagrangian description of rubber maerial characterized by nearly incompressible hyper-elastic behavior under a large deformation assumption. Validity of the model as well as effects of different algorithms corresponding to incompresibility constraints and time integration methods on numerical solution responses are also demonstrated.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1995년도 추계 학술대회논문집
• /
• pp.175-181
• /
• 1995

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• 동력기계공학회지
• /
• 제14권2호
• /
• pp.48-54
• /
• 2010

An improved numerical scheme, to which the hydroelastic method is adapted, is introduced for predicting the motion and structural responses of tension leg platforms(TLPs) in regular waves. The numerical approach in this work is based on a combination of the three dimensional source distribution method and the finite element method. The hydrodynamic interactions among TLP members, such as columns and pontoons, are included in the motion and structural response analysis. The drag forces on the submerged slender members, which are proportional to the square of relative velocity, are included in order to estimate the responses of members with better accuracy. Comparisons with other results verify the works in this paper.

• Steel and Composite Structures
• /
• 제24권5호
• /
• pp.523-536
• /
• 2017

A layerwise (LW) formulation based on the Galerkin method is presented to investigate the three-dimensional stress state in long sandwich plate which is subjected to tension force and pure bending moment. Based on the Galerkin method and the LW discretization approach, the equilibrium equations of elasticity for the long plate are written in the weak form and discretized through the thickness of the plate. The discretized equations are written in terms of displacement components of the numerical layers. The governing equations of the plate are solved analytically for the free edge boundary conditions. The distribution of stress state especially the 3D stress state in the vicinity of the edges of the sandwich plate which is subjected to tension and pure bending is studied. In order to increase the accuracy, the out of plane stresses are obtained by integrating the equilibrium equations of elasticity. The convergence and accuracy of the predictions are studied and various numerical results are presented for distribution of the in-plane and out of plane stresses in symmetric and un-symmetric sandwich plates.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2001년도 봄 학술발표회 논문집
• /
• pp.97-104
• /
• 2001

The magnitude and the orientation of in-situ stresses contribute to ground displacement and stresses in the field of underground space. This paper investigates in-situ stresses at various depth on the basis of 392 data which were determined by over-coring and hydro-fracturing test methods in the Korea peninsula. The result shows that in-situ stress distribution are more or less non-uniform through the Granite and Gneiss sub-area, and that the K-value in the Volcanic sub-area are below 1 at the deep depth. Also, the result of three dimensional numerical analyses of tunnel shows that the direction and magnitude of displacement around tunnel are much effected by the stress difference between the maximum and the minimum horizontal stress.