Numerical simulations of two-dimensional (2D) turbulent thermal convection for inhomogeneous boundary condition are investigated using the lattice Boltzmann method (LBM). This study mainly appraises the temporal evolution and the scaling behavior of global quantities and of small-scale turbulence properties. The research results show that the flow is dominated by large-scale structures in the turbulence regime. Mushroom plumes emerge at both ends of each heat source, and smaller plumes increasingly rise. It is found that the gradient of root mean-square (rms) vertical velocities and the gradient of the rms temperature in the bottom boundary layer decreases with time evolution. It is further observed that the temporal evolution of the Kolmogorov scale, the kinetic-energy dissipation rates and thermal dissipation rates agree well with the theoretical predictions. It is also observed that there is a range of linear scaling in the 2nd-order structure functions of the velocity and temperature fluctuations and mixed velocity-temperature structure function.

The motion and rotation of an ellipsoidal particle inside square tubes and rectangular tubes with the confinement ratio R/a∈(1.0,4.0) are studied by the lattice Boltzmann method (LBM), where R and a are the radius of the tube and the semi-major axis length of the ellipsoid, respectively. The Reynolds numbers (Re) up to 50 are considered. For the prolate ellipsoid inside square and rectangular tubes, three typical stable motion modes which depend on R/a are identified, namely, the kayaking mode, the tumbling mode, and the log-rolling mode are identified for the prolate spheroid. The diagonal plane strongly attracts the particle in square tubes with 1.2≤R/a<3.0. To explore the mechanism, some constrained cases are simulated. It is found that the tumbling mode in the diagonal plane is stable because the fluid force acting on the particle tends to diminish the small displacement and will bring it back to the plane. Inside rectangular tubes the particle will migrate to a middle plane between short walls instead of the diagonal plane. Through the comparisons between the initial unstable equilibrium motion state and terminal stable mode, it is seems that the particle tend to adopt the mode with smaller kinetic energy.

In this paper, a new formulation is proposed to evaluate the origin intensity factors (OIFs) in the singular boundary method (SBM) for solving 3D potential problems with Dirichlet boundary condition. The SBM is a strong-form boundary discretization collocation technique and is mathematically simple, easy-to-program, and free of mesh. The crucial step in the implementation of the SBM is to determine the OIFs which isolate the singularities of the fundamental solutions. Traditionally, the inverse interpolation technique (IIT) is adopted to calculate the OIFs on Dirichlet boundary, which is time consuming for large-scale simulation. In recent years, the new methodology has been developed to efficiently calculate the OIFs on Neumann boundary, but the Dirichlet problem remains an open issue. This study employs the subtracting and adding-back technique based on the integration of the fundamental solution over the whole boundary to develop a new formulation of the OIFs on 3D Dirichlet boundary. Several problems with varied domain shapes and boundary conditions are carried out to validate the effectiveness and feasibility of the proposed scheme in comparison with the SBM based on inverse interpolation technique, the method of fundamental solutions, and the boundary element method.

We consider a scalar conservation law with zero-flux boundary conditions imposed on the boundary of a rectangular multidimensional domain. We study monotone schemes applied to this problem. For the Godunov version of the scheme, we simply set the boundary flux equal to zero. For other monotone schemes, we additionally apply a simple modification to the numerical flux. We show that the approximate solutions produced by these schemes converge to the unique entropy solution, in the sense of [7], of the conservation law. Our convergence result relies on a BV bound on the approximate numerical solution. In addition, we show that a certain functional that is closely related to the total variation is nonincreasing from one time level to the next. We extend our scheme to handle degenerate convection-diffusion equations and for the one-dimensional case we prove convergence to the unique entropy solution.

In this article, by applying the Stokes projection and an orthogonal projection with respect to curl and div operators, some new error estimates of finite element method (FEM) for the stationary incompressible magnetohydrodynamics (MHD) are obtained. To our knowledge, it is the first time to establish the error bounds which are explicitly dependent on the Reynolds numbers, coupling number and mesh size. On the other hand, The uniform stability and convergence of an Oseen type finite element iterative method for MHD with respect to high hydrodynamic Reynolds number Re and magnetic Reynolds number Rm, or small δ=1–σ with

(C0, C1 are constants depending only on Ω and F is related to the source terms of equations) are analyzed under the condition that . Finally, some numerical tests are presented to demonstrate the effectiveness of this algorithm.

In this paper, a simplified lattice Boltzmann method (SLBM) without evolution of the distribution function is developed for simulating incompressible viscous flows. This method is developed from the application of fractional step technique to the macroscopic Navier-Stokes (N-S) equations recovered from lattice Boltzmann equation by using Chapman-Enskog expansion analysis. In SLBM, the equilibrium distribution function is calculated from the macroscopic variables, while the non-equilibrium distribution function is simply evaluated from the difference of two equilibrium distribution functions. Therefore, SLBM tracks the evolution of the macroscopic variables rather than the distribution function. As a result, lower virtual memories are required and physical boundary conditions could be directly implemented. Through numerical test at high Reynolds number, the method shows very nice performance in numerical stability. An accuracy test for the 2D Taylor-Green flow shows that SLBM has the second-order of accuracy in space. More benchmark tests, including the Couette flow, the Poiseuille flow as well as the 2D lid-driven cavity flow, are conducted to further validate the present method; and the simulation results are in good agreement with available data in literatures.

