Velocity–vorticity–helicity formulation and a solver for the Navier–Stokes equations

For the three-dimensional incompressible Navier–Stokes equations, we present a formulation featuring velocity, vorticity and helical density as independent variables. We find the helical density can be observed as a Lagrange multiplier corresponding to the divergence-free constraint on the vorticity variable, similar to the pressure in the case of the incompressibility condition for velocity. As one possible practical application of this new formulation, we consider a time-splitting numerical scheme based on an alternating procedure between vorticity–helical density and velocity–Bernoulli pressure systems of equations. Results of numerical experiments include a comparison with some well-known schemes based on pressure–velocity formulation and illustrate the competitiveness on the new scheme as well as the soundness of the new formulation.     

Stability and bifurcation analysis in the delay-coupled van der Pol oscillators

In this paper, the dynamics of a system of two van der Pol equations with a finite delay are investigated. We show that there exist the stability switches and a sequence of Hopf bifurcations occur at the zero equilibrium when the delay varies. Using the theory of normal form and the center manifold theorem, the explicit expression for determining the direction of the Hopf bifurcations and the stability of the bifurcating periodic solutions are derived.     

求解二维Poisson方程的重心插值配点法

采用重心Lagrange插值配点法计算了二维Poisson方程.采用重心Lagrange插值法构造近似函数,由配点法离散Poisson方程及其边界条件.数值算例表明方法具有理论简单、计算精度高的特点.     

Dissipative Time Dependent Density Functional Theory

The simplest density functional theory due to Thomas, Fermi, Dirac and Weizsäcker is employed to describe the non-equilibrium thermodynamic evolution of an electron gas. The temperature effect is introduced via the Fermi-Dirac entropy, while the irreversible dynamics is described by a non-linear diffusion equation. A dissipative Kohn-Sham equation is also proposed, which improves the Thomas-Fermi-Weizsäcker kinetic functional.     

One-dimensional photonic crystal design

In this article we present a method to determine the band spectrum, band gaps, and discrete energy levels, of a one-dimensional photonic crystal with localized impurities. For one-dimensional crystals with piecewise constant refractive indices we develop an algorithm to recover the refractive index distribution from the period map. Finally, we derive the relationship between the period map and the scattering matrix containing the information on the localized modes.     

Traveling waves for a reaction–diffusion–advection predator–prey model

In this paper we study a reaction–diffusion–advection predator–prey model in a river. The existence of predator-invasion traveling wave solutions and prey-spread traveling wave solutions in the upstream and downstream directions is established and the corresponding minimal wave speeds are obtained. While some crucial improvements in theoretical methods have been established, the proofs of the existence and nonexistence of such traveling waves are based on Schauder’s fixed-point theorem, LaSalle’s invariance principle and Laplace transform. Based on theoretical results, we investigate the effect of the hydrological and biological factors on minimal wave speeds and hence on the spread of the prey and the invasion of the predator in the river. The linear determinacy of the predator–prey Lotka–Volterra system is compared with nonlinear determinacy of the competitive Lotka–Volterra system to investigate the mechanics of linear and nonlinear determinacy.     

Development of a multi‐scale ocean model by using particle Laplacian method for anisotropic mass transfer

The direct injection of CO₂ into the deep ocean is one of the feasible ways for the mitigation of the global warming, although there is a concern about its environmental impact near the injection point. To minimize its biological impact, it is necessary to make CO₂ disperse as quickly as possible, and it is said that injection with a pipe towed by a moving ship is effective for this purpose. Because the injection ship moves over a spatial scale of O(10²km), a mesoscale model is necessary to analyse the dispersion of CO₂. At the same time, since it is important to investigate high CO₂ concentration near the injection point, a small‐scale model is also required. Therefore, in this study, a numerical model was developed to analyse CO₂ dispersion in the deep ocean by using a fixed mesoscale and a moving small‐scale grid systems, the latter of which is nested and moves in the former along the trajectory of the moving ship. To overcome the artificial diffusion of mass concentration at the interface of the two different grid systems and to keep its spatial accuracy almost the same as that in the small‐scale, a particle Laplacian method was adopted and newly modified for anisotropic diffusion in the ocean. Copyright © 2010 John Wiley & Sons, Ltd.     

The multi-concentration and two-dimensional capillary electrophoresis method for the analysis of drugs in urine samples

A novel method has been developed by integration of multi-concentration and two-dimensional(2D) capillary electrophoresis(CE) for simultaneous enhancement of detection sensitivity and separation power in complex samples.Capillary zone electrophoresis(CZE) was used as the first dimension separation according to mobilities,from which the effluent fractions were further analyzed by micellar electrokinetic capillary chromatography(MEKC) acting as the second dimension.Cation-selective exhaustive injection(CSEI) ...     

A well‐balanced upstream flux‐splitting finite‐volume scheme for shallow‐water flow simulations with irregular bed topography

This study extends the upstream flux‐splitting finite‐volume (UFF) scheme to shallow water equations with source terms. Coupling the hydrostatic reconstruction method (HRM) with the UFF scheme achieves a resultant numerical scheme that adequately balances flux gradients and source terms. The proposed scheme is validated in three benchmark problems and applied to flood flows in the natural/irregular river with bridge pier obstructions. The results of the simulations are in satisfactory agreement with the available analytical solutions, experimental data and field measurements. Comparisons of the present results with those obtained by the surface gradient method (SGM) demonstrate the superior stability and higher accuracy of the HRM. The stability test results also show that the HRM requires less CPU time (up to 60%) than the SGM. The proposed well‐balanced UFF scheme is accurate, stable and efficient to solve flow problems involving irregular bed topography. Copyright © 2009 John Wiley & Sons, Ltd.     

Penalty finite element method for Navier–Stokes equations with nonlinear slip boundary conditions

The penalty finite element method for Navier–Stokes equations with nonlinear slip boundary conditions is investigated in this paper. Since this class of nonlinear slip boundary conditions include the subdifferential property, the weak variational formulation is a variational inequality problem of the second kind. Using the penalty finite element approximation, we obtain optimal error estimates between the exact solution and the finite element approximation solution. Finally, we show the numerical results which are in full agreement with the theoretical results. Copyright © 2011 John Wiley & Sons, Ltd.