MOC/SN耦合三维中子输运程序KYCORE开发与初步验证

KYCORE程序是中国核动力研究设计院开发的径向MOC（特征线方法）与轴向SN耦合三维中子输运程序。KYCORE将二维MOC与一维SN通过角通量实现高精度耦合,并通过粗网有限差分实现快速收敛,是目前可工程化应用于三维中子输运计算中精度最高的方法之一。介绍了2D/1D计算与加速理论,并通过与蒙特卡罗程序的计算对比,数值验证了KYCORE三维中子计算的准确性与高效性。     

An Exner-based coupled model for two-dimensional transient flow over erodible bed

Transient flow over erodible bed is solved in this work assuming that the dynamics of the bed load problem is described by two mathematical models: the hydrodynamic model, assumed to be well formulated by means of the depth averaged shallow water equations, and the Exner equation. The Exner equation is written assuming that bed load transport is governed by a power law of the flow velocity and by a flow/sediment interaction parameter variable in time and space. The complete system is formed by four coupled partial differential equations and a genuinely Roe-type first order scheme has been used to solve it on triangular unstructured meshes. Exact solutions have been derived for the particular case of initial value Riemann problems with variable bed level and depending on particular forms of the solid discharge formula. The model, supplied with the corresponding solid transport formulae, is tested by comparing with the exact solutions. The model is validated against laboratory experimental data of different unsteady problems over erodible bed.     

Two-dimensional coupled mathematical modeling of fluvial processes with intense sediment transport and rapid bed evolution

Alluvial rivers may experience intense sediment transport and rapid bed evolution under a high flow regime, for which traditional decoupled mathematical river models based on simplified conservation equations are not applicable. A two-dimensional coupled mathematical model is presented, which is generally applicable to the fluvial processes with either intense or weak sediment transport. The governing equations of the model comprise the complete shallow water hydrodynamic equations closed with Manning roughness for boundary resistance and empirical relationships for sediment exchange with the erodible bed. The second-order Total-Variation-Diminishing version of the Weighted-Average-Flux method, along with the HLLC approximate Riemann Solver, is adapted to solve the governing equations, which can properly resolve shock waves and contact discontinuities. The model is applied to the pilot study of the flooding due to a sudden outburst of a real glacial-lake. Supported by the National Basic Research and Development Program of China (973 Program) (Grant No. 2007CB14106), the National Natural Science Foundation of China (Grant No. 50459001), and the Key Project of Chinese Academy of Sciences (Grant No. KZCX3-SW-357-02)     

MOC/SN耦合三维中子输运程序KYCORE开发与初步验证

KYCORE程序是中国核动力研究设计院开发的径向MOC（特征线方法）与轴向SN耦合三维中子输运程序。KYCORE将二维MOC与一维SN通过角通量实现高精度耦合,并通过粗网有限差分实现快速收敛,是目前可工程化应用于三维中子输运计算中精度最高的方法之一。介绍了2D/1D计算与加速理论,并通过与蒙特卡罗程序的计算对比,数值验证了KYCORE三维中子计算的准确性与高效性。     

基于格子Boltzmann方法的超声速非平衡流模拟

基于D1Q4可压缩格子Boltzmann模型,按照流通矢量分裂方法的思路,采用坐标旋转技术构造求解三维带化学反应Navier-Stokes方程对流通量求解器.结合有限体积法求解三维化学非平衡流Navier-Stokes方程,采用时间算子分裂算法解决化学反应刚性问题,数值模拟超声速化学非平衡流的三个经典算例.数值结果表明:在高马赫数下,采用D1Q4可压缩格子Boltzmann模型构造的三维对流通量求解器数值模拟中没有出现非物理解,同时在超声速化学非平衡流场中正确分辨激波、燃烧波等物理现象,精度和分辨率均较高,验证了本文构造的三维对流通量求解器的可靠性,拓宽了D1Q4可压缩格子Boltzmann模型的应用范围,为计算超声速化学非平衡流提供一种新方法.     

