• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.4
• /
• pp.76-82
• /
• 2013

Motion reduction of an offshore structure at resonant frequency is essential for avoiding critical damage to the topside and mooring system. A damping plate has a distinct advantage in reducing the motion of a floating structure by increasing the added mass and the damping coefficient. In this study, the heave motion responses of a circular cylinder with an impermeable and a permeable damping plate attached at the bottom of the cylinder were investigated thru a model test. The viscous damping coefficients for various combinations of porosity were obtained from a free-decay test by determining the ratio between any pair of successive amplitudes. Maximum energy dissipation occurred at a porous plate with a porosity P = 0.1008. Experimental results for regular and irregular waves were compared with an analytical solution by Cho (2011). The measured heave RAO and spectrum reasonably followed the trends of the predicted values. A significant motion reduction at resonant frequency was pronounced and the heaving-motion energy calculated by the integration of the area under the heave motion spectrum was reduced by more than 75% by the damping plate. However, additional energy dissipation by eddies of strong vorticity and flow separation inside a porous damping plate was not found in the present experiments.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.6
• /
• pp.374-385
• /
• 2013

Sediment particle motions have been numerically simulated over a sinusoidal ripple. Turbulent boundary layer flows are generated by Large Eddy Simulation, and the sediment particle motions are simulated using Lagrangian particle tracking method. Two unsteady flow conditions are used in the experiment by employing two different wave amplitudes while keeping other conditions such as wave period same. As expected, the amount of suspended sediment particles is clearly dependent on the wave amplitude as it is increasing with increasing flow intensity. However, it is also observed that the pattern of suspension may be different as well due to the only different condition caused by wave amplitude. Specially, the time of maximum sediment suspension within the wave period is not coincident between the two cases because sediment suspension is strongly affected by the existence of turbulent eddies that are formed at different times over the ripple between the two cases as well. The role of these turbulent eddies on sediment suspension is important as it is also confirmed in previous researches. However, it is also found the time of these eddies' formation may also dependent on the wave amplitude over rippled beds. Therefore, it has been proved that various flow as well as geometric conditions under waves has to be considered in order to have better understanding on the sediment suspension process over ripples. In addition, it is found that high turbulent energy and strong upward flow velocities occur during the time of eddy formation, which also supports high suspension rate at these time steps. The results indicate that the relationship between the structure of flows and bedforms has to be carefully examined in studying sediment suspension at coastal regions.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.6
• /
• pp.732-742
• /
• 1985

Measurements of mean velocity and turbulence characteristics are obtained with a linearized constant temperature hot-wire anemometer in a two-dimensional turbulent jet discharging parallel to a flate. Wall static pressure distribution is also measure. The Reynolds number based on the jet nozzle width (D) is about 42,000 and the step height is 2.5D. The reattachment length is found to be 7.5D by using both wool tuft and oil methods. Upstream of the reattachment point, there exist double coherent structures and mean velocity, Reynolds stresses and triple product profiles are asymmetric about jet center line due to the influence of streamline curvature and recirculating flow region. Near the reattachment point, wall static pressure and turbulence quantities change its shape rapidly because of the large eddies by the solid wall. Especially, turbulence intensity has a maximum value in the reattachment regin, then decreases slowly in the redeveloping wall jet ragion. Downstream of X/D=14, a single large scale eddy structure is formed. Far downstream affer the reattachment(X/D.geq.18) mean velocity profile, the decay of maximum velocity and the variation of jet half width are nearly similar to those of plane wall jet, but the Reynolds stresses are higher than those of the latter.

• Ocean and Polar Research
• /
• v.23 no.2
• /
• pp.141-160
• /
• 2001

The challenge of predicting the Japanese coastal ocean motivated Frontier Observational Research System for Global Change (FORSGC) and the Japan Marine Science and Technology Center (JAMSTEC) to start a multiyear observational programme in the upstream Kuroshio in November 2000. This field effort, the Kuroshio Observation Program (KOP), should enable us to determine the barotropic and baroclinic components of the western boundary current system, thus, to better understand interactions of the currents with mesoscale eddies, the Kuroshio instabilities, and path bimodality. We, then, will be able to improve modeling predictability of the mesoscale, seasonal, and inter-annual processes in the midstream Kuroshio near the Japanese main islands by using this knowledge. The KOP is focused on an enhanced regional coverage of the sea surface height variability and the baroclinic structure of the mainstream Kuroshio in the East China Sea, the Ryukyu Current east of the Ryukyu's, and the Kuroshio recirculation. An attractive approach of the KOP is a development of a new data acquisition system via acoustic telemetry of the observational data. The monitoring system will provide observations for assimilation into extensive numerical models of the ocean circulation, targeting the real-time monitoring of the Japanese coastal waters.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.5 no.3
• /
• pp.186-194
• /
• 2000

