ALE寬頻元素法於不可壓縮流之發展與應用

Development and Applications of ALE Spectral Element Method for Simulation of Incompressible

蘇炯彰(Jiong-Zhang Su)

指導教授 ： 顏瑞和

國立臺灣大學/工學院/機械工程學研究所/碩士(2006年)

https://doi.org/10.6342/NTU.2006.01991

摘要

本文的主要目的在於發展ALE法的寬頻元素法用於處理具移動邊界的不可壓縮流問題，並應用此工具於流-固相互作用的研究。統御方程式的運動描述法採用ALE法是為了有效處理邊界的移動與減少網格的扭曲，而空間離散所採用的寬頻元素法具有高精確度的特性，在同一個數值精確度下需要計算的自由度比低階的方法相對來的少，程式上要處理的資料也就比較少，有助於提高流場的程式效率也有利於動態網格的應用，且使用的網格也可較大，可減輕網格的扭曲問題；本文以一維、二維的伯格方程式及二維不可壓縮之Navier-Stokes方程式驗證、討論程式的正確性與移動網格對數值收斂的影響；在實例的應用上，將會測試、討論三個例子，包括凹陷移動的管內流、流體流經一橫向振動圓柱與流體流經兩並肩排列橫向振動圓柱，前兩個例子在趨勢與大小上都與文獻的模擬結果符合，確認本文所發展的工具於實際問題的可靠性，最後用於觀察兩並肩排列橫向振動圓柱的現象，除數值結果的討論外，並與實驗結果比較同異之處。

關鍵字

移動邊界 ； ALE ； 寬頻元素法 ； 渦漩剝離 ； 振動圓柱 ； 共振

並列摘要

The main purpose of this thesis is to develop a spectral element method based on the ALE formulation for incompressible flow with moving boundary problems and present the application for fluid-structure interaction problem. Arbitrary Lagrangian Eulerian(ALE) formulation can offer robust treatment of the moving boundary and handle distortions of the computational mesh. Spectral element method, delivering high-order convergence characteristics, allow utilization of relatively fewer degrees of freedom than low-order methods for a desired accuracy. This is an advantage in reduction of data,efficiency of the algorithms and the application of dynamic mesh. Also, relatively larger elements than low-order methods enables spectral element ALE algorithms to reduce the problem of distortions. One dimensional and two dimensional Burger’s equation, incompressible two dimensional Navier-Stokes equation are applied to test the numerical accuracy and the effect of moving mesh on numerical convergence. The following two examples, flow in a channel with a moving indentation and flow past an transverse oscillating cylinder, show quantitative agreement with the numerical results and demonstrate the effectiveness of the method in some of these application. Detailed results are presented for last case, flow past an transverse oscillating cylinder, and preliminarily compared with the experimental results.

並列關鍵字

moving boundary ； ALE ； spectral element method ； vortex shedding ； oscillating cylinder ； lock-in

參考文獻

[3] C. W. Hirt, A. A. Amsden, and J. L. Cook, "Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds," Journal of Computational Physics, vol. 14, pp. 227-253, 1974.

[6] T. J. R. Hughes, W. K. Liu, and T. K. Zimmermann, "Lagrangian-Eulerian Finite-Element Formulation for Incompressible Viscous Flows," Computer Methods in Applied Mechanics and Engineering, vol. 29, pp. 329-349, 1981.

[7] W. K. Liu, T. Belytschko, and H. Chang, "An Arbitrary Lagrangian-Eulerian Finite-Element Method for Path-Dependent Materials," Computer Methods in Applied Mechanics and Engineering, vol. 58, pp. 227-245, 1986.

[8] F. P. T. Baaijems, "An U-Ale Formulation of 3-D Unsteady Viscoelastic Flow," International Journal for Numerical Methods in Engineering, vol. 36, pp. 1115-1143, 1993.

[9] E. Kuhl, H. Askes, and P. Steinmann, "An ALE formulation based on spatial and material settings of continuum mechanics. Part 1: Generic hyperelastic formulation," Computer Methods in Applied Mechanics and Engineering, vol. 193, pp. 4207-4222, 2004.