Unstructured finite volume discretisation of bed friction and convective flux in solute transport models linked to the shallow water equations

The finite volume discretisation of the shallow water equations has been the subject of many previous studies, most of which deal with a well-balanced conservative discretisation of the convective flux and bathymetry. However, the bed friction discretisation has not been so profusely analysed in previous works, while it may play a leading role in certain applications of shallow water models. In this paper we analyse the numerical discretisation of the bed friction term in the two-dimensional shallow water equations, and we propose a new unstructured upwind finite volume discretisation for this term. The new discretisation proposed improves the accuracy of the model in problems in which the bed friction is a relevant force in the momentum equation, and it guarantees a perfect balance between gravity and bed friction under uniform flow conditions. The relation between the numerical scheme used to solve the hydrodynamic equations and the scheme used to solve a scalar transport model linked to the shallow water equations, is also analysed in the paper. It is shown that the scheme used in the scalar transport model must take into consideration the scheme used to solve the hydrodynamic equations. The most important implication is that a well-balanced and conservative scheme for the scalar transport equation cannot be formulated just from the water depth and velocity fields, but has to consider also the way in which the hydrodynamic equations have been solved.     

A novel simulation theory and model system for multi-field coupling pipe-flow system

Due to the lack of a theoretical basis for multi-field coupling in many system-level models, a novel set of system-level basic equations for flow/heat transfer/combustion coupling is put forward. Then a finite volume model of quasi-1D transient flow field for multi-species compressible variable-cross-section pipe flow is established by discretising the basic equations on spatially staggered grids. Combining with the 2D axisymmetric model for pipe-wall temperature field and specific chemical reaction mechanisms, a finite volume model system is established; a set of specific calculation methods suitable for multi-field coupling system-level research is structured for various parameters in this model; specific modularisation simulation models can be further derived in accordance with specific structures of various typical components in a liquid propulsion system. This novel system can also be used to derive two sub-systems: a flow/heat transfer two-field coupling pipe-flow model system without chemical reaction and species diffusion; and a chemical equilibrium thermodynamic calculation-based multi-field coupling system. The applicability and accuracy of two sub-systems have been verified through a series of dynamic modelling and simulations in earlier studies. The validity of this system is verified in an air–hydrogen combustion sample system. The basic equations and the model system provide a unified universal theory and numerical system for modelling and simulation and even virtual testing of various pipeline systems.     

A primitive-variable Riemann method for solution of the shallow water equations with wetting and drying

A Riemann flux that uses primitive variables rather than conserved variables is developed for the shallow water equations with nonuniform bathymetry. This primitive-variable flux is both conservative and well behaved at zero depth. The unstructured finite-volume discretization used is suitable for highly nonuniform grids that provide resolution of complex geometries and localized flow structures. A source-term discretization is derived for nonuniform bottom that balances the discrete flux integral both for still water and in dry regions. This primitive-variable formulation is uniformly valid in wet and dry regions with embedded wetting and drying fronts. A fully nonlinear implicit scheme and both nonlinear and time-linearized explicit schemes are developed for the time integration. The implicit scheme is solved by a parallel Newton-iterative algorithm with numerically computed flux Jacobians. A concise treatment of characteristic-variable boundary conditions with source terms is also given. Computed results obtained for the one-dimensional dam break on wet and dry beds and for normal-mode oscillations in a circular parabolic basin are in very close agreement with the analytical solutions. Other results for a forced breaking wave with friction interacting with a sloped bottom demonstrate a complex wave motion with wetting, drying and multiple interacting wave fronts. Finally, a highly nonuniform, coastline-conforming unstructured grid is used to demonstrate an unsteady simulation that models an artificial coastal flooding due to a forced wave entering the Gulf of Mexico.     

New Symmetry Reductions,Dromions—Like and Compacton Solutions for a 2D BS（m,n） Equations Hierarchy with Fully Nonlinear Dispersion

