In this paper, a switch function-based gas-kinetic scheme (SF-GKS) is presented for the simulation of inviscid and viscous compressible flows. With the finite volume discretization, Euler and Navier-Stokes equations are solved and the SF-GKS is applied to evaluate the inviscid flux at cell interface. The viscous flux is obtained by the conventional smooth function approximation. Unlike the traditional gas-kinetic scheme in the calculation of inviscid flux such as Kinetic Flux Vector Splitting (KFVS), the numerical dissipation is controlled with a switch function in the present scheme. That is, the numerical dissipation is only introduced in the region around strong shock waves. As a consequence, the present SF-GKS can well capture strong shock waves and thin boundary layers simultaneously. The present SF-GKS is firstly validated by its application to the inviscid flow problems, including 1-D Euler shock tube, regular shock reflection and double Mach reflection. Then, SF-GKS is extended to solve viscous transonic and hypersonic flow problems. Good agreement between the present results and those in the literature verifies the accuracy and robustness of SF-GKS.

A counterexample is constructed. It confirms that the error of conforming finite element solution is proportional to the coefficient jump, when solving interface elliptic equations. The Scott-Zhang operator is applied to a nonconforming finite element. It is shown that the nonconforming finite element provides the optimal order approximation in interpolation, in L2-projection, and in solving elliptic differential equation, independent of the coefficient jump in the elliptic differential equation. Numerical tests confirm the theoretical finding.

A mathematical multi-zone ice accretion model used in the numerical simulation of icing on airfoil surface based on three water states, namely, continuous film, rivulets and beads is studied in this paper. An improved multi-zone roughness model is proposed. According to the flow state of liquid water and film flow, rivulets flow governing equations are established to calculate film mass distribution, film velocity, rivulet wetness factor and rivulet mass distribution. Force equilibrium equations of droplet are used to establish the critical conditions of water film broken into rivulets and rivulets broken into beads. The temperature conduction inside the water layer and ice layer is considered. Using the proposed model ice accretion on a NACA0012 airfoil profile with a 4° angle of attack under different icing conditions is simulated. Different ice shapes like glaze ice, mixed ice and rime ice are obtained, and the results agree well with icing wind tunnel experiment data. It can be seen that, water films are formed on the surface, and heights of the films vary with icing time and locations. This results in spatially-temporally varying surface roughness and heat transfer process, ultimately affects the ice prediction. Model simulations indicate that the process of water film formation and evolution cannot be ignored, especially under glaze ice condition.

In this article, Sinc collocation method is considered to obtain the numerical solution of integral algebraic equation of index-1 by reducing it to an explicit system of algebraic equation. It is shown that Sinc collocation solution can produce an error of order . Moreover, Sinc method is applied to several examples to illustrate the accuracy and implementation of the method.

In this paper, we discuss the dependence of the solutions on the parameters (order, initial function, right-hand function) of fractional neutral delay differential equations (FNDDEs). The corresponding theoretical results are given respectively. Furthermore, we present some numerical results that support our theoretical analysis.

This paper presents an analysis of unsteady flow of incompressible fractional Maxwell fluid filled in the annular region between two infinite coaxial circular cylinders. The fluid motion is created by the inner cylinder that applies a longitudinal time-dependent shear stress and the outer cylinder that is moving at a constant velocity. The velocity field and shear stress are determined using the Laplace and finite Hankel transforms. Obtained solutions are presented in terms of the generalized G and R functions. We also obtain the solutions for ordinary Maxwell fluid and Newtonian fluid as special cases of generalized solutions. The influence of different parameters on the velocity field and shear stress are also presented using graphical illustration. Finally, a comparison is drawn between motions of fractional Maxwell fluid, ordinary Maxwell fluid and Newtonian fluid.

This paper presents a lattice Boltzmann (LB) method based study aimed at numerical simulation of aeroacoustic phenomenon in flows around a symmetric obstacle. To simulate the compressible flow accurately, a potential energy double-distribution-function (DDF) lattice Boltzmann method is used over the entire computational domain from the near to far fields. The buffer zone and absorbing boundary condition is employed to eliminate the non-physical reflecting. Through the direct numerical simulation, the flow around a circular cylinder at Re=150, M=0.2 and the flow around a NACA0012 airfoil at Re=10000, M=0.8, α=0° are investigated. The generation and propagation of the sound produced by the vortex shedding are reappeared clearly. The obtained results increase our understanding of the characteristic features of the aeroacoustic sound.

This article addresses the resolution of the inverse problem for the parameter identification in orthotropic materials with a number of measurements merely on the boundaries. The inverse problem is formulated as an optimization problem of a residual functional which evaluates the differences between the experimental and predicted displacements. The singular boundary method, an integration-free, mathematically simple and boundary-only meshless method, is employed to numerically determine the predicted displacements. The residual functional is minimized by the Levenberg-Marquardt method. Three numerical examples are carried out to illustrate the robustness, efficiency, and accuracy of the proposed scheme. In addition, different levels of noise are added into the boundary conditions to verify the stability of the present methodology.

In this paper, we study carbon emission trading whose market is gaining in popularity as a policy instrument for global climate change. The mathematical model is presented for pricing options on CO2 emission allowance futures with jump diffusion processes, and a so-called fitted finite volume method is proposed to solve the pricing model for the spatial discretization, in which the Crank-Nicolson is employed for time stepping. In addition, the stability and the convergence of the fully discrete scheme are given, and some numerical results, which are compared with the closed form solution and the Monte Carlo simulation solution, are provided to demonstrate the rates of convergence and the robustness of the numerical method.

This paper extends an algorithm of P1-conservative interpolation on triangular meshes to tetrahedral meshes and thus constructs an approach of solution reconstruction for three-dimensional problems. The conservation property is achieved by local mesh intersection and the mass of a tetrahedron of the current mesh is calculated by the integral on its intersection with the background mesh. For each current tetrahedron, the overlapped background tetrahedrons are detected efficiently. A mesh intersection algorithm is proposed to construct the intersection of a current tetrahedron with the overlapped background tetrahedron and mesh the intersection region by tetrahedrons. A localization algorithm is employed to search the host units in background mesh for each vertex of the current mesh. In order to enforce the maximum principle and avoid the loss of monotonicity, correction of nodal interpolated solution on tetrahedral meshes is given. The performance of the present solution reconstruction method is verified by numerical experiments on several analytic functions and the solution of the flow around a sphere.

In this paper, we consider a least squares nonconforming finite element of low order for solving the transport equations. We give a detailed overview on the stability and the convergence properties of our considered methods in the stability norm. Moreover, we derive residual type a posteriori error estimates for the least squares nonconforming finite element methods under H–1-norm, which can be used as the error indicators to guide the mesh refinement procedure in the adaptive finite element method. The theoretical results are supported by a series of numerical experiments.