• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.1
• /
• pp.1-11
• /
• 2016

The objective of this paper is to present an analytical solution in deep water waves and verify the validity of the theory (Shin, 2015). Hence this is a follow-up to Shin (2015). Instead of a variational approach, another approach was considered for a more accurate assessment in this study. The products of two coefficients were not neglected in this study. The two wave profiles from the KFSBC and DFSBC were evaluated at N discrete points on the free-surface, and the combination coefficients were determined for when the two curves pass the discrete points. Thus, the solution satisfies the differential equation (DE), bottom boundary condition (BBC), and the kinematic free surface boundary condition (KFSBC) exactly. The error in the dynamic free surface boundary condition (DFSBC) is less than 0.003%. The wave theory was simplified based on the assumption tanh $D{\approx}1$ in this paper. Unlike the perturbation method, the results are possible for steep waves and can be calculated without iteration. The result is very simple compared to the 5th Stokes' theory. Stokes' breaking-wave criterion has been checked in this study.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.11
• /
• pp.1907-1916
• /
• 2003

During decades, there has been much progress in understanding of the inelastic behavior of the materials and numerous inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. To obtain the increment of state variable, its evolution laws are linearized by several approximation methods, such as general midpoint rule(GMR) or general trapezoidal rule(GTR). In this investigation, semi-implicit integration schemes using GTR and GMR were developed and implemented into ABAQUS by means of UMAT subroutine. The comparison of integration schemes was conducted on the simple tension case, and simple shear case and nonproportional loading case. The fully implicit integration(FI) was the most stable but amplified the truncation error when the nonlinearity of state variable is strong. The semi-implicit integration using GTR gave the most accurate results at tension and shear problem. The numerical solutions with refined time increment were always placed between results of GTR and those of FI. GTR integration with adjusting midpoint parameter can be recommended as the best integration method for viscoplastic equation considering nonlinear kinematic hardening.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Smart Structures and Systems
• /
• v.18 no.1
• /
• pp.155-181
• /
• 2016

In the present paper we discuss the stability and the post-buckling behaviour of thin piezoelastic plates. The first part of the paper is concerned with the modelling of such plates. We discuss the constitutive modelling, starting with the three-dimensional constitutive relations within Voigt's linearized theory of piezoelasticity. Assuming a plane state of stress and a linear distribution of the strains with respect to the thickness of the thin plate, two-dimensional constitutive relations are obtained. The specific form of the linear thickness distribution of the strain is first derived within a fully geometrically nonlinear formulation, for which a Finite Element implementation is introduced. Then, a simplified theory based on the von Karman and Tsien kinematic assumption and the Berger approximation is introduced for simply supported plates with polygonal planform. The governing equations of this theory are solved using a Galerkin procedure and cast into a non-dimensional formulation. In the second part of the paper we discuss the stability and the post-buckling behaviour for single term and multi term solutions of the non-dimensional equations. Finally, numerical results are presented using the Finite Element implementation for the fully geometrically nonlinear theory. The results from the simplified von Karman and Tsien theory are then verified by a comparison with the numerical solutions.

• Korean Journal of Hydrosciences
• /
• v.5
• /
• pp.1-16
• /
• 1994

In this paper distributed models for simulating spatially and temporally varied moving storm in a watershed were developed. The complete simulation in a watershed is achieved through two sequential flow simulations which are overland flow simulation and channel network flow simulation. Two dimensional continuity equation and momentum equation of kinematic approximation were used in the overland flow simulation. On the other hand, in the channel network simulation two types of governing equations which are one dimensional continuity and momentum equations between two adjacent sections in a channel, and continuity and energy equations at a channel junction were applied. The finite difference formulations were used in the channel network model. Macks Creek Experimental Watershed in Idaho, USA was selected as a target watershed and the moving storm on August 23, 1965, which continued from 3:30 P.M. to 5:30 P.M., was utilized. The rainfall intensity fo the moving storm in the watershed was temporally varied and the storm was continuously moved from one place to the other place in a watershed. Furthermore, runoff parameters, which are soil types, vegetation coverages, overland plane slopes, channel bed slopes and so on, are spatially varied. The good agreement between the hydrograph simulated using distributed models and the hydrograph observed by ARS are Shown. Also, the conservations of mass between upstreams and downstreams at channel junctions are well indicated and the wpatial and temporal vaiability in a watershed is well simulated using suggested distributed models.

