Problems of orthogonal properties become more serious when the traditional transfinite interpolation (TFI) dynamic mesh method is employed for large deformations. Based on an analysis of the geometric relationship and interpolation features, an improvement for the present TFI method is proposed with a rotation correction. A new spring-TFI hybrid dynamic mesh method is developed for a structured mesh. First, each block of the computation domain is divided into several sub-blocks. Then, a spring network which connects the corners of all sub-blocks is established to smooth the mesh by means of spring analysis. Finally, a modified TFI method is used for calculating the inner deformations of the sub-blocks. Computational results of typical two and three dimensional viscous grids indicate that good orthogonal and smoothing properties can be achieved by rotation correction for large mesh deformations. In addition, the computational efficiency is slightly decreased than the traditional TFI method, but improved by 1 or 2 orders of magnitude when compared with the spring analogy method.
ZHANG Bing, HAN Jinglong
. Spring-TFI Hybrid Dynamic Mesh Method with Rotation Correction[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2011
, 32(10)
: 1815
-1823
.
DOI: CNKI:11-1929/V.20110419.1703.006
[1] Gordon W N, Hall C A. Construction of curvilinear coordinate systems and application to mesh generation[J]. International Journal for Numerical Methods in Engineering, 1973, 7(4): 461-477.
[2] Thompson J F, Thames F C, Mastin C W. Automatic numerical generation of body-fitted curvilinear coordinate system for field containing any number of arbitrary two-dimensional bodies[J]. Journal of Computational Physics, 1974, 15(3): 299-319.
[3] Reuther J J. Aerodynamics shape optimization of complex aircraft configurations via an adjoint formulation. AIAA-1996-20094, 1996.
[4] Byun C, Guruswamy G P. A parallel multi-block moving grid method for aeroelastic applications on full aircraft. AIAA-1998-24782, 1998.
[5] Jones W T, Samareh-Abolhassani J. A grid generation system for multi-disciplinary design optimization. AIAA-1995-1689, 1995.
[6] Tsai H M, Wong A S F, Cai J, et al. Unsteady flow calculations with a parallel multiblock moving mesh algorithm[J]. AIAA Journal, 2001, 39(6): 1021-1029.
[7] Gaitonde A L, Fiddes S P. Three-dimensional moving mesh method for the calculation of unsteady transonic flows[J]. The Aeronautical Journal, 1995, 99(984): 150-160.
[8] Farhat C, Degand C, Koobus B, et al. Torsional springs for two dimensional dynamic unstructured fluid meshes[J]. Computer Methods in Applied Mechanics and Engineering, 1998, 1(63): 231-245.
[9] Blom F J. Considerations on the spring analogy[J]. International Journal of Numerical Methods in Fluids, 2000, 32(6): 647-668.
[10] Degand C, Farhat C. A three-dimensional torsional spring analogy method for unstructured dynamic meshes[J]. Computers and Structures, 2002, 80(3-4): 305-316.
[11] Zeng D H, Ethier C R. A semi-torsional spring analogy model for updating unstructured meshes in 3d moving domains[J]. Finite Elements in Analysis and Design, 2005, 41(11-14): 1118-1139.
[12] Chew L P. Constrained Delaunay triangulations[J]. Algorithmica, 1989, 4(1): 97-108.
[13] Liu X Q, Qin N, Xia H. Fast dynamic grid deformation based on Delaunay graph mapping[J]. Journal of Computational Physics, 2006, 211(2): 405-423.
[14] 刘学强, 李青, 柴建忠, 等. 一种新的动网格方法及其应用[J]. 航空学报, 2008, 29(4): 817-822. Liu Xueqiang, Li Qing, Chai Jianzhong, et al. A new dynamic grid algrithm and its application[J]. Acta Aeronautica et Astronautica Sinca, 2008, 29(4): 817-822. (in Chinese).
[15] Chen P C, Hill L R. A three-dimensional boundary element method for CFD/CSD grid interfacing. AIAA-1999-1213, 1999.
[16] Johnson A A, Tezduyar T E. Simulation of multiple sphere falling in a liquid-filled tube[J]. Computer Methods in Applied Mechanics and Engineering, 1996, 134(4): 351-373.
CJA