NUMERICAL SOLUTION OF QPNS EQUATIONS USING SUPG FEM

A numerical formulation is proposed to integrate the set of parabolized Navier-Stokes equations using SUPG FEM. This method introduces a time-dependent relaxation effect that suppresses all divergent or departure solutions. The new scheme is termed QPNS / SUPG, which leads to a well-posed initial value problem and therefore is stable for forward integration. In order to test the capabilities concerning the shock capture and calculating the boundary layer flow, the supersonic laminar flows over a wedge and a flat plate have been evaluated using the present method.

NUMERICAL SIMULATION OF THE FLOW PAST WING-BODY USING EULER EQUATIONS

AbstractThis paper gives H-O type grid generation and applies finite-volume method forwing-body flow.A 3-D Euler code for wing-body has been developed. We improve in thiscode the distribution of the number of grid points over wing and body surfaces,and keepthe training edge of the wing being one of the grid lines even when the edge has sweepbackor sweep forward.The code developed can provide not only the spanwise pressure distribu-tion as usual ou cross-flow planes,but also the chordwise Pressure distributions.The com-putation for NASA TN D-712 wing-body model with the present code shows that thecomputed pressure distributions are in very good agreement with the experiment.

A numerical study for some efficient time-integration schemes using barotropic equations

Three kinds of methods, i. e., explicit, semi-implicit, and semi-implicit and semi-Lagrangian method, are tested in the time-integration of shallow-water equations on rotating sphere. Helpful results are available from experiments, especially about the accuracy and efficiency of different semi-implicit and semi-Lagrangian schemes.

A Radial Basis Function Method with Improved Accuracy for Fourth Order Boundary Value Problems

Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.

NUMERICAL SOLUTION OF NAVIER-STOKES EQUATION OF UNSTEADY SEPARATED FLOWS DUE TO A SPOILER'S OSCILLATION

The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for ψ respectively. The upwind scheme is used for the convective terms. The moving boundary conditions are specially treated, and the effects of outlet conditions on the flow field are abo examined. Numerical results obtained show that the spoiler's oscillation induces forming, growing and shedding of the vortices. The shedding frequency of vortices is equal to that of the spoiler's oscillation. The forced unsteady separated flows under the present investigation depend mainly on the reduced frequency. At low reduced frequency, the vortices shed from the spoiler interact weakly with each other, and move downstream at an almost uniform speed of 038 V∞. At high reduced frequency, the interaction between the adjacent vortices strengthens. They close up to and rotate around each other, and eventually, merge into one vortex.

NUMERICAL METHOD FOR SOLVING THE EULER EQUATION FOR UNSTEADY TRANSONIC FLOWS OVER OSCILLATING AIRFOILS

Through transformations, the time-dependent boundary condition on the airfoil contour and the boundary condition at infinity are brought fixed to the boundaries of a finite domain. The boundary conditions can thus be satisfied exactly without increasing the computational time. The novel scheme is useful for computing transonic, strong disturbance, unsteady flows with high reduced frequencies. The scheme makes use of curvefitted orthogonal meshes and the lattice control technique to obtain the optimal grid distribution. The numerical results are satisfactory.

A High-Order Conservative Numerical Method for Gross-Pitaevskii Equation with Time-Varying Coefficients in Modeling BEC

We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation （BEC）. This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.

SEVERAL NEW TYPES OF FINITE-DIFFERENCE SCHEMES FOR SHALLOW-WATER EQUATION WITH INITIAL-BOUNDARY VALUE AND THEIR NUMERICAL EXPERIMENT

This paper proposed several new types of finite-difference methods for the shallow water equation in absolute coordinate system and put forward an effective two-step predictor-corrector method, a compact and iterative algorithm for five diagonal matrix. Then the iterative method was used for a multi-grid procedure for shallow water equation. A t last, an initial-boundary value problem was considered, and the numerical results show that the linear sinusoidal wave would successively evolve into conoidal wave.

THE STABILITY AND CONVERGENCE OF THE FINITE ANALYTIC METHOD FOR THE NUMERICAL SOLUTION OF CONVECTIVE DIFFUSION EQUATION

In this paper we make a close study of the finite analytic method by means of the maximum principles in differential equations and give the proof of the stability and convergence of the finite analytic method.

NUMERICAL SOLUTION OF QUASILINEAR SINGULARLY PERTURBED ORDINARY DIFFERENTIAL EQUATION WITHOUT TURNING POINTS

In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an algorithm whose accuracy is good for arbitrary e>0 .

NUMERICAL SOLUTION OF THE SINGULARLY PERTURBED PROBLEM FOR THE HYPERBOLIC EQUATION WITH INITIAL JUMP

In this paper we consider the initial-boundary value problem for a second order hyperbolic equation with initial jump. The bounds on the derivatives of the exact solution are given. Then a difference scheme is constructed on a non-uniform grid. Finally, uniform convergence of the difference solution is proved in the sense of the discrete energy norm.

A New Modified EWMA Control Chart for Monitoring Processes Involving Autocorrelated Data

