• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.741-750
• /
• 1995

The turbulent buoyancy-driven flow in 2-dimensional enclosed cavities heated from the vertical side is numerically calculated for both cases of a Rayleigh number of 5*10$^{10}$ for air and 2.5*10$^{10}$ for water. Three different turbulence models are considered : standard k-.epsilon. model of Ozoe and low-Reynolds-number model of Lam and Bremhorst, and another low-Reynolds-number model of Davidson. The results indicate that the use of low-Reynolds number models is recommended for the indoor airflow computation, and the results from Davidson model are reasonably close to the reported experimental data. A sensitivity study shows that the amounts of wall-heat transfer and the velocity profiles with the Lam and Bremhorst model largely depend on the choice of the wall function for .epsilon..

• Tunnel and Underground Space
• /
• v.22 no.2
• /
• pp.93-105
• /
• 2012

This paper compares the two- and three-dimensional (2D and 3D) approaches for the numerical determination of the equivalent mechanical properties of fractured rock masses. Both orthogonally-fractured model and discrete fracture networks (DFN) were used for the geometry and 2D models were cut in various directions from 3D model to compare their mechanical properties. Geological data were loosely based on the data available from Sellafield, UK. Analytical method based on compliance tensor transformation was used for investigation in orthogonally fractured rock and numerical experiments were conducted on fractured rock mass with DFN geometry. It is shown that 2D approach always overestimates the elastic modulus of fractured rock masses by a factor of up to around two because fractures are assumed to be perpendicular to the model plane in 2D problems. Poisson ratios tend to have larger values in 2D analysis while there is opposite trend in some sections. The study quantitatively demonstrates the limitation of the 2D approach that uses the simplified model from true 3D geometry.

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2617-2629
• /
• 1995

The flow through a centrifugal compressor rotor was calculated using the quasi-3-dimensional and fully 3-dimensional Navier-Stokes solution methods. The calculated results, obtained during the development of the computer codes for both methods are discussed. In the inviscid quasi 3-dimensional analysis, stream function formulation was used for the blade to blade (B-B) plane calculations, and the streamline curvature method was used for the meridional (H-S) plane calculations. In the viscous 3-dimensional flow analysis, a control volume method based on a general rotating curvilinear coordinate system was used to solve the time-averaged Navier-Stokes equations, and a standard k-.epsilon. model was used to obtain eddy viscosity. The quasi-3-dimensional analysis reasonably predicts the pressure distributions and requires much less computation time in the region where viscous effects are not strong; however, it fails to predict velocity field and loss mechanism through the impeller passage. The viscous 3-dimensional flow analysis shows reasonable pressure distributions and typical jet-wake flow field through the impeller passage. Secondary flow and total pressure distributions on cross-sectional planes explain the loss mechanisms through the impeller.

• Journal of the Korea Concrete Institute
• /
• v.13 no.4
• /
• pp.379-388
• /
• 2001

Existing microplane models for concrete use three-dimensional spherical microplanes in the analysis of two-dimensional planar members as well as three-dimensional members. Also, they do not accurately describe the post-cracking behavior of reinforced concrete in tension-compression. In this study, a new microplane model is developed to overcome the disadvantages of the existing models. Instead of the spherical microplanes, the proposed microplane model uses disk microplanes involving a less number of microplanes and two-dimensional stresses and strains. As the result, the proposed model is more effective in numerical calculations. Also, the concept of the strain boundary is introduced to describe accurately the compressive behavior of reinforced concrete with tensile cracks in tension-compression. The validity of the proposed model is verified by comparison with existing experiments. In this paper, the microplane model and the numerical techniques involved in the finite element analysis are described in detail.

• Transactions of the Korean hydrogen and new energy society
• /
• v.26 no.3
• /
• pp.260-270
• /
• 2015

One-dimensional packed bed reactor model accounting for interfacial and intra-particle gradients was developed and based on it numerical analyses were performed to investigate the dynamic behavior of a commercial scale methanation reactor. Methanation reaction was almost complete near the reactor inlet and gases with equilibrated composition were discharged from the reactor. Both the intra-particle temperature gradient and differential surface temperature rise were found to be severe near the reactor inlet. To reduce the possible degradation or fracture of catalyst particles and prevent local overheating on the catalyst, addition of inert material can be an effective way.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.44-53
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed for a 3D component of a nuclear system code and a component-scale analysis tool. A two-fluid three-field model is used for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. An unstructured grid is adopted for realistic simulations of the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been applied to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the modified numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing very well.

• Computers and Concrete
• /
• v.15 no.3
• /
• pp.391-410
• /
• 2015

One of the most common applications of ferrocement is the manufacturing of thin stiffened plates which are prone to buckling. This study focuses on the investigation of the behavior of a ferrocement plate, stiffened in both directions by means of an appropriate grid of ribs. In the present paper detailed three-dimensional numerical Finite Element models are formulated for the simulation of the behavior of the structure under study, which are able to take into account both the geometric and material non-linearities that are present in the subject at hand (plasticity, cracking, large displacements). The difference among the formulated models lies on the use of different types of finite elements. The numerical results obtained by each model are compared and the most efficient model is determined. Finally, this model is in the sequel used for the further investigation of the effect of different parameters on the ultimate load capacity, such as the initial out-of-plane imperfection of the plate and the interaction between the axial loads in both directions.

• Journal of information and communication convergence engineering
• /
• v.7 no.3
• /
• pp.361-365
• /
• 2009

This paper has been presented the transport characteristics of FinFET using the analytical potential model based on the Poisson's equation in subthreshold and threshold region. The threshold voltage is the most important factor of device design since threshold voltage decides ON/OFF of transistor. We have investigated the variations of threshold voltage and drain induced barrier lowing according to the variation of geometry such as the length, width and thickness of channel. The analytical potential model derived from the three dimensional Poisson's equation has been used since the channel electrostatics under threshold and subthreshold region is governed by the Poisson's equation. The appropriate boundary conditions for source/drain and gates has been also used to solve analytically the three dimensional Poisson's equation. Since the model is validated by comparing with the three dimensional numerical simulation, the subthreshold current is derived from this potential model. The threshold voltage is obtained from calculating the front gate bias when the drain current is $10^{-6}A$.

• Journal of applied mathematics & informatics
• /
• v.26 no.3_4
• /
• pp.519-530
• /
• 2008

This paper studies the multi-domain coupled system of one dimensional Arctic temperature field and establishes identification model about the thermodynamic parameters of sea ice (heat storage capacity, density and conductivity) by the so-called output least-square estimate according to the temperature data acquired by a monitor buoy installed in the Arctic ocean. By the optimal control theory, the existence and dependability of weak solution and the identifiability of identification model have been given. Moreover, necessary optimality condition is proposed. Furthermore, the optimal algorithm for the identification model is constructed. By using the optimal thermodynamic parameters of Arctic sea ice, the numerical simulation is implemented, and the numerical results of temperature distribution of Arctic sea ice are demonstrated.