高超声速流动计算中LU-SGS隐式算法的应用

在高超声速条件下,对原始LU-SGS格式及其改进方法的收敛速度做了深入的比较分析,目的是进一步更好地将LU-SGS算法用于工程上复杂外形的计算模拟中。二维圆柱,三维钝锥及空天飞机算例的结果表明:对于高超声速粘性流动的计算,粘性项应进行隐式处理;BLU-SGS方法给出的内迭代方式的收敛性优于DP-LUR方法所给出的内迭代方式;LU-SGS算法中雅克比系数矩阵的计算方式对计算量及收敛性影响较大,若采用精确的矩阵形式,则在流动无分离情况下能取得快速收敛的效果,而在含有流动分离的情况,因受稳定性的影响,精确的矩阵形式的收敛表现不及对角近似形式。

全隐LU-SGS算法在高超声速热化学非平衡流刚性问题中的应用

在高超声速热化学非平衡流动计算中,当地气体能量松弛时间、化学反应特征时间与流动时间推进步长量级差异过大会带来严重数值刚性问题,且在高雷诺数条件下,壁面、拐角等强干扰区网格加密使得该问题加剧,导致初始最大CFL数极小,收敛速度缓慢。原始LU-SGS算法仅考虑化学反应源项和对流项的隐式处理,通过推导黏性项Jacobian矩阵谱半径并采用对角近似处理,发展了热化学非平衡FLU-SGS和BLU-SGS两种全隐LU-SGS算法;针对高焓二维圆柱和轴对称返回舱算例,对比改进前后三种算法的收敛特性。结果表明,FLU-SGS及BLU-SGS算法能够快速建立强黏性干扰和大分离流场、解决热化学非平衡复杂流计算中的刚性问题,实现初始最大CFL数3至5个量级的提升,加速收敛效果明显。

CHARACTERISTIC GALERKIN METHOD FOR CONVECTION-DIFFUSION EQUATIONS AND IMPLICIT ALGORITHM USING PRECISE INTEGRATION

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.

PREDICTOR-CORRECTOR ALGORITHMS FOR SOLVING GENERALIZED MIXED IMPLICIT QUASI-EQUILIBRIUM PROBLEMS

A new class of generalized mixed implicit quasi-equilibrium problems （GMIQEP） with four-functions is introduced and studied. The new class of equilibrium problems includes many known generalized equilibrium problems and generalized mixed implicit quasi-variational inequality problems as many special cases. By employing the auxiliary principle technique, some predictor-corrector iterative algorithms for solving the GMIQEP are suggested and analyzed. The convergence of the suggested algorithm only requires the continuity and the partially relaxed implicit strong monotonicity of the mappings

EXISTENCE AND ALGORITHM OF SOLUTIONS FOR GENERAL MULTIVALUED MIXED IMPLICIT QUASI-VARIATIONAL INEQUALITIES

A new class of general multivalued mixed implicit quasi-variational inequalities in a real Hilbert space was introduced, which includes the known class of generalized mixed implicit quasi-variational inequalities as a special case , introduced and studied by Ding Xie-ping . The auxiliary variational principle technique was applied to solve this class of general multivalued mixed implicit quasi-variational inequalities. Firstly, a new auxiliary variational inequality with a proper convex , lower semicontinuous , binary functional was defined and a suitable functional was chosen so that its unique minimum point is equivalent to the solution of such an auxiliary variational inequality . Secondly , this auxiliary variational inequality was utilized to construct a new iterative algorithm for computing approximate solutions to general multivalued mixed implicit quasi-variational inequalities . Here , the equivalence guarantees that the algorithm can generate a sequence of approximate solutions. Finally, the existence of solutions and convergence of approximate solutions for general multivalued mixed implicit quasi-variational inequalities are proved. Moreover, the new convergerce criteria for the algorithm were provided. Therefore, the results give an affirmative answer to the open question raised by M. A . Noor, and extend and improve the earlier and recent results for various variational inequalities and complementarity problems including the corresponding results for mixed variational inequalities, mixed quasi-variational inequalities and quasi-complementarity problems involving the single-valued and set- valued mappings in the recent literature .

Existence and algorithm of solutions for a system of generalized mixed implicit equilibrium problems in Banach spaces

A new system of generalized mixed implicit equilibrium problems is introduced and studied in Banach spaces. First, the notion of the Yosida proximal mapping for generalized mixed implicit equilibrium problems is introduced. By using the notion, a system of generalized equation problems is considered, and its equivalence with the system of generalized mixed implicit equilibrium problems is also proved. Next, by applying the system of generalized equation problems, we suggest and analyze an iterative algorithm to compute the approximate solutions of the system of generalized mixed implicit equilibrium problems. The strong convergence of the iterative sequences generated by the algorithm is proved under quite mild conditions. The results are new and unify and generalize some recent results in this field.