• Structural Engineering and Mechanics
• /
• v.13 no.2
• /
• pp.135-154
• /
• 2002

This paper presents the convected material frame approach to study the nonlinear behavior of inelastic frame structures. The convected material frame approach is a modification of the co-rotational approximation by incorporating an adaptive convected material frame in the basic definition of the displacement vector and strain tensor. In the formulation, each discrete element is associated with a local coordinate system that rotates and translates with the element. For each load increment, the corresponding strain-displacement and nodal force-stress relationships are defined in the updated local coordinates, and based on the updated element geometry. The rigid body motion and deformation displacements are decoupled for each increment. This modified approach incorporates the geometrical nonlinearities through the continuous updating of the material frame geometry. A generalized nonlinear function is used to derive the inelastic constitutive relation and the kinematic hardening is considered. The equation of motion is integrated by an explicit procedure and it involves only vector assemblage and vector storage in the analysis by assuming a lumped mass matrix of diagonal form. Several numerical examples are demonstrated in close agreement with the solutions obtained by the ANSYS code. Numerical studies show that the proposed approach is capable of investigating large deflection of inelastic planar structures and providing an excellent numerical performance.

• Journal of Korea Water Resources Association
• /
• v.49 no.6
• /
• pp.459-467
• /
• 2016

LISFLOOD-FP model in which channel flows are resolved separately from the floodplain flows using either a kinematic or diffusive wave approximation has been used to analyze flooding behavior on the lower Mankyung River influenced by backwater. A calibration and validation process was applied using the previous flood events to assess the model performance. Sensitivity analysis was conducted for main calibrated parameters, such as Manning roughness coefficient and downstream boundary condition. Also, we examined the effect of warm-up for the initial conditions. The results show that the computed hydrograph is in good agreement with measured data on the study reach, even though it was a hydrologic kinematic wave model. The sensitive analysis show that the difference between the computed results may be greater depending on the used calibrated parameters and that the sufficient calibration/validation process against various flood events is necessary. If the flood inundation simulation is performed using the validated model, it is expected to be able to contribute about river planning and policy decision-making for flood damage reduction.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.3
• /
• pp.157-167
• /
• 2015

In order to provide simple and accurate wave theory in design of offshore structure, an analytical approximation is introduced in this paper. The solution is limited to flat bottom having a constant water depth. Water is considered as inviscid, incompressible and irrotational. The solution satisfies the continuity equation, bottom boundary condition and non-linear kinematic free surface boundary condition exactly. Error for dynamic condition is quite small. The solution is suitable in description of breaking waves. The solution is presented with closed form and dispersion relation is also presented with closed form. In the last century, there have been two main approaches to the nonlinear problems. One of these is perturbation method. Stokes wave and Cnoidal wave are based on the method. The other is numerical method. Dean's stream function theory is based on the method. In this paper, power series method was considered. The power series method can be applied to certain nonlinear differential equations (initial value problems). The series coefficients are specified by a nonlinear recurrence inherited from the differential equation. Because the non-linear wave problem is a boundary value problem, the power series method cannot be applied to the problem in general. But finite number of coefficients is necessary to describe the wave profile, truncated power series is enough. Therefore the power series method can be applied to the problem. In this case, the series coefficients are specified by a set of equations instead of recurrence. By using the set of equations, the nonlinear wave problem has been solved in this paper.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.112-122
• /
• 1993

A potential-based panel method is presented for the analysis of a partially or supercavitating two-dimensional hydrofoil at a finite submergence beneath a free surface, treating without approximation the effects of the finite Froude number and the hydrostatic pressure. Free surface sources and normal dipoles are distributed on the foil and cavity surfaces, their strength being determined by satisfying the kinematic and dynamic boundary conditions on the foil-cavity boundary. The cavity surface is determined iteratively as a part of the solution. Numerical results show that the wave profile is altered significantly due to the presence of the cavity. The buoyancy effect due to the hydrostatic pressure, which has usually been neglected in most of the cavitating flow analysis, is found playing an important role, especially for the supercavitating hydrofoil; the gravity field increases the cavity size in shallow submergence, but decreases it when deeply submerged, while the lift reduces at all submergence depth.