An implicit method using contravariant velocity components and its application to calculations in a harbour-channel area

The key problem in the computation of fluid dynamics using fine boundary-fitted grids is how to improve the numerical stability and decrease the calculating quantity. To solve this problem, implicit schemes should be adopted since explicit schemes may bring about a great increase in computation quantity according to the Courant-FrledrichsLewy condition. Whereas the adoption of implicit schemes is difficult to be realized because of the existence of two partial derivatives of surface elevations with respect to variables of alternative direction coordinates in each momentum equation in non-rectangular coordinates. With an aim to design an implicit scheme in non-reetangular ccordinates in the present paper, new momentum equations with the contravariant components of velocity vector are derived based on the shallow water dynamic equations in generalized curvilinear coordinates. In each equation, the coefficients before the two detivatives of surface elevations have different orders of magnitude, i. e., the derivative with the larger ceefficient rnay play a more important role than that with the smaller one. With this advantage, the ADI scheme can then be easily employed to improve the numerical stability and decrease the calculating quantity. The calculation in a harbour and a channel in Macau nearshore area shows that the implicit model is effective in calculating current fields in small size areas.

Implicit numerical scheme based on SMAC method for unsteady incompressible Navier-Stokes equations

An implicit numerical scheme is developed based on the simplified marker and cell （SMAC） method to solve Reynolds-averaged equations in general curvilinear coordinates for three-dimensional （3-D） unsteady incompressible turbulent flow. The governing equations include the Reynolds-averaged momentum equations, in which contravariant velocities are unknown variables, pressure-correction Poisson equation and k- s turbulent equations. The governing equations are discretized in a 3-D MAC staggered grid system. To improve the numerical stability of the implicit SMAC scheme, the higherorder high-resolution Chakravarthy-Osher total variation diminishing （TVD） scheme is used to discretize the convective terms in momentum equations and k- e equations. The discretized algebraic momentum equations and k- s equations are solved by the time-diversion multiple access （CTDMA） method. The algebraic Poisson equations are solved by the Tschebyscheff SLOR （successive linear over relaxation） method with alternating computational directions. At the end of the paper, the unsteady flow at high Reynolds numbers through a simplified cascade made up of NACA65-410 blade are simulated with the program written according to the implicit numerical scheme. The reliability and accuracy of the implicit numerical scheme are verified through the satisfactory agreement between the numerical results of the surface pressure coefficient and experimental data. The numerical results indicate that Reynolds number and angle of attack are two primary factors affecting the characteristics of unsteady flow.

Complex Spacetime Frame: Four-Vector Identities and Tensors

This paper provides derivation of some basic identities for complex four-component vectors defined in a complex four-dimensional spacetime frame specified by an imaginary temporal axis. The resulting four-vector identities take exactly the same forms of the standard vector identities established in the familiar three-dimensional space, thereby confirming the consistency of the definition of the complex four-vectors and their mathematical operations in the general complex spacetime frame. Contravariant and covariant forms have been defined, providing appropriate definitions of complex tensors, which point to the possibility of reformulating differential geometry within a spacetime frame.

Coordinates over Complex Terrain in Atmospheric Model

In the terrain following coordinate,Gal-Chen and Somerville[1]and other proposed a vertical coordinate z∝(z−z_(bottom))/(ztop−z_(bottom))and constant spatial intervals ofδx andδy along the other directions.Because the variation ofδx andδy was ignored,their coordinate does not really follow the terrain.It fails to reproduce the divergence and curl over a complex terrain.Aligning the coordinate with real terrain,the divergence and curl we obtained from the curvilinear coordinate are consistent with the Cartesian coordinate.With a modification,the simulated total mass,energy,and momentum from the Navier-Stokes equations are conserved and in agreement with those calculated from Cartesian coordinate.

W^t-approximation representations over quasi k-Gorenstein algebras

The notions of quasi k-Gorenstein algebras and W^t-approximation representations are introduced. The existence and uniqueness (up to projective equivalences) of W^t-approximation representations over quasi k-Gorenstein algebras are established. Some applications of W^t-approximation representations to homologically finite subcategories are given.

Numerical simulation of wave transformation,breaking and runup by a contravariant fully non-linear Boussinesq equations model

In this paper we propose a new model based on a contravariant integral form of the fully non-linear Boussinesq equations （FNBE） in order to simulate wave transformation phenomena, wave breaking, runup and nearshore currents in computational domains representing the complex morphology of real coastal regions. The above-mentioned contravariant integral form, in which Christoffel symbols are absent, is characterized by the fact that the continuity equation does not include any dispersive term. The Boussinesq equation system is numerically solved by a hybrid finite volume-f＇mite difference scheme. A high-order upwind weighted essentially non-oscillatory （WENO） finite volume scheme that involves an exact Riemann solver is implemented. The wave breaking is represented by discontinuities of the weak solution of the integral form of the non-linear shallow water equations （NSWE）. On the basis of the shock-capturing high order WENO scheme a new procedure, for the computation of the structure of the solution of a Riemann problem associated with a wet/dry front, is proposed in order to simulate the run up hydrodynamics in swash zone. The capacity of the proposed model to correctly represent wave propagation, wave breaking, run up and wave induced currents is verified against test cases present in literature. The results obtained are compared with experimental measures, analytical solutions or alternative numerical solutions. The proposed model is applied to a real case regarding the simulation of wave fields and nearshore currents in the coastal region opposite San Mauro Cilento （Italy）.

Spectrum of L-fuzzy Prime Ideals of a Distributive Lattice

A topological space denoted by Fspec(D),called L-fuzzy prime spectrum of a bounded distributive lattice D is introduced.This space Fspec(D) is compact and it contains a subspace homeomorphic with the prime spectrum of D which is dense in it.The correspondence associating D to the topological space Fspec(D) is shown to define a contravarient functor from the category of bounded distributive lattices into the category of compact topological spaces.