• Ocean Systems Engineering
• /
• 제1권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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• 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.

• 대한기계학회논문집B
• /
• 제27권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
• /
• 제28권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.

• 한국해안·해양공학회논문집
• /
• 제33권3호
• /
• pp.93-100
• /
• 2021

고정된 격자시스템에서 임의형상의 불투과 경계를 갖는 물체와 유체와 연성해석이 가능한 IB(Immersed Boundary)법이 개발 된 이후로 다양한 CFD 모델에서 IB법의 활용이 증가하고 있다. 기존의 IB법의 대부분은 구조물의 경계면에서 산정되는 유체력으로부터 수치적으로 경계조건을 만족시키는 directing-forcing법이나 구조물 내부에 가상셀을 위치시켜 보간을 통해 경계조건을 만족시키는 ghost-cell법들로 알고리즘이 복잡하다. 본 연구에서는 고정된 격자시스템에서 가동물체형 구조물 해석이 가능함과 더불어 3차원으로의 확장도 용이한 SIB(Simplified Immersed Boundary)법을 제안하였다. 본 연구에서 제안한 SIB법은 각 상(phase)의 밀도함수가 국소질량의 중심과 함께 이동하는 것으로 가정한 단일유체모델(one-field model for immiscible two-phase fluid)을 기초로 하였다. 또한 이동하는 고체상태의 구조물을 취급하기 위해 고체의 밀도함수를 이용한 체적가중평균법을 적용하고, 수치확산을 방지하기 위해 이류계산에는 CIP법을 적용하였다. 제안된 SIB법의 해석성능을 검토하기 위해 자유수면으로 낙하하는 물체에 대한 수치모의를 수행하였다. 수치해석결과는 자유수면으로 낙하하는 물체를 양호하게 재현하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2005년도 창립60주년 기념 추계 학술대회
• /
• pp.65.1-65.1
• /
• 2005

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• 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.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2001년도 춘계학술대회 논문집
• /
• 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

• 대한기계학회논문집
• /
• 제13권6호
• /
• pp.1310-1320
• /
• 1989

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

• Smart Structures and Systems
• /
• 제4권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.