In order to study the local refinement issue of the horizontal resolution for a global model with Spherical Centroidal Voronoi Tessellations (SCVTs), the SCVTs are set to 10242 cells and 40962 cells respectively using the density function. The ratio between the grid resolutions in the high and low resolution regions (hereafter RHL) is set to 1:2, 1:3 and 1:4 for 10242 cells and 40962 cells, and the width of the grid transition zone (for simplicity, WTZ) is set to 18° and 9° to investigate their impacts on the model simulation. The ideal test cases, i.e. the cosine bell and global steady-state nonlinear zonal geostrophic flow, are carried out with the above settings. Simulation results showthat the larger the RHL is, the larger the resulting error is. It is obvious that the 1:4 ratio gives rise to much larger errors than the 1:2 or 1:3 ratio; the errors resulting from the WTZ is much smaller than that from the RHL. No significant wave distortion or reflected waves are found when the fluctuation passes through the refinement region, and the error is significantly small in the refinement region. Therefore,when designing a local refinement scheme in the global model with SCVT, the RHL should be less than 1:4, i.e., the error is acceptable when the RHL is 1:2 or 1:3.

In this paper, an alternative lattice Boltzmann (LB)model for incompressible flows is proposed. By modifying directly the moments of the equilibrium distribution function (EDF), the continuous expression of the EDF in tensor Hermite polynomials is derived using the moment expansion and then discretizedwith the discrete velocity vectors of the D2Q9 lattice. The present model as well as its counterpart, the incompressible LB model proposed by Guo, reproduces the incompressible Navier-Stokes (N-S) equations for both steady and unsteady flows. Besides, an alternative pressure formula, which represents the pressure as the diagonal part of the stress tensor, is adopted in the present model. Furthermore, in order to enhance the stability of the present LB model, an additional relaxation time pertaining to the non-hydrodynamic mode is added to the BGK collision operator. The present LB model is validated by two benchmark tests: the cavity flow with different Reynolds number (Re) and the flow past an impulsively started cylinder at Re=40 and 550.

The focus of this article is to present the projected finite element method for solving systems of reaction-diffusion equations on evolving closed spheroidal surfaces with applications to pattern formation. The advantages of the projected finite element method are that it is easy to implement and that it provides a conforming finite element discretization which is “logically” rectangular. Furthermore, the surface is not approximated but described exactly through the projection. The surface evolution law is incorporated into the projection operator resulting in a time-dependent operator. The time-dependent projection operator is composed of the radial projection with a Lipschitz continuous mapping. The projection operator is used to generate the surface mesh whose connectivity remains constant during the evolution of the surface. To illustrate the methodology several numerical experiments are exhibited for different surface evolution laws such as uniform isotropic (linear, logistic and exponential), anisotropic, and concentration-driven. This numerical methodology allows us to study new reaction-kinetics that only give rise to patterning in the presence of surface evolution such as the activator-activator and short-range inhibition; long-range activation.

In this paper, we propose a compact scheme to numerically study the coupled stochastic nonlinear Schrödinger equations. We prove that the compact scheme preserves the discrete stochastic multi-symplectic conservation law, discrete charge conservation law and discrete energy evolution law almost surely. Numerical experiments confirm well the theoretical analysis results. Furthermore, we present a detailed numerical investigation of the optical phenomena based on the compact scheme. By numerical experiments for various amplitudes of noise, we find that the noise accelerates the oscillation of the soliton and leads to the decay of the solution amplitudes with respect to time. In particular, if the noise is relatively strong, the soliton will be totally destroyed. Meanwhile, we observe that the phase shift is sensibly modified by the noise. Moreover, the numerical results present inelastic interaction which is different from the deterministic case.

A cheapest stable nonconforming finite element method is presented for solving the incompressible flow in a square cavity without smoothing the corner singularities. The stable cheapest nonconforming finite element pair based on P1×P0 on rectangularmeshes [29] is employed with a minimal modification of the discontinuous Dirichlet data on the top boundary, where is the finite element space of piecewise constant pressures with the globally one-dimensional checker-board pattern subspace eliminated. The proposed Stokes elements have the least number of degrees of freedom compared to those of known stable Stokes elements. Three accuracy indications for our elements are analyzed and numerically verified. Also, various numerous computational results obtained by using our proposed element show excellent accuracy.

In this paper, numerical simulations of solid particle erosion phenomena on H-type finned circular/elliptic tube surface, which is of great significance to the antiwear design of heat exchanger, are presented. The Eulerian-Lagrangian approach is applied to simulate the dilute gas-solid flow through H-type finned circular/elliptic tubes. A semi-empirical model is adopted to predict the erosion rate. The dynamics behavior of the entrained solid particles in the flow is presented. The geometry of eroded tube surface is changed with the predicted erosion which is taken into account by a UDF and the flow field is re-solved for the eroded tube surface at every time step. The influences of ten parameters (the tube bundle arrangement, particle size, particle concentration, fluid Reynolds number, fin thickness, fin pitch, fin length, fin width, slit width and the transverse tube pitch) on the maximum erosion depth of the H-type circular/elliptic finned tube surface are investigated. Using H-type finned elliptic tube surface can effectively reduce the erosion rate of tube surface comparedwith that using H-type finned circular tube surface. The erosion in in-line arrangement is less severe than that in staggered arrangement. With the increase of particle size, particle concentration and the fluid Reynolds number, the erosion rate of the tube surface rises. The numerically predicted effect of Reynolds number is in good agreement with previous test data. Among the six geometry parameters, the most influential parameter is the transverse tube pitch.

