NOVEL IMMERSED BOUNDARY-LATTICE BOLTZMANN METHOD BASED ON FEEDBACK LAW

The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.

A local pseudo arc-length method for hyperbolic conservation laws

A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.

Identifying influential spreaders in complex networks based on entropy weight method and gravity law

In complex networks,identifying influential spreader is of great significance for improving the reliability of networks and ensuring the safe and effective operation of networks.Nowadays,it is widely used in power networks,aviation networks,computer networks,and social networks,and so on.Traditional centrality methods mainly include degree centrality,closeness centrality,betweenness centrality,eigenvector centrality,k-shell,etc.However,single centrality method is onesided and inaccurate,and sometimes many nodes have the same centrality value,namely the same ranking result,which makes it difficult to distinguish between nodes.According to several classical methods of identifying influential nodes,in this paper we propose a novel method that is more full-scaled and universally applicable.Taken into account in this method are several aspects of node’s properties,including local topological characteristics,central location of nodes,propagation characteristics,and properties of neighbor nodes.In view of the idea of the multi-attribute decision-making,we regard the basic centrality method as node’s attribute and use the entropy weight method to weigh different attributes,and obtain node’s combined centrality.Then,the combined centrality is applied to the gravity law to comprehensively identify influential nodes in networks.Finally,the classical susceptible-infected-recovered(SIR)model is used to simulate the epidemic spreading in six real-society networks.Our proposed method not only considers the four topological properties of nodes,but also emphasizes the influence of neighbor nodes from the aspect of gravity.It is proved that the new method can effectively overcome the disadvantages of single centrality method and increase the accuracy of identifying influential nodes,which is of great significance for monitoring and controlling the complex networks.

SOLUTION OF THE RAYLEIGH PROBLEM FOR A POWER-LAW NON-NEWTONIAN CONDUCTING FLUID VIA GROUP METHOD

An investigation is made of the magnetic Rayleigh problem where a semi_infinite plate is given an impulsive motion and thereafter moves with constant velocity in a non_Newtonian power law fluid of infinite extent. The solution of this highly non_linear problem is obtained by means of the transformation group theoretic approach. The one_parameter group transformation reduces the number of independent variables by one and the governing partial differential equation with the boundary conditions reduce to an ordinary differential equation with the appropriate boundary conditions. Effect of the some parameters on the velocity u(y,t) has been studied and the results are plotted.

SPECTRAL/HP ELEMENT METHOD WITH HIERARCHICAL RECONSTRUCTION FOR SOLVING NONLINEAR HYPERBOLIC CONSERVATION LAWS

The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.

Godunov＇s method for initial-boundary value problem of scalar conservation laws

This paper is concerned with Godunov＇s scheme for the initial-boundary value problem of scalar conservation laws. A kind of Godunov＇s scheme, which satisfies the boundary entropy condition, was given. By use of the scheme, numerical simulation for the weak entropy solution to the initial-boundary value problem of scalar conservation laws is conducted.

Conservative and Easily Implemented Finite Volume Semi-Lagrangian WENO Methods for 1D and 2D Hyperbolic Conservation Laws

The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux in time into the integration in space. Compared with the traditional semi-Lagrange scheme, the scheme devised here tries to directly evaluate the average fluxes along cell edges. It is this difference that makes the scheme in this paper simple to implement and easily extend to nonlinear cases. The procedure of evaluation of the average fluxes only depends on the high-order spatial interpolation. Hence the scheme can be implemented as long as the spatial interpolation is available, and no additional temporal discretization is needed. In this paper, the high-order spatial discretization is chosen to be the classical 5th-order weighted essentially non-oscillatory spatial interpolation. In the end, 1D and 2D numerical results show that this method is rather robust. In addition, to exhibit the numerical resolution and efficiency of the proposed scheme, the numerical solutions of the classical 5th-order WENO scheme combined with the 3rd-order Runge-Kutta temporal discretization (WENOJS) are chosen as the reference. We find that the scheme proposed in the paper generates comparable solutions with that of WENOJS, but with less CPU time.

