非结构网格上的三维浅水流动数值模型

针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。     

三角形网格下求解二维浅水方程的和谐Godunov格式

为保证计算格式的和谐性,通过特殊的底坡源项处理技术,在三角形网格上建立了求解二维浅水流动方程的具有空间二阶精度的Godunov格式。应用准确Riemann解求解法向数值通量,用改正的干底Riemann解处理动边界问题。经典型算例和钱塘江河口涌潮计算验证,表明模型健全,分辨率高,具有较大的推广应用价值。     

Instability of planetary flows using Riemann curvature: a numerical study

ABSTRACT

The instability of ideal non-divergent zonal flows on the sphere is determined numerically by the instability criterion of Arnold (Ann. Inst. Fourier 1966, 16, 319) for the sectional curvature. Zonal flows are unstable for all perturbations besides for a small set which are in approximate resonance. The planetary rotation is stable and the presence of rotation reduces the instability of perturbations.     

Propagation and run-up of nearshore tsunamis with HLLC approximate Riemann solver

A numerical model describing the propagation and run-up process of nearshore tsunamis in the vicinity of shorelines is developed based on an approximate Riemann solver. The governing equations of the model are the nonlinear shallow-water equations. The governing equations are discretized explicitly by using a finite volume method. The nonlinear terms in the momentum equations are solved with the Harten-Lax-van Leer-Contact (HLLC) approximate Riemann solver. The developed model is first applied to prediction of water motions in a parabolic basin, and propagation and subsequent run-up process of nearshore tsunamis around a circular island. Computed results are then compared with available analytical solutions and laboratory measurements. Very reasonable agreements are observed.     

浅水流动计算中—阶有限体积法Osher格式的实现

近年来,一阶有限体积法Osher格式已在二维浅水明流的一批模型问题和应用实例中获得成功。本文首次讨论其算法实现的种种问题。核心是建立单元水力模型-阶梯流,在数学上可用一类特殊的黎曼问题来描述。将该问题化作气体动力学中的黎曼问题近似求解,然后对结果加以校正。还在理论分析和数值试验的基础上详细讨论了各种外部边界条件、内部边界和动边界的处理,构成完整的算法。     

A Godunov-Type Scheme for Atmospheric Flows on Unstructured Grids: Euler and Navier-Stokes Equations

In recent years there has been a growing interest in using Godunov-type methods for atmospheric flow problems. Godunov's unique approach to numerical modeling of fluid flow is characterized by introducing physical reasoning in the development of the numerical scheme (van Leer, 1999). The construction of the scheme itself is based upon the physical phenomenon described by the equation sets. These finite volume discretizations are conservative and have the ability to resolve regions of steep gradients accurately, thus avoiding dispersion errors in the solution. Positivity of scalars (an important factor when considering the transport of microphysical quantities) is also guaranteed by applying the total variation diminishing condition appropriately. This paper describes the implementation of a Godunov-type finite volume scheme based on unstructured adaptive grids for simulating flows on the meso-, micro- and urban-scales. The Harten-Lax-van Leer-Contact (HLLC) approximate Riemann solver used to calculate the Godunov fluxes is described in detail. The higher-order spatial accuracy is achieved via gradient reconstruction techniques after van Leer and the total variation diminishing condition is enforced with the aid of slope-limiters. A multi-stage explicit Runge-Kutta time marching scheme is used for maintaining higher-order accuracy in time. The scheme is conservative and exhibits minimal numerical dispersion and diffusion. The subgrid scale diffusion in the model is parameterized via the Smagorinsky-Lilly turbulence closure. The scheme uses a non-staggered mesh arrangement of variables (all quantities are cell-centered) and requires no explicit filtering for stability. A comparison with exact solutions shows that the scheme can resolve the different types of wave structures admitted by the atmospheric flow equation set. A qualitative evaluation for an idealized test case of convection in a neutral atmosphere is also presented. The scheme was able to simulate the onset of Kelvin-Helmholtz type instability and shows promise in simulating atmospheric flows characterized by sharp gradients without using explicit filtering for numerical stability.     

