物理学报
物理學報
물이학보
2013年
6期
331-337
,共7页
苏铁熊%马理强?%刘谋斌%常建忠
囌鐵熊%馬理彊?%劉謀斌%常建忠
소철웅%마리강?%류모빈%상건충
液滴%固壁面%光滑粒子动力学%表面张力
液滴%固壁麵%光滑粒子動力學%錶麵張力
액적%고벽면%광활입자동역학%표면장력
liquid drop%solid surfaces%smoothed particle hydrodynamics%surface tension
采用改进的光滑粒子动力学(SPH)方法对液滴冲击固壁面问题进行了数值模拟.为了提高传统SPH方法的计算精度和数值稳定性,在传统的SPH方法的基础上对粒子方法中的密度和核梯度进行了修正,采用了考虑黎曼解法的SPH流体控制方程,构造了一种新型的粒子间相互作用力(IIF)模型来模拟表面张力的影响.应用改进的SPH方法对液滴冲击固壁面问题进行了数值模拟.计算结果表明:新型的IIF模型能够较好地模拟表面张力的影响,改进的SPH方法能够精细地描述液滴与固壁面相互作用过程中液滴的内部压力场演变和自由面形态变化,液滴的铺展因子随初始韦伯数的增大而增大,数值模拟结果与实验得到的结果基本一致.
採用改進的光滑粒子動力學(SPH)方法對液滴遲擊固壁麵問題進行瞭數值模擬.為瞭提高傳統SPH方法的計算精度和數值穩定性,在傳統的SPH方法的基礎上對粒子方法中的密度和覈梯度進行瞭脩正,採用瞭攷慮黎曼解法的SPH流體控製方程,構造瞭一種新型的粒子間相互作用力(IIF)模型來模擬錶麵張力的影響.應用改進的SPH方法對液滴遲擊固壁麵問題進行瞭數值模擬.計算結果錶明:新型的IIF模型能夠較好地模擬錶麵張力的影響,改進的SPH方法能夠精細地描述液滴與固壁麵相互作用過程中液滴的內部壓力場縯變和自由麵形態變化,液滴的鋪展因子隨初始韋伯數的增大而增大,數值模擬結果與實驗得到的結果基本一緻.
채용개진적광활입자동역학(SPH)방법대액적충격고벽면문제진행료수치모의.위료제고전통SPH방법적계산정도화수치은정성,재전통적SPH방법적기출상대입자방법중적밀도화핵제도진행료수정,채용료고필려만해법적SPH류체공제방정,구조료일충신형적입자간상호작용력(IIF)모형래모의표면장력적영향.응용개진적SPH방법대액적충격고벽면문제진행료수치모의.계산결과표명:신형적IIF모형능구교호지모의표면장력적영향,개진적SPH방법능구정세지묘술액적여고벽면상호작용과정중액적적내부압력장연변화자유면형태변화,액적적포전인자수초시위백수적증대이증대,수치모의결과여실험득도적결과기본일치.
@@@@In this paper, we present a numerical simulation of a single liquid drop impacting onto solid surface with smoothed particle hydrodynamics (SPH). SPH is a Lagrangian, meshfree particle method, and it is attractive in dealing with free surfaces, moving interfaces and deformable boundaries. The SPH model includes an improved approximation scheme with corrections to kernel gradient and density to improve computational accuracy. Riemann solver is adopted to solve equations of fluid motion. An new inter-particle interaction force is used for modeling the surface tension effects, and the modified SPH method is used to investigate liquid drop impacting onto solid surfaces. It is demonstrated that the inter-particle interaction force can effectively simulate the effect of surface tension. It can well describe the dynamic process of morphology evolution and the pressure field evolution with accurate and stable results. The spread factor increases with the increase of the initial Weber number. The numerical results are in good agreement with the theoretical and experimental results in the literature.