• Structural Engineering and Mechanics
• /
• v.74 no.5
• /
• pp.679-698
• /
• 2020

In this paper, the meshless local Petrov-Galerkin (MLPG) method is developed for dynamic analysis of non-symmetric nanocomposite cylindrical shell equations of elastic wave motion with nonlinear grading patterns under shock loading. The mechanical properties of the nanocomposite cylinder are obtained based on a micro-mechanical model. In this study, four kinds of grading patterns are assumed for carbon nanotube mechanical properties. The displacements can be approximated using shape function so, the multiquadrics (MQ) Radial Basis Functions (RBF) are used as the shape function. In order to discretize the derived equations in time domains, the Newmark time approximation scheme with suitable time step is used. To demonstrate the accuracy of the present method for dynamic analysis, at the first a problem verifies with analytical solution and then the present method compares with the finite element method (FEM), finally, the present method verifies by using the element free Galerkin (EFG) method. The comparison shows the high capacity and accuracy of the present method in the dynamic analysis of cylindrical shells. The capability of the present method to dynamic analysis of non-symmetric nanocomposite cylindrical shell is demonstrated by dynamic analysis of the cylinder with different kinds of grading patterns and angle of nanocomposite reinforcements. The present method shows high accuracy, efficiency and capability to dynamic analysis of non-symmetric nanocomposite cylindrical shell, which it furnishes a ground for a more flexible design.

• Journal of the Society of Naval Architects of Korea
• /
• v.52 no.6
• /
• pp.471-477
• /
• 2015

When a ship sails in a seaway, the resistance on a ship increases due to incident waves and winds. The magnitude of added resistance amounts to about 15–30% of a calm-water resistance. An accurate prediction of added resistance in waves, therefore, is essential to evaluate the performance of a ship in a real sea state and to design an optimum hull form from the viewpoint of the International Maritime Organization (IMO) regulations such as Energy Efficiency Design Index (EEDI) and Energy Efficiency Operational Indicator (EEOI). The present study considers added resistance problem of AFRAMAX-class tankers with the conventional bow and Ax-bow shapes. Added resistance due to waves is successfully calculated using 1) a three-dimensional time-domain seakeeping computations based on a Rankine panel method (three-dimensional panel) and 2) a commercial CFD program (STAR-CCM+). In the hydrodynamic computations of a three-dimensional panel method, geometric nonlinearity is accounted for in Froude-Krylov and restoring forces using simple wave corrections over exact wet hull surface of the tankers. Furthermore, a CFD program is applied by performing fully nonlinear computation without using an analytical formula for added resistance or empirical values for the viscous effect. Numerical computations are validated through four degree-of-freedom model-scale seakeeping experiments in regular head waves at the deep towing tank of Hyundai Heavy Industries.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.7 s.100
• /
• pp.788-798
• /
• 2005

The monopole theory has long been used to model air-pumped effect from the elastic cavities in car tire. This approach models the change of an air as a Piston moving backward and forward on a spring and equates local air movements exactly with the volume changes of the system. Thus, the monopole theory has a restricted domain of applicability due to the usual assumption of a small amplitude acoustic wave equation and acoustic monopole theory This paper describes an approach to predict the air-pumping noise of a car tyre with CFD/Kirchhoff integral method. The tyre groove is simply modeled as piston-cavity-sliding door geometry and with the aid of CFD technique flow properties in the groove of rolling car tyre are acquired.'rhese unsteady flow data are used as a air-pumping source in the next CFD calculation of full tyre-road geometry. Acoustic far field is predicted from Kirchhoff integral method by using unsteady flow data in space and time which is provided by the CFD calculation of full tyre-road domain. This approach can cover the non-linearity of acoustic monopole theory with the aid of Non-linear governing equation in CFD calculation. The method proposed in this paper is applied to the prediction of air-pumping noise of simply modeled car tyre and through the predicted results, the influence of nonlinear effect on air-pumping noise propagation is investigated.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.6
• /
• pp.461-468
• /
• 2012