The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state is the superconducting phase for which the Cooper pairs have a non-zero total momentum, depending on the splitting of the Fermi surface sheets for electrons with opposite spin. In infinite systems the momentum is a continuous function of the temperature. In this paper, we have shown how the finite size of the system, through the discretized geometry of the Fermi surface, affects the physical properties of the FFLO state by introducing discontinuities in the Cooper pair momentum. Our calculation in an isotropic system show that the superconducting state with two opposite Cooper pair momenta is more stable than state with one momentum also in nano-size systems, where finite size effects play a crucial role.

This paper systematically compares the numerical implementation and computational cost between the Fourier spectral iterative perturbation method (FSIPM) and the finite element method (FEM) in solving partial differential equilibrium equations with inhomogeneous material coefficients and eigen-fields (e.g., stress-free strain and spontaneous electric polarization) involved in phase-field models. Four benchmark numerical examples, including inhomogeneous elastic, electrostatic, and steady-state heat conduction problems demonstrate that (1) the FSIPM rigorously requires uniform hexahedral (3D) and quadrilateral (2D) mesh and periodic boundary conditions for numerical implementation while the FEM permits arbitrary mesh and boundary conditions; (2) the FSIPM solutions are comparable to their FEM counterparts, and both of them agree with the analytic solutions, (3) the FSIPM is much faster in solving equilibrium equations than the FEM to achieve the accurate solutions, thus exhibiting a greater potential for large-scale 3D computations.

In this paper, we propose a Static Condensation Reduced Basis Element (SCRBE) approach for the Reynolds Lubrication Equation (RLE). The SCRBE method is a computational tool that allows to efficiently analyze parametrized structures which can be decomposed into a large number of similar components. Here, we extend the methodology to allow for a more general domain decomposition, a typical example being a checkerboard-pattern assembled from similar components. To this end, we extend the formulation and associated a posteriori error bound procedure. Our motivation comes from the analysis of the pressure distribution in plain journal bearings governed by the RLE. However, the SCRBE approach presented is not limited to bearings and the RLE, but directly extends to other component-based systems. We show numerical results for plain bearings to demonstrate the validity of the proposed approach.

In this study an explicit Finite Difference Method (FDM) based scheme is developed to solve the Maxwell's equations in time domain for a lossless medium. This manuscript focuses on two unique aspects – the three dimensional time-accurate discretization of the hyperbolic system of Maxwell equations in three-point non-staggered grid stencil and it's application to parallel computing through the use of Graphics Processing Units (GPU). The proposed temporal scheme is symplectic, thus permitting conservation of all Hamiltonians in the Maxwell equation. Moreover, to enable accurate predictions over large time frames, a phase velocity preserving scheme is developed for treatment of the spatial derivative terms. As a result, the chosen time increment and grid spacing can be optimally coupled. An additional theoretical investigation into this pairing is also shown. Finally, the application of the proposed scheme to parallel computing using one Nvidia K20 Tesla GPU card is demonstrated. For the benchmarks performed, the parallel speedup when compared to a single core of an Intel i7-4820K CPU is approximately 190x.

In this work we discuss the numerical discretization of the time-dependent Maxwell's equations using a fully explicit leap-frog type discontinuous Galerkin method. We present a sufficient condition for the stability and error estimates, for cases of typical boundary conditions, either perfect electric, perfect magnetic or first order Silver-Müller. The bounds of the stability region point out the influence of not only the mesh size but also the dependence on the choice of the numerical flux and the degree of the polynomials used in the construction of the finite element space, making possible to balance accuracy and computational efficiency. In the model we consider heterogeneous anisotropic permittivity tensors which arise naturally in many applications of interest. Numerical results supporting the analysis are provided.

We present a high-order discontinuous Galerkin (DG) method for the time domain Maxwell's equations in three-dimensional heterogeneous media. New hierarchical orthonormal basis functions on unstructured tetrahedral meshes are used for spatial discretization while Runge-Kutta methods for time discretization. A uniaxial perfectly matched layer (UPML) is employed to terminate the computational domain. Exponential convergence with respect to the order of the basis functions is observed and large parallel speedup is obtained for a plane-wave scattering model. The rapid decay of the out-going wave in the UPML is shown in a dipole radiation simulation. Moreover, the low frequency electromagnetic fields excited by a horizontal electric dipole (HED) and a vertical magnetic dipole (VMD) over a layered conductive half-space and a high frequency ground penetrating radar (GPR) detection for an underground structure are investigated, showing the high accuracy and broadband simulation capability of the proposed method.

