• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.4
• /
• pp.204-210
• /
• 1998

This paper presents the development of the probability density function applicable for wave heights (peak-to-trough excursions) in finite water depth including shallow water depth. The probability distribution applicable to wave heights of a non-Gaussian random process is derived based on the concept of the maximum entropy method. When wave heights are limited by breaking wave heights (or water depth) and only first and second moments of wave heights are given, the probability density function developed is closed form and expressed in terms of wave parameters such as $H_m$(mean wave height), $H_{rms}$(root-mean-square wave height), $H_b$(breaking wave height). When higher than third moment of wave heights are given, it is necessary to solve the system of nonlinear integral equations numerically using Newton-Raphson method to obtain the parameters of probability density function which is maximizing the entropy function. The probability density function thusly derived agrees very well with the histogram of wave heights in finite water depth obtained during storm. The probability density function of wave heights developed using maximum entropy method appears to be useful in estimating extreme values and statistical properties of wave heights for the design of coastal structures.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.4
• /
• pp.525-532
• /
• 2010

In this paper, wave breakers which occur in two dimensional coasts are simulated using a SPH(Smoothed Particle Hydrodynamics) method which represents the movement of fluidic physical volume with particles. As continuative fluid is approximated to the particles, the simulations are performed using fully Lagrangian method without any grid system. Two-dimensional Navier-Stokes equations and continuity equation are used for the numerical simulations. To generate incident waves, a piston type wavemaker is employed. The accuracy of the wave which is numerically generated by the wavemaker is verified by comparing with analytical results. The computations are carried out with various wave heights and slopes. The wave patterns generated through the numerical simulations are compared with several existing experimental and computational results. Agreement between the experimental data and the computation results is comparatively good. Also, the breaker depth index and the breaker height index from the present calculations are compared with the existing experimental results, and the tendency is very similar.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.10a
• /
• pp.198-203
• /
• 2003

Unsteady two-dimensional nonlinear waves which are generated by the numerical wavemaker of plunging type are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier-Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Marker-density function method is adopted for the simulation of wave breaking phenomena, and the computations are carried out with various wave amplitudes and two section shapes of wavemaker. The computation results are compared with other existing computational and experimental results, and the agreement between the experimental data and the computation results is good.

• Nuclear Engineering and Technology
• /
• v.40 no.6
• /
• pp.477-488
• /
• 2008

The objective of the paper is to analyze the thermally induced density wave oscillations in water cooled boiling water reactors. A transient thermal hydraulic model is developed with a characteristics-based implicit finite-difference scheme to solve the nonlinear mass, momentum and energy conservation equations in a time-domain. A two-phase flow was simulated with a one-dimensional homogeneous equilibrium model. The model treats the boundary conditions naturally and takes into account the compressibility effect of the two-phase flow. The axial variation of the heat flux profile can also be handled with the model. Unlike the method of characteristics analysis, the present numerical model is computationally inexpensive in terms of time and works in a Eulerian coordinate system without the loss of accuracy. The model was validated against available benchmarks. The model was extended for the purpose of studying the flow-induced density wave oscillations in forced circulation and natural circulation boiling water reactors. Various parametric studies were undertaken to evaluate the model's performance under different operating conditions. Marginal stability boundaries were drawn for type-I and type-II instabilities in a dimensionless parameter space. The significance of adiabatic riser sections in different boiling reactors was analyzed in detail. The effect of the axial heat flux profile was also investigated for different boiling reactors.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.10a
• /
• pp.232-236
• /
• 2003

The anchor is laid on seabed and the main engine is worked to against incident environmental loads in typoon. As the main engine is broken down in the storm, the anchor chain is cutted and the vessel is drifted. Although a ship is moored by two point mooring lines to keep the her position, a ship is crashed into a rock because of typoon and the accident of oil spilling may be occured. In this paper, we studied the position-keeping of a ship which is analyized based on the slow motion maneuvering equations considering wave, current and wind. The direct integration method is employed to estimate wave loads. The current forces are calculated by using mathematical of MMG. The two point mooring forces are quasisatatically evaluated by using the catenary equation. The coefficeints of wind forces are modeled from Isherwood’s emperical data and the variation of wind speed is estimated by wind spectrum. The nonlinear motions of a two point moored ship are simulated considering wave, current, wind load in time domain.

