• Journal of the Society of Naval Architects of Korea
• /
• v.46 no.5
• /
• pp.488-497
• /
• 2009

The freak wave, also known as New-Year-Wave in the north Atlantic, is relatively large and spontaneous ocean surface wave that can sink even large ships and destroy maritime structures. To understand oceanic conditions that develop freak waves, we simulated and generated two versions of scale-downed waves (1:64 and 1:42) in a numerical wave tank and compared the results with the experiment in wave flume. Both of the breaking and non-breaking waves were generated in the simulation. The numerical simulation was implemented based on the finite volume method and a genetic optimization algorithm. Random values were assigned as the initial values for the parameter in the control function, which produced signals representing the motion of wave-maker. The same signal obtained from the optimization process was used for both of the simulation and the experiment. By varying the object function and restrictions of the simulation, a best profile of design wave was selected based on the characteristics, height and period of simulated waves. Results showed that the simulation and experiment with the scale of 1:42 agreed better with freak waves in the natural condition. The presented simulation method will contribute to saving the time and cost for conducting subsequent response analyses of motion under freak waves in the course of the model test for ship and maritime structure.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.6
• /
• pp.1355-1364
• /
• 1994

This paper analyzes the running mechanism of flexible and thin foil above rotating protrusion through a numerical simulation. The scope of analysis is confined to the phenomena of elastohy-drodynamic lubrication between the stationary and rotary drums with a running protrusion and thin foil. This mathematical model is based on the modified Reynolds equation and the equation of plate, considering the geometry of protrusion, running direction of protrusion, and the effect of geometric nonlinearity. Finite element method is adopted as a numerical simulation technique to solve the avobe coupled nonlinear equations. In numerical analysis, the effects of the scanning angle in Reynolds equation and the nonlinear term in plate equation are evaluated. Furthermore, the simulation is applied to the situation that thin foil is located in the entire drums (stationary and rotary drums).

• Journal of Environmental Science International
• /
• v.2 no.2
• /
• pp.103-113
• /
• 1993

The land and sea breeze over the Pusan coastal area is studied by three dimensional mesoscale numerical model. According to the results of the simulation experiments, both Pusan areas and Kimhae areas, the sea breeze began at 0800LST and the strongest at 1500LST and then at 1800LST. After midnight, the sea breeze changed about the land breeze and become weaker than that of the sea breeze in the daytime. Comparisons between calculations and observations showed that the characteristics of diurnal variation and v-component of the wind velocity relatively is similar to the Pusan areas. On the Kimhae areas, however, observations showed time lag which compared to the results of simulation experiments in the velocity of sea breeze and diurnal variation. From the above results, comparisons between calculations and observations is much more similar to the coastal areas than on the inland area.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1059-1069
• /
• 2023

Dam-break flow occurs when an elevated dam suddenly collapses, resulting in the catastrophic release of rapid and uncontrolled impounded water. This study compares laminar and turbulent closure models for simulating three-dimensional dam-break flows using OpenFOAM. The Reynolds-Averaged Navier-Stokes (RANS) model, specifically the k-ε model, is employed to capture turbulent dissipation. Two scenarios are evaluated based on a laboratory experiment and a modified multi-layered block obstacle scenario. Both models effectively represent dam-break flows, with the turbulent closure model reducing oscillations. However, excessive dissipation in turbulent models can underestimate water surface profiles. Improving numerical schemes and grid resolution enhances flow recreation, particularly near structures and during turbulence. Model stability is more significantly influenced by numerical schemes and grid refinement than the use of turbulence closure. The k-ε model's reliance on time-averaging processes poses challenges in representing dam-break profiles with pronounced discontinuities and unsteadiness. While simulating turbulence models requires extensive computational efforts, the performance improvement compared to laminar models is marginal. To achieve better representation, more advanced turbulence models like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are recommended, necessitating small spatial and time scales. This research provides insights into the applicability of different modeling approaches for simulating dam-break flows, emphasizing the importance of accurate representation near structures and during turbulence.

• Investigative Magnetic Resonance Imaging
• /
• v.23 no.1
• /
• pp.38-45
• /
• 2019

Purpose: To demonstrate the high-resolution numerical simulation of the respiration-induced dynamic $B_0$ shift in the head using generalized susceptibility voxel convolution (gSVC). Materials and Methods: Previous dynamic $B_0$ simulation research has been limited to low-resolution numerical models due to the large computational demands of conventional Fourier-based $B_0$ calculation methods. Here, we show that a recently-proposed gSVC method can simulate dynamic $B_0$ maps from a realistic breathing human body model with high spatiotemporal resolution in a time-efficient manner. For a human body model, we used the Extended Cardiac And Torso (XCAT) phantom originally developed for computed tomography. The spatial resolution (voxel size) was kept isotropic and varied from 1 to 10 mm. We calculated $B_0$ maps in the brain of the model at 10 equally spaced points in a respiration cycle and analyzed the spatial gradients of each of them. The results were compared with experimental measurements in the literature. Results: The simulation predicted a maximum temporal variation of the $B_0$ shift in the brain of about 7 Hz at 7T. The magnitudes of the respiration-induced $B_0$ gradient in the x (right/left), y (anterior/posterior), and z (head/feet) directions determined by volumetric linear fitting, were < 0.01 Hz/cm, 0.18 Hz/cm, and 0.26 Hz/cm, respectively. These compared favorably with previous reports. We found that simulation voxel sizes greater than 5 mm can produce unreliable results. Conclusion: We have presented an efficient simulation framework for respiration-induced $B_0$ variation in the head. The method can be used to predict $B_0$ shifts with high spatiotemporal resolution under different breathing conditions and aid in the design of dynamic $B_0$ compensation strategies.