AN IMPLICIT ALGORITHM OF THIN LAYER EQUATIONS IN VISCOUS，TRANSONIC，TWO－PHASE NOZZLE FLOW

Omitting viscosity along flow direction, we have simplified the dimensionless N-Sequations in arbitrary curved coordinate system as the thin layer equations. Using theimplicit approximate-factorization algorithm to solve the gas-phase governing equ-ations and the characteristic method to follow the tracks of particles, we then obtainedthe full coupled numerical method of two-phase.transonic, turbulent flow. Here, par- ticle size may be grouped, the subsonic boundary condition at entry of nozzle is ireatedby quasi-characteristic method in reference plane and the algebraic model is used forturbulent flow. These methods are applied in viscous two-phase flow. calculation of ro-cket nozzle and in the prediciton of thrust and specific impulse for solid propellant ro-cket motor. The calculation results are in good agreement with the measurerment va-lues. Moreover, the influences of different particle radius, different particle mass frac-tion and particle size grouped on flow field have been discussed, and the influences of particle two-dimensional radial velosity component and viscosity on specific impulse ofrocket motor have been analysed.The method of this paper possesses the advantage of saving computer time. More important, the effect is more obvious for the calculation of particle size being grouped.

Existence and Algorithm of Solutions for Generalized Mixed Implicit Quasi-Variational-Like Inequalities

By applying the auxiliary variational principle technique, the existence of solutions for a new class of generalized mixed implicit quasi-variational-like inequalities and the convergence criteria of a new iterative algorithm to compute approximate solutions are proved in Hilbert spaces. The obtained result is a improvement over and generalization of the main theorem proposed by Ding.

Iterative algorithm for solutions to new system of generalized mixed implicit equilibrium

A new system of generalized mixed implicit equilibrium problems （SGMIEP） involving nonmonotone set-valued mappings is introduced and studied in real reflexive Banach spaces. First, an auxiliary mixed equilibrium problem （AMEP） is introduced. The existence and the uniqueness of the solutions to the AMEP are proved under quite mild assumptions without any coercive conditions. Next, by using the solution mapping of the AMEP, a system of generalized equation problems （SGEP） is considered, and its equivalence with the SGMIEP is shown. By using the SGEP, a new iterative algorithm for solving the SGMIEP is proposed and analyzed. The strong convergence of the iterative sequences generated by the algorithm is proved under suitable conditions. These results are new, which unify and generalize some recent results in this field.

A New Eighth Order Implicit Block Algorithms for the Direct Solution of Second Order Ordinary Differential Equations

This paper focuses on derivation of a uniform order 8 implicit block method for the direct solution of general second order differential equations through continuous coefficients of Linear Multi-step Method (LMM). The continuous formulation and its first derivatives were evaluated at some selected grid and off grid points to obtain our proposed method. The superiority of the method over the existing methods is established numerically.

Full-vectorial finite-difference beam propagation method based on the modified alternating direction implicit method

A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.

ON A CLASS OF GENERALIZED NONLINEAR IMPLICIT QUASIVARIATIONAL INCLUSIONS

In this paper, a new class of generalized nonlinear implicit quasivariational inclusions involving a set-valued maximal monotone wrapping are studied. A existence theorem of solutions for this class of generalized nonlinear implicit quasivariational inclusions is Proved without compactness assumptions. A new iterative algorithm for finding approximate solutions of the generalized nonlinear implicit quasivariational inclusions is suggested and analysed and the convergence of iterative sequence generated by the new algorithm is also given, As special cases, some known results in this field are also discussed.

An implicit upwind parabolized Navier-Stokes code for chemically nonequilibrium flows

The previously developed single-sweep parabolized Navier-Stokes （SSPNS） space marching code for ideal gas flows has been extended to compute chemically nonequilibrium flows. In the code, the strongly coupled set of gas dynamics, species conservation, and turbulence equations is integrated with the implicit lower-upper symmetric GaussSeidel （LU-SGS） method in the streamwise direction in a space marching manner. The AUSMPW＋ scheme is used to calculate the inviscid fluxes in the crossflow direction, while the conventional central scheme for the viscous fluxes. The k-g two-equation turbulence model is used. The revised SSPNS code is validated by computing the Burrows-Kurkov non-premixed H2/air supersonic combustion flows, premixed H2/air hypersonic combustion flows in a three-dimensional duct with a 15° compression ramp, as well as the hypersonic laminar chemically nonequilibrium air flows around two 10° half-angle cones. The results of these calculations are in good agreement with those of experiments, NASA UPS or Prabhu＇s PNS codes. It can be concluded that the SSPNS code is highly efficient for steady supersonic/ hypersonic chemically reaction flows when there is no large streamwise separation.

