• International Journal of Naval Architecture and Ocean Engineering
• /
• 제7권3호
• /
• pp.610-625
• /
• 2015

In this paper the wave run-up around a circular column in regular waves is numerically calculated to investigate the applicability of the Modified Marker-Density (MMD) method to prediction of wave run-up around an offshore platform. The MMD method is one of the methods to define the highly nonlinear free surface. The governing equations are the Navier-Stokes equations and the continuity equation which are computed in Cartesian grid system. To validate incident waves generated by numerical simulation, those are compared with the solutions of the Stokes $5^{th}$ order wave theory. The wave run-up simulations are performed varying the steepness and period of incident waves as referred experimental data. The numerical results are compared to the experimental data and the results show good agreements.

• 대한설비공학회:학술대회논문집
• /
• 대한설비공학회 2006년도 하계학술발표대회 논문집
• /
• pp.116-121
• /
• 2006

Numerical simulations for dendritic growth of crystals are conducted in this study by the level set method. The effect of order of difference is tested for reinitialization error in simple problems and authors founded in case of 1st order of difference that very fine grids have to be used to minimize the error and higher order of difference is desirable to minimize the reinitialization error The 2nd and 4th order Runge-Kutta scheme in time and 3rd and 5th order of WENO schemes with Godunov scheme are applied for space discretization. Numerical results are compared with the analytical theory, phase-field method and other researcher's level set method.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• 제2권1호
• /
• pp.51-62
• /
• 2008

This paper proposes a second-order implicit constraint enforcement method which yields enhanced controllability compared to a first-order implicit constraints enforcement method. Although the proposed method requires solving a linear system twice, it yields superior accuracy from the constraints error perspective and guarantees the precise and natural movement of objects, in contrast to the first-order method. Thus, the proposed method is the most suitable for exact prediction simulations. This paper describes the numerical formulation of second-order implicit constraints enforcement. To prove its superiority, the proposed method is compared with the firstorder method using a simple two-link simulation. In this paper, there is a reasonable discussion about the comparison of constraints error and the analysis of dynamic behavior using kinetic energy and potential energy.

• Geomechanics and Engineering
• /
• 제2권3호
• /
• pp.161-176
• /
• 2010

Discrete Element Method (DEM) simulations were developed to investigate the effectiveness of geosynthetic reinforcement and the effectiveness of maintenance techniques performed on a simulated ballast bed subjected to dynamic loading. The results from four samples subjected each one to a total of 425 load cycles are presented: one unreinforced and unmaintained sample, one unmaintained but reinforced sample, one unreinforced sample subjected to maintenance in the form of stoneblowing after 200 load cycles, and one unreinforced sample subjected to maintenance in the form of tamping after 200 load cycles. The obtained values of permanent deformation as a function of the applied number of load cycles for the four cases are presented together allowing a comparison of the effectiveness of each technique. Moreover, snapshots of the simulated track sections are presented at different moments of the simulations. The simulations indicated that the geosynthetic reinforcement may not be beneficial for the analyzed case while stoneblowing was the most effective maintenance technique.

• Ocean Systems Engineering
• /
• 제9권3호
• /
• pp.241-259
• /
• 2019

This paper presents a fully three-dimensional numerical approach for analyzing deepwater drilling riser-conductor system vortex-induced vibrations (VIV) including nonlinear soil-structure interactions (SSI). The drilling riser-conductor system is modeled as a tensioned beam with linearly distributed tension and is solved by a fully implicit discretization scheme. The fluid field around the riser-conductor system is obtained by Finite-Analytic Navier-Stokes (FANS) code, which numerically solves the unsteady Navier-Stokes equations. The SSI is considered by modeling the lateral soil resistance force according to nonlinear p-y curves. Overset grid method is adopted to mesh the fluid domain. A partitioned fluid-structure interaction (FSI) method is achieved by communication between the fluid solver and riser motion solver. A riser-conductor system VIV simulation without SSI is firstly presented and served as a benchmark case for the subsequent simulations. Two SSI models based on a nonlinear p-y curve are then applied to the VIV simulations. Also, the effects of two key soil properties on the VIV simulations of riser-conductor systems are studied.

• Computers and Concrete
• /
• 제25권2호
• /
• pp.181-192
• /
• 2020

