计算力学学报
計算力學學報
계산역학학보
CHINESE JOURNAL OF COMPUTATIONAL MECHANICS
2004年
4期
412-418
,共7页
冯其京%何长江%张敏%梁仙红%张树道
馮其京%何長江%張敏%樑仙紅%張樹道
풍기경%하장강%장민%량선홍%장수도
欧拉%界面%二维%三维
歐拉%界麵%二維%三維
구랍%계면%이유%삼유
Euler%interface%2 dimension%3 dimension
本文所提算法适用于二维和三维多介质流体力学两步欧拉数值方法中输运计算的混合网格(包括自由面网格)界面处理.在一个混合网格中,界面被近似地看作直线(二维)或平面(三维).整个方法分为三步:(1)第一步,用混合网格周围的八个网格的介质面积份额(二维)或二十六个网格的介质体积份额(三维)确定界面的法线方向;第二步,用混合网格的本身的介质面积份额(二维)或体积份额(三维)确定界面的方程(位置);第三步,用此直线方程求出通过网格边界的流以及下一时刻网格的面积份额(二维)或体积份额(三维).最后给出了用此方法所做的一些数值计算及与SLIC算法的比较.
本文所提算法適用于二維和三維多介質流體力學兩步歐拉數值方法中輸運計算的混閤網格(包括自由麵網格)界麵處理.在一箇混閤網格中,界麵被近似地看作直線(二維)或平麵(三維).整箇方法分為三步:(1)第一步,用混閤網格週圍的八箇網格的介質麵積份額(二維)或二十六箇網格的介質體積份額(三維)確定界麵的法線方嚮;第二步,用混閤網格的本身的介質麵積份額(二維)或體積份額(三維)確定界麵的方程(位置);第三步,用此直線方程求齣通過網格邊界的流以及下一時刻網格的麵積份額(二維)或體積份額(三維).最後給齣瞭用此方法所做的一些數值計算及與SLIC算法的比較.
본문소제산법괄용우이유화삼유다개질류체역학량보구랍수치방법중수운계산적혼합망격(포괄자유면망격)계면처리.재일개혼합망격중,계면피근사지간작직선(이유)혹평면(삼유).정개방법분위삼보:(1)제일보,용혼합망격주위적팔개망격적개질면적빈액(이유)혹이십륙개망격적개질체적빈액(삼유)학정계면적법선방향;제이보,용혼합망격적본신적개질면적빈액(이유)혹체적빈액(삼유)학정계면적방정(위치);제삼보,용차직선방정구출통과망격변계적류이급하일시각망격적면적빈액(이유)혹체적빈액(삼유).최후급출료용차방법소주적일사수치계산급여SLIC산법적비교.
The algorithm proposed in this paper can be used to deal with interface of mixed cell (including free surface cell) in remapping step of twostep 2-d and 3-d multi-component Eulerian hydrodynamics numerical methods. In a mixed cell the interface is regarded as a straight line in 2-d problems or plane in 3-d problems.The whole method is composed of three steps: (1) Using area fractions of eight cells in 2-d problems or volume fractions of twenty six cells in 3-d problems surrounding the mixed cell to determine the normal orientation of the interface. (2) Determining the equation (position) of the interface with the area fraction (2-d) or volume fraction (3-d) of the mixed cell. (3) Calculating the flux passing across the boundary of cell, and area fraction (2-d) or volume fraction (3-d) of cell at next time. In the last part of this paper, some 2-d and 3-d results of calculation in comparison with results of MEPH (simplified SLIC) are presented.