• Journal of Ocean Engineering and Technology
• /
• v.18 no.3
• /
• pp.8-12
• /
• 2004

The anchor is laid on the seabed, and the main engine is working against incident environmental loads in a typhoon. As the main engine is broken Mum in the storm, the anchor chain is cut and the vessel drifts. Although a ship is moored by two-point mooring lines to maintain her position, it has crashed into a rock because of a typhoon, resulting in a possible accidental oil spillage. In this paper, we studied maintenance of a ship's position, which is analyzed based on the slow motion maneuvering equations considering wave, current, and wind. To estimate wave loads, the direct integration method is employed. The current forces are calculated, using MMG (Mathematical Modeling Group). Th two-point mooring forces are quasi-statistically evaluated, using the catenary equation. Th coefficients of wind forces are modeled from Isherwood's empirical data, and the variation of wind speed is estimated by wind spectrum. The nonlinear motions of a two-point moored ship are simulated, considering wave, current, and wind load, in specific domain of time.

• International Journal of Aeronautical and Space Sciences
• /
• v.4 no.1
• /
• pp.34-44
• /
• 2003

An advanced aeroelastic analysis using a computational structural dynamics (CSD), finite element method (FEM) and computational fluid dynamics (CFD) is presented in this Paper. A general aeroelastic analysis system is originally developed and applied to realistic design problems in the transonic flow region, where strong shock wave interactions exist. The present computational approach is based on the modal-based coupled nonlinear analysis with the matched-point concept and adopts the high-speed parallel processing technique on the low-cost network based PC-clustered machines. It can give very accurate and useful engineering data on the structural dynamic design of advanced flight vehicles. For the nonlinear unsteady aerodynamics in high transonic flow region, Euler equations using the unstructured grid system have been applied to easily consider complex configurations. It is typically shown that the advanced numerical approach can give very realistic and practical results for design engineers and safe flight tests. One can find that the present study conducts a virtual flutter flight test which are usually very dangerous in reality.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.2
• /
• pp.1-12
• /
• 1996

Hydrodynamic characteristics of shipping water on deck are investigated by using a simplified two-dimensional model. Formulation of the shipping water on deck leads to a nonlinear hyperbolic system of equations based on the shallow-water wave theory. Time-domain solution of these equations are obtained numerically using a finite difference method which adopts predictor-corrector method for time-marching and 2nd upwind differencing method for convection term calculation. To confirm the validity of the present numerical method, we calculated some shallow-water wave problems accompanying a bore and compared the obtained results with the analytic solutions. We found good agreements between them. Though the calculation results of shipping water on deck, we show that the shipping water flows into the deck as a rarefying wave arid grows into a bore after colliding with a deck structure. Also we examined the effects of acceleration and slope of deck and found that they have significant influences on the flow of shipping water.

• Journal of Ocean Engineering and Technology
• /
• v.7 no.1
• /
• pp.13-24
• /
• 1993

The dispersion in the oscillatory flow generated by gravitational waves above the spatially periodic repples is studied. The steady parts of equations describing the orbit of the passive particle in a two dimensional field are assumed to be simply trigonometric functions. From the view point of nonlinear dynamics, the motion of the particle is chaotic under externally time-periodic perturbations which come from the wave motion. Two cases considered here are; (i) shallow water, and (ii) deep water approximation.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.80-86
• /
• 1994

A numerical method for simulations of inviscid incompressible flow fields around a ship advancing on the free surface is developed. A body fitted coordinate system, generated by numerically solving elliptic type partial differential equations is used to conform the ship and free surface configurations. Three dimensional Euler equations transformed to the non-staggered body fitted coordinate system are discretised by finite difference method. Time and spatial derivatives are discretised by forward and centered differencings, respectively, and artificial dissipations are added to discretised convection terms for improvements of numerical stability. At each time steps, free surface elevations are recomputed to satisfy nonlinear free surface conditions. Poisson equations for pressure field are solved iteratively and the velocity field for next time step is extrapolated. To verify the developed numerical method, flow fields around a Wigley model are simulated(Fn=0.250-0.408) and compared with experimental data to show good agreements.