Numerical Method for Non-Linear Conservation Laws: Inviscid Burgers Equation

This paper deals with the Burgers equation which is the most common model used in the nonlinear conservation laws. Here the theoretical aspect of conservation law is discussed by using inviscid Burgers equation. At first, we introduce the general non-linear conservation law as a partial differential equation and its solution procedure by the method of characteristic. Next, we present the weak solution of the problem with entropy condition. Taking into account shock wave and rarefaction wave, the Riemann problem has also been discussed. Finally, the finite volume method is considered to approximate the numerical solution of the inviscid Burgers equation with continuous and discontinuous initial data. An illustration of the problem is provided by some examples. Moreover, the Godunov method provides a good approximation for the problem.

A Riemann-Solver Free Spacetime Discontinuous Galerkin Method for General Conservation Laws

This paper summarizes a Riemann-solver-free spacetime discontinuous Galerkin method developed for general conservation laws. The method integrates the best features of the spacetime Conservation Element/Solution Element (CE/SE) method and the discontinuous Galerkin (DG) method. The core idea is to construct a staggered spacetime mesh through alternate cell-centered CEs and vertex-centered CEs within each time step. Inside each SE, the solution is approximated using high-order spacetime DG basis polynomials. The spacetime flux conservation is enforced inside each CE using the DG concept. The unknowns are stored at both vertices and cell centroids of the spatial mesh. However, the solutions at vertices and cell centroids are updated at different time levels within each time step in an alternate fashion. Thanks to the staggered spacetime formulation, there are no left and right states for the solution at the spacetime interface. Instead, the solution available to evaluate the flux is continuous across the interface. Therefore, no (approximate) Riemann solvers are needed to provide a unique numerical flux. The current method can be used to solve arbitrary conservation laws including the compressible Euler equations, shallow water equations and magnetohydrodynamics (MHD) equations without the need of any form of Riemann solvers. A set of benchmark problems of various conservation laws are presented to demonstrate the accuracy of the method.

Research on the Method of Balance Sheet Quality Inspection based on the Ford Law

This paper based on Benford＇s Law （Benford Law）, the Shanghai stock exchange 15 listed bank balance sheets for the past 5 years data for the study sample, results show that the model in the inspection report can quantify the quality in the process of accounting information quality comparison and has the outstanding advantages of simple, clear, this paper focuses on the Benford＇s law of balance sheet quality inspection method based on in our country and how to test the quality of balance sheets in the introduction and promotion.

New High-Order Numerical Methods for Hyperbolic Systems of Nonlinear PDEs with Uncertainties

In this paper,we develop new high-order numerical methods for hyperbolic systems of nonlinear partial differential equations(PDEs)with uncertainties.The new approach is realized in the semi-discrete finite-volume framework and is based on fifth-order weighted essentially non-oscillatory(WENO)interpolations in(multidimensional)random space combined with second-order piecewise linear reconstruction in physical space.Compared with spectral approximations in the random space,the presented methods are essentially non-oscillatory as they do not suffer from the Gibbs phenomenon while still achieving high-order accuracy.The new methods are tested on a number of numerical examples for both the Euler equations of gas dynamics and the Saint-Venant system of shallow-water equations.In the latter case,the methods are also proven to be well-balanced and positivity-preserving.

Spectral Element Viscosity Methods for Nonlinear Conservation Laws on the Semi-Infinite Interval

In this paper we propose a spectral element: vanishing viscosity (SEW) method for the conservation laws on the semi-infinite interval. By using a suitable mapping, the problem is first transformed into a modified conservation law in a bounded interval, then the well-known spectral vanishing viscosity technique is generalized to the multi-domain case in order to approximate this trarsformed equation more efficiently. The construction details and convergence analysis are presented. Under a usual assumption of boundedness of the approximation solutions, it is proven that the solution of the SEW approximation converges to the uniciue entropy solution of the conservation laws. A number of numerical tests is carried out to confirm the theoretical results.