In this study, a Cartesian-grid method based on finite volume approach is applied to simulate the ship motions in large amplitude waves. Fractional step method is applied for pressure-velocity coupling and TVD limiter is used to interpolate the cell face value for the discretization of convective term. Water, air, and solid phases are identified by using the concept of volume-fraction function for each phase. In order to capture the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with weighed line interface calculation (WLIC) method which considers multidimensional information. The volume fraction of solid body embedded in the Cartesian grid system is calculated using a level-set based algorithm, and the body boundary condition is imposed by a volume weighted formula. Numerical simulations for the two-dimensional barge type model and Wigley hull in linear waves have been carried out to validate the newly developed code. To demonstrate the applicability for highly nonlinear wave-body interactions such as green water on the deck, numerical analysis on the large-amplitude motion of S175 containership is conducted and all computational results are compared with experimental data.

• 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.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2292-2314
• /
• 2006

Efficient numerical method is developed for the prediction of aerodynamic noise generation and propagation in low Mach number flows such as aeolian tone noise. The proposed numerical method is based on acoustic/viscous splitting techniques of which acoustic solvers use simplified linearised Euler equations, full linearised Euler equations and nonlinear perturbation equations as acoustic governing equations. All of acoustic equations are forced with immersed surface dipole model which is developed for the efficient computation of aerodynamic noise generation and propagation in low Mach number flows in which dipole source, originating from unsteady pressure fluctuation on a solid surface, is known to be more efficient than quadrupole sources. Multi-scale overset grid technique is also utilized to resolve the complex geometries. Initially, aeolian tone from single cylinder is considered to examine the effects that the immersed surface dipole models combined with the different acoustic governing equations have on the overall accuracy of the method. Then, the current numerical method is applied to the simulation of the aeolian tones from twin cylinders aligned perpendicularly to the mean flow and separated 3 diameters between their centers. In this configuration, symmetric vortices are shed from twin cylinders, which leads to the anti-phase of the lift dipoles and the in-phase of the drag dipoles. Due to these phase differences, the directivity of the fluctuating pressure from the lift dipoles shows the comparable magnitude with that from the drag dipoles at 10 diameters apart from the origin. However, the directivity at 100 diameters shows that the lift-dipole originated noise has larger magnitude than, but still comparable to, that of the drag-dipole one. Comparison of the numerical results with and without mean flow effects on the acoustic wave emphasizes the effects of the sheared background flows around the cylinders on the propagating acoustic waves, which is not generally considered by the classic acoustic analogy methods. Through the comparison of the results using the immersed surface dipole models with those using point sources, it is demonstrated that the current methods can allow for the complex interactions between the acoustic wave and the solid wall and the effects of the mean flow on the acoustic waves.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.18 no.1
• /
• pp.1-12
• /
• 2013

Acoustic backscatter profiles were measured by Eco-sounder along an east-west section in the mid-eastern Yellow Sea and at an anchoring station in the low salinity region off the Keum River estuary in September 2012, with observing physical water property structure by CTD. Tidal front was established around the sand ridge developed in 50 m depth region. Internal waves measured by Eco-sounder during low tide period in the eastern side of the sand ridge were nonlinear depression waves with wave height of 15 m and mean wavelength of 500 m. These waves were interpreted into tidal internal waves that were produced by tidal current flowing over the sand ridge to the southeast. When weakly non-linear soliton model was applied, propagation speed and period of these internal depression wave were 50 m/s and 16~18 min. Red tides by Dinoflagelates Cochlodinium were observed in the sea surface where strong acoustic scattering layer was raised up to 7 m. Hourly CTD profiles taken at the anchoring station off the Keum River estuary showed the halocline depth change by tidal current and land-sea breeze. When tidal current flowed strongly to the northeast during flood period and land-breeze of 7 m/s blew to the west, the halocline was temporally raised up as much as 2 m and acoustic profile images showed a complex structure in the surface layer within 5-m depth: in tens of seconds the declined acoustic structure of strong and weak scattering signals alternatively appeared with entrainment and intrusion shape. These acoustic profile structures in the surface mixed layer were observed for the first time in the coastal sea of the mid-eastern Yellow Sea. The acoustic profile images and turbidity data suggest that relatively transparent low-layer water be intruded or entrained into the turbid upper-layer water by vertical shear between flood current and land breeze-induced surface current.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.6
• /
• pp.306-318
• /
• 2018

