1、自然環(huán)境水體中的分層流動(dòng)通常是由溫度、主動(dòng)輸運(yùn)物質(zhì)(如鹽度、泥沙等)或者它們的組合作用誘導(dǎo)產(chǎn)生的。由于各組分對(duì)水體密度的貢獻(xiàn)以及其自身的擴(kuò)散特性的差異使得該類分層流動(dòng)更為復(fù)雜。本文從兩流體的角度，采用。Level Set方法追蹤兩流體分層界面的運(yùn)動(dòng)，對(duì)每相流體分別求解，并考慮界面處兩流體間的相互作用，模擬了污水排海等分層流問題，詳細(xì)闡述了海水入侵過程和機(jī)理。本文研究的主要成果包括： (1)采用Level Set方法追蹤界面的時(shí)空

2、變化。本文基于快速行進(jìn)法和源點(diǎn)掃描法的思想，提出了非結(jié)構(gòu)同位網(wǎng)格下Level Set的實(shí)現(xiàn)策略。并通過實(shí)例，對(duì)Level Set法與VOF法作了比較分析。 (2)界面附近網(wǎng)格點(diǎn)密度、擴(kuò)散系數(shù)等物理參數(shù)處理，二維三維采用了不同的策略。對(duì)二維情形，根據(jù)符號(hào)距離函數(shù)精確計(jì)算各網(wǎng)格單元的相體積分布，而后可計(jì)算網(wǎng)格單元的密度和擴(kuò)散系數(shù)；對(duì)三維情形，直接采用正弦函數(shù)光滑處理。對(duì)于界面物理量的擴(kuò)展采用一階的isobaric fix ghost

3、處理。 (3)將密度方程和動(dòng)量方程統(tǒng)一成一個(gè)積分方程，采用梯度重構(gòu)策略，結(jié)合限制器，改進(jìn)了兩流體的SIMPLE(SIMPLEC)算法和人工壓縮算法，建立了非結(jié)構(gòu)網(wǎng)格下N-S方程的求解系統(tǒng)。改進(jìn)了出口邊界的處理策略，能有效地處理“無反射”出口邊界問題。對(duì)不同的密度比的兩流體系統(tǒng)進(jìn)行了驗(yàn)證分析。 (4)建立了k-ε紊流模型，壁面采用對(duì)數(shù)率壁函數(shù)處理技術(shù)。 (5)從兩流體角度對(duì)污水排海問題建立了二維和三維數(shù)學(xué)模型，詳細(xì)