Comprehensive two-dimensional river ice model based on boundary-fitted coordinate transformation method

River ice is a natural phenomenon in cold regions, influenced by meteorology, geomorphology, and hydraulic conditions. River ice processes involve complex interactions between hydrodynamic, mechanical, and thermal processes, and they are also influenced by weather and hydrologic conditions. Because natural rivers are serpentine, with bends, narrows, and straight reaches, the commonly-used one-dimensional river ice models and two-dimensional models based on the rectangular Cartesian coordinates are incapable of simulating the physical phenomena accurately. In order to accurately simulate the complicated river geometry and overcome the difficulties of numerical simulation resulting from both complex boundaries and differences between length and width scales, a two-dimensional river ice numerical model based on a boundary-fitted coordinate transformation method was developed. The presented model considers the influence of the frazil ice accumulation under ice cover and the shape of the leading edge of ice cover during the freezing process. The model is capable of determining the velocity field, the distribution of water temperature, the concentration distribution of frazil ice, the transport of floating ice, the progression, stability, and thawing of ice cover, and the transport, accumulation, and erosion of ice under ice cover. A MacCormack scheme was used to solve the equations numerically. The model was validated with field observations from the Hequ Reach of the Yellow River. Comparison of simulation results with field data indicates that the model is capable of simulating the river ice process with high accuracy.

Numerical Modeling of the Hyperbolic Mild-Slope Equation in Curvilinear Coordinates

The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccurate in areas with irregular shorelines, such as estuaries and harbors. Based on the hyperbolic mild-slope equation in Cartesian coordinates, the numerical model in orthogonal curvilinear coordinates is developed. The transformed model is discretized by the finite difference method and solved by the ADI method with space-staggered grids. The numerical predictions in curvilinear co- ordinates show good agreemenl with the data obtained in three typical physical expedments, which demonstrates that the present model can be used to simulate wave propagation, for normal incidence and oblique incidence, in domains with complicated topography and boundary conditions.

Numerical Simulation of Wave-Induced Currents Combined with Parabolic Mild-Slope Equation in Curvilinear Coordinates

Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.

3D Numerical Simulation of Overbank Flow in Non-Orthogonal Curvilinear Coordinates

The velocity field in meandering compound channels with overhank flow is highly three dimensional. To date, its features have been investigated experimentally and little research has been undertaken to investigate the feasibility of reproducing these velocity fields with computer models. If computer modeling were to prove successful in this context, it could become a useful prediction technique and research tool to enhance our understanding of natural river dynamics. A 3-D k-E turbulence hydrodynamic model in curvilinear coordinates is established to simulate the overhank flow. The bodyfitted coordinate is adopted in the horizontal plane, the part grid is adopted in the vertical direction, and the wall-function method is employed to simulate the bed resistance. The model is applied to the simulation of the meandering channel with straight flood plain banks, and the main velocities and secondary velocities for both the longitudinal and cross sections are presented. Comparison and analysis show that the results of simulation are fit to reflect the results of experiment. These results show the application value of the model to 3D overhank flow.

A Numerical Model of Nonlinear Wave Propagation in Curvilinear Coordinates

For the simulation of the nonlinear wave propagation in coastal areas with complex boundaries, a numerical model is developed in curvilinear coordinates. In the model, the Boussinesq-type equations including the dissipation terms are em- ployed as the governing equations. In the present model, the dependent variables of the transformed equations are the free surface elevation and the utility velocity variables, instead of the usual primitive velocity variables. The introduction of utility velocity variables which are the products of the contravariant components of the velocity vector and the Jacobi ma- trix can make the transformed equations relatively concise, the treatment of lateral boundary conditions easier and the de- velopment of the program simpler. The predictor-corrector method and five-point finite-difference scheme are employed to discretize the time derivatives and the spatial ones, respectively. The numerical model is tested for three cases. It is found that the numerical results are in good agreement with the analytical results and experimental data.

SPACE-TIME CONSERVATION SCHEME IN NONORTHOGONAL CURVILINEAR COORDINATES

A scheme of space time conservation (STC) based on the method of space time conservation element and solution element (CE/SE) is represented in a nonorthogonal curvilinear coordinate system. The corresponding initial and boundary conditions are discussed. It is seen that in the nonorthogonal coordinates the scheme maintains the advantages of the STC method, and is noted for its simple structure, clear physical meaning, rapid calculation and high accuracy. It is easy to extend to the multidimensional flow. The numerical results for a 2D Euler equation show good agreement with those from other computational methods and the experiment.

Numerical model of flow and pollutant transportation in non-orthogonal curvilinear coordinates

