• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.6
• /
• pp.436-449
• /
• 2013

In this paper, a numerical analysis is carried out to study the characteristics of supercavitating flows and the drag of relatively simple two-dimensional and axisymmetric bodies which can be used for supercavity generation device, cavitator, of a high-speed underwater vehicle. In order to investigate the suitability of numerical models, cavity flows around the hemispherical head form and two-dimensional wedge are calculated with combinations of three turbulence models(standard $k-{\epsilon}$, realizable $k-{\epsilon}$, Reynolds stress) and two cavitation models(Schnerr-Sauer, Zwart-Gerber-Belamri). From the results, it is confirmed that the calculated cavity flow is more affected by the turbulence model than the cavitation model. For the calculation of steady state cavity flows, the convergence in case of the realizable $k-{\epsilon}$ model is better than the other turbulence models. The numerical result of the Schnerr-Sauer cavitation model is changed less by turbulence model and more robust than the Zwart-Gerber-Belamri model. Thus the realizable $k-{\epsilon}$ turbulence model and the Schnerr-Sauer cavitation model are applied to calculate supercavitating flows around disks, two dimensional $10^{\circ}$ and $30^{\circ}$ wedges. In case of the disk, the cavitation number dependences of the cavity size and the drag coefficient predicted are similar to either experimental data or Reichardt's semi-empirical equations, but the drag coefficient is overestimated about 3% higher than the Reichardt's equation. In case of the wedges, the cavitation number dependences of the cavity size are similar to experimental data and Newman's linear theory, and the agreement of the cavity length predicted and Newman's linear theory becomes better as decreasing cavitation number. However, the drag coefficients of wedges agree more with experimental data than those of Newman's analytic solution. The cavitation number dependences of the drag coefficients of both the disk and the wedge appear linear and simple formula for estimating the drag of supercavitating disks and wedges are suggested. Consequently, the CFD scheme of this study can be applied for numerical analysis of supercavitating flows of the cavitator and the cavitator design.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.3
• /
• pp.180-201
• /
• 2020

It has been widely known that the effect of diffracted waves at the tip of composite breakwater with finite length causes the change of standing wave height along the length of breakwater, the spatial change of wave pressure on caisson, and the occurrence of meandering damage on the different sliding distance in sequence. It is hard to deal with the spatial change of wave force on trunk of breakwater through the two-dimensional experiment and/or numerical analysis. In this study, two and three-dimensional numerical techniques with olaFlow model are used to approach the spatial change of wave force including the impulsive breaking wave pressure applied to trunk of breakwater, the effect of rear region, and the occurrence of diffracted waves at the tip of caisson located on the high crested rubble mound. In addition, it is thoroughly studied the mean wave height, mean horizontal velocity, and mean turbulent kinetic energy through the numerical analysis. In conclusion, it is confirmed that the larger wave pressure occurs at the front wall of caisson around the still water level than the original design conditions when it generates the shock-crushing wave pressure checked by not two-dimensional analysis, but three-dimensional analysis through the change of wave pressure applied to the caisson along the length of breakwater.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.387-390
• /
• 2007

The mixing characteristics of heated water discharged from different types of effluents are simulated using a commercial software, Flow-3D model. In order to verify the model, the comparison of numerical results with the experimental data are conducted for each type of effluent, the submerged and surface overfall. It is observed that the numerical results show a reasonable agreement with the experimental data. Based on this study the application to the heated water discharge problem in the field can be expected as a further study and it can be the fundamental data when determining the type of effluent in a powerplant.

• Geotechnical Engineering
• /
• v.12 no.6
• /
• pp.127-138
• /
• 1996

The in-situ consolidation behavior of drainage system-installed deposits has three dimensional characteristics. Therefore, for an approximate 2-D plane strain consolidation analysis, it is necessary to convert the 3-D spatial flow of actual cases into the laminar flow simulated by the 2-D plane strain model. . In this paper, in order to properly model the effect of three dimensional characteristics, an equivalent and efficient model has been applied in a finite element technique for the analysis of the drainage system-installed soil deposits. The equivalent two dimensional model involves equivalent permeabilities and drainage widths. To validate the equivalent two dimensional model, three dimensional analyses were per formed by using the ABAQUS program and the results of 3-D analyses were compared with those of the 2-D analyses. By using the proposed equivalent model, one may be able to appropriately predict the consolidation behavior of drainage system-installed soft deposits.