In order to evaluate the seismic capacity of massive vertical type breakwaters which have intensively been deployed along the coast of South Korea over the last two decades, we carry out the preliminary numerical simulation against the PoHang, GyeongJu, Hachinohe 1, Hachinohe 2, Ofunato, and artificial seismic waves based on the measured time series of ground acceleration. Numerical result shows that significant sliding can be resulted in once non-negligible portion of seismic energy is shifted toward the longer period during its propagation process toward the ground surface in a form of shear wave. It is well known that during these propagation process, shear waves due to the seismic activity would be amplified, and non-negligible portion of seismic energy be shifted toward the longer period. Among these, the shift of seismic energy toward the longer period is induced by the viscosity and internal friction intrinsic in the soil. On the other hand, the amplification of shear waves can be attributed to the fact that the shear modulus is getting smaller toward the ground surface following the descending effective stress toward the ground surface. And the weakened intensity of soil as the number of attacking shear waves are accumulated can also contribute these phenomenon (Das, 1993). In this rationale, we constitute the numerical model using the model by Hardin and Drnevich (1972) for the weakened shear modulus as shear waves go on, and shear wave equation, in the numerical integration of which $Newmark-{\beta}$ method and Modified Newton-Raphson method are evoked to take nonlinear stress-strain relationship into account. It is shown that the numerical model proposed in this study could duplicate the well known features of seismic shear waves such as that a great deal of probability mass is shifted toward the larger amplitude and longer period when shear waves propagate toward the ground surface.

• Atmosphere
• /
• v.23 no.3
• /
• pp.293-305
• /
• 2013

Nonhydrostatic effects on convectively forced mesoscale flows in two dimensions are numerically investigated using a nondimensional model. An elevated heating that represents convective heating due to deep cumulus convection is specified in a uniform basic flow with constant stability, and numerical experiments are performed with different values of the nonlinearity factor and nonhydrostaticity factor. The simulation result in a linear system is first compared to the analytic solution. The simulated vertical velocity field is very similar to the analytic one, confirming the high accuracy of nondimensional model's solutions. When the nonhydrostaticity factor is small, alternating regions of upward and downward motion above the heating top appear. On the other hand, when the nonhydrostaticity factor is relatively large, alternating updraft and downdraft cells appear downwind of the main updraft region. These features according to the nonhydrostaticity factor appear in both linear and nonlinear flow systems. The location of the maximum vertical velocity in the main updraft region differs depending on the degrees of nonlinearity and nonhydrostaticity. Using the Taylor-Goldstein equation in a linear, steady-state, invscid system, it is analyzed that evanescent waves exist for a given nonhydrostaticity factor. The critical wavelength of an evanescent wave is given by ${\lambda}_c=2{\pi}{\beta}$, where ${\beta}$ is the nonhydrostaticity factor. Waves whose wavelengths are smaller than the critical wavelength become evanescent. The alternating updraft and downdraft cells are formed by the superposition of evanescent waves and horizontally propagating parts of propagating waves. Simulation results show that the horizontal length of the updraft and downdraft cells is the half of the critical wavelength (${\pi}{\beta}$) in a linear flow system and larger than ${\pi}{\beta}$ in a weakly nonlinear flow system.

• Sleep Medicine and Psychophysiology
• /
• v.14 no.1
• /
• pp.42-48
• /
• 2007

Introduction: Detrended fluctuation analysis (DFA) is used as a way of studying nonlinearity of EEG. In this study, DFA is applied on sleep EEG of normal subjects to look into its nonlinearity in terms of EEG channels and sleep stages. Method: Twelve healthy young subjects (age：$23.8{\pm}2.5$ years old, male：female=7：5) have undergone nocturnal polysomnography (nPSG). EEG from nPSG was classified in terms of its channels and sleep stages and was analyzed by DFA. Scaling exponents (SEs) yielded by DFA were compared using linear mixed model analysis. Results: Scaling exponents (SEs) of sleep EEG were distributed around 1 showing long term temporal correlation and self-similarity. SE of C3 channel was bigger than that of O1 channel. As sleep stage progressed from stage 1 to slow wave sleep, SE increased accordingly. SE of stage REM sleep did not show significant difference when compared with that of stage 1 sleep. Conclusion: SEs of Normal sleep EEG showed nonlinear characteristic with scale-free fluctuation, long-range temporal correlation, self-similarity and self-organized criticality. SE from DFA differentiated sleep stages and EEG channels. It can be a useful tool in the research with sleep EEG.