• Journal of Internet Computing and Services
• /
• v.11 no.2
• /
• pp.99-107
• /
• 2010

In this paper a blind watermarking algorithm for digital image is presented. The proposed method operates in wavelet domain. The watermark is decomposed into 2D-IMFs using BEMD which is the 2-dimensional extension of 1 dimensional empirical mode decomposition. The CDMA based on SS technique is applied to watermark embedding and detection process. In the watermark embedding process, each IMF of watermark is embedded into middle frequency subimages in wavelet domain, so subimages just include partial information about embedded watermark. By characteristics of BEMD, when the partial information of watermark is synthesized, the original watermark is reconstructed. The experimental results show that the proposed watermarking algorithm is imperceptible and moreover is robust against JPEG compression, common image processing distortions.

• Structural Engineering and Mechanics
• /
• v.14 no.2
• /
• pp.151-170
• /
• 2002

Moving force identification is a very important inverse problem in structural dynamics. Most of the identification methods are eventually converted to a linear algebraic equation set. Different ways to solve the equation set may lead to solutions with completely different levels of accuracy. Based on the measured bending moment responses of the bridge made in laboratory, this paper presented the time domain method (TDM) and frequency-time domain method (FTDM) for identifying the two moving wheel loads of a vehicle moving across a bridge. Directly calculating pseudo-inverse (PI) matrix and using the singular value decomposition (SVD) technique are adopted as means for solving the over-determined system equation in the TDM and FTDM. The effects of bridge and vehicle parameters on the TDM and FTDM are also investigated. Assessment results show that the SVD technique can effectively improve identification accuracy when using the TDM and FTDM, particularly in the case of the FTDM. This improved accuracy makes the TDM and FTDM more feasible and acceptable as methods for moving force identification.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.379-392
• /
• 2002

Numerical flow computations around an aeroelastic 3D square cylinder immersed in the turbulent boundary layer are shown. Present computational code can be characterized by three numerical aspects which are 1) the method of artificial compressibility is adopted for the incompressible flow computations, 2) the domain decomposition technique is used to get better grid point distributions, and 3) to achieve the conservation law both in time and space when the flow is computed a with moving and transformed grid, the time derivatives of metrics are evaluated using the time-and-space volume. To provide time-dependant inflow boundary conditions satisfying prescribed time-averaged velocity profiles, a convenient way for generating inflow turbulence is proposed. The square cylinder is modeled as a 4-lumped-mass system and it vibrates with two-degree of freedom of heaving motion. Those blocks which surround the cylinder are deformed according to the cylinder's motion. Vigorous oscillations occur as the vortex shedding frequency approaches cylinder's natural frequencies.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.469-476
• /
• 2007

A domain/boundary decomposition technique is applied to carry out efficient finite element analyses of 3-D contact problems. Appropriate penalty functions are selected for connecting an interface and contact interfaces with neighboring subdomains that satisfy continuity constraints. As a consequence, all the effective stiffness matrices have positive definiteness, and computational efficiency can be improved to a considerable degree. If necessary, any complex-shaped 3-D domain can be divided into several simple-shaped subdomains without considering the conformity of meshes along the interface. With a set of numerical examples, the basic characteristics of computational efficiency are investigated carefully.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.21-28
• /
• 2024

In this study, a polyhedral domain decomposition method for Smoothed Particle Hydrodynamics (SPH) analysis is introduced. SPH which is one of meshless methods is a numerical analysis method for fluid flow simulation. It can be useful for analyzing fluidic soil or fluid-structure interaction problems. SPH is a particle-based method, where increased particle count generally improves accuracy but diminishes numerical efficiency. To enhance numerical efficiency, parallel processing algorithms are commonly employed with the Cartesian coordinate-based domain decomposition method. However, for parallel analysis of complex geometric shapes or fluidic problems under dynamic boundary conditions, the Cartesian coordinate-based domain decomposition method may not be suitable. The introduced polyhedral domain decomposition technique offers advantages in enhancing parallel efficiency in such problems. It allows partitioning into various forms of 3D polyhedral elements to better fit the problem. Physical properties of SPH particles are calculated using information from neighboring particles within the smoothing length. Methods for sharing particle information physically separable at partitioning and sharing information at cross-points where parallel efficiency might diminish are presented. Through numerical analysis examples, the proposed method's parallel efficiency approached 95% for up to 12 cores. However, as the number of cores is increased, parallel efficiency is decreased due to increased information sharing among cores.