• Proceedings of the KSR Conference
• /
• 2007.05a
• /
• pp.1319-1324
• /
• 2007

Many researches have been performed to analyze the smoke movement in tunnel fires by using field model. Recently, FDS(Fire Dynamics Simulator) v.4, which is one of the field model and developed from NIST(National Institute of Standards and Technology), is widely used. In tunnel fires, FDS can show detail results in local point, but it has difficulties in boundary condition and taking long computing time as the number of grid increases. So, there is a need to use alternative method for tunnel fire simulation. A zone model is different kind of CFD method and solves ordinary differential equation based on conservation and auxiliary equations. It shows good macroscopic view in less computing time compared to field model. In this study, therefore, to confirm the applicability of CFAST in tunnel fire analysis, numerical simulations using CFAST are conducted to analyze smoke movement in longitudinal ventilation reduced-scale tunnel fires. Then the results are compared with experimental results. The differences of temperature and critical velocity between numerical results and experimental data are over $30^{\circ}C$ and 0.9m/s, respectively. These values are out of error range. It shows that CFAST 6.0 is hard to be used for tunnel fire simulation.

• Tunnel and Underground Space
• /
• v.21 no.3
• /
• pp.216-224
• /
• 2011

In this study, Class II behavior of rock failure process under uniaxial and biaxial compression has been numerically simulated using bonded particle model. Class II behavior of rock was simulated by radial strain controlled uniaxial and biaxial compression tests using a suggested method of ISRM. Micro-parameters used in the simulation were determined based on the laboratory uniaxial compression tests carried out at ${\"{A}}sp{\"{o}}$ Hard Rock Laboratory, Sweden. Class II behavior of ${\"{A}}sp{\"{o}}$ rock was effectively simulated using newly proposed numerical technique in this study, and the results of numerical simulations show good similarity with the complete stress-strain curves for Class II behavior obtained from the laboratory tests.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.1
• /
• pp.8-18
• /
• 2015

A numerical simulation for a nonslender BWB UCAV configuration with a rounded leading edge and span of 1.0 m was performed to analyze its aerodynamic characteristics. Numerical results were compared with experimental data obtained at a free stream velocity of 50 m/s and at angles of attack from -4 to $26^{\circ}$. The Reynolds number, based on the mean chord length, is $1.25{\times}106$. 3D multi-block hexahedral grids are used to guarantee good grid quality and to efficiently resolve the boundary layer. Menter's shear stress transport model and two transition models (${\gamma}-Re_{\theta}$ model and ${\gamma}$ model) were used to assess the effect of the laminar/turbulent transition on the flow characteristics. Aerodynamic coefficients, such as drag, lift, and the pitching moment, were compared with experimental data. Drag and lift coefficients of the UCAV were predicted well while the pitching moment coefficient was underpredicted at high angles of attack and influenced strongly by the selected turbulent models. After assessing the pressure distribution, skin friction lines and velocity field around UCAV configuration, it was found that the transition effect should be considered in the prediction of aerodynamic characteristics of vortical flow fields.

• Transactions of the Korean hydrogen and new energy society
• /
• v.22 no.3
• /
• pp.363-371
• /
• 2011

In this paper, a three-dimensional hydrogen desorption model is developed to precisely study the hydrogen desorption kinetics and resultant heat and mass transport phenomena in metal hydride hydrogen storage vessels. The metal hydride hydrogen desorption model, i.e. governed by the conservation of mass, momentum, and thermal energy is first experimentally validated against the temperature evolution data measured on a cylindrical $LaNi_5$ metal hydride vessel. The equilibrium pressure used for hydrogen desorption simulations is derived as a function of H/M atomic ratio and temperature based on the experimental data in the literature. The numerical simulation results agree well with experimental data and the 3D desorption model successfully captures key experimental trends during hydrogen desorption process. Both the simulation and experiment display an initial sharp decrease in the temperature mainly caused by relatively slow heat supply rate from the vessel external wall. On the other hand, the effect of heat supply becomes influential at the latter stages, leading to smooth increase in the vessel temperature in both simulation and experiment. This numerical study provides the fundamental understanding of detailed heat and mass transfer phenomena during hydrogen desorption process and further indicates that efficient design of storage vessel and heating system is critical to achieve fast hydrogen discharging performance.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.13 no.1
• /
• pp.15-18
• /
• 2003

High purity single crystalline calcium fluoride ($CaF_2$) has excellent optical transmission characteristics down to deep UV and is therefore selected as the main optical material for the next generation of lithography apparatus operating at wavelength of 157 nm. The growth of large sized $CaF_2$ single crystals with the required properties for this optical application can be achieved only by optimizing the crystal growth process by the aid of numerical simulation. This needs especially a precise calculation of the heat transport and temperature distribution in the solid and liquid $CaF_2$ under crystal growth conditions. As $CaF_2$ is considered to be semitransparent, the internal radiative heat transfer in $CaF_2$ plays an decisive role in the simulation of the heat transport. On the other hand it is very difficult to obtain quantitative experimental data for evaluating numerical models as $CaF_2$ is extremely corrosive at high temperatures. In this work we present a newly developed experimental technique to perform temperature measurements in $CaF_2$-crystal as well as in the melt under conditions of crystal growth process. These experimental results are compared to calculated temperature data, which were obtained by using different numerical models concerning the internal heat transfer in semitransparent $CaF_2$. It will be shown, that an advanced model, which was developed by the authors, gives a much better agreement with experimental data as a standard model, which was taken from the literature.