• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.6
• /
• pp.332-349
• /
• 2016

OLAFOAM is the powerful CFD code and is an expanded version of $OpenFOAM^{(R)}$, for wave mechanics simulation. The $OpenFOAM^{(R)}$ does provide many solvers to correspond to each object of the numerical calculation in a variety of fields. OLAFOAM's governing equation bases on VARANS (Volume-Averaged Reynolds-Averaged Navier-Stokes) equation, and the finite volume method is applied to numerical techniques. The program is coded in C++ and run on the Linux operating system. First of all, in this study, OLAFOAM was validated for 1) wave transformation inside porous structure under bore and regular wave conditions, 2) wave transformation by submerged breakwater under regular wave condition, and 3) regular wave transformation and resultant vertical velocity distribution under current by comparison with existing laboratory measurements. Hereafter, this study, which is almost no examination carried out until now, analyzed closely variation characteristics of water surface level, wave height, frequency spectrum, breaking waves, averaged velocity and turbulent kinetic energy around porous submerged breakwater in the wave and current coexisting field for the case of permeable or impermeable rear beach. It was revealed that the wave height fluctuation according to current direction(following or opposing) was closely related to the turbulent kinetic energy, and others.

• Ocean Systems Engineering
• /
• v.6 no.1
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.7
• /
• pp.922-930
• /
• 2020

Marine accidents involving fishing boats, caused by a loss of stability, have been increasing over the last decade. One of the main reasons for these accidents is a sudden wind attacks. In this regard, the wind loads acting on the ship hull need to be estimated accurately for safety assessments of the motion and maneuverability of the ship. Therefore, this study aims to develop a computational model for the inlet boundary condition and to numerically estimate the wind load acting on a fishing boat. In particular, wind loads acting on a fishing boat at the wind speed profile boundary condition were compared with the numerical results obtained under uniform wind speed. The wind loads were estimated at intervals of 15° over the range of 0° to 180°, and i.e., a total of 13 cases. Furthermore, a numerical mesh model was developed based on the results of the mesh dependency test. The numerical analysis was performed using the RANS-based commercial solver STAR-CCM+ (ver. 13.06) with the k-ω turbulent model in the steady state. The wind loads for surge, sway, and heave motions were reduced by 39.5 %, 41.6 %, and 46.1 % and roll, pitch, and yaw motions were 48.2 %, 50.6 %, and 36.5 %, respectively, as compared with the values under uniform wind speed. It was confirmed that the developed inlet boundary condition describing the wind speed gradient with respect to height features higher accuracy than the boundary condition of uniform wind speed. The insights obtained in this study can be useful for the development of a numerical computation method for ships.

• Journal of Navigation and Port Research
• /
• v.46 no.6
• /
• pp.491-500
• /
• 2022

It has been noted that an accurate estimation of wind loads on offshore structures such as an FLNG (Liquefied Natural Gas Floating P roduction Storage Offloading Units, LNG FPSOs) with a large topside plays an important role in the safety design of hull and mooring system. Therefore, the present study aims to develop a computational model for estimating the wind load acting on an FLNG. In particular, it is the sequel to the previous research by the author. The numerical computation model in the present study was modified based on the previous research. Numerical analysis for estimating wind loads was performed in two conditions for an interval of wind direction (α), 15° over the range of 0° to 360°. One condition is uniform wind speed and the other is the NPD model reflecting the wind speed profile. At first, the effect of sand-grain roughness on the speed profile of the NPD model was studied. Based on the developed NPD model, mesh convergence tests were carried out for 3 wind headings, i.e. head, quartering, and beam. Finally, wind loads on 6-degrees of freedom were numerically estimated and compared by two boundary conditions, uniform speed, and the NPD model. In the present study, a commercial RANS-based viscous solver, STAR-CCM+ (ver. 17.02) was adopted. In summary, wind loads in surge and yaw from the wind speed profile boundary condition were increased by 20.35% and 34.27% at most. Particularly, the interval mean of sway (45° < α <135°, 225° < α < 315°) and roll (60° < α < 135°, 225° < α < 270°) increased by 15.60% and 10.89% against the uniform wind speed (10m/s) boundary condition.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.