The planar 2D k-ε double equations' turbulence model was adopted and transformed into non-orthogonal curvilinear coordinates. The concentration convection-diffusion was introduced to planar 2D SIMPLEC algorithm of flow in non-orthogonal curvilinear coordinates. The numerical model of pollutant transportation in non-orthogonal curvilinear coordinates was constructed. The model was applied to simulate the flow and pollutant concentration fields. In the testing concentration field, two optimal operations of contamination discharging both along bank and in the centerline at the first bend of the meandering channel were adopted. Comparison with available data showed the model developed was successful, was valuable to engineering application.

Curve integral with path independent in orthogonal curvilinear coordinate system

It is explored that the line integral is a path independent in two or three arbitrary dimensional orthogonal curvilinear coordinate systems, which is based on the integral condition with the path independent in two or three dimensional rectangular coordinate systems. Firstly, according to the coordinate transformation, the condition that the line integral is the path independent in the polar coordinate system is obtained easily from the Green's theorem in two-dimensional rectangular coordinate system and the condition is extended to arbitrary two-dimension orthogonal curvilinear coordinates. Secondly, through the coordinate transformation relationship and the area projection method, the Stokes formula in three-dimensional rectangular coordinate system is promoted to the spherical coordinate system and cylindrical coordinate system, and the condition that the line integral is a path independent is obtained. Furthermore, the condition is extended to arbitrary three-dimension orthogonal curvilinear coordinates. Lastly, the conclusions are made.

Expression of strain tensor in orthogonal curvilinear coordinates

Based on an analysis of connotation and extension of the concept of the orthogonal curvilinear coordinates, we have deduced a platform of strain tensor expression of Cartesian coordinates, which turns out to be a function of Lame coefficient and unit vector. By using transform matrix between Cartesian coordinates and orthogonal eurvilinear coordinates, we have deduced a mathematical expression for correcting displacement vector differential in orthogonal curvilinear coordinates, and given a general expression of strain tensor in orthogonal curvilinear coordinates.

Generalization of the Second Order Vector Potential Formulation for Arbitrary Non-Orthogonal Curvilinear Coordinates Systems from the Covariant Form of Maxwell's Equations

A great number of semi-analytical models, notably the representation of electromagnetic fields by integral equations are based on the second order vector potential (SOVP) formalism which introduces two scalar potentials in order to obtain analytical expressions of the electromagnetic fields from the two potentials. However, the scalar decomposition is often known for canonical coordinate systems. This paper aims in introducing a specific SOVP formulation dedicated to arbitrary non-orthogonal curvilinear coordinates systems. The electromagnetic field representation which is derived in this paper constitutes the key stone for the development of semi-analytical models for solving some eddy currents moelling problems and electromagnetic radiation problems considering at least two homogeneous media separated by a rough interface. This SOVP formulation is derived from the tensor formalism and Maxwell’s equations written in a non-orthogonal coordinates system adapted to a surface characterized by a 2D arbitrary aperiodic profile.

Electrodynamics in Curvilinear Coordinates and the Equation of a Geodesic Line

The creation of the theory of relativity, which discovered the equivalence of mass and energy, showed that the concept of a point charge, used in the formulation of Coulomb’s law, one of the basic laws of classical electrodynamics, contradicts the famous formula establishing the equivalence of mass and energy. But the discovery of quarks makes it possible to present classical electrodynamics in a form free from the indicated contradiction. In the article, having considered the electromagnetic field in a curvilinear coordinate system, a theory has been created that expands our understanding of the electromagnetic field, the nature of quarks, the nature of strong interaction, and the connection between strong interaction and electromagnetic interaction. This theory is based on the principle of equivalence of an electromagnetic field to a free material particle formulated in the article and the law of formation of elementary particles from an electromagnetic field that follows from it.

Staggered Calculation of Turbulent Flows Over a 3D Bluff Body Using Curvilinear Nonorthogonal Coordinates

A finite difference method is developed to predict turbulent flows over 3D bluffbodies. The K-ε turbulence model with Launder and Spalding's wall treatment isemployed. The solution alsorithm is based on a body fitted nonorthogonalcurvilinear eourdinate system and a stagsered grid arrangement. The covariantvelocity components are chosen as dependent variables. Convective fluxes aredescribed by the Power haw Scheme. The grids are generated with an ellipticgrid generator using control functions. Results obtained are compared withexporiment measurements and other calculations.

Algebraic-Stress Turbulent Model of Planar 2-D Pollutant Transportation in Curvilinear Coordinates

The shallow-water equations and pollutant convective-diffusive equation are transformed into curvilinear coordinate system. The anisotropic algebraic-stress turbulent model is introduced to simulate the turbulence items, and algebraic-stress turbulent model of planar 2-D pollutant convective-diffusive in curvilinear coordinates is build. The meandering channel with measured data of concentration in lab is adopted to validate the model, and the distribution figure of pollutant concentration field calculated through this model and that of the k-ε model, which show the model is superior to k-ε turbulent model in dealing with anisotropy character of flow.