• Proceedings of the KSRS Conference
• /
• 2008.10a
• /
• pp.7-10
• /
• 2008

We present a technique for constructing a digital elevation model (DEM) from contours. The elevation of each ground point in DEM is computed by interpolating the heights of the two adjacent contours of the point. The technique decomposes each sub-domain between adjacent contours into a set of sub-regions. The decomposition is accomplished by constructing a medial axis of the sub-domain. Each sub-region in the decomposition is classified into a variety of terrain features like hillsides, valleys, ridges, etc. The elevations of points are interpolated with different methods according to terrain features they belong to. For a given point in hillside, an approximate gradient line passing through the point is determined and the elevation of the point is interpolated from the known elevations of the two adjacent contours along the approximate gradient line. The univariate monotone rational Hermite spline is used for the interpolation. The DEM constructed by the technique is to be used to orthorectify the high-resolution KOMPSAT3 imagery.

• Bulletin of the Korean Mathematical Society
• /
• v.53 no.3
• /
• pp.853-874
• /
• 2016

Based on a particular overlapping domain decomposition technique, a parallel finite element discretization algorithm for the generalized Stokes equations is proposed and investigated. In this algorithm, each processor computes a local approximate solution in its own subdomain by solving a global problem on a mesh that is fine around its own subdomain and coarse elsewhere, and hence avoids communication with other processors in the process of computations. This algorithm has low communication complexity. It only requires the application of an existing sequential solver on the global meshes associated with each subdomain, and hence can reuse existing sequential software. Numerical results are given to demonstrate the effectiveness of the parallel algorithm.

• Earthquakes and Structures
• /
• v.2 no.4
• /
• pp.323-336
• /
• 2011

The problem of reconstruction of complete building response from a limited number of response measurements is considered. The response at the intermediate degrees of freedom is reconstructed by using piecewise cubic Hermite polynomial interpolation in time domain. The piecewise cubic Hermite polynomial interpolation is preferred over the spline interpolation due to its trend preserving character. It has been shown that factorization of response data in variable separable form via singular value decomposition can be used to derive the complete set of normal modes of the structural system. The time domain principal components can be used to derive empirical transfer functions from which the natural frequencies of the structural system can be identified by peak-picking technique. A reduced-rank approximation for the system flexibility matrix can be readily constructed from the identified mass-orthonormal mode shapes and natural frequencies.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.599-604
• /
• 2000

In order to study unsteady aerodynamic loads on high speed trains passing by each other at the speed of 350km/h, three-dimensional flow fields around trains during the crossing event are numerically simulated using the three-dimensional Euler equations. The Roe's FDS with MUSCL interpolation is employed to simulate wave phenomena properly. An efficient moving grid system based on domain decomposition techniques is developed to analyze the unsteady flow field induced by the restricted motion of a train on a rail. The numerical simulations of the trains passing by on the double-track are carried out to study the effect of the train nose-shape, the train length and the existence of tunnel when the crossing event occur. Unsteady aerodynamic loads side force and drag force-acting on the train during the crossing are numerically predicted and anlayzed. It is found that the strength of the side force mainly depends on the nose-shape, and that of drag force on tunnel existence. And it is observed that the push-pull like impulsive force successively acts on each car and acts in different directions between the neighborhood cars. The maximum change of the impulsive force reaches about 3 tons. These aerodynamic force data are absolutely necessary for the evaluation of the stability of the high speed multi-car train. The results also indicate the effectiveness of the present numerical method for the simulation of unsteady flow field induced by the bodies in the relative motion.