We study the numerical identification of an unknown portion of the boundary on which either the Dirichlet or the Neumann condition is provided from the knowledge of Cauchy data on the remaining, accessible and known part of the boundary of a two-dimensional domain, for problems governed by Helmholtz-type equations. This inverse geometric problem is solved using the plane waves method (PWM) in conjunction with the Tikhonov regularization method. The value for the regularization parameter is chosen according to Hansen's L-curve criterion. The stability, convergence, accuracy and efficiency of the proposed method are investigated by considering several examples.

The paper presents the size-dependant behaviors of the carbon nanotubes under electrostatic actuation using the modified couple stress theory and homotopy perturbation method. Due to the less accuracy of the classical elasticity theorems, the modified couple stress theory is applied in order to capture the size-dependant properties of the carbon nanotubes. Both of the static and dynamic behaviors under static DC and step DC voltages are discussed. The effects of various dimensions and boundary conditions on the deflection and pull-in voltages of the carbon nanotubes are to be investigated in detail via application of the homotopy perturbation method to solve the nonlinear governing equations semi-analytically.

A new sharp interface method with the combination of Ghost Fluid Method (GFM) and Cut Cell scheme is developed to study compressible multi-phase flows with clear interfaces. Straight-line cutting is applied on the cells passed by the interface. A new real-ghost mixing method is presented and applied around the cut cells to deal with very small cut cells. A cut face reconstruction method similar to volume of fluid is applied to deal with topological change problems. A high order Level Set (LS) method is applied to evolve the free interface, with the Level Set velocities from exact Riemann solver on the cut faces. Various 1D and 2D numerical examples are tested to show the robustness and ability of the present method in wide flow variable domains. This method benefits from cut cell on the sharp interface description, shows good conservation performance, and does not have the topological change difficulty of the full cut cell method presented in Chang, Deng & Theofanous, J. Comput. Phys., 242 (2013), pp. 946–990.

In this paper, we consider the transform magnetic (TM) model of electromagnetic scattering in the cavity. By the Polynomial Preserving Recovery technique, we present superconvergence analysis for the vertex-edge-face type finite element. From the numerical example, we can see that the provided method is efficient and stable.

In this paper, we propose a new energy-preserving scheme and a new momentum-preserving scheme for the modified regularized long wave equation. The proposed schemes are designed by using the discrete variational derivative method and the finite volume element method. For comparison, we also propose a finite volume element scheme. The conservation properties of the proposed schemes are analyzed and we find that the energy-preserving scheme can precisely conserve the discrete total mass and total energy, the momentum-preserving scheme can precisely conserve the discrete total mass and total momentum, while the finite volume element scheme merely conserve the discrete total mass. We also analyze their linear stability property using the Von Neumann theory and find that the proposed schemes are unconditionally linear stable. Finally, we present some numerical examples to illustrate the effectiveness of the proposed schemes.

In the paper, a new time second-order mass-conserved implicit/explicit domain decomposition method (DDM) for the diffusion equations is proposed. In the scheme, firstly, we calculate the interface fluxes of sub-domains from the obtained solutions and fluxes at the previous time level, for which we apply high-order Taylor’s expansion and transfer the time derivatives to spatial derivatives to improve the accuracy. Secondly, the interior solutions and fluxes in sub-domains are computed by the implicit scheme and using the relations between solutions and fluxes, without any correction step. The mass conservation is proved and the convergence order of the numerical solutions is proved to be second-order in both time and space steps. The super-convergence of numerical fluxes is also proved to be second-order in both time and space steps. The scheme is stable under the stable condition r≤3/5. The important feature is that the proposed domain decomposition scheme is mass-conserved and is of second order convergence in time. Numerical experiments confirm the theoretical results.

Numerical investigation of the underexpanded sonic coaxial jets is carried out using large eddy simulation for three typical inner nozzle lip-thicknesses. Various fundamental mechanisms dictating the flow phenomena including shock structure, shear layer evolution and sound production are investigated. It is found that the inner nozzle lip induces a recirculation zone between inner and outer jets, which significantly influences the behaviors of shock structures and shear layers. The sound properties of the coaxial jets are further analyzed in detail. As the inner lip-thickness increases, the helical screech mode switches to an axisymmetric one and high-frequency screech also occurs with an oscillation frequency of recirculation zone. Based on the temporal Fourier transform and correlation analysis, the primary sources of low- and high-frequency screeches are associated with the downstream shock cells in the jet column and the secondary shock structures in the outer annular jet, respectively. The proper orthogonal decomposition analysis reveals that the dominant structures constructed by the most energetic modes shift from the downstream shock cells region to the upstream secondary shock region as the lip-thickness increases. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the coherent structures and sound properties in sonic coaxial jets.

