• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.355-372
• /
• 2011

Based on the fully nonlinear velocity potential theory, the liquid sloshing in a three dimensional tank under random excitation is studied. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing scheme, B-spline curve, is applied to both the longitudinal and transverse directions of the tank to eliminate the possible saw-tooth instabilities. When the tank is undergoing one dimensional regular motion of small amplitude, the calculated results are found to be in very good agreement with linear analytical solution. In the simulation, the normal standing waves, travelling waves and bores are observed. The extensive calculation has been made for the tank undergoing specified random oscillation. The nonlinear effect of random sloshing wave is studied and the effect of peak frequency used for the generation of random oscillation is investigated. It is found that, even as the peak value of spectrum for oscillation becomes smaller, the maximum wave elevation on the side wall becomes bigger when the peak frequency is closer to the natural frequency.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.299-306
• /
• 2006

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The 'Hamilton-Jacobi (H-J)' equation and computationally robust numerical technique of 'up-wind scheme' lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H -J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes is not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.701-713
• /
• 2003

The BEM is a highly efficient method in the sense of economical computation. However, boundary integration is not easy for the complex geometry and moving surface, e.g. a rotating blade. Thus, Kirchhoff surface is designed in an effort to overcome the difficulty resulting from complex boundary conditions. A Kirchhoff surface is a fictitious surface which envelopes acoustic sources of main concern. Acoustic sources may be distributed on each Kirchhoff surface element according to their acoustic characteristics. In this study, an axial fan is assumed to have unsteady loading noise as a dominant source. Dipole sources can be modeled to solve the FW-H equation. Acoustic field is then computed by determining Kirchhoff surface on which near-field is implemented, to analyze the effect of Kirchhoff surface on it. The optimal shape and the location of Kirchhoff surface are discussed by comparing with experimental data acquired in an anechoic chamber.

• Geomechanics and Engineering
• /
• v.28 no.2
• /
• pp.107-116
• /
• 2022

The steady state of unsaturated soil takes a long time to achieve. The soil seepage behaviours and hydraulic properties depend highly on the wetting/drying rate. It is observed that the soil-water characteristic curve (SWCC) is dependent on the wetting/drying rate, which is known as the dynamic effect. The dynamic effect apparently influences the scanning curves and will substantially affect the seepage behavior. However, the previous models commonly ignore the dynamic effect and cannot quantitatively describe the hysteresis scanning loops under dynamic conditions. In this study, a dynamic hysteresis model for SWCC is proposed considering the dynamic change of contact angle and the moving of the contact line. The drying contact angle under dynamic condition is smaller than that under static condition, while the wetting contact angle under dynamic condition is larger than that under static condition. The dynamic contact angle is expressed as a function of the saturation rate according to the Laplace equation. The model is given by a differential equation, in which the slope of the scanning curve is related to the slope of the boundary curve by means of contact angle. Empirical models can simulate the boundary curves. Given the two boundary curves, the scanning curve can be well predicted. In this model, only two parameters are introduced to describe the dynamic effect. They can be easily obtained from the experiment, which facilitates the calibration of the model. The proposed model is verified by the experimental data recorded in the literature and is proved to be more convenient and effective.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.33 no.3
• /
• pp.93-100
• /
• 2021

Since the IB (Immersed Boundary) method, which can perform coupling analysis with objects and fluids having an impermeable boundary of arbitrary shape on a fixed grid system, has been developed, the IB method in various CFD models is increasing. The representative IB methods are the directing-forcing method and the ghost cell method. The directing-forcing type method numerically satisfies the boundary condition from the fluid force calculated at the boundary surface of the structure, and the ghost-cell type method is a computational method that satisfies the boundary condition through interpolation by placing a virtual cell inside the obstacle. These IB methods have a disadvantage in that the computational algorithm is complex. In this study, the simplified immersed boundary (SIB) method enables the analysis of temporary structures on a fixed grid system and is easy to expand to three proposed dimensions. The SIB method proposed in this study is based on a one-field model for immiscible two-phase fluid that assumes that the density function of each phase moves with the center of local mass. In addition, the volume-weighted average method using the density function of the solid was applied to handle moving solid structures, and the CIP method was applied to the advection calculation to prevent numerical diffusion. To examine the analysis performance of the proposed SIB method, a numerical simulation was performed on an object falling to the free water surface. The numerical analysis result reproduced the object falling to the free water surface well.

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.65.1-65.1
• /
• 2005

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.1174-1178
• /
• 2005

In this study, flood inundations have been simulated by using the numerical model FLUMEN solving the shallow-water equations with a finite volume method. Before applying to a real problem, the numerical model is first applied to simplified problems. Obtained numerical results are verified by comparing to available analytical solutions and laboratory measurements. Reasonable agreements are observed. The model is then applied to a simulation of flood events with real geometries. The results of the present study provide basic informations for a flood inundation map.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2001.04a
• /
• pp.435-442
• /
• 2001

A substructural model for a part of a long continuous bridge is identified using measured acceleration at limited locations within the part. Boundary spring constants and structural properties are idenfitied using a system identification method. The proposed method has been examined through simulation studies for static and dynamic responses and the results are discussed in the paper. The method is applied to an actual plate-girder gerber bridge with modal response obtained from a moving truck and construction blasting vibration

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.6
• /
• pp.1310-1320
• /
• 1989

본 논문에서는 비교적 많은 계산이 수행되어온 수평원관 주위의 상변화 열전달 문제를 예로 들어, 상변화 물질이 유한한 열전도율을 갖는 일정 두께의 원관 주위에 존재할 때의 상변화 문제에 대하여 초기 액상의 영역을 임의로 가정하는 불합리점을 제거하고 문제를 풀 수 있는 수치해석 방법을 제시하고자 한다.

• Smart Structures and Systems
• /
• v.4 no.2
• /
• pp.123-135
• /
• 2008

Damping capacity of SMA damping devices is simulated numerically under distinct geometry and loading conditions. Two-dimensional numerical simulations are performed on the basis of a phenomenological model of dynamics of martensite-austenite phase boundaries. Results of the simulations predict the time delay and the value of the stress transferred to other parts of a construction by a damper device.