• Journal of Biosystems Engineering
• /
• v.22 no.4
• /
• pp.421-432
• /
• 1997

This study was conducted to investigate dynamic characteristics of agricultural tractor with a particular interest in ride vibrations when it is subjected to various excitation forces. As the first part of it this paper describes development of dynamic model of a tractor-trailer system and its equations of motions. An 3 dimensional 16-degree-of-freedom dynamic model for a tractor-trailer system was developed and its equations of motions were derived, which will be used to investigate the effects of irregular ground surface and excitation forces due to the engine mounted on the tractor. And the excitation forces were also formulated analytically. The transition matrix method and QR algorithm were proposed for numerical solution of the equation of motions fur the developed model. The later parts of the study will include a proof of the model and optimization from which tractors can be designed to minimize the ride vibrations. This will be presented in the second and third papers to be followed shortly.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.8
• /
• pp.855-861
• /
• 2011

We predict the rolled-strip profile for a 6-high mill using thin rolling theory and a numerical model. In the numerical model, we calculate the distributions of the contact pressures between the rolls and the rolling pressure between the strip and the work roll in the transverse direction using the geometric structure of the 6-high mill and the boundary conditions. We determine the distribution of the rolling pressure in the rolling direction via a thin-foil rolling model using Fleck's theory. We calculate the three-dimensional elastic deformation of the work roll using the pressures of the width and rolling directions. We then obtain the three-dimensional strip profile via the elastic deformation of the work roll during the rolling process. The profile is verified by a thin cold-rolling test and FE simulation.

• Architectural research
• /
• v.11 no.1
• /
• pp.15-23
• /
• 2009

An application of the constrained adaptive topology optimization (CATO) algorithm is described for three-dimensional topology optimization of engineering structures. The enhanced assumed strain lower order solid finite element (FE) is used to evaluate the values of objective and constraint functions required in optimization process. The strain energy (SE) terms such as elastic and modal SEs are employed as the objective function to be minimized and the initial volume of structures is introduced as the constraint function. The SIMP model is adopted to facilitate the material redistribution and also to produce clearer and more distinct structural topologies. The linearly weighted objective function is introduced to consider both static and dynamic characteristics of structures. Several numerical tests are tackled and it is used to investigate the performance of the proposed three-dimensional topology optimization process. From numerical results, it is found to be that the CATO algorithm is easy to implement and extremely applicable to produce the reasonable optimum topologies for three dimensional optimization problems.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.1B
• /
• pp.115-123
• /
• 2008

In the present study, a three dimensional hydrodynamic characteristics around the fully submerged dual breakwaters with a sand beach has been investigated numerically using a 3-D numerical scheme, which can determine the eddy viscosity with LES turbulence model and is able to consider wave-structure-seabed interaction in 3-dimensional wave field (LES-WASS-3D), recently developed by Hur and Lee (2007). Based on the numerical experiments, strong vortices can be generated fore and aft edges of the structures, and propagate lee sides. Thus relatively large circulation flows are occurred around submerged breakwaters. The 3-D flow hydrodynamic characteristics have been examined by mean flows and mean vortices for various x-y, x-z sections and y-z layers. Wave height distribution and wave set-up around and over submerged breakwaters, and breaking point migration toward shore side is discussed in detail.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.243-248
• /
• 2007

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

• Computers and Concrete
• /
• v.8 no.3
• /
• pp.327-341
• /
• 2011

Three-dimensional graphic objects created by MATLAB are exported to the AUTOCAD program through the MATLAB handle functions. The imported SAT format files are used to produce the finite element mesh for MSC.PATRAN. Based on the Monte-Carlo random sample principle, the material heterogeneity of cement composites with randomly distributed fibers is described by the WEIBULL distribution function. In this paper, a concept called "soft region" including micro-defects, micro-voids, etc. is put forward for the simulation of crack propagation in fiber-reinforced cement composites. The performance of the numerical model is demonstrated by several examples involving crack initiation and growth in the composites under three-dimensional stress conditions: tensile loading; compressive loading and crack growth along a bimaterial interface.