Solving diffusion equation using wavelet method

A wavelet-based numerical method based on wavelets of Hermite cubic splines is presented for computing singularly perturbed convection-dominated diffusion equation. The advantages of the method are explained. To improve the accuracy of singular areas, wavelets are configured hierarchically for solving algebraic equations. Numerical examples show that the proposed method has good efficiency and precision.     

APPLICATION OF TWO OPERATORS TRANSFORM FROM CLASS A OPERATOR TO THE CLASS OF HYPONORMAL OPERATOR

In this article, we give an operator transform T (*) from class A operator to the class of hyponormal operators. It is different from the operator transform T defined by M. Ch and T. Yamazaki. Then, we show that σ(T ) = σ( T (*)) and σa(T )\{0} = σa( T (*))\{0}, in case T belongs to class A. Next, we obtain some relations between T and T (9).     

L~p CONTINUITY OF HRMANDER SYMBOL OPERATORS OpS_（0,0） ~m AND NUMERICAL ALGORITHM

If we use Littlewood-Paley decomposition, there is no pseudo-orthogonality for Ho¨rmander symbol operators OpS m 0 , 0 , which is different to the case S m ρ,δ (0 ≤δ < ρ≤ 1). In this paper, we use a special numerical algorithm based on wavelets to study the L p continuity of non infinite smooth operators OpS m 0 , 0 ; in fact, we apply first special wavelets to symbol to get special basic operators, then we regroup all the special basic operators at given scale and prove that such scale operator’s continuity decreases very fast, we sum such scale operators and a symbol operator can be approached by very good compact operators. By correlation of basic operators, we get very exact pseudo-orthogonality and also L 2 → L 2 continuity for scale operators. By considering the influence region of scale operator, we get H 1 (= F 0 , 2 1 ) → L 1 continuity and L ∞→ BMO continuity. By interpolation theorem, we get also L p (= F 0 , 2 p ) → L p continuity for 1 < p < ∞ . Our results are sharp for F 0 , 2 p → L p continuity when 1 ≤ p ≤ 2, that is to say, we find out the exact order of derivations for which the symbols can ensure the resulting operators to be bounded on these spaces.     

STABILITY OF EQUILIBRIA OF NEMATIC LIQUID CRYSTALLINE POLYMERS

We provide an analytical study on the stability of equilibria of rigid rodlike nematic liquid crystalline polymers (LCPs) governed by the Smoluchowski equation with the Maier-Saupe intermolecular potential. We simplify the expression of the free energy of an orientational distribution function of rodlike LCP molecules by properly selecting a coordinate system and then investigate its stability with respect to perturbations of orientational probability density. By computing the Hessian matrix explicitly, we are able to prove the hysteresis phenomenon of nematic LCPs: when the normalized polymer concentration b is below a critical value b*(6.7314863965), the only equilibrium state is isotropic and it is stable; when b*< b < 15/2, two anisotropic (prolate) equilibrium states occur together with a stable isotropic equilibrium state. Here the more aligned prolate state is stable whereas the less aligned prolate state is unstable. When b > 15/2, there are three equilibrium states: a stable prolate state, an unstable isotropic state and an unstable oblate state.     

Iterative operator‐splitting methods for nonlinear differential equations and applications

In this article, we consider iterative operator‐splitting methods for nonlinear differential equations with bounded and unbounded operators. The main feature of the proposed idea is the embedding of Newton's method for solving the split parts of the nonlinear equation at each step. The convergence properties of such a mixed method are studied and demonstrated. We confirm with numerical applications the effectiveness of the proposed scheme in comparison with the standard operator‐splitting methods by providing improved results and convergence rates. We apply our results to deposition processes. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1026–1054, 2011     

A nonstandard finite difference scheme for contaminant transport with kinetic langmuir sorption

This study focuses on a contaminant transport model with Langmuir sorption under nonequilibrium conditions. The numerical instabilities of the standard finite difference schemes including the upwind method are investigated. By using the nonstandard finite difference method, a better finite difference model is constructed. The numerical simulation on a specific system configuration proves the advantages of the new finite difference model. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 767–785, 2011     

甲流病毒传播的年龄结构模型的稳定性

研究甲型H1N1流感病毒的传播规律,建立年龄结构具有接种措施的SEIR流行病模型,给出了疾病流行的阈值并证明了地方病平衡点的稳定性问题.最后根据一些实际数据,进行数值模拟进而对模型的合理性加以完善,借助模型预测下一阶段甲流爆发的可能性并提出相关应对措施.     

Accurate 3-D finite element simulation of elastic wave propagation with the combination of explicit and implicit time-integration methods

One of the big issues in finite element solutions of wave propagation problems is the presence of spurious high-frequency oscillations that may lead to divergent results at mesh refinement. The paper deals with the extension of the new two-stage time-integration technique developed in our previous papers to the solution of wave propagation problems with explicit time-integration methods.The explicit central difference method is used for accurate time-integration of the semi-discrete system of elastodynamics at the stage of basic computations and allows spurious high-frequency oscillations. To filter these oscillations, pre- or/and post-processing (the filtering stage) is applied using a few time increments of the implicit time-continuous Galerkin method with large numerical dissipation.A special calibration procedure is used for the selection of the minimum necessary amount of numerical dissipation (in terms of a time increment) at the filtering stage. In contrast to existing approaches that use a time-integration method with the same dissipation (or artificial viscosity) for all time increments, the new technique yields accurate and non-oscillatory results for wave propagation problems without interaction between user and computer code. The solutions of 3-D wave propagation and impact problems show the effectiveness of the new approach.     

磁流体压缩管道热电磁流动特性研究

利用磁流体五波模型对低磁雷诺数下压缩管道中磁流体流动进行数值模拟。该模型由带有电磁作用强制项的Navier-Stokes方程组与电势Poisson方程组成,数值格式分别采用严格保证熵条件的熵条件格式和中心差分格式。数值模拟对不同磁作用数下的不同几何外形管道进行数值模拟研究,结果表明在磁流体压缩管道中,由于发生器模式提取...     

Coquet-type formulas for the rarefied weighted Thue–Morse sequence

Newman proved for the classical Thue–Morse sequence, ((−1)^s(n))_n≥0, that for all N∈N with real constants satisfying c₂>c₁>0 and λ=log3/log4. Coquet improved this result and deduced , where F(x) is a nowhere-differentiable, continuous function with period 1 and η(N)∈{−1,0,1}. In this paper we obtain for the weighted version of the Thue–Morse sequence that for the sum a Coquet-type formula exists for every r∈{0,1,2} if and only if the sequence of weights is eventually periodic. From the specific Coquet-type formulas we derive parts of the weak Newman-type results that were recently obtained by Larcher and Zellinger.