• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.5
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.1
• /
• pp.11-22
• /
• 2019

In this paper, a Rankine source method is applied and validated to analyze the hydrodynamic response of a three-dimensional floating structure in the frequency domain. The boundary value problems for radiation and diffraction problem are solved by using a desingularized indirect boundary integral equation method (DIBIEM). The DIBIEM is simpler and faster than conventional methods based on the numerical surface integration of Green's function because the singularities of Green's function are located outside of fluid regions. In case of floating structure with complex geometry, it is difficult to desingularize the singularities of Green's function consistently. Therefore a mixed approach is carried out in this study. The mixed approach is partially desingularized except singularities of the body. Wave drift loads are calculated by the middle-field formulation method that is mathematically simple and has fast convergence. In order to validate the accuracy of the developed program, various numerical simulations are carried out and these results are analyzed and compared with previously published calculations and experiments.

• Journal of the Korea Society of Computer and Information
• /
• v.3 no.2
• /
• pp.85-97
• /
• 1998

The existing edge detection methods can not represent the real edge of object at fitting point or detect the edge which has unsufficient connecting trait. Especially, the two-fold thick edge detected by these methods cannot coincide real boundary of subject and it's location. To overcome these problems, we introduce the Genetic Algorithm(GA) in edge detection. The energy function is the value of fixel's satisfaction degree to edge condition. And it consists of the fitness value to image formation type, fitness value to connecting trait to it's neighboring edge and evalulation function which can represents the edge at fitting point as one fixel. This method is superior to remove the noise in edge detection than the existing methods. And it also detects the clear and exact edge because it can find the one fixel which is located at fitting point and has strong connecting trait.

• 한국신재생에너지학회:학술대회논문집
• /
• 2005.06a
• /
• pp.534-537
• /
• 2005

Free Vibration Analysis using Non-dimensional Dynamic Influence Function (NDIF) is extended to arbitrarily shaped plates including polygonal plates. Since the corners of polygonal plates have indefinite normal directions and additional boundary conditions related to a twisting moment at a corner along with moment and shear force zero conditions, it is not easy to apply the NDIF method to polygonal plates wi th the free boundary condition. Moreover, owing to the fact that the local polar coordinate system, which has been introduced for free plates with smoothly varying edges, cannot be employed for the straight edges of the polygonal plates, a new coordinate system is required for the polygonal plates. These problems are solved by developing the new method of modifying a corner into a circular arc and setting the normal direction at the corner to an average value of normal direct ions of two edges adjacent to the corner. Some case studies for plates with various shapes show that the proposed method gives credible natural frequencies and mode shapes for various polygons that agree well with those by an exact method or FEM (ANSYS).

• Journal of applied mathematics & informatics
• /
• v.28 no.1_2
• /
• pp.163-173
• /
• 2010

We consider the system $$\{{{-{\Delta}_pu\;=\;{\lambda}f(\upsilon),\;\;\;x\;{\in}\;{\Omega}, \atop -{\Delta}_q{\upsilon}\;=\;{\mu}g(u),\;\;\;x\;{\in}\;{\Omega},} \atop u\;=\;\upsilon\;=\;0,\;\;\;x\;{\in}\;{\partial\Omega},}$$ where ${\Delta}_pu\;=\;div(|{\nabla}_u|^{p-2}{\nabla}_u)$, ${\Delta}_{q{\upsilon}}\;=\;div(|{\nabla}_{\upsilon}|^{q-2}{\nabla}_{\upsilon})$, p, $q\;{\geq}\;2$, $\Omega$ is a ball in $\mathbf{R}^N$ with a smooth boundary $\partial\Omega$, $N\;{\geq}\;1$, $\lambda$, $\mu$ are positive parameters, and f, g are smooth functions that are negative at the origin and f(x) ~ $x^m$ g(x) ~ $x^n$ for x large for some m, $n\;{\geq}\;0$ with mn < (p - 1)(q - 1). We establish the existence and uniqueness of positive radial solutions when the parameters $\lambda$ and $\mu$ are large.

• Smart Structures and Systems
• /
• v.2 no.3
• /
• pp.209-224
• /
• 2006

The use of smart dampers to optimally control the response of structures is on the increase. To maximize the potential use of such damper systems, their accurate modeling and assessment of their performance is of vital interest. In this study, the performance of a controllable fluid dashpot damper, in terms of damper forces, damper dynamic range and damping force hysteretic loops, respectively, is studied mathematically. The study employs a damper Bingham-Maxwell (BingMax) model whose mathematical formulation is developed using a Fourier series technique. The technique treats this one-dimensional Navier-Stokes's momentum equation as a linear superposition of initial-boundary value problems (IBVPs): boundary conditions, viscous term, constant Direct Current (DC) induced fluid plug and fluid inertial term. To hold the formulation applicable, the DC current level to the damper is supplied as discrete constants. The formulation and subsequent simulation are validated with experimental results of a commercially available magneto rheological (MR) dashpot damper (Lord model No's RD-1005-3) subjected to a sinusoidal stroke motion using a 'SCHENK' material testing machine in the Materials Laboratory at the University of Technology, Sydney.