In this study, the oceanic responses to given atmospheric boundary conditions are investigated using a mesoscale ocean circulation model. The numerical experiments are divided into two parts: One is, so called, spin-up experiment and the other is reproduction experiment. The spin-up experiment simulates climatic state of ocean by integrating the ocean model with upper boundary conditions of the monthly mean atmospheric climate data. In the reproduction experiment, for the reproduction of major oceanic changes around Korean Peninsula during the period of 1980-1998 (19 years), the model has been integrated under the boundary condition of the 19year monthly mean atmosphere data. The spined-up state of ocean generated from the spin-up experiment is assigned to the initial boundary condition of the reproduction experiment. In the spin-up experiment, the model properly simulates the major features of circulation structure around Korean Peninsula; such as separation of East Korean Warm Current (EKWC), formation of the polar front, cold water band associated with the small scale eddies in the East Sea, the formation of front along west coast, and the seasonal variation of circulation pattern caused by changing upwind current in the West Sea. In the reproduction experiment, the model has shown the interannual sea surface temperature variations and a warming trend of about 0.5$^{\circ}$C during the period around Korean Peninsula, as in the case of the observation. Therefore, it is concluded that the model is capable of simulating not only the mean states but also the variabilities of ocean under the given atmosphere boundary conditions.

• Fisheries and Aquatic Sciences
• /
• v.13 no.2
• /
• pp.167-181
• /
• 2010

To predict changes in the marine environment of the Beolgyo Stream Estuary in Jeonnam Province, South Korea, where cohesive tidal flats cover a broad area and a large bridge is under construction, this study conducted numerical simulations involving tidal flow and cohesive sediment transport. A wetting and drying (WAD) technique for tidal flats from the Princeton Ocean Model (POM) was applied to a large-scale-grid hydrodynamic module capable of evaluating the flow resistance of structures. Derivation of the eddy viscosity coefficient for wakes created by structures was accomplished through the explicit use of shear velocity and Chezy's average velocity. Furthermore, various field observations, including of tide, tidal flow, suspended sediment concentrations, bottom sediments, and water depth, were performed to verify the model and obtain input data for it. In particular, geologic parameters related to the evaluation of settling velocity and critical shear stresses for erosion and deposition were observed, and numerical tests for the representation of suspended sediment concentrations were performed to determine proper values for the empirical coefficients in the sediment transport module. According to the simulation results, the velocity variation was particularly prominent around the piers in the tidal channel. Erosion occurred mainly along the tidal channels near the piers, where bridge structures reduced the flow cross section, creating strong flow. In contrast, in the rear area of the structure, where the flow was relatively weak due to the formation of eddies, deposition and moderated erosion were predicted. In estuaries and coastal waters, changes in the flow environment caused by artificial structures can produce changes in the sedimentary environment, which in turn can affect the local marine ecosystem. The numerical model proposed in this study will enable systematic prediction of changes to flow and sedimentary environments caused by the construction of artificial structures.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.8 no.2
• /
• pp.94-110
• /
• 2003

This study utilizes the dataset of Topex/Poseidon(T/P) altimeter sea surface height (1992-2000 yr., 286 cycles)to investigate the tempore-spatial variability in the East (Japan) Sea. Optimal interpolation (Ol) technique was applied to the pre-processed T/P dataset (level 2) to produce sea surface height anomaly (SSHA) map on regular grids. Spectral analyses of the timeseries of the SSHA at chosen stations and empirical orthogonal function (EOF) analysis of the SSHA in the entire East Sea were made. Distribution of the SSHA can be divided by the southern and northern regions sharply by the polar front situated in the middle of the East Sea. The southern region under the direct influence of the Tsushima Current exhibits higher amplitude of the SSHA fluctuation, while the northern region does relatively smaller one. The spatio-temporal variability of the SSHA in the East Sea can be characterized by the five modes of the EOFs accounting for more than 85% of the total variance. The first mode dominates the SSHA variation in the entire domain with strong seasonal and inter-annual periods accounting for the 72.3% of the total variance. The other modes (up to 5th account for 14%) are responsible for the SSHA variation associated with the local current system, meandering of the polar frontal axis, and mesoscale eddies. Spectral peaks with significant confluence level show semi-annual, annual and interannual (2, 3-4 years) periods.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.24 no.2
• /
• pp.298-317
• /
• 2019

The physical characteristics of the Ulleung Warm Eddy (UWE) and its relationship with the East Korea Warm Current (EKWC) were analyzed using the CMEMS (Copernicus Marine Environment Monitoring Service) satellite altimetry data and the CTD data of the National Institute of Fisheries Science (NIFS) near the Ulleung Basin from 1993 to 2017. The distribution of the UWEs coupled with EKWC accounts for 81% of the total number of the UWEs. Only 7% of the total eddies are completely separated from the EKWC. The UWE has the characteristics of high temperature and high salinity water inside of it when it is formed from the EKWC. However, when the UWE is wintering, its internal structure changes greatly. In the winter, surface homogeneous layer of $10^{\circ}C$ and 34.2 psu inside of the UWE is produced by vertical convection from sea-surface cooling, and deepened to a maximum depth of approximately 250 m in early spring. In summer, the UWE changes into a structure with a stratified structure in the upper layer within a depth of 100 m and a homogeneous layer made in winter in the lower layer. 62 UWEs were produced for 25 years from 1993 to 2017. on average, 2.5 UWEs were formed annually, and the average life span was 259 days (approximately 8.6 months). The average size of the UWEs is 98 km in the east-west direction and 109 km in the north-south direction. The average size of UWE using satellite altimetric data is estimated to be 1~25 km smaller than that using water temperature cross-sectional data.