In order to overcome the possible singularity associated with the Point Interpolation Method (PIM), the Radial Point Interpolation Method (RPIM) was proposed by G. R. Liu. Radial basis functions (RBF) was used in RPIM as basis functions for interpolation. All these radial basis functions include shape parameters. The choice of these shape parameters has been and stays a problematic theme in RBF approximation and interpolation theory. The object of this study is to contribute to the analysis of how these shape parameters affect the accuracy of the radial PIM. The RPIM is studied based on the global Galerkin weak form performed using two integration technics: classical Gaussian integration and the strain smoothing integration scheme. The numerical performance of this method is tested on their behavior on curve fitting, and on three elastic mechanical problems with regular or irregular nodes distributions. A range of recommended shape parameters is obtained from the analysis of different error indexes and also the condition number of the matrix system. All resulting RPIM methods perform very well in term of numerical computation. The Smoothed Radial Point Interpolation Method (SRPIM) shows a higher accuracy, especially in a situation of distorted node scheme.

A new simple mathematical method has been proposed to predict rock stress around a noncircular tunnel and the method is calibrated and validated with a numerical model. It can be found that the tunnel shapes and polar angles affect the applicable zone of the theoretical model significantly and the applicable zone of a rectangular tunnel was obtained using this method. The method can be used to predict the values of the concentrated stress, and to analyze the change rate of rock stress and back to calculate the mechanical boundary condition in the applicable zone. The results of the stress change rate indicate that the horizontal stress is negatively related to the vertical boundary load and positively related to the horizontal boundary load. The vertical stress is negatively related to the horizontal boundary load and positively related to the vertical boundary load. These findings can be used to explain the evolution of the vertical increment in stress obtained with field-based borehole stress monitoring.

The Stokes axisymmetric flow of an incompressible viscous fluid past a micropolar fluid spheroid whose shape deviates slightly from that of a sphere is studied analytically. The boundary conditions used are the vanishing of the normal velocities, the continuity of the tangential velocities, continuity of shear stresses and spin-vorticity relation at the surface of the spheroid. The hydrodynamic drag force acting on the fluid spheroid is calculated. An exact solution of the problem is obtained to the first order in the small parameter characterizing the deformation. It is observed that due to increase spin parameter value, the drag coefficient decreases. Well known results are deduced and comparisons are made with classical viscous fluid and micropolar fluid.

In this work, the a-posteriori error indicator with an explicit formula for p-version finite element methods in square is investigated, and its reliable and efficient properties are deduced. Especially, this a-posteriori error indicator is determined by the right hand itemof themodel. We reformulate this a-posteriori error indicator with finite coefficients, which can be easily calculated during applications.

In this paper, we first discuss the well-posedness of linearizing equations, and then study the stability and unstability of the 3-D compressible Euler Equation, by analysing the existence of saddle point. In addition, we give the existence of local solutions of the compressible Euler equation.

This paper proposes an extrapolation cascadic multigrid (EXCMG) method to solve elliptic problems in domains with reentrant corners. On a class of λ-graded meshes, we derive some new extrapolation formulas to construct a high-order approximation to the finite element solution on the next finer mesh using the numerical solutions on two-level of grids (current and previous grids). Then, this high-order approximation is used as the initial guess to reduce computational cost of the conjugate gradient method. Recursive application of this idea results in the EXCMG method proposed in this paper. Finally, numerical results for a crack problem and an L-shaped problem are presented to verify the efficiency and effectiveness of the proposed EXCMG method.

In this paper, a double-passage shape correction (DPSC) method is presented for simulation of unsteady flows around vibrating blades and aeroelastic prediction. Based on the idea of phase-lagged boundary conditions, the shape correction method was proposed aimed at efficiently dealing with unsteady flow problems in turbomachinery. However, the original single-passage shape correction (SPSC) may show the disadvantage of slow convergence of unsteady solutions and even produce nonphysical oscillation. The reason is found to be related with the disturbances on the circumferential boundaries that can not be damped by numerical schemes. To overcome these difficulties, the DPSC method is adopted here, in which the Fourier coefficients are computed from flow variables at implicit boundaries instead of circumferential boundaries in the SPSC method. This treatment actually reduces the interaction between the calculation of Fourier coefficients and the update of flow variables. Therefore a faster convergence speed could be achieved and also the solution stability is improved. The present method is developed to be suitable for viscous and turbulent flows. And for real three-dimensional (3D) problems, the rotating effects are also considered. For validation, a 2D oscillating turbine cascade, a 3D oscillating flat plate cascade and a 3D practical transonic fan rotor are investigated. Comparisons with experimental data or other solutions and relevant discussions are presented in detail. Numerical results show that the solution accuracy of DPSC method is favorable and at least comparable to the SPSC method. However, fewer iteration cycles are needed to get a converged and stable unsteady solution, which greatly improves the computational efficiency.