Design and Analysis of Terminal Guidance Law for Kinetic Interceptors based on Pulse Engine

A new terminal guidance law is proposed based on a solid propellant pulse engine and an improved proportional navigation method to address the terminal guidance issue for kinetic interceptors.On this basis,the start-stop curve of the pulse motor during the terminal guidance process is designed,along with its start-up logic.The effectiveness of the proposed guidance strategy is verified through simulation.

A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.

Symmetry Reductions, Exact Solutions and Conservation Laws of Asymmetric Nizhnik-Novikov-Veselov Equation

By applying a direct symmetry method, we get the symmetry of the asymmetric Nizhnik-Novikov-Veselov equation （ANNV）. Taking the special case, we have a finite-dimensional symmetry. By using the equivalent vector of the symmetry, we construct an eight-dimensional symmetry algebra and get the optimal system of group-invariant solutions. To every case of the optimal system, we reduce the ANNV equation and obtain some solutions to the reduced equations. Furthermore, we find some new explicit solutions of the ANNV equation. At last, we give the conservation laws of the ANNV equation.

Criteria for Beverloo's scaling law

Beverloo's scaling law can describe the flow rate of grains discharging from hoppers. In this paper, we show that the Beverloo's scaling law is valid for varying material parameters. The flow rates from a hopper with different hopper and orifice sizes(D, D_0) are studied by running large-scale simulations. When the hopper size is fixed, the numerical results show that Beverloo's law is valid even if the orifice diameter is very large and then the criteria for this law are discussed.To eliminate the effect of walls, it is found that the criteria can be suggested as D-D_0≥ 40 d or D/D_0≥ 2. Interestingly,it is found that there is still a scaling relation between the flow rate and orifice diameter if D/D_0 is fixed and less than 2.When the orifice diameter is close to the hopper size, the velocity field changes and the vertical velocities of grains above the free fall region are much larger. Then, the free fall arch assumption is invalid and Beverloo's law is inapplicable.

Multi-symplectic Runge-Kutta methods for Landau-Ginzburg-Higgs equation

Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations （PDEs） derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.

LOCAL AND PARALLEL FINITE ELEMENT METHOD FOR THE MIXED NAVIER-STOKES/DARCY MODEL WITH BEAVERS-JOSEPH INTERFACE CONDITIONS

In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed and analyzed. Optimal errors are obtained and numerical experiments are presented to show the efficiency and effectiveness of the local and parallel finite element algorithm.

ITERATIVE ALGORITHM FOR AXIALLY ACCELERATING STRINGS WITH INTEGRAL CONSTITUTIVE LAW

A numerical method is proposed to simulate the transverse vibrations of a viscoelastic moving string constituted by an integral law. In the numerical computation, the Galerkin method based on the Hermite functions is applied to discretize the state variables, and the Runge- Kutta method is applied to solve the resulting differential-integral equation system. A linear iterative process is designed to compute the integral terms at each time step, which makes the numerical method more efficient and accurate. As examples, nonlinear parametric vibrations of an axially moving viscoelastic string are analyzed.

Vibrations of FGM thin cylindrical shells with exponential volume fraction law

In this paper, the influence of an exponential volume fraction law on the vibration frequencies of thin functionally graded cylindrical shells is studied. Material properties in the shell thickness direction are graded in accordance with the exponential law. Expressions for the strain-displacement and curvature-displacement relationships are taken from Love＇s thin shell theory. The Rayleigh-Ritz approach is used to derive the shell eigenfrequency equation. Axial modal dependence is assumed in the characteristic beam functions. Natural frequencies of the shells are observed to be dependent on the constituent volume fractions. The results are compared with those available in the literature for the validity of the present methodology.