• Coupled systems mechanics
• /
• v.4 no.3
• /
• pp.263-277
• /
• 2015

This paper presents a coupled FEM-BEM strategy for the numerical analysis of elastodynamic problems where infinite-domain models and complex heterogeneous media are involved, rendering a configuration in which neither the Finite Element Method (FEM) nor the Boundary Element Method (BEM) is most appropriate for the numerical analysis. In this case, the coupling of these methodologies is recommended, allowing exploring their respective advantages. Here, frequency domain analyses are focused and an iterative FEM-BEM coupling technique is considered. In this iterative coupling, each sub-domain of the model is solved separately, and the variables at the common interfaces are iteratively updated, until convergence is achieved. A relaxation parameter is introduced into the coupling algorithm and an expression for its optimal value is deduced. The iterative FEM-BEM coupling technique allows independent discretizations to be efficiently employed for both finite and boundary element methods, without any requirement of matching nodes at the common interfaces. In addition, it leads to smaller and better-conditioned systems of equations (different solvers, suitable for each sub-domain, may be employed), which do not need to be treated (inverted, triangularized etc.) at each iterative step, providing an accurate and efficient methodology.

• Advances in Computational Design
• /
• v.3 no.2
• /
• pp.165-190
• /
• 2018

In this article, the static behavior of non-prismatic sandwich beams composed of functionally graded (FG) materials is investigated for the first time. Two types of beams in which the variation of elastic modulus follows a power-law form are studied. The principle of minimum total potential energy is applied along with the Ritz method to derive and solve the governing equations. Considering conventional boundary conditions, Chebyshev polynomials of the first kind are used as auxiliary shape functions. The formulation is developed within the framework of well-known Timoshenko and Reddy beam theories (TBT, RBT). Since the beams are simultaneously tapered and functionally graded, bending and shear stress pushover curves are presented to get a profound insight into the variation of stresses along the beam. The proposed formulations and solution scheme are verified through benchmark problems. In this context, excellent agreement is observed. Numerical results are included considering beams with various cross sectional types to inspect the effects of taper ratio and gradient index on deflections and stresses. It is observed that the boundary conditions, taper ratio, gradient index value and core to the thickness ratio significantly influence the stress and deflection responses.

• Structural Engineering and Mechanics
• /
• v.4 no.5
• /
• pp.513-527
• /
• 1996

In this paper, both an approximate expression and an exact expression for the contribution factor of an element to the natural frequency of the finite element discretized system of a structure in general and a membrane in particular have been derived from the energy conservation principle and the finite element formulation of structural eigenvalue problems. The approximate expression for the contribution factor of an element is used to predict and determine the elements to be removed in an iteration since it depends only on the quantities associated with the old system in the iteration. The exact expression for the contribution factor of an element makes it possible to check whether the element is correctly removed at the end of an iteration because it depends on both the old system and the new system in the iteration. Thus, the combined use of the approximate expression and the exact expression allows a considerable number of elements to be removed in a single iteration so that the efficiency of the evolutionary structural optimization method can be greatly improved for solving the natural frequency optimization problem of a structure. A square membrane with different boundary supports has been chosen to investigate the general evolutionary path for the fundamental natural frequency of the structure. The related results indicated that if the objective of a structural optimization is to raise the fundamental natural frequency of the structure to an optimal value, the general evolutionary path during its optimization is that the elements are gradually removed along the direction from the area surrounded by the contour of the highest value to that surrounded by the contour of the lowest value.

• International Journal of Ocean System Engineering
• /
• v.2 no.2
• /
• pp.116-138
• /
• 2012

As in the other fields of mechanics, analytical methods represent an important analysis tool in marine hydrodynamics. The analytical approach is interesting for different reasons : it gives reference results for numerical codes verification, it gives physical insight into some complicated problems, it can be used as a simplified predesign tool, etc. This approach is of course limited to some simplified geometries (cylinders, spheres, ...), and only the case of one or more cylinders, truncated or not, will be considered here. Presented methods are basically eigenfunction expansions whose complexity depends on the boundary conditions. The hydrodynamic boundary value problem (BVP) is formulated within the usual assumptions of potential flow and is additionally simplified by the perturbation method. By using this approach, the highly nonlinear problem decomposes into its linear part and the higher order (second, third, ...) corrections. Also, periodicity is assumed so that the time dependence can be factorized i.e. the frequency domain formulation is adopted. As far as free surface flows are concerned, only cases without or with small forward speed are sufficiently simple to be solved semi-analytically. The problem of the floating body advancing in waves with arbitrary forward speed is far more complicated. These remarks are also valid for the general numerical methods where the case of arbitrary forward speed, even linearized, is still too difficult from numerical point of view, and "it is fair to say that there exists at present no general practical numerical method for the wave resistance problem" [9], and even less for the general seakeeping problem. We note also that, in the case of bluff bodies like cylinders, the assumptions of the potential flow are justified only if the forward speed is less than the product of wave amplitude with wave frequency.