• Steel and Composite Structures
• /
• v.27 no.3
• /
• pp.255-271
• /
• 2018

This paper deals with the transient dynamic analysis and elastic wave propagation in a functionally graded graphene platelets (FGGPLs)-reinforced composite thick hollow cylinder, which is subjected to shock loading. A micromechanical model based on the Halpin-Tsai model and rule of mixture is modified for nonlinear functionally graded distributions of graphene platelets (GPLs) in polymer matrix of composites. The governing equations are derived for an axisymmetric FGGPLs-reinforced composite cylinder with a finite length and then solved using a hybrid meshless method based on the generalized finite difference (GFD) and Newmark finite difference methods. A numerical time discretization is performed for the dynamic problem using the Newmark method. The dynamic behaviors of the displacements and stresses are obtained and discussed in detail using the modified micromechanical model and meshless GFD method. The effects of the reinforcement of the composite cylinder by GPLs on the elastic wave propagations in both displacement and stress fields are obtained for various parameters. It is concluded that the proposed micromechanical model and also the meshless GFD method have a high capability to simulate the composite structures under shock loadings, which are reinforced by FGGPLs. It is shown that the modified micromechanical model and solution technique based on the meshless GFD method are accurate. Also, the time histories of the field variables are shown for various parameters.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.119-125
• /
• 2004

Unsteady nonlinear wave motions on the free surface over a plane beach of constant slope 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. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. The second-order Stokes wave theory and solitary wave theory are employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with the plane beach of constant slope in surf zone. The results are compared with each other. The marker-density method is better then the hybrid method. Also they are compared with other existing experimental results. The Agreement between the experimental data and the computation results is good.

• Korean Journal of Optics and Photonics
• /
• v.15 no.6
• /
• pp.555-562
• /
• 2004

We propose a hybrid method combining the implicit with the explicit methods in order to reduce the calculation time and improve the convergence problem of the 3-dimensional finite-difference beam propagation method. The numerical simulation of a directional coupler is carried out by the proposed method. It is found from the simulation results that the calculation speed of our method is 10 times faster than that of direct solving techniques.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.7-13
• /
• 2011

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.

• Journal of Magnetics
• /
• v.21 no.3
• /
• pp.442-449
• /
• 2016

In this paper, an effective numerical analysis method is proposed for calculating dosimetry of the wireless power transfer system operating low-frequency ranges. The finite-difference time-domain (FDTD) method is widely used to analyze bio-electromagnetic field problems, which require high resolution, such as a heterogeneous whole-body voxel human model. However, applying the standard method in the low-frequency band incurs an inordinate number of time steps. We overcome this problem by proposing a modified finite-difference time-domain method which utilizes a quasi-static approximation with the surface equivalence theorem. The analysis results of the simple model by using proposed method are in good agreement with those from a commercial electromagnetic simulator. A simulation of the induced electric fields in a human head voxel model exposed to a wireless power transmission system provides a realistic example of an application of the proposed method. The simulation results of the realistic human model with the proposed method are verified by comparing it with the conventional FDTD method.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.1-7
• /
• 1997

A FA-FD hybrid method, developed for solving three-dimensional incompressible Navier-Stokes equations, is applied to calculate three-dimensional laminar flows through a square duct with a 90° bend. The method discretizes the convective terms in the primary flow direction with 3rd-order upwind finite-differences and the convective and diffusive terms in the transverse directions with the two-dimensional finite analytic method. The non-staggered grid system is used and the pressure-velocity coupling is achieved by a global iteration procedure based on the PISO algorithm. Detailed comparisons between the computed solutions and the available experimental data are given mainly for the velocity distributions at cross-sections in a 90° bend of a square duct with both fully developed and developing entry flows. Although the computational result shows generally a good agreement with the experimental data, there are some significant discrepancies underlining the necessity of more accurate numerical methods as well as reliable experimental data for their validation.

• Earthquakes and Structures
• /
• v.5 no.2
• /
• pp.161-205
• /
• 2013

We study seismically induced, anti-plane strain wave motion in a non-homogeneous geological region containing tunnels. Two different scenarios are considered: (a) The first models two tunnels in a finite geological region embedded within a laterally inhomogeneous, layered geological profile containing a seismic source. For this case, labelled as the first boundary-value problem (BVP 1), an efficient hybrid technique comprising the finite difference method (FDM) and the boundary element method (BEM) is developed and applied. Since the later method is based on the frequency-dependent fundamental solution of elastodynamics, the hybrid technique is defined in the frequency domain. Then, an inverse fast Fourier transformation (FFT) is used to recover time histories; (b) The second models a finite region with two tunnels, is embedded in a homogeneous half-plane, and is subjected to incident, time-harmonic SH-waves. This case, labelled as the second boundary-value problem (BVP 2), considers complex soil properties such as anisotropy, continuous inhomogeneity and poroelasticity. The computational approach is now the BEM alone, since solution of the surrounding half plane by the FDM is unnecessary. In sum, the hybrid FDM-BEM technique is able to quantify dependence of the signals that develop at the free surface to the following key parameters: seismic source properties and heterogeneous structure of the wave path (the FDM component) and near-surface geological deposits containing discontinuities in the form of tunnels (the BEM component). Finally, the hybrid technique is used for evaluating the seismic wave field that develops within a key geological cross-section of the Metro construction project in Thessaloniki, Greece, which includes the important Roman-era historical monument of Rotunda dating from the 3rd century A.D.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.36-42
• /
• 2006

The existing numerical approximations of convection flux, especially the spatial higher-order difference schemes, in unstructured cell-centered finite volume methods are examined in detail with each other and evaluated with respect to the accuracy through their application to a 2-D benchmark problem. Six higher-order schemes are examined, which include two second-order upwind schemes, two central difference schemes and two hybrid schemes. It is found that the 2nd-order upwind scheme by Mathur and Murthy(1997) and the central difference scheme by Demirdzic and Muzaferija(1995) have more accurate prediction performance than the other higher-order schemes used in unstructured cell-centered finite volume methods.

• Korean Journal of Hydrosciences
• /
• v.5
• /
• pp.85-97
• /
• 1994

A hybrid finite difference method for the longitudinal dispersion equation, which is based on combining the Holly-Preissmann scheme with fifth-degree Hermite interpolating polynomial and the generalized Crank-Nicholson scheme, is described and comparatively evaluated with other characteristics-based numerical methods. Longitudinal dispersion of an instantaneously-loaded pollutant source is simulated, and computational results are compared with the exact solution. The present method is free from wiggles regardless of the Courant number, and exactly reproduces the location of the peak concentration. Overall accuracy of the computation increases for smaller value of the weighting factor, $\theta$of the model. Larger values of $\theta$ overestimates the peak concentration. Smaller Courant number yields better accuracy, in general, but the sensitivity is very low, especially when the value of $\theta$ is small. From comparisons with the hybrid method using cubic interpolating polynomial and with splitoperator methods, the present method shows the best performance in reproducing the exact solution as the advection becomes more dominant.