Solution to the Fokker–Planck Equation with Piecewise-Constant Drift

We study the solution to the Fokker–Planck equation with piecewise-constant drift, taking the case with two jumps in the drift as an example. The solution in Laplace space can be expressed in closed analytic form, and its inverse can be obtained conveniently using some numerical inversion methods. The results obtained by numerical inversion can be regarded as exact solutions, enabling us to demonstrate the validity of some numerical methods for solving the Fokker–Planck equation. In particular, we use the solved problem as a benchmark example for demonstrating the fifth-order convergence rate of the finite difference scheme proposed previously [Chen Y and Deng X Phys. Rev. E 100(2019) 053303].     

An approach to estimating and extrapolating model error based on inverse problem methods:towards accurate numerical weather prediction

Model error is one of the key factors restricting the accuracy of numerical weather prediction(NWP). Considering the continuous evolution of the atmosphere, the observed data(ignoring the measurement error) can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers’ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.     

Remarks on interpretations of the Eotvos experiment and misinterpretation of E=mc^2

The Eotvos experiment on the verification of equivalence between inertial mass and gravitational mass of a body is famous for its accuracy. A question is, however, can these experimental results be applied to the case of a physical space in general relativity, where the space coordinates could be arbitrary？ It is pointed out that it can be validly applied because it has been proven that Einstein＇s equivalence principle for a physical space must have a frame of reference with the Euclidean-like structure. Will claimed further that such an overall accuracy can be translated into an accuracy of the equivalence between inertial mass and each type of energy. It is shown that, according to general relativity, such a claim is incorrect. The root of this problem is due to an inadequate understanding of special relativity that produced the famous equation E=mc^2, which must be understood in terms of energy conservation. Concurrently, it is pointed out that this error is a problem in Will＇s book, ‘Theory and Experiment in Gravitational Physics＇.     

Exponential B-spline collocation method for numerical solution of the generalized regularized long wave equation

The aim of the present paper is to present a numerical algorithm for the time-dependent generalized regularized long wave equation with boundary conditions. We semi-discretize the continuous problem by means of the Crank–Nicolson finite difference method in the temporal direction and exponential B-spline collocation method in the spatial direction. The method is shown to be unconditionally stable. It is shown that the method is convergent with an order of O(k2+ h2).Our scheme leads to a tri-diagonal nonlinear system. This new method has lower computational cost in comparison to the Sinc-collocation method. Finally, numerical examples demonstrate the stability and accuracy of this method.     

Moment searching algorithm for bioluminescence tomography

To avoid the ill-posedness in the inverse problem of bioluminescence tomography, a moment searching algorithm fusing the finite element method （FEM） with the moment concept in theoretical mechanics is developed. In the algorithm, the source＇s information is mapped to the surface photon flux density by FEM, and the source＇s position is modified with the feedback through the algorithm of barycenter searching, which makes full use of the position information of the photon flux density on surface. The position is modified in every iterative step and will finally converge to the real source＇s value theoretically.     

~ Transform Demonstration of Dark Soliton Solutions Found by Inverse Scattering

One of the basic problems about the inverse scattering transform for solving a completely integrable nonlinear evolutions equation is to demonstrate that the Jost solutions obtained from the inverse scattering equations of Cauchy integral satisfy the Lax equations. Such a basic problem still exists in the procedure of deriving the dark soliton solutions of the NLS equation in normal dispersion with non-vanishing boundary conditions through the inverse scattering transform. In this paper, a pair of Jost solutions with same analytic properties are composed to be a 2 × 2 matrix and then another pair are introduced to be its right inverse confirmed by the Liouville theorem. As they are both 2 × 2 matrices, the right inverse should be the left inverse too, based upon which it is not difficult to show that these Jost solutions satisfy both the first and second Lax equations. As a result of compatibility condition, the dark soliton solutions definitely satisfy the NLS equation in normal dispersion with non-vanishing boundary conditions.     

An application of radiation boundary condition to scattering problem of acoustic waves in fluids

A numerical method of solving acoustic wave scattering pnblemin fluids is described.Radiation boundary condition(RBC)obtained byfactorization method of Helmholtz equation is applied to transforming theexterior boundary value problem in unbounded region into one in a finiteregion.Combined with RBC and scatterer surface boundary condition,Helmholtz equation is solved numerically by the finite difference method.Computational results for sphere and prolate spheroidal scatterers are inexcellent agreement with eigenfunction solutions and much better than theresults of OSRC method.     