Control charts are one of the tools in statistical process control widely used for monitoring,measuring,controlling,improving the quality,and detecting problems in processes in variousfields.The average run length(ARL)can be used to determine the efficacy of a control chart.In this study,we develop a new modified exponentially weighted moving average(EWMA)control chart and derive explicit formulas for both one and the two-sided ARLs for a p-order autoregressive(AR(p))process with exponential white noise on the new modified EWMA control chart.The accuracy of the explicit formulas was compared to that of the well-known numerical integral equation(NIE)method.Although both methods were highly consistent with an absolute percentage difference of less than 0.00001%,the ARL using the explicit formulas method could be computed much more quickly.Moreover,the performance of the explicit formulas for the ARL on the new modified EWMA control chart was better than on the modified and standard EWMA control charts based on the relative mean index(RMI).In addition,to illustrate the applicability of using the proposed explicit formulas for the ARL on the new modified EWMA control chart in practice,the explicit formulas for the ARL were also applied to a process with real data from the energy and agriculturalfields.

An inverse problem to estimate an unknown order of a Riemann-Liouville fractional derivative for a fractional Stokes' first problem for a heated generalized second grade fluid

In this paper,we propose a numerical method to estimate the unknown order of a Riemann-Liouville fractional derivative for a fractional Stokes＇ first problem for a heated generalized second grade fluid.The implicit numerical method is employed to solve the direct problem.For the inverse problem,we first obtain the fractional sensitivity equation by means of the digamma function,and then we propose an efficient numerical method,that is,the Levenberg-Marquardt algorithm based on a fractional derivative,to estimate the unknown order of a Riemann-Liouville fractional derivative.In order to demonstrate the effectiveness of the proposed numerical method,two cases in which the measurement values contain random measurement error or not are considered.The computational results demonstrate that the proposed numerical method could efficiently obtain the optimal estimation of the unknown order of a RiemannLiouville fractional derivative for a fractional Stokes＇ first problem for a heated generalized second grade fluid.

Benchmarks of 3D Laplace Equation Solvers in a Cubic Configuration for Streamer Simulation

The aim of this paper is to test a developed SOR R＆B method using the Chebyshev accelerator algorithm to solve the Laplace equation in a cubic 3D configuration. Comparisons are made in terms of precision and computing time with other elliptic equation solvers proposed in the open source LIS library. The first results, obtained by using a single core on a HPC, show that the developed SOR R＆B method is efficient when the spectral radius needed for the Chebyshev acceleration is carefully pre-estimated. Preliminary results obtained with a parallelized code using the MPI library are also discussed when the calculation is distributed over one hundred cores.

THE APPLICATION OF INTEGRAL EQUATIONS TO THE NUMERICAL SOLUTION OF NONLINEAR SINGULAR PERTURBATION PROBLEMS

The nonlinear singular perturbation problem is solved numerically on nonequidistant meshes which are dense in the boundary layers. The method presented is based on the numerical solution of integral equations [1]. The fourth order uniform accuracy of the scheme is proved. A numerical experiment demonstrates the effectiveness of the method.

Monte-Carlo simulation of a stochastic differential equation

For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.

A key factor to the spin parameter of uniformly rotating compact stars: crust structure

We study the dimensionless spin parameter j ≡ cJ/（GM2） of different kinds of uniformly rotating compact stars, including traditional neutron stars, hyperonic neutron stars and hybrid stars, based on relativistic mean field theory and the MIT bag model. It is found that j ~ 0.7, which had been suggested in traditional neutron stars, is sustained for hyperonic neutron stars and hybrid stars with M 〉 0.5 MG. Not the interior but rather the crust structure of the stars is a key factor to determine jmax for three kinds of selected compact stars. Furthermore, a universal formula j = 0.63（f/fK） -- 0.42（f/fK）2 ＋ 0.48（f/fK）z is suggested to determine the spin parameter at any rotational frequency f smaller than the Keplerian frequency fK.

NUMERICAL ANALYSIS FOR A MEAN-FIELD EQUATION FOR THE ISING MODEL WITH GLAUBER DYNAMICS

In this paper, a mean-field equation of motion which is derived by Penrose (1991) for the dynamic Ising model with Glauber dynamics is considered. Various finite difference schemes such as explicit Euler scheme, predictor-corrector scheme and some implicit schemes are given and their convergence, stabilities and dynamical properties are discussed. Moreover, a Lyapunov functional for the discrete semigroup {S}(n>0) is constructed. Finally, numerical examples are computed and analyzed. it shows that the model is a better approximation to Cahn-Allen equation which is mentioned in Penrose (1991).

EXPONENTIAL INTEGRATORS FOR STOCHASTIC SCHRODINGER EQUATIONS DRIVEN BY ITO NOISE

We study an explicit exponential scheme for the time discretisation of stochastic SchrS- dinger Equations Driven by additive or Multiplicative It6 Noise. The numerical scheme is shown to converge with strong order 1 if the noise is additive and with strong order 1/2 for multiplicative noise. In addition, if the noise is additive, we show that the exact solutions of the linear stochastic Sehr6dinger equations satisfy trace formulas for the expected mass, energy, and momentum （i. e., linear drifts in these quantities）. Furthermore, we inspect the behaviour of the numerical solutions with respect to these trace formulas. Several numerical simulations are presented and confirm our theoretical results.