Sensitivity equations for hyperbolic conservation law-based flow models

The present paper focuses on the governing equations for the sensitivity of the variables to the parameters in flow models that can be described by one-dimensional scalar, hyperbolic conservation laws. The sensitivity is shown to obey a hyperbolic, scalar conservation law. The sensitivity is a conserved scalar except in the case of discontinuous flow solutions, where an extra, point source term must be added to the equations in order to enforce conservation. The propagation speed of the sensitivity waves being identical to that of the conserved variable in the original conservation law, the system of conservation laws formed by the original hyperbolic equation and the equation satisfied by the sensitivity is linearly degenerate. A consequence on the solution of the Riemann problem is that rarefaction waves for the variable of the original equation result in vacuum regions for the sensitivity. The numerical solution of the hyperbolic conservation law for the sensitivity by finite volume methods requires the implementation of a specific shock detection procedure. A set of necessary conditions is defined for the discretisation of the source term in the sensitivity equation. An application to the one-dimensional kinematic wave equation shows that the proposed numerical technique allows analytical solutions to be reproduced correctly. The computational examples show that first-order numerical schemes do not yield satisfactory numerical solutions in the neighbourhood of moving shocks and that higher-order schemes, such as the MUSCL scheme, should be used for sharp transients.     

ON LONG-WAVE BREAKING Ⅳ. NUMERICAL RESULTS

ImODUcrIONBreakingwavesaretheagentSformanyimPohantupperoceanproassesinvolvingtransferofhobontalmornentumfromwindwivestosurfaceimtS.AfaIniliarandspancularpropertyoflongwivespropagatinginaninfiniteorhaif4nfinitechannelwithbottomsl0peisthemeCanismofwavebrmking.Alltheusualtheories(srnallamPlitudeapproxthetionsofAiryandStokes,noTilinearshallowwitertheory,K0rteweg-DeVireSequattonsforsolitaryandcnoidalWhves)areessentiallyapproxirnations,validonlyWhenthefluidadetionissuffidenhysmallcomPatalto…     

基于精确Riemann求解器的明满流过渡过程模拟

Preissmann窄缝法模拟明满流过渡过程方法简单,但存在明显的非物理振荡,抑制非物理振荡是该方法应用的关键。基于Godunov格式和精确Riemann求解器对明满流过渡过程进行模拟,针对Riemann问题代数恒等式在明满流交界处不光滑问题,提出了三阶收敛方法与二分法结合的迭代求解方法,保证迭代收敛至真实解;针对由于变量空间重构方法不能准确表达变量在空间中真实物理状态而导致的非物理振荡,提出了基于精确Riemann解的变量空间重构方法,准确表达激波间断在单元内的空间分布状态,从机理上抑制了非物理振荡。实例研究表明,数值计算结果与解析解或实测值吻合良好,研究成果为明满流过渡过程的高精度数值模拟提供了新的方法。     

An approximate Riemann solver for sensitivity equations with discontinuous solutions

The sensitivity of a model output (called a variable) to a parameter can be defined as the partial derivative of the variable with respect to the parameter. When the governing equations are not differentiable with respect to this parameter, problems arise in the numerical solution of the sensitivity equations, such as locally infinite values or instability. An approximate Riemann solver is thus proposed for direct sensitivity calculation for hyperbolic systems of conservation laws in the presence of discontinuous solutions. The proposed approach uses an extra source term in the form of a Dirac function to restore sensitivity balance across the shocks. It is valid for systems such as the Euler equations for gas dynamics or the shallow water equations for free surface flow. The method is first detailed and its application to the shallow water equations is proposed, with some test cases such as dike- or dam-break problems with or without source terms. An application to a two-dimensional flow problem illustrates the superiority of direct sensitivity calculation over the classical empirical approach.