High-order fully implicit SIMPLE-based model for fully implicit simulation of upward two-phase flow

The drift-flux model has a practical importance in two-phase flow analysis.In this study,a finite volume solution is developed for a transient four-equation drift-flux model through the staggered mesh,leading to the development of a fully implicit discretization method.The main advantage of the fully implicit method is its unconditional stability.Newton's scheme is a popular method of choice for the solution of a nonlinear system of equations arising from fully implicit discretization of field equations.However,the lack of convergence robustness and the construction of Jacobian matrix have created several difficulties for the researchers.In this paper,a fully implicit model is developed based on the SIMPLE algorithm for two-phase flow simulations.The drawbacks of Newton's method are avoided in the developed model.Different limiter functions are considered,and the stabilized method is developed under steady and transient conditions.The results obtained by the numerical modeling are in good agreement with the experimental data.As expected,the results prove that the developed model is not restricted by any stability limit.

g-η-Monotone mapping and resolvent operator technique for solving generalized implicit variational-like inclusions

A new class of g-η-monotone mappings and a class of generalized implicit variational-like inclusions involving g-η-monotone mappings are introduced. The resolvent operator of g-η-monotone mappings is defined and its Lipschitz continuity is presented, An iterative algorithm for approximating the solutions of generalized implicit wriational- like inclusions is suggested and analyzed. The convergence of iterative sequences generated by the algorithm is also proved,

A Class of Parallel Implicit Runge-Kutta Formulas

A class of parallel implicit Runge-Kutta formulas is constructed for multiprocessor system. A family of parallel implicit two-stage fourth order Runge-Kutta formulas is given. For these formulas, the convergence is proved and the stability analysis is given. The numerical examples demonstrate that these formulas can solve an extensive class of initial value problems for the ordinary differential equations.

Set-Valued Non-Linear Random Implicit Quasivariational Inclusions

In this paper, we propose iterative algorithms for set valued nonlinear random implicit quasivariational inclusions. We define the related random implicit proximal operator equations and establish an equivalence between them. Finally, we prove the existence and convergence of random iterative sequences generated by random iterative algorithms.

An Implicit Block LU-SGS Algorithm-Based Lattice Boltzmann Flux Solver for Simulation of Hypersonic Flows

This paper proposes a stable and efficient implicit block Lower-Upper Symmetric-Gauss-Seidel(LU-SGS)algorithm-based lattice Boltzmann flux solver(LBFS)for simulation of hypersonic flows.In this method,the finite volume method(FVM)is applied to discretize the Navier-Stokes equations,and the LBFS is utilized to evaluate the numerical flux at the cell interface.In LBFS,the local solution of discrete velocity Boltzmann equation(DVBE)with the non-free parameter D1Q4 lattice Boltzmann model is adopted to reconstruct the inviscid flux across the cell interface,and the viscous flux is approximated by conventional smooth function approach.In order to improve the robustness and convergence rate of the simulation for hypersonic flows,especially for problems with complex geometry,the implicit block LU-SGS algorithm is introduced to solve resultant discrete governing equations.A double cone model at Mach number of Ma=9.86 is firstly simulated to validate the proposed scheme,and a hypersonic flight vehicle with wings and rudders at Mach number of Ma=5.56 is then calculated to extend the application in practical engineering problems.Numerical results show that the proposed scheme could offer a more accurate and effective prediction for hypersonic flows.

GENETIC ALGORITHM IN REDUCTION OF NUMERICAL DISPERSION OF 3-D ADI-FDTD METHOD

A new method to reduce the numerical dispersion of the three-dimensional Alternating Di-rection Implicit Finite-Difference Time-Domain (3-D ADI-FDTD) method is proposed. Firstly,the numerical formulations of the 3-D ADI-FDTD method are modified with the artificial anisotropy,and the new numerical dispersion relation is derived. Secondly,the relative permittivity tensor of the artificial anisotropy can be obtained by the Adaptive Genetic Algorithm (AGA). In order to demon-strate the accuracy and efficiency of this new method,a monopole antenna is simulated as an exam-ple. And the numerical results and the computational requirements of the proposed method are com-pared with those of the conventional ADI-FDTD method and the measured data. In addition the re-duction of the numerical dispersion is investigated as the objective function of the AGA. It is found that this new method is accurate and efficient by choosing proper objective function.

An Implicitly Restarted Block Arnoldi Method in a Vector-Wise Fashion

In this paper, we develop an implicitly restarted block Arnoldi algorithm in a vector-wise fashion. The vector-wise construction greatly simplifies both the detection of necessary deflation and the actual deflation itself, so it is preferable to the block-wise construction. The numerical experiment shows that our algorithm is effective.