Nowadays, steel fiber reinforced concretes (SFRCs) are widely used in practical applications. Significant experimental research has thus been carried out to determine the constitutive equations that represent the behavior of SFRCs under multiaxial loadings. However, numerical modelling of SFRCs via FEM has been challenging due to the complexities of the implementation of these constitutive equations. In this study, following the literature, a plasticity model is constructed for the behavior of SFRCs that involves the Willam-Warnke failure surface with the relevant evolution laws and a non-associated flow rule for determining the plastic deformations. For the precise (yet rapid) integration of the constitutive equations, an explicit substepping scheme consisting of yield intersection and drift correction algorithms is employed and thus implemented in ABAQUS via UMAT. The FEM model includes various material parameters that are determined from the experimental data. Three sets of parameters are used in the numerical simulations. While the first set is from the experiments that are conducted in this study on SFRC specimens with various contents of steel fibers, the other two sets are from the experiments reported in the literature. The response of SFRCs under multiaxial compression obtained from various numerical simulations are compared with the experimental data. The good agreement between numerical results and the experimental data indicates that not only the adopted plasticity model represents the behavior of SFRCs very well but also the implemented integration scheme can be employed in practical applications of SFRCs.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 춘계 학술대회논문집
• /
• pp.153-158
• /
• 2005

Numerical simulations of two-dimensional steady incompressible lid-driven flow in a square cavity are presented to verify the validity of a new solution code(PowerCFD) with unstructured grids. The code uses the non-staggered(collocated) grid approach which is very popular for incompressible flow analysis because of its numerical efficiency on the curvilinear or unstructured grids. Solutions are obtained for configurations with a Reynolds number as high as 10,000 with both rectangular and hybrid types of unstructured grid mesh. Interesting features of the flow are presented in detail and comparisons are made with benchmark solutions found in the literature. It is found that the code is capable of producing accurately the nature of the lid-driven cavity flow at high Reynolds numbers.

• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2006년도 춘계학술대회
• /
• pp.269-272
• /
• 2006

Numerical analysis of wind turbine scale effect was performed by using computational fluid dynamics. For the numerical analysis of wind turbine. Three dimensional Navier-Stokes solver with various turbulence models was tested and realizable k-e turbulence model was selected for the simulation of wind turbines. To validate the present method, performance of NREL (National Renewable Energy Laboratory) Phase VI wind turbine model was analyzed and compared with experiment and blind test data. Using the present method, numerical simulations for various size of wind tunnel model were carried out and characteristics were observed in detail. The power loss due to the interference between wind turbine and nacelle was also computed for relatively larger nacelle installation in wind tunnel test. The present results showed good correlations with experimental data and reasonable trends of scale effect of wind turbine.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제10권2호
• /
• pp.119-128
• /
• 2018

A numerical method is proposed to simulate level ice interaction with ship in transverse and longitudinal directions in time domain. A novel method is proposed to simulate non-symmetric transverse force in a stochastic way. On the basis of observations from the model tests, the simulation of longitudinal force combines the ice bending force acting on the waterline, submersion force below the waterline and ice friction forces caused by transverse force and ice floes rotation amidships. In the simulations the ship was fixed and towed through an intact ice sheet at a certain speed. The setup of the numerical simulation is similar to the ice tank setup as much as possible. The simulated results are compared with model tests data and the results show good agreement with the measurement.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 1994년도 Proceedings of the Korea Automatic Control Conference, 9th (KACC) ; Taejeon, Korea; 17-20 Oct. 1994
• /
• pp.71-75
• /
• 1994

This paper presents a method of optimal design of an automatic transfer system which is controlled by the electro-pneumatic servo scheme. The electro-pneumatic automatic transfer system can move parts to desired points or displace defective parts. The dynamic performance of the system can be examined by observing the behavior of the output. The output of the servo control system is the motion of the cylinder, pneumatic actuator. The dynamic performance of the cylinder is governed by the parameters of the components of the entire system. The optimal design can be accomplished by selecting of the parameters such that the desired dynamic performance of the cylinder is obtained. The optimal set of parameters might be obtained through the repeated simulations. Repeated simulations, however, is not effective to determine the optimal set of parameters since the set of parameters is large. This paper presents modeling, application of an optimization method, and the numerical results. The optimization algorithm utilizes the concept of the conjugate gradient method. The results show that the suggested optimization scheme can render faster convergence of iteration compared to other method based on an algebraic optimization method and can reduce the design efforts.