We have found two types of important exact solutions,compacton solutions,which are solitary waves with the property that after colliding with their own kind,they re-emerge with the same coherent shape very much as the solitons do during a completely elastic interaction,in the (1 1)D,(1 2)D and even (1 3)D models,and dromion solutions (exponentially decaying solutions in all direction) in many (1 2)D and (1 3)D models.In this paper,symmetry reductions in (1 2)D are considered for the break soliton-type equation with fully nonlinear dispersion (called BS(m,n) equation)ut b(u^m)xxy 4b(u^n δx^-1uy)x=0,which is a generalized model of (1 2)D break soliton equation ut buxxy 4buuy 4buxδx^-1uy=0,by using the extended direct reduction method.As a result,six types of symmetry reductions are obtained.Starting from the reduction equations and some simple transformations,we obtain the solitary wavke solutions of BS(1,n) equations,compacton solutions of BS(m,m-1) equations and the compacton-like solution of the potential form (called PBS(3,2)) ωxt b(ux^m)xxy 4b(ωx^nωy)x=0.In addition,we show that the variable ∫^x uy dx admits dromion solutions rather than the field u itself in BS(1,n) equation.     

A Riemann solver for unsteady computation of 2D shallow flows with variable density

A novel 2D numerical model for vertically homogeneous shallow flows with variable horizontal density is presented. Density varies according to the volumetric concentration of different components or species that can represent suspended material or dissolved solutes. The system of equations is formed by the 2D equations for mass and momentum of the mixture, supplemented by equations for the mass or volume fraction of the mixture constituents. A new formulation of the Roe-type scheme including density variation is defined to solve the system on two-dimensional meshes. By using an augmented Riemann solver, the numerical scheme is defined properly including the presence of source terms involving reaction. The numerical scheme is validated using analytical steady-state solutions of variable-density flows and exact solutions for the particular case of initial value Riemann problems with variable bed level and reaction terms. Also, a 2D case that includes interaction with obstacles illustrates the stability and robustness of the numerical scheme in presence of non-uniform bed topography and wetting/drying fronts. The obtained results point out that the new method is able to predict faithfully the overall behavior of the solution and of any type of waves.     

Formation of ripples over a sand bed submitted to a turbulent shear flow

We investigate the process of ripple formation when a sand bed is submitted to a steady and turbulent liquid flow. The sand transport dynamics is described in terms of a simple relaxation law which accounts for the fact that the transport rate does not adapt instantaneously to its equilibrium value. The equilibrium sand flux is evaluated using a standard law based on the estimation of the flow shear stress calculated at the sand bed surface. The latter is estimated from an analytical resolution of the flow over a deformed sand bed which is based on the Jackson and Hunt calculation [J.C.R. Hunt, Quart. J. R. Met. Soc. 101, 929 (1975)]. Within this model, we investigate the stability of the sand bed and are able to derive analytical scaling laws for the wavelength and phase velocity of the most dangerous mode. In the deep flow limit, the model predicts the occurrence of a single mode of instability corresponding to the formation of ripples. Predictions of our model are compared with previous models and available experimental data.     

Modeling and numerical approximation of a 2.5D set of equations for mesoscale atmospheric processes

The set of 3D inviscid primitive equations for the atmosphere is dimensionally reduced by a Discontinuous Galerkin discretization in one horizontal direction. The resulting model is a 2D system of balance laws with a source term depending on the layering procedure and the choice of coupling fluxes, which is established in terms of upwind considerations. The “2.5D” system is discretized via a WENO–TVD scheme based in a flux-limiter approach. We study four tests cases related to atmospheric phenomena to analyze the physical validity of the model.     

LTE computation of axisymmetric arc-flow interaction incircuit-breakers

