Symmetries of the hyperbolic shallow water equations and the Green–Naghdi model in Lagrangian coordinates

The observation that the hyperbolic shallow water equations and the Green–Naghdi equations in Lagrangian coordinates have the form of an Euler–Lagrange equation with a natural Lagrangian allows us to apply Noether's theorem for constructing conservation laws for these equations. In this study the complete group analysis of these equations is given: admitted Lie groups of point and contact transformations, classification of the point symmetries and all invariant solutions are studied. For the hyperbolic shallow water equations new conservation laws which have no analog in Eulerian coordinates are obtained. Using Noether's theorem a new conservation law of the Green–Naghdi equations is found. The dependence of solutions on the parameter is illustrated by self-similar solutions which are invariant solutions of both models.     

A FINITE VOLUME SOLVER FOR THE SHALLOW WATER EQUATIONS IN VARYING BED ENVIRONMENTS

Abstract

A numerical scheme for solving the shallow-water equations is presented. An analogy is made between flows governed by shallow-water equations and the Euler system of equations used in gas dynamics. An emphasis is placed on the difference presented by the bathymetry in hydraulic systems. The discretization of the governing equations is based on Roe's flux difference-splitting solver, initially developed for solving inviscid compressible flows. The spatial discretization is handled within a finite-volume context by using triangles or quadrilaterals as the basic control-volume cells. This approach enables an easy and flexible treatment of general geometries. A development of the boundary conditions tailored for the current scheme is given. Fundamental validation tests are presented.     

A weighted least squares method for first-order hyperbolic systems

The paper presents a generalization of the classical L²-norm weighted least squares method for the numerical solution of a first-order hyperbolic system. This alternative least squares method consists of the minimization of the weighted sum of the L² residuals for each equation of the system. The order of accuracy of global conservation of each equation of the system is shown to be inversely proportional to the weight associated with the equation. The optimal relative weights between the equations are then determined in order to satisfy global conservation of the energy of the physical system. As an application of the algorithm, the shallow water equations on an irregular domain are first discretized in time and then solved using Laplace modification and the proposed least squares method.     

Roll waves structure in two-layer Hele–Shaw flows

The mathematical model of inhomogeneous fluid motion in a Hele–Shaw cell is proposed. Based on this model the equations for describing two-layer flows and development of roll waves at the interface are derived. Conditions of roll waves existence are formulated in terms of Whitham criterion. Numerical calculations of the interface position are provided. It is shown that small perturbations of the interface in the inlet section of the channel lead to the roll waves for certain parameters of the flow. Two-parametric class of exact solutions corresponding to the roll waves regime is obtained. Diagrams of critical depths of roll waves development are constructed.     

Water Vapor Diffusion Through Soil as Affected by Temperature and Aggregate Size

Water vapor diffusion through the soil is an important part in the total water flux in the unsaturated zone of arid or semiarid regions and has several significant agricultural and engineering applications because soil moisture contents near the surface are relatively low. Water vapor diffusing through dry soil is absorbed for both long and short terms. Long-term absorption allows more water to enter than exit the soil, as reflected in the concentration gradient. Short-term absorption leads to an apparent reduction in the diffusion rate, as reflected in the diffusion coefficient. This investigation studied the effects of soil temperature and porosity on the isothermal diffusion of water vapor through soil. The diffusion model consisted of 25.4 cm × 8.9 cm × 20.3 cm Plexiglas box divided into two compartments by a partition holding a soil reservoir. Water vapor moved from a container suspended by a spring in one compartment, through the porous medium in the center of the model, to calcium chloride in a container suspended by a spring in the other compartment. The porous materials consisted of aggregates of varying size (2–2.8, 1–2, and 0.5–1 mm) of a Fayatte silty clay loam (a fine-silty, mixed mesic Typic Hapludalf). The flow rates of water vapor were measured at temperatures of 10, 20, 30, and 40°C. Warmer temperatures increased the rate of diffusion through dry soil while reduced the amount of water absorbed by that soil. Reducing porosity slowed the rate of diffusion and increased the amount of water absorbed. The dry soil in this study absorbed from 1/8 to 2/3 of the diffusing water. Maximum absorption rates occurred with the most compact soil samples at the highest temperature, though the maximum absorption as a percentage of the diffusing water was in the compact samples at the lowest temperature. The diffusivity equation D/D ₀ = [(S – 0.1)/0.9]² fit the D/D ₀ values obtained from these data if a coefficient of 1/3 or 1/3.5 is added to correct for the time delays caused by temporary sorption of the diffusing water vapor. The data, influenced by the interaction of water vapor and soil materials, represent a diffusion rate lower than the diffusion rate that would have resulted without this interaction. Mention of trade names, proprietary products, or specific equipment is intended for reader information only and does not constitute a guarantee or warranty by the USDA-ARS nor does it imply approval of the product named to the exclusion of other products. An erratum to this article can be found at     

Diffusion models with microstructure

Models of diffusion are presented which recognize the local geometry of individual cells or storage sites and the exchange of flux on the micro-scale of these cells. Such models have been obtained by homogenization, but here we indicate stronger existence-uniqueness results of parabolic type can be obtained directly. Connections between these models and their historical development will be described.     

Numerical solution of the generalized Serre equations with the MacCormack finite-difference scheme