The Jiles-Atherton (J-A) model is a commonly used physics-based model in describing the hysteresis characteristics of ferromagnetic materials. However, citations and interpretation of this model in literature have been non-uniform. Solution methods for solving numerically this model has not been studied adequately. In this paper, through analyzing the mathematical properties of equations and the physical mechanism of energy conservation, we point out some unreasonable descriptions of this model and develop a relatively more accurate, modified J-A model together with its numerical solution method. Our method employs a fixed point method to compute anhysteretic magnetization. We obtain the susceptibility value of the anhysteretic magnetization analytically and apply the 4th order Runge-Kutta method to the solution of total magnetization. Computational errors are estimated and then precisions of the solving method in describing various materials are verified. At last, through analyzing the effects of the accelerating method, iterative error and step size on the computational errors, we optimize the numerical method to achieve the effects of high precision and short computing time. From analysis, we determine the range of best values of some key parameters for fast and accurate computation.

This paper presents an approach using the method of separation of variables applied to 2D Helmholtz equations in the Cartesian coordinate. The solution is then computed by a series solutions resulted from solving a sequence of 1D problems, in which the 1D solutions are computed using pollution free difference schemes. Moreover, non-polluted numerical integration formulae are constructed to handle the integration due to the forcing term in the inhomogeneous 1D problems. Consequently, the computed solution does not suffer the pollution effect. Another attractive feature of this approach is that a direct method can be effectively applied to solve the tridiagonal matrix resulted from numerical discretization of the 1D Helmholtz equation. The method has been tested to compute 2D Helmholtz solutions simulating electromagnetic scattering from an open large cavity and rectangular waveguide.

Scale similarity is found in many natural phenomena in the universe, from fluid dynamics to astrophysics. In large eddy simulations of turbulent flows, some sub-grid scale (SGS) models are based on scale similarity. The earliest scale similarity SGS model was developed by Bardina et al., which produced SGS stresses with good correlation to the true stresses. In the present study, we perform a mathematical analysis of scale similarity. The analysis has revealed that the ratio of the resolved stress to the SGS stress is γ2, where γ is the ratio of the second filter width to the first filter width, under the assumption of small filter width. The implications of this analysis are discussed in the context of large eddy simulation.

We construct a nonlinear monotone finite volume scheme for three-dimensional diffusion equation on tetrahedral meshes. Since it is crucial important to eliminate the vertex unknowns in the construction of the scheme, we present a new efficient eliminating method. The scheme has only cell-centered unknowns and can deal with discontinuous or tensor diffusion coefficient problems on distorted meshes rigorously. The numerical results illustrate that the resulting scheme can preserve positivity on distorted tetrahedral meshes, and also show that our scheme appears to be approximate second-order accuracy for solution.

The 2D Maxwell eigenproblems are studied from a new point of view. An electromagnetic problem is cast from the Lagrangian system with single variable into the Hamiltonian system with dual variables. The electric and magnetic components transverse to the wave propagation direction are treated as dual variables to each other. Higher order curl-conforming and divergence-conforming vector basis functions are used to construct dual vector spectral elements. Numerical examples demonstrate some unique advantages of the proposed method.

We present a numerical method to solve the Vlasov-Poisson-Fokker-Planck (VPFP) system using the NRxx method proposed in [4, 7, 9]. A globally hyperbolic moment system similar to that in [23] is derived. In this system, the Fokker-Planck (FP) operator term is reduced into the linear combination of the moment coefficients, which can be solved analytically under proper truncation. The non-splitting method, which can keep mass conservation and the balance law of the total momentum, is used to solve the whole system. A numerical problem for the VPFP system with an analytic solution is presented to indicate the spectral convergence with the moment number and the linear convergence with the grid size. Two more numerical experiments are tested to demonstrate the stability and accuracy of the NRxx method when applied to the VPFP system.

Pressure-correction projection finite element methods (FEMs) are proposed to solve nonstationary natural convection problems in this paper. The first-order and second-order backward difference formulas are applied for time derivative, the stability analysis and error estimates of the semi-discrete schemes are presented using energy method. Compared with characteristic variational multiscale FEM, pressure-correction projection FEMs are more efficient and unconditionally energy stable. Ample numerical results are presented to demonstrate the effectiveness of the pressure-correction projection FEMs for solving these problems.

An idea of designing oscillation-less and high-resolution hybrid schemes is proposed and several types of hybrid schemes based on this idea are presented on block-structured grids. The general framework, for designing various types of hybrid schemes, is established using a Multi-dimensional Optimal Order Detection (MOOD) method proposed by Clain, Diot and Loubère [1]. The methodology utilizes low dissipation or dispersion but less robust schemes to update the solution and then implements robust and high resolution schemes to deal with problematic situations. A wide range of computational methods including central scheme, MUSCL scheme, linear upwind scheme and Weighted Essentially Non Oscillatory (WENO) scheme have been applied in the current hybrid schemes framework. Detailed numerical studies on classical test cases for the Euler system are performed, addressing the issues of the resolution and non-oscillatory property around the discontinuities.