• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.121-122
• /
• 2003

In this study a microchannel with various arrangement of blocks is newly proposed. This design comprises periodically arranged simple blocks. In this configuration, the stirring is greatly enhanced at a certain geometric parameter set. To characterize the flow field and the stirring effect both the numerical and experimental methods were employed. To obtain the velocity field, three-dimensional numerical computation to the Navier Stokes equations are performed by using a commercial code, FLUENT 6.0. The fluid-flow solutions are then cast into studying the characteristics of stirring with the aid of Lyapunov exponent. The numerical results show that the particles' trajectories in the microchannel heavily depend on the block arrangement. It was shown that the stirring is significantly enhanced at larger block-height and it reaches maximum when the height is 0.8 times the channel width. We also studied the effect of the block stagger angle, and it turns out that the stirring performance is the best at the block angel ${45^\circ}$.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.28-38
• /
• 2001

A numerical method for the assessment and correction of tunnel wall interference effects on forced-oscillation testing is presented. The method is based on the wall pressure signature method using computed wall pressure distributions. The wall pressure field is computed using unsteady three-dimensional full Navier-Stokes solver for a 70-degree pitching delta wing in a wind tunnel. Approximately-factorized alternate direction implicit (AF-ADI) scheme is advanced in time by solving block tri-diagonal matrices. The algebraic Baldwin-Lomax turbulence, model is included to simulate the turbulent flow effect. Also, dual time sub-iteration with, local, time stepping is implemented to improve the convergence. The computed wall pressure field is then imposed as boundary conditions for Euler re-simulation to obtain the interference flow field. The static computation shows good agreement with experiments. The dynamic computation demonstrates reasonable physical phenomena with a good convergence history. The effects of the tunnel wall in upwash and blockage are analyzed using the computed interference flow field for several reduced frequencies and amplitudes. The corrected results by pressure signature method agree well with the results of free air conditions.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.825-830
• /
• 2004

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc (DCA) compressor cascades. Two types of double-circular-arc cascades were used in this analysis. The appropriate turbulence model for compressor analysis was selected among the conventional turbulence models such as Baldwin-Lomax, k-$\varepsilon$ and k-$\varepsilon$ models. The results of current study were compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE/PWIM algorithm for collocated grid and hybrid scheme for the convective terms were the main features of numerical tools. As commonly known, turbulence modeling is very important for the prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. For selection of turbulence model, 2-D analysis was performed. And then, k-$\varepsilon$ turbulence model with wall function chosen as the reasonable turbulence model for 3-D calculation was used to increase the efficiency of computation times. A reasonable result of 3-D flow pattern passing through the double-circular-arc cascade was obtained.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.1
• /
• pp.45-58
• /
• 1992

A three-dimensional hydrodynamic numerical model is used to compute the annual and seasonal meteorologically-induced residual circulation on the Yellow Sea and the East China Sea continental shelf. The model is formulated having irregular coastal boundaires and non-uniform depth distribution representative of nature. The previous three-dimensional model of the East China Sea (Choi. 19U) has been further refined to resolve the flow over the continental shelf in more detail. The mesh resolution of the present finite difference grid system used is 4 minutes latitude by 5 minutes longitude over the entire shelf. The circulation pattern showing depth and spatial distribution of currents over the Yellow Sea and the East China Sea is presented. Meteorologically-induced currents are subsequently used to compute turn-over times for the three depths (surface. mid-depth. bottom) and the total water column of various regions of the Yellow Sea and the East China Sea.

• Economic and Environmental Geology
• /
• v.27 no.1
• /
• pp.41-47
• /
• 1994

A fast Hankel transform (FHT) algorithm (Anderson, 1982) is applied to numerical evaluation of many Green's tensor integrals encountered in three-dimensional electromagnetic modeling using integral equations. Efficient computation of Hankel transforms is obtained by a combination of related and lagged convolutions which are available in the FHT. We express Green's tensor integrals for a layered half-space, and rewrite those to a form of related functions so that the FHT can be applied in an efficient manner. By use of the FHT, a complete or full matrix of the related Hankel transform can be rapidly and accurately calculated for about the same computation time as would be required for a single direct convolution. Computing time for a five-layer half-space shows that the FHT is about 117 and 4 times faster than conventional direct and multiple lagged convolution methods, respectively.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2005.06a
• /
• pp.812-815
• /
• 2005

The incremental forming process employs several tens or hundreds of continuous local strokes, so the entire process is difficult to analyze due to much computation time and large computer memory. The objective of this work is to propose a new numerical scheme of the finite element method, automatic expansion of domain (AED), and to reduce computation time and computer memory. In the AED scheme, an effective analysis domain in each local forming step is defined and then the domain is automatically expanded in accordance with the repeated process. In order to verify the validity of the criterion for the AED scheme and the applicability of the AED scheme, two-dimensional incremental plane-strain forging process is first analyzed using the proposed scheme with various criteria and full domain. In addition, three-dimensional incremental radial forging process is analyzed to verify the applicability of the proposed scheme to a practical incremental forging process.

• Water for future
• /
• v.16 no.4
• /
• pp.221-236
• /
• 1983

In semi-enclosed shallow sea areas typified by the Yellow sea and the East China Sea, currents and sea surface variations are predominantly tidal. During the recent years two-dimensional numerical hydrodynamic model of the Yellow Sea and the East China Sea has been developed, based on the vertically-integrated equations of motion and continuity, capable of reproducing amplitudes and phases of the principal components of tides to satisfiable accuracy. As a subsequent development a three-dimensional hydrodynamical nymerical model covering the Yellow Sea and the East China Sea has been formulated to investigate the vertical distribution of horizontal tidal current and the response of the continented to investigate the vertical distribution of horizontal tidal current and the response of the continental shelf sea to steady uniform wind stress field imposed over the surface. Features of the M2 tidal current and the wind-induced three-dimensional current structure determined from the computation have been examined and discussed.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2005.04a
• /
• pp.327-334
• /
• 2005

In this paper, using an adjoint variable method, we develop a design sensitivity analysis (DSA) method applicable to 3-Dimensional heat conduction problems in steady state. Also, a topology design optimization method is developed using the developed DSA method. Design sensitivity expressions with respect to the thermal conductivity are derived. Since the already factorized system matrix is utilized to obtain the adjoint solution, the cost for the sensitivity computation is trivial. For the topology design optimization, the design variables are parameterized into normalized bulk material densities. The objective function and constraint are the thermal compliance of structures and allowable material volume, respectively, Through several numerical examples, the developed DSA method is verified to yield efficiency and accurate sensitivity results compared with finite difference ones. Also, the topology optimization yields physical meaningful results.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.3
• /
• pp.114-119
• /
• 2001

Goda formula (Goda, 1973) has been used in the determination of wave pressures acting on a large size caisson such as the pier of the cable stayed bridge at sea. Goda formula, however, is to evaluate the wave pressures acting the infinite vertical caisson of composite breakwater so that it can`t be applied to a large caisson with finite width and length because of diffraction effects. In the present study, three dimensional nonlinear frequence domain method based on perturbation method and boundary integral method is applied to the computation of the linear and nonlinear wave pressures acting on the front of a large size caisson under the variation of its width and length, and angle of incident wave. The numerical results are compared to Goda\`s ones, and then the characteristics of wave pressure distributions acting on a large size caisson are discussed.