This paper describes a two-dimensional numerical model to solve the generalized Serre equations. In order to solve the system equations, written in the conservative form, we use an explicit finite-difference method based on the MacCormack time-splitting scheme. The numerical method and the computational model are validated by comparing one- and two-dimensional numerical solutions with theoretical and experimental results. Finally, the two-dimensional model (in a horizontal plane) is tested in a domain with complicated boundary conditions.     

Adaptive Q‐tree Godunov‐type scheme for shallow water equations

This paper presents details of a second‐order accurate, Godunov‐type numerical model of the two‐dimensional shallow water equations (SWEs) written in matrix form and discretized using finite volumes. Roe's flux function is used for the convection terms and a non‐linear limiter is applied to prevent unwanted spurious oscillations. A new mathematical formulation is presented, which inherently balances flux gradient and source terms. It is, therefore, suitable for cases where the bathymetry is non‐uniform, unlike other formulations given in the literature based on Roe's approximate Riemann solver. The model is based on hierarchical quadtree (Q‐tree) grids, which adapt to inherent flow parameters, such as magnitude of the free surface gradient and depth‐averaged vorticity. Validation tests include wind‐induced circulation in a dish‐shaped basin, two‐dimensional frictionless rectangular and circular dam‐breaks, an oblique hydraulic jump, and jet‐forced flow in a circular reservoir. Copyright © 2001 John Wiley & Sons, Ltd.     

An unsymmetric constitutive equation for anisotropic viscoelastic fluid

A continuum constitutive theory of corotational derivative type is developed for the anisotropic viscoelastic fluid–liquid crystalline (LC) polymers. A concept of anisotropic viscoelastic simple fluid is introduced. The stress tensor instead of the velocity gradient tensor D in the classic Leslie–Ericksen theory is described by the first Rivlin–Ericksen tensor A and a spin tensor W measured with respect to a co-rotational coordinate system. A model LCP-H on this theory is proposed and the characteristic unsymmetric behaviour of the shear stress is predicted for LC polymer liquids. Two shear stresses thereby in shear flow of LC polymer liquids lead to internal vortex flow and rotational flow. The conclusion could be of theoretical meaning for the modern liquid crystalline display technology. By using the equation, extrusion–extensional flows of the fluid are studied for fiber spinning of LC polymer melts, the elongational viscosity vs. extension rate with variation of shear rate is given in figures. A considerable increase of elongational viscosity and bifurcation behaviour are observed when the orientational motion of the director vector is considered. The contraction of extrudate of LC polymer melts is caused by the high elongational viscosity. For anisotropic viscoelastic fluids, an important advance has been made in the investigation on the constitutive equation on the basis of which a series of new anisotropic non-Newtonian fluid problems can be addressed. The project supported by the National Natural Science Foundation of China (10372100, 19832050) (Key project). The English text was polished by Yunming Chen.     

Thermodynamic analysis and numerical resolution of a turbulent - fully ionized plasma flow model

This paper is devoted to the modeling and numerical resolution of a non-dissipative compressible turbulent plasma flow model involving three temperatures (turbulence, ions and electrons). The first step is to derive such a model. To do this, an analysis of the Reynolds averaged Euler equations (the k-model) is carried out. It is shown that thermodynamic requirements enable the derivation of an equation of state for turbulent variables. This equation of state is of the same type as those of an ideal gas. In this context, the various thermodynamic variables of turbulence can be obtained (energy, pressure, temperature etc.). This hyperbolic conservative model has exactly the same structure as the two temperatures plasma model of Zeldovich. Thanks to the clear structure of these two models, the turbulent plasma model is derived and involves three temperatures. The second step is to derive an accurate numerical scheme for its solution. A linearized Riemann solver and a positive HLLC type solver are derived and embedded into a conventional Godunov scheme. It is shown that this method requires important corrections to preserve contact discontinuities and temperatures monotonicity. The corrections are based upon a non-conservative formulation of the turbulence and electrons energy equations, while total energy conservation is preserved. The modified method behaves correctly with contacts and shocks.Received: 11 February 2002, Revised: 19 June 2003, Accepted: 7 August 2003, Published online: 14 October 2003 Correspondence to: R. Saurel     

海岸波浪场模型研究进展

从建模原理、波浪在近岸区域传播的众多机制、模型的类别、优势、局限性以及模型在未来的发展趋势等方面,综述了在海岸工程实践中广泛运用的以下两大类海岸波浪场预测模型的最新研究进展:(1)能量平衡模型.它一般用来预测海洋深水波候,已发展到相当完善的阶段,例如,最为著名的WAM3G模型.这种模型在海岸工程中的作用就在于可以模拟施加在波浪上的随时间变化的风场效应.(2)质量、动量守恒模型.它在海岸工程中应用最为普遍,并且内容丰富,数值技巧多样.目前包含了以下代表性的模型:缓坡方程、抛物型方程、非线性浅水方程、高阶Boussinesq型方程、Green-Naghdi理论.      

Traffic Flow Models with Phase Transitions

The theory of hyperbolic conservation laws has been successfully applied to the study of vehicular traffic flows. We present here some models showing phase transitions, that in terms of traffic flows correspond to two distinct behaviors, free or congested.