THE MACCORMACK SCHEME AND THE BOUNDARY-FITTED COORDINATE SYSTEM FOR NUMERICAL SIMULATION OF SUPER CRITICAL FREE SURFACE FLOWS

This paper is concerned with the mathematical model for the numerical simulation of supercritical surface flows. A boundary-fitted coordinate system was used to overcome the difficulties and inaccuracy associated with the determination of flow characteristics near the fl ow boundaries. The MacCormack scheme was applied for the solution of the transfo rmed system of equations. Comparisions between computed results and experimental data show a satisfactory agreement.

Computing Depth-averaged Flows Using Boundary-fitted Coordinates and Staggered Grids

A depth averaged nonlinear k ε model for turbulent flows in complex geometries has been developed in a boundary fitted coordinate system. The SIMPLEC procedure is used to develop an economical discrete method for staggered grids to analyze flows in a 90° bend. This paper describes how to change a program in rectangular coordinate into a boundary fitted coordinate. The results compare well with experimental data for flow in a meandering channel showing the efficiency of the model and the discrete method.

仿射空间拉普拉斯算子探讨

为了进一步研究提升惯性导航系统的精度,对拉普拉斯算子进行讨论:拉普拉斯算子是个微分算子,拉普拉斯方程又名调和方程、位势方程,求解拉普拉斯方程是物理学和力学等领域经常遇到的一类重要数学问题;基于曲线坐标系和张量理论,分别给出n维仿射空间中的拉普拉斯算子表达式,并证明不同表述的等价性;最后,基于格林定理和变方方法,分别给出3维仿射空间拉普拉斯算子的表达式。

NUMERICAL SIMULATION OF 2D HORIZONTAL COOLING WATER DISCHARGE IN GENERALIZED CURVILINEAR COORDINATE

2D horizontal model is one of the major mathematical methods for the research on cooling water discharge from the power plant. In this paper, the shallow water equations are transformed under the generalized curvilinear coordinate system and the elliptic differential equations are used to generate curvilinear grids, so a model in generalized curviline ar coordinate is presented to simulate 2D horizontal cooling water, Governing equations of the model are discretized by finite volume method, and non-staggered grids and SIMPLE algorithm are introduced to simplify the program during the discretization. This model is used to simulate the movement of cooling water in a simplified meandering channel and a natural channel, calculating results indicate this model can correctly reflect the movement rules of cooling water, which verifies the model can be applied in engineering practice.

MATHEMATICAL MODEL FOR 2-D TIDAL FLOW AND WATER QUALITY WITH ORTHOGONAL CURVILINEAR COORDINATES

This paper presents a numerical method for simulating the 2-D tidal flow andwater quality with the orthogonal curvilinear coordinates. In order to overcome the computationaldifficulties in natural rivers, such as the complicated boundary figures, the great disparitybetween length and width of computational domain, etc. , orthogonal boundary-filled grid was used.The irregular domain in physical plane was transformed into a rectangular domain in a transformedplane, and the depth-averaged momentum equations and mass equation were given and discretized basedon the alternating direction implicit finite difference scheme in curvilinear coordinates. Theapplication of the presented method was illustrated by an example of analyzing the Yangtze River inthe vicinity of Nanjing city. A fair agreement between the measured data and computed resultsdemonstrates the validity of the developed method.

NUMERICAL SIMULATION OF 2-D TIDAL FLOW AND WATER QUALITY UNDER THE CURVILINEAR COORDINATES

This paper presents a numerical method to simulate the 2-D tidal flow and water quality under the curvilinear coordinates. In order to overcome the computational difficulties in natural rivers, such as the complicated boundary figures, the great disparity between length and width of computational domain, etc. , boundary-fitted grid is used, the irregular domain in physical plane is transformed into a rectangular domain in transformed plane, and the depth-averaged momentum equations and mass equation are rewritten and discretized based on the finite volume techniques in curvilinear coordinates. Practical application of the method is illustrated by an example for the Dachangzhen Section of the Yangtze River. A fair agreement between the values measured and computed demonstrates the validity of the method developed.

A NON-UNIFORM SEDIMENT TRANSPORT MODEL WITH THE BOUNDARY-FITTING ORTHOGONAL COORDINATE SYSTEM

A 2-D non-uniform sediment mat hmatical model in the boundary-fitting orthogonal coordinate system was develop ed in this paper. The governing equations, the numerical scheme, the boundary con ditions, the movable boundary technique and the numerical solutions were present ed. The model was verified by the data of the reach 25km upstream the Jialingjia ng estuary and the 44km long main stream of the Chongqing reach of the Yangtze r iver. The calculated results show that, the water elevation, the velocity distri bution and the river bed deformation are in agreement with the measured data.