A method of solving the stiffness problem in Biot＇s poroelastic equations using a staggered high-order finite-difference

In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.     

Microscopic Theory of Blue Phases Ⅰ and Ⅱ of Liquid Crystal

The microscopic theory of the blue phases of chiral liquid crystal is proposed. Beginning with the potential between two molecules, by using the cell model of liquid, applying statistical physical method, the distribution function and the free energy of the system are obtained. By using variational approach and zero-order approximation, the differential equation that the order parameter tensor of the blue phase can satisfy is obtained. Then we change the differential equation to the eigenequation problem in quantum mechanics. Considering the symmetry of the blue phases, the order parameter tensors of blue phases Ⅰ and Ⅱ can be made up of the eigenvectors. Our results are the same as the results of Ginzberg-Landau‘s phenomenological theory. The parameters in the order parameter tensors that we calculate in the located system are close to the predecessors‘ results.     

Microscopic Theory of Blue Phases Ⅰ and Ⅱ of Liquid Crystal

The microscopic theory of the blue phases of chiral liquid crystal is proposed. Beginning with the potential between two molecules, by using the cell model of liquid, applying statistical physical method, the distribution function and the free energy of the system are obtained. By using variational approach and zero-order approximation, the differential equation that the order parameter tensor of the blue phase can satisfy is obtained. Then we change the differential equation to the eigenequation problem in quantum mechanics. Considering the symmetry of the blue phases,the order parameter tensors of blue phases Ⅰ and Ⅱ can be made up of the eigenvectors. Our results are the same as the results of Ginzberg-Landau‘s phenomenological theory. The parameters in the order parameter tensors that we calculate in the located system are close to the predecessors‘ results.     

An approximation for the boundary optimal control problem of a heat equation defined in a variable domain

In this paper, we consider a numerical approximation for the boundary optimal control problem with the control constraint governed by a heat equation defined in a variable domain. For this variable domain problem, the boundary of the domain is moving and the shape of theboundary is defined by a known time-dependent function. By making use of the Galerkin finite element method, we first project the original optimal control problem into a semi-discrete optimal control problem governed by a system of ordinary differential equations. Then, based on the aforementioned semi-discrete problem, we apply the control parameterization method to obtain an optimal parameter selection problem governed by a lumped parameter system, which can be solved as a nonlinear optimization problem by a Sequential Quadratic Programming （SQP） algorithm. The numerical simulation is given to illustrate the effectiveness of our numerical approximation for the variable domain problem with the finite element method and the control parameterization method.     

Capillary Adhesion of Microbeams： Finite Deformation Analysis

Capillary force may cause adhesion of devices at micro- and nano-scales. Considering the fact that large deformation is often involved in adhesion of microbeams, we analysed the capillary adhesion of two beams using finite deformation elasticity theory. The critical adhesion condition can be obtained from the present method as a function of the bending stiffness, Young＇s contact angle, the spacing of the two beams as well as the surface tensions of the solid and liquid phases. The solution for the capillary adhesion of a beam with a rigid substrate is also given. The results from the finite deformation analysis are compared with that of infinitesimal deformation method in order to show the necessity of accounting for the nonlinear effect associated with large deflection. The method adopted in this study can also be used to solve other adhesion problems associated with van der Waals force or electrostatic force.     

An improved boundary element-free method （IBEFM） for two-dimensional potential problems

The interpolating moving least-squares (IMLS) method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method (BEFM), combining the boundary integral equation (BIE) method with the IMLS method, the improved boundary element-free method (IBEFM) for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker δ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker δ function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.     

Robust H∞ observer-based of a class of complex control for synchronization dynamical networks

This paper is concerned with the robust Hoo synchronization problem for a class of complex dynamical networks by applying the observer-based control. The proposed feedback control scheme is developed to ensure the asymptotic stability of the augmented system, to reconstruct the non-measurable state variables of each node and to improve the H∞ performance related to the synchronization error and observation error despite the external disturbance. Based on the Lyapunov stability theory, a synchronization criterion is obtained under which the controlled network can be robustly stabilized onto a desired state with a guaranteed H∞ performance. The controller and the observer gains can be given by the feasible solutions of a set of linear matrix inequalities （LMIs）. The effectiveness of the proposed control scheme is demonstrated by a numerical example through simulation.     