A numerical procedure for the computation of arc-flow interaction in gas-blast circuit breakers is presented. In the proposed approach the flow is obtained via the solution of the Euler equations and the coupling with the arc by a source term in the energy equation. This source term is composed of an ohmic energy production part and a radiative transport part. The equations are solved by a time-marching procedure using a finite-volume discretization. The model is applied to the computation of an axisymmetric arc in a model breaker. The results reproduce qualitatively the major features of arc-flow interaction and show the presence of a localized form of choking of the flow around the arc boundary     

Onsager relations and hydrodynamic balance equations in 2D quantum wells

We clarify the role of heat flux in the hydrodynamic balance equations in 2D quantum wells, facilitating the formulation of an Onsager relation within the framework of this theory. We find that the Onsager relation is satisfied within the framework of the 2D hydrodynamic balance equation transport theory at sufficiently high density. The condition of high density is consonant with the requirement of strong electron-electron interactions for the validity of our balance equation formulation.     

Stability of an erodible bed in various shear flows

The 2D laminar quasi-steady asymptotically simplified and linearized flow with a simplified mass transport of sediments is solved over a slowly erodible bed in various laminar basic shear flow (steady, oscillating or decelerating). The simplified mass transport equation includes the two following phenomena: flux of erosion when the skin friction goes over a threshold value, and a non local effect coming either from an inertial effect or from a slope effect. It is shown that the bed is always unstable for small wave numbers. Examples of long time evolution in various shear régimes are presented, wave trains of ripples are created and merge into a unique bump. This coarsening process is such that the maximum wave length obeys a power law with time.     

非正交曲线坐标下三维粘性流动数值分析

本文基于非正交曲线坐标与相应的非正交速度分量下导得的守恒型N—S方程，讨论了求解三维粘性流动的数值方法，计算中显式时间推进算法与Baldwin—Lomax湍流模型被采用，应用本工作发展的程序，作为算例计算了一个沿径向非等截面环形叶栅的三维粘性流场，得到了诸如三维压力分布，总压损失分布以及十分清晰的二次流动图景等丰富的流场信息。     

GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method

Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates (S_n) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.     

Lattice Boltzmann Simulation of Free-Surface Temperature Dispersion in Shallow Water Flows

We develop a lattice Boltzmann method for modeling free-surface temperature dispersion in the shallow water flows. The governing equations are derived from the incompressible Navier-Stokes equations with assumptions of shallow water flows including bed frictions, eddy viscosity, wind shear stresses and Coriolis forces. The thermal effects are incorporated in the momentum equation by using a Boussinesq approximation. The dispersion of free-surface temperature is modelled by an advection-diffusion equation. Two distribution functions are used in the lattice Boltzmann method to recover the flow and temperature variables using the same lattice structure. Neither upwind discretization procedures nor Riemann problem solvers are needed in discretizing the shallow water equations. In addition, the source terms are straightforwardly included in the model without relying on well-balanced techniques to treat flux gradients and source terms. We validate the model for a class of problems with known analytical solutions and we also present numerical results for sea-surface temperature distribution in the Strait of Gibraltar.     

Simulation of the relationship between porosity and tortuosity in porous media with cubic particles

Tortuosity is an important parameter used in areas such as vascular medicine,neurobiology,and the field of soil permeability and diffusion to express the mass transport in porous media.It is a function of the porosity and the shape and distribution of particles.In this paper,the tortuosity of cubic particles is calculated.With the assumption that the porous medium is homogeneous,the problem is converted to the micro-level over a unit cell,and geometry models of flow paths are proposed.In three-dimensional(3D) cells,the flow paths are too complicated to define.Hence,the 3D models are converted to two-dimensional(2D) models to simplify the calculation process.It is noticed that the path in the 2D model is shorter than that in the 3D model.As a result,triangular particles and the interaction are also taken into consideration to account for the longer distance respectively.We have proposed quadrate particle and interaction(QI) and quadrate and triangular particle(QT) models with cubic particles.Both models have shown good agreement with the experimental data.It is also found that they can predict the toruosities of some kinds of porous media,like freshwater sediment and Negev chalk.