• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.119-127
• /
• 1993

Recently, ADI scheme has been a most common tool for solving shallow water equation numerically. But ADI models of tidal flow is likely to cause so called ADI effect in such a region of the Yellow Sea which shows complex topography and has submarine canyons especially. To overcome this, a finite difference algorithm is developed which adopts fully implicit method and preconditioned conjugate gradient squared method. Applying the algorithm including simulation of intertidal zone to Sae-Man-Keum. velocity fields and flooding/drying phenomena are simulated well in spite of complex topography.

• Nuclear Engineering and Technology
• /
• v.40 no.6
• /
• pp.477-488
• /
• 2008

The objective of the paper is to analyze the thermally induced density wave oscillations in water cooled boiling water reactors. A transient thermal hydraulic model is developed with a characteristics-based implicit finite-difference scheme to solve the nonlinear mass, momentum and energy conservation equations in a time-domain. A two-phase flow was simulated with a one-dimensional homogeneous equilibrium model. The model treats the boundary conditions naturally and takes into account the compressibility effect of the two-phase flow. The axial variation of the heat flux profile can also be handled with the model. Unlike the method of characteristics analysis, the present numerical model is computationally inexpensive in terms of time and works in a Eulerian coordinate system without the loss of accuracy. The model was validated against available benchmarks. The model was extended for the purpose of studying the flow-induced density wave oscillations in forced circulation and natural circulation boiling water reactors. Various parametric studies were undertaken to evaluate the model's performance under different operating conditions. Marginal stability boundaries were drawn for type-I and type-II instabilities in a dimensionless parameter space. The significance of adiabatic riser sections in different boiling reactors was analyzed in detail. The effect of the axial heat flux profile was also investigated for different boiling reactors.

• Solar Energy
• /
• v.5 no.2
• /
• pp.77-83
• /
• 1985

An analysis is performed to investigate the influence of the buoyancy force and the thickness variation of melting layer in the containment that is filled with phase-change Material surrounding a cylindrical heating tube during melting process. The phase-change material is assumed to be initially solid at its phase-change temperature and the remaining solid at any given time is still at the phase-change temperature and neglecting the effect of heat transfer occuring within the solid. At the start of melting process, the thickness of melting layer is assumed to be a stefan-problem and after the starting process, the change of temperature and velocity is calculated using a two dimensional finite difference method. The governing equations for velocity and temperature are solved by a finite difference method which used SIMPLE (Semi Implicit Method Pressure linked Equations) algorithm. Results are presented for a wide range of Granshof number and in accordance with the time increment and it is founded that two dimensional fluid flow occurred by natural convection decreases the velocity of melting process at the bottom of container. The larger the radius of heating tube, the higher heat transfer is occurred in the melting layer.

• Steel and Composite Structures
• /
• v.49 no.5
• /
• pp.587-602
• /
• 2023

This article presents the numerical modelling of transient heat transfer in highly heterogeneous composite materials where the thermal conductivity, specific heat and density are assumed to be directional-dependent. This article uses a coupled finite element-finite difference scheme to perform the transient heat transfer analysis of unidirectional (1D) and multidirectional (2D/3D) functionally graded composite panels. Here, 1D/2D/3D functionally graded structures are subjected to nonuniform heat source and inhomogeneous boundary conditions. Here, the multidirectional functionally graded materials are modelled by varying material properties in individual or in-combination of spatial directions. Here, fully spatial-dependent material properties are evaluated using Voigt's micromechanics scheme via multivariable power-law functions. The weak form is obtained through the Galerkin method and solved further via the element-space and time-step discretisation through the 2D-isoparametric finite element and the implicit backward finite difference schemes, respectively. The present model is verified by comparing it with the previously reported results and the commercially available finite element tool. The numerous illustrations confirm the significance of boundary conditions and material heterogeneity on the transient temperature responses of 1D/2D/3D functionally graded panels.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.51 no.11
• /
• pp.559-566
• /
• 2002

In this paper, we present an accurate and stable method for the solution of the transient electromagnetic response from the conducting wire structures using the time domain integral equation. By using an implicit scheme with the central finite difference approximation for the time domain electric field integral equation, we obtain the transient response from a wire scatterer illuminated by a plane wave and a conducting wire antenna with an impressed voltage source. Also, we consider a wire above a 3-dimensional conducting object. Numerical results are presented, which show the validity of the presented methodology, and compared with a conventional method using backward finite difference approximation.

• Solar Energy
• /
• v.2 no.1
• /
• pp.1-8
• /
• 1982

A theoretical study by numerical method has been performed on the natural convection of an air contained in enclosures. The enclosures have rectangular cross section with one vertical wall heated and the other cooled, and with two horizontal partition plates of finite thermal conductivity. Steady two-dimensional flow was assumed. The computation was executed by means of the Implicit Alternating Direction (I.A.D) finite-difference method. Two partition plates of Aluminium whose thickness were 0.05mm was employed in computation. Isothemals, streamlines, local Nusselt numbers and mean Nusselt numbers were obtained for various Grashof numbers and aspect ratio and these results were compared with those in the case of the enclosure with two horizontal insulated plates. From the present results, the heat transfer in the case of partition plates was greater than that in the case of insulation. This study suggests a method to measure the overall heat-transfer of coefficient in double walls which supported by partition plates for insulative construction.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1243-1251
• /
• 1994

The implicit finite difference program was developed to evaluate the relationship between time and consolidation ratio within the zone of vertical drain effective radius. In the evaluation, the excess pore water pressure was considered to dissipate in two directions, namely, vertical and radial flow direction. To calculate subsoil stress increments in the soil due to multi-step embanking, the foundation soil was assumed to be an isotropic and homogeneous elastic medium and the initial excess pore water pressure was estimated by using Skempton's parameters whose condition is plane strain and elastic phase of pore pressure response within the soft ground. Regarding to the settlement estimation, immediate and primary consolidation settlements were calculated. The secondary or delayed consolidation settlement was not considered. Numerically calculated excess pore water pressure and settlements were similar to the measured data in situ. Thus, this method can be used to predict the time-consolidation ratio of each layer treated by vertical drain method.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.155-174
• /
• 2006

Parabolized stability equations for compressible flows in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Compressible and incompressible flat plate flow stability under two-dimensional and three¬dimensional disturbances has been investigated to test the present code. Results of the present computation are found to be in good agreement with the multiple scale analysis and DNS data. Stability calculation results by the present PSE code for compressible boundary layer at Mach numbers ranging from 0.02 to 1.5 are also presented and are again seen to be as accurate as the spectral method.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.35-46
• /
• 1996

원심형 임펠러 내부 유로등 큰 곡률을 수반하는 터보기계 요소의 유동해석을 위한 계산코드를 개발하였다. 이 코드에서는 곡선좌표계에 유도된 3차원 비압축성 Navier-Stokes의 운동 방정식을 SMAC 음해법으로 푼다. 이 코드를 이용하여 유로의 단면이 정사각형이고 90도로 굽은 덕트내부의 층류 입구유동을 해석하고, 굽은 관 특유의 유동현상을 수치모사하였다 또한 곡관부 입구에서 충분히 발달한 유동, 또는 발달중인 유동이 유입될 경우에 이것이 곡관부 내부의 유동에 미치는 영향을 상·하류의 계산영역이 서로 다를 몇몇 유동장에 대하여 조사하고, 본 계산에서 얻어진 결과와 실형결과와의 비교로 본 3차원 유동해석 코드의 유효성을 검토 하였다.