Ion-acoustic waves in plasma of warm ions and isothermal electrons using time-fractional KdV equation

The ion-acoustic solitary wave in collisionless unmagnetized plasma consisting of warm ions-fluid and isothermal electrons is studied using the time fractional KdV equation. The reductive perturbation method has been employed to derive the Korteweg-de Vries equation for small but finite amplitude ion-acoustic wave in warm plasma. The Lagrangian of the time fractional KdV equation is used in a similar form to the Lagrangian of the regular KdV equation with fractional derivative for the time differentiation. The variation of the functional of this Lagrangian leads to the Euler-Lagrange equation that gives the time fractional KdV equation. The variational-iteration method is used to solve the derived time fractional KdV equation. The calculations of the solution are carried out for different values of the time fractional order. These calculations show that the time fractional can be used to modulate the electrostatic potential wave instead of adding a higher order dissipation term to the KdV equation. The results of the present investigation may be applicable to some plasma environments,such as the ionosphere plasma.     

The inverse problem based on a full dispersive wave equation

The inverse problem for harmonic waves and wave packets was studied based on a full dispersive wave equation.First,a full dispersive wave equation which describes wave propagation in nondissipative microstructured linear solids is established based on the Mindlin theory,and the dispersion characteristics are discussed.Second,based on the full dispersive wave equation,an inverse problem for determining the four unknown coefficients of wave equation is posed in terms of the frequencies and corresponding wave numbers of four different harmonic waves,and the inverse problem is demonstrated with rigorous mathematical theory. Research proves that the coefficients of wave equation related to material properties can be uniquely determined in cases of normal and anomalous dispersions by measuring the frequencies and corresponding wave numbers of four different harmonic waves which propagate in a nondissipative microstructured linear solids.     

Robust H_∞ observer-based control for synchronization of a class of complex dynamical networks

This paper is concerned with the robust H∞ synchronization problem for a class of complex dynamical networks by applying the observer-based control. The proposed feedback control scheme is developed to ensure the asymptotic stability of the augmented system, to reconstruct the non-measurable state variables of each node and to improve the H∞ performance related to the synchronization error and observation error despite the external disturbance. Based on the Lyapunov stability theory, a synchronization criterion is obtained under which the controlled network can be robustly stabilized onto a desired state with a guaranteed H∞ performance. The controller and the observer gains can be given by the feasible solutions of a set of linear matrix inequalities (LMIs). The effectiveness of the proposed control scheme is demonstrated by a numerical example through simulation.     

Application of Tikhonov regularization method to wind retrieval from scatterometer data Ⅱ: cyclone wind retrieval with consideration of rain

According to the conclusion of the simulation experiments in paper I, the Tikhonov regularization method is applied to cyclone wind retrieval with a rain-effect-considering geophysical model function (called GMF+Rain). The GMF+Rain model which is based on the NASA scatterometer-2 (NSCAT2) GMF is presented to compensate for the effects of rain on cyclone wind retrieval. With the multiple solution scheme (MSS), the noise of wind retrieval is effectively suppressed, but the influence of the background increases. It will cause a large wind direction error in ambiguity removal when the background error is large. However, this can be mitigated by the new ambiguity removal method of Tikhonov regularization as proved in the simulation experiments. A case study on an extratropical cyclone of hurricane observed with SeaWinds at 25-km resolution shows that the retrieved wind speed for areas with rain is in better agreement with that derived from the best track analysis for the GMF+Rain model, but the wind direction obtained with the two-dimensional variational (2DVAR) ambiguity removal is incorrect. The new method of Tikhonov regularization effectively improves the performance of wind direction ambiguity removal through choosing appropriate regularization parameters and the retrieved wind speed is almost the same as that obtained from the 2DVAR.     

Chaos Control and Anti-control via a Fuzzy Neural Network Inverse System Method

We propose a new method for chaos control and anti-control,which is referred to as the fuzzy-neural network inverse system method (FNNIS).The Sugeno-type fuzzy-neural network (FNN) is employed to learn the kinetics of the system to be controlled.Then the FNN model is used with the inverse system method to make the system to be controlled to track the reference input.If the system to be controlled is chaotic and the reference input is non-chaotic,chaos control can be implemented via the FNNIS method.If the system to be controlled is nonchaotic and the reference imput is chaotic,chaos anti-control can be implemented.Theorems about the effect of the FNN model error upon control are established.The simulation results show that this method is feasible and effective for chaos control and anti-control.