Immersed boundary method for unsteady kinetic model equations

Predicting unsteady flows and aerodynamic forces for large displacement motion of microstructures requires transient solution of Boltzmann equation with moving boundaries. For the inclusion of moving complex boundaries for these problems, three immersed boundary method flux formulations (interpolation, relaxation, and interrelaxation) are presented. These formulations are implemented in a 2‐D finite volume method solver for ellipsoidal‐statistical (ES)‐Bhatnagar‐Gross‐Krook (BGK) equations using unstructured meshes. For the verification, a transient analytical solution for free molecular 1‐D flow is derived, and results are compared with the immersed boundary (IB)‐ES‐BGK methods. In 2‐D, methods are verified with the conformal, non‐moving finite volume method, and it is shown that the interrelaxation flux formulation gives an error less than the interpolation and relaxation methods for a given mesh size. Furthermore, formulations applied to a thermally induced flow for a heated beam near a cold substrate show that interrelaxation formulation gives more accurate solution in terms of heat flux. As a 2‐D unsteady application, IB/ES‐BGK methods are used to determine flow properties and damping forces for impulsive motion of microbeam due to high inertial forces. IB/ES‐BGK methods are compared with Navier–Stokes solution at low Knudsen numbers, and it is shown that velocity slip in the transitional rarefied regime reduces the unsteady damping force. Copyright © 2015 John Wiley & Sons, Ltd.     

Binding numbers for fractional ID-k-factor-critical graphs

LetG be a graph,and k≥2 be a positive integer.A graph G is fractional independentset-deletable k-factor-critical(in short,fractional ID-k-factor-critical),if G I has a fractional k-factor for every independent set I of G.The binding number bind(G)of a graph G is defined as bind(G)=min|NG(X)||X|:=X V(G),NG(X)=V(G).In this paper,it is proved that a graph G is fractional ID-k-factor-critical if n≥6k 9 and bind(G)(3k 1)(n 1)kn 2k+2.     

Maximum likelihood estimation for second level fMRI data analysis with expectation trust region algorithm

The trust region method which originated from the Levenberg–Marquardt (LM) algorithm for mixed effect model estimation are considered in the context of second level functional magnetic resonance imaging (fMRI) data analysis. We first present the mathematical and optimization details of the method for the mixed effect model analysis, then we compare the proposed methods with the conventional expectation-maximization (EM) algorithm based on a series of datasets (synthetic and real human fMRI datasets). From simulation studies, we found a higher damping factor for the LM algorithm is better than lower damping factor for the fMRI data analysis. More importantly, in most cases, the expectation trust region algorithm is superior to the EM algorithm in terms of accuracy if the random effect variance is large. We also compare these algorithms on real human datasets which comprise repeated measures of fMRI in phased-encoded and random block experiment designs. We observed that the proposed method is faster in computation and robust to Gaussian noise for the fMRI analysis. The advantages and limitations of the suggested methods are discussed.     

Cantor-type cylindrical-coordinate method for differential equations with local fractional derivatives

In this Letter, we propose to use the Cantor-type cylindrical-coordinate method in order to investigate a family of local fractional differential operators on Cantor sets. Some testing examples are given to illustrate the capability of the proposed method for the heat-conduction equation on a Cantor set and the damped wave equation in fractal strings. It is seen to be a powerful tool to convert differential equations on Cantor sets from Cantorian-coordinate systems to Cantor-type cylindrical-coordinate systems.     

Resonant control of solitons

It is shown that the effect of “scattering on resonance” can be used to control envelope solitons in the driven nonlinear Schrödinger equation. The control occurs by the frequency modulated driving with multiple crossing of the resonant frequency of the soliton.     

Interbasis expansions for isotropic harmonic oscillator

The exact solutions of the isotropic harmonic oscillator are reviewed in Cartesian, cylindrical polar and spherical coordinates. The problem of interbasis expansions of the eigenfunctions is solved completely. The explicit expansion coefficients of the basis for given coordinates in terms of other two coordinates are presented for lower excited states. Such a property is occurred only for those degenerated states for given principal quantum number n.     

Lattice Boltzmann method for the fractional sub‐diffusion equation

A lattice Boltzmann model for the fractional sub‐diffusion equation is presented. By using the Chapman–Enskog expansion and the multiscale time expansion, several higher‐order moments of equilibrium distribution functions and a series of partial differential equations in different time scales are obtained. Furthermore, the modified partial differential equation of the fractional sub‐diffusion equation with the second‐order truncation error is obtained. In the numerical simulations, comparisons between numerical results of the lattice Boltzmann models and exact solutions are given. The numerical results agree well with the classical ones. Copyright © 2015 John Wiley & Sons, Ltd.     

Equilibrium analysis for a mass-conserving model in presence of cavitation

We study the existence of equilibrium positions for the load problem in Lubrication Theory. The problem consists of two surfaces in relative motion separated by a small distance filled by a lubricant. The system is described by the modified Reynolds equation (Elrod–Adams model) which describes the behavior of the lubricant and an extra integral equation given the balance of forces. The balance of forces allows to obtain the unknown position of the surfaces, defined with one degree of freedom.     

Properties of ε-value functions for a two person differential game

We consider the problem of approximate minimax for the Bolza problem of optimal control. Starting from the method of dynamic programming (Bellman) we define the ?-value function to be the approximation for the value function being a solution to the Hamilton–Jacobi equation.     

变形谐振子势的能谱方程

讨论了3种变形谐振子势:左右两边不同参数的谐振子势、左边方形势加右边谐振子势和谐振子势中间加δ势中的能量本征态函数.这些函数都可以由厄米函数表示.由波函数及其一次导数在原点的衔接条件,得到了能谱方程.