Compact implicit integration factor methods for some complex-valued nonlinear equations

The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.

Calculation methods of lubricant film pressure distribution of radial grooved thrust bearings

In order to calculate the pressure distribution of radial grooved thrust bearing, analytical and numerical methods were applied respectively. Grooved region and land region were linked by u- sing the mass conservations principle at the groove/land boundary in each method. The block-weight approach was implemented to deal with the non-coincidence of mesh and radial groove pattern in nu- merical method. It was observed that the numerical solutions had higher precision as mesh number exceed 70 x 70, and the relaxation iteration of differential scheme presented the fastest convergence speed when relaxation factor was close to 1.94.

Research Progress and Identification Method of Apple Stress Resistance

In the research, changes of apple chemistry, and molecule, under stresses, are n terms of morphology, physiology, bio- illustrated and research and identifica- tion methods of apple resistance are explored involving drought-resistance, flood-re- sistance, salt-stress resistance, cold-hardiness and heat-resistance. In addition prospects of apple resistance research are proposed, as well.

Applications of the first integral method to nonlinear evolution equations

In this paper, we establish travelling wave solutions for some nonlinear evolution equations. The first integral method is used to construct the travelling wave solutions of the modified Benjamin-Bona-Mahony and the coupled Klein-Gordon equations. The obtained results include periodic and solitary wave solutions. The first integral method presents a wider applicability to handling nonlinear wave equations.

