主权项 |
一种非结构化网格嵌套波浪数值模拟方法,其特征是,包括如下步骤:1)建立大区域的非结构化网格波浪数学模型,其控制方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mi>N</mi><mo>+</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac><msub><mi>C</mi><mi>x</mi></msub><mi>N</mi><mo>+</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac><msub><mi>C</mi><mi>y</mi></msub><mi>N</mi><mo>+</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>σ</mi></mrow></mfrac><msub><mi>C</mi><mi>σ</mi></msub><mi>N</mi><mo>+</mo><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>θ</mi></mrow></mfrac><msub><mi>C</mi><mi>θ</mi></msub><mi>N</mi><mo>=</mo><mfrac><mi>s</mi><mi>σ</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581361430000011.GIF" wi="1440" he="109" /></maths>上式左边第一项为波作用量N随时间的变化率,N=N(σ,θ)=E(σ,θ)/σ,E(σ,θ)为能谱密度,σ为波浪频率,θ为波浪方向;上式左边第二项和第三项分别表示N(σ,θ)在空间x、y方向上的传播;上式左边第四项表示N(σ,θ)在σ空间因为流场和水深所引起的变化;上式左边第五项表示N(σ,θ)在θ空间的传播,亦即水深及流场所引起的折射;方程右边的S代表以谱密度表示的源汇项,包括风能输入、波与波之间非线性相互作用和由于底摩擦、白浪、深度诱导破碎所引起的能量损耗;所述C<sub>x</sub>、C<sub>y</sub>、C<sub>σ</sub>、C<sub>θ</sub>分别代表在x、y、σ、θ空间的波浪传播速度;其中,C<sub>x</sub>、C<sub>y</sub>、C<sub>σ</sub>、C<sub>θ</sub>分别为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>C</mi><mi>x</mi></msub><mo>=</mo><mfrac><mi>dx</mi><mi>dt</mi></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><mn>1</mn><mo>+</mo><mfrac><mrow><mn>2</mn><mi>kd</mi></mrow><mrow><mi>sinh</mi><mrow><mo>(</mo><mi>skd</mi><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mfrac><mrow><mi>σ</mi><msub><mi>k</mi><mi>x</mi></msub></mrow><msup><mi>k</mi><mn>2</mn></msup></mfrac><mo>+</mo><msub><mi>U</mi><mi>x</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581361430000012.GIF" wi="1397" he="126" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>C</mi><mi>y</mi></msub><mo>=</mo><mfrac><mi>dy</mi><mi>dt</mi></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><mn>1</mn><mo>+</mo><mfrac><mrow><mn>2</mn><mi>kd</mi></mrow><mrow><mi>sinh</mi><mrow><mo>(</mo><mn>2</mn><mi>kd</mi><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mfrac><mrow><mi>σ</mi><msub><mi>k</mi><mi>y</mi></msub></mrow><msup><mi>k</mi><mn>2</mn></msup></mfrac><mo>+</mo><msub><mi>U</mi><mi>y</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581361430000013.GIF" wi="1396" he="147" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>C</mi><mi>σ</mi></msub><mo>=</mo><mfrac><mi>dσ</mi><mi>dt</mi></mfrac><mo>=</mo><mfrac><mrow><mo>∂</mo><mi>σ</mi></mrow><mrow><mo>∂</mo><mi>d</mi></mrow></mfrac><mo>[</mo><mfrac><mrow><mo>∂</mo><mi>d</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mover><mi>U</mi><mo>→</mo></mover><mo>·</mo><mo>▿</mo><mi>d</mi><mo>]</mo><mo>-</mo><msub><mi>C</mi><mi>g</mi></msub><mover><mi>k</mi><mo>→</mo></mover><mo>·</mo><mfrac><mrow><mo>∂</mo><mover><mi>U</mi><mo>→</mo></mover></mrow><mrow><mo>∂</mo><mi>s</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581361430000014.GIF" wi="1397" he="163" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>C</mi><mi>θ</mi></msub><mo>=</mo><mfrac><mi>dθ</mi><mi>dt</mi></mfrac><mo>=</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><mo>[</mo><mfrac><mrow><mo>∂</mo><mi>σ</mi></mrow><mrow><mo>∂</mo><mi>d</mi></mrow></mfrac><mfrac><mrow><mo>∂</mo><mi>d</mi></mrow><mrow><mo>∂</mo><mi>m</mi></mrow></mfrac><mo>+</mo><mover><mi>k</mi><mo>→</mo></mover><mo>·</mo><mfrac><mover><mi>U</mi><mo>→</mo></mover><mrow><mo>∂</mo><mi>m</mi></mrow></mfrac><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581361430000015.GIF" wi="1394" he="182" /></maths>其中,<img file="FDA0000581361430000016.GIF" wi="262" he="85" />为波数,k<sub>x</sub>,k<sub>y</sub>分别为x、y空间的分量,k为波数,d为水深,<img file="FDA0000581361430000017.GIF" wi="272" he="89" />为流速,U<sub>x</sub>,u<sub>y</sub>分别为x、y空间的2分量,s为沿θ方向空间坐标,m为垂直于s的坐标,相对频率<img file="FDA0000581361430000018.GIF" wi="327" he="82" />ω为波浪的固有频率;算子<img file="FDA0000581361430000019.GIF" wi="50" he="79" />定义为<img file="FDA0000581361430000021.GIF" wi="396" he="96" /><img file="FDA0000581361430000022.GIF" wi="48" he="85" />为波速;t为时间。2)建立非结构化网格嵌套的浅海工程海域波浪数学模型;根据所述大区域非结构化网格波浪数学模型提供给所述浅海工程海域波浪数学模型的完全嵌套的频谱和方向谱边界条件(包括将要嵌套的计算域边界处节点的坐标、频率、方向、能谱或变密度等要素),然后应用所述步骤1)中的控制方程和有关的辅助方程,建立浅海工程海域波浪数学模型,进行工程海域波浪数值计算。3)求解所述非结构化网格嵌套波浪数学模型。 |