Exact solutions of nonlinear fractional differential equations by (G'/G)-expansion method

In this paper, we use the fractional complex transform and the （G＇/G）-expansion method to study the nonlinear fractional differential equations and find the exact solutions. The fractional complex transform is proposed to convert a partial fractional differential equation with Jumarie＇s modified Riemann-Liouville derivative into its ordinary differential equation. It is shown that the considered transform and method are very efficient and powerful in solving wide classes of nonlinear fractional order equations.

A novel simulation method for positive corona current pulses

A novel two-dimensional （2D） simulation method of positive corona current pulses is proposed. A control-volume- based finite element method （CV-FEM） is used to solve continuity equations, and the Galerkin finite element method （FEM） is used to solve Poisson＇s equation. In the proposed method, photoionization is considered by adopting an exact Helmholtz photoionization model. Furthermore, fully implicit discretization and variable time step are used to ensure the time-efficiency of the present method. Finally, the method is applied to a positive rod-plane corona problem. The numerical results are in agreement with the experimental results, and the validity of the proposed method is verified.

Interval finite difference method for steady-state temperature field prediction with interval parameters

A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.

A GENERALIZED GAUSS-TYPE QUADRATURE FORMULA AND ITS APPLICATIONS TO PSEUDOSPECTRAL METHOD

A generalized Gauss-type quadrature formula is introduced, which assists in selection of collocation points in pseudospectral method for differential equations with two-point derivative boundary conditions. Some results on the related Jacobi interpolation are established. A pseudospectral scheme is proposed for the Kuramoto-Sivashisky equation. A skew symmetric decomposition is used for dealing with the nonlinear convection term. The stability and convergence of the proposed scheme are proved. The error estimates are obtained. Numerical results show the efficiency of this approach.

Variational iteration method for solving compressible Euler equations

This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.

Investigation of third-grade non-Newtonian blood flow in arteries under periodic body acceleration using multi-step differential transformation method

In this paper, a non-Newtonian third-grade blood in coronary and femoral arteries is simulated analytically and numerically. The blood is considered as the third- grade non-Newtonian fluid under the periodic body acceleration motion and the pulsatile pressure gradient. The hybrid multi-step differential transformation method （Hybrid-Ms- DTM） and the Crank-Nicholson method （CNM） are used to solve the partial differential equation （PDE）, and a good agreement between them is observed in the results. The effects of the some physical parameters such as the amplitude, the lead angle, and the body acceleration frequency on the velocity and shear stress profiles are considered. The results show that increasing the amplitude, Ag, and reducing the lead angle of body acceleration, 9, make higher velocity profiles on the center line of both arteries. Also, the maximum wall shear stress increases when Ag increases.

High order multiplication perturbation method for singular perturbation problems

This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.

Optimization of multibody systems based on the generalized-α projection method for DAEs

Efficient optimization strategy of multibody systems is developed in this paper. Aug- mented Lagrange method is used to transform constrained optimal problem into unconstrained form firstly. Then methods based on second order sensitivity are used to solve the unconstrained problem, where the sensitivity is solved by hybrid method. Generalized-α method and generalized-α projection method for the differential-algebraic equation, which shows more efficient properties with the lager time step, are presented to get state variables and adjoint variables during the optimization procedure. Numerical results validate the accuracy and efficiency of the methods is presented.

Discrete time transfer matrix method for dynamics of multibody system with flexible beams moving in space

In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space. Formulations and numerical example of a rigid- flexible-body three pendulums system moving in space are given to validate the method. Using the new method to study the dynamics of multi-rigid-flexible-body system mov- ing in space, the global dynamics equations of system are not needed, the orders of involved matrices of the system are very low and the computational speed is high, irrespec- tive of the size of the system. The new method is simple, straightforward, practical, and provides a powerful tool for multi-rigid-flexible-body system dynamics.

Infinitesimal dividing modeling method for dual suppliers inventory model with random lead times

As one of the basic inventory cost models, the （Q, τ）inventory cost model of dual suppliers with random procurement lead time is mostly formulated by using the concepts of ＂effective lead time＂ and ＂lead time demand＂, which may lead to an imprecise inventory cost. Through the real-time statistic of the inventory quantities, this paper considers the precise （Q, τ） inventory cost model of dual supplier procurement by using an infinitesimal dividing method. The traditional modeling method of the inventory cost for dual supplier procurement includes complex procedures. To reduce the complexity effectively, the presented method investigates the statistics properties in real-time of the inventory quantities with the application of the infinitesimal dividing method. It is proved that the optimal holding and shortage costs of dual supplier procurement are less than those of single supplier procurement respectively. With the assumption that both suppliers have the same distribution of lead times, the convexity of the cost function per unit time is proved. So the optimal solution can be easily obtained by applying the classical convex optimization methods. The numerical examples are given to verify the main conclusions.

Chirped Waves for a Generalized (2 + 1)-Dimensional Nonlinear Schrdinger Equation

The exact chirped soliton-like and quasi-periodic wave solutions of （2 ＋ 1）-dimensional generalized nonlinear Schr6dinger equation including linear and nonlinear gain （loss） with variable coefficients are obtained detalledly in this paper. The form and the behavior of solutions are strongly affected by the modulation of both the dispersion coefficient and the nonlinearity coefficient. In addition, self-similar soliton-like waves precisely piloted from our obtained solutions by tailoring the dispersion and linear gain （loss）.

Validation and application of three-dimensional discontinuous deformation analysis with tetrahedron finite element meshed block

In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.

Coiflet solution of strongly nonlinear p-Laplacian equations

A new boundary extension technique based on the Lagrange interpolat- ing polynomial is proposed and used to solve the function approximation defined on an interval by a series of scaling Coiflet functions, where the coefficients are used as the single-point samplings. The obtained approximation formula can exactly represent any polynomials defined on the interval with the order up to one third of the length of the compact support of the adopted Coiflet function. Based on the Galerkin method, a Coifiet-based solution procedure is established for general two-dimensional p^Laplacian equations, following which the equations can be discretized into a concise matrix form. As examples of applications, the proposed modified wavelet Galerkin method is applied to three typical p-Laplacian equations with strong nonlinearity. The numerical results justify the efficiency and accuracy of the method.

Bright and dark soliton solutions for some nonlinear fractional differential equations

In this work, we propose a new approach, namely ansatz method, for solving fractional differential equations based on a fractional complex transform and apply it to the nonlinear partial space-time fractional modified Benjamin-Bona- Mahoney （mBBM） equation, the time fractional mKdV equation and the nonlinear fractional Zoomeron equation which gives rise to some new exact solutions. The physical parameters in the soliton solutions： amplitude, inverse width, free parameters and velocity are obtained as functions of the dependent model coefficients. This method is suitable and more powerful for solving other kinds of nonlinear fractional PDEs arising in mathematical physics. Since the fractional deriva- tives are described in the modified Riemann-Liouville sense.

Modeling of Diffusion Transport through Oral Biofilms with the Inverse Problem Method

Aim The purpose of this study was to develop a mathe-matical model to quantitatively describe the passive trans-port of macromolecules within dental biofilms. Methodology Fluorescently labeled dextrans with different molecular mass （3 kD,10 kD,40 kD,70 kD,2 000 kD） were used as a series of diffusion probes. Streptococcus mutans,Streptococcus sanguinis,Actinomyces naeslundii and Fusobacterium nucleatum were used as inocula for biofilm formation. The diffusion processes of different probes through the in vitro biofilm were recorded with a confocal laser microscope. Results Mathematical function of biofilm penetration was constructed on the basis of the inverse problem method. Based on this function,not only the relationship between average concentration of steady-state and molecule weights can be analyzed,but also that between penetrative time and molecule weights. Conclusion This can be used to predict the effective concentration and the penetrative time of anti-biofilm medicines that can diffuse through oral biofilm. Further-more,an improved model for large molecule is proposed by considering the exchange time at the upper boundary of the dental biofilm.

THREE SOLUTIONS FOR A FRACTIONAL ELLIPTIC PROBLEMS WITH CRITICAL AND SUPERCRITICAL GROWTH

In this paper, we deal with the existence and multiplicity of solutions to the frac- tional elliptic problems involving critical and supercritical Sobolev exponent via variational arguments. By means of the truncation combining with the Moser iteration, we prove that our problem has at least three solutions.