The effect of variable viscosity on MHD viscoelastic fluid flow and heat transfer over a stretching sheet

An analysis has been carried out to study the momentum and heat transfer characteristics in an incompressible electrically conducting non-Newtonian boundary layer flow of a viscoelastic fluid over a stretching sheet. The partial differential equations governing the flow and heat transfer characteristics are converted into highly non-linear coupled ordinary differential equations by similarity transformations. The effect of variable fluid viscosity, Magnetic parameter, Prandtl number, variable thermal conductivity, heat source/sink parameter and thermal radiation parameter are analyzed for velocity, temperature fields, and wall temperature gradient. The resultant coupled highly non-linear ordinary differential equations are solved numerically by employing a shooting technique with fourth order Runge–Kutta integration scheme. The fluid viscosity and thermal conductivity, respectively, assumed to vary as an inverse and linear function of temperature. The analysis reveals that the wall temperature profile decreases significantly due to increase in magnetic field parameter. Further, it is noticed that the skin friction of the sheet decreases due to increase in the Magnetic parameter of the flow characteristics.     

A Two-Grid Finite Element Approximation for Nonlinear Time Fractional Two-Term Mixed Sub-Diffusion and Diffusion Wave Equations

In this paper, we develop a two-grid method (TGM) based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations. A two-grid algorithm is proposed for solving the nonlinear system, which consists of two steps: a nonlinear FE system is solved on a coarse grid, then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution. The fully discrete numerical approximation is analyzed, where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with order $\alpha\in(1,2)$ and $\alpha_{1}\in(0,1)$. Numerical stability and optimal error estimate $O(h^{r+1}+H^{2r+2}+\tau^{\min\{3-\alpha,2-\alpha_{1}\}})$ in $L^{2}$-norm are presented for two-grid scheme, where $t,$ $H$ and $h$ are the time step size, coarse grid mesh size and fine grid mesh size, respectively. Finally, numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.     

Fully Discrete $A$-$ϕ$ Finite Element Method for Maxwell's Equations with a Nonlinear Boundary Condition

In this paper we present a fully discrete $A$-$ϕ$ finite element method to solve Maxwell's equations with a nonlinear degenerate boundary condition, which represents a generalization of the classical Silver-Müller condition for a non-perfect conductor. The relationship between the normal components of the electric field $E$ and the magnetic field $H$ obeys a power-law nonlinearity of the type $H× n = n × (|E× n|^{α-1}E × n)$ with $α ∈ (0,1]$. We prove the existence and uniqueness of the solutions of the proposed $A$-$ϕ$ scheme and derive the error estimates. Finally, we present some numerical experiments to verify the theoretical result.     

Fully discrete T‐ψ finite element method to solve a nonlinear induction hardening problem

We study an induction hardening model described by Maxwell's equations coupled with a heat equation. The magnetic induction field is assumed a nonlinear constitutional relation and the electric conductivity is temperature‐dependent. The T ‐ψ method is to transform Maxwell's equations to the vector–scalar potential formulations and to solve the potentials by means of the finite element method. In this article, we present a fully discrete T ‐ψ finite element scheme for this nonlinear coupled problem and discuss its solvability. We prove that the discrete solution converges to a weak solution of the continuous problem. Finally, we conclude with two numerical experiments for the coupled system.     

Banach空间一类非线性微分积分方程的整体解

该文研究 Banach空间中一类非线性 Volterra型微分积分方程在无穷区间 R 上的耦合最小最大拟解及解的整体存在性 .利用单调迭代方法及 Monch不动点定理 ,给出了该类方程耦合最小最大拟解及解的整体存在性定理 ,改进、推广了 [1 - 2 ]中的相应结果     

Evans functions and bifurcations of standing wave fronts of a nonlinear system of reaction diffusion equations

Consider the following nonlinear system of reaction diffusion equations arising from mathematical neuroscience $\frac{\partial u}{\partial t}=\frac{\partial^2u}{\partial x^2}+\alpha[\beta H(u-\theta)-u]-w,~ \frac{\partial w}{\partial t}=\varepsilon(u-\gamma w).$ Also consider the nonlinear scalar reaction diffusion equation $\frac{\partial u}{\partial t}=\frac{\partial^2u}{\partial x^2}+\alpha[\beta H(u-\theta)-u].$ In these model equations, $\alpha>0$, $\beta>0$, $\gamma>0$, $\varepsilon>0$ and $\theta>0$ are positive constants, such that $0<2\theta<\beta$. In the model equations, $u=u(x,t)$ represents the membrane potential of a neuron at position $x$ and time $t$, $w=w(x,t)$ represents the leaking current, a slow process that controls the excitation.\\indent The main purpose of this paper is to couple together linearized stability criterion (the equivalence of the nonlinear stability, the linear stability and the spectral stability of the standing wave fronts) and Evans functions (complex analytic functions) to establish the existence, stability, instability and bifurcations of standing wave fronts of the nonlinear system of reaction diffusion equations and to establish the existence and stability of the standing wave fronts of the nonlinear scalar reaction diffusion equation.     

非线性耦合Klein-Gordon方程组和耦合Schr\"{o}dinger-Boussinesq方程组的双行波解

本文提出了一种全新复合$(\frac{G''}{G})$展开方法,运用这种新方法并借助符号计算软件构造了非线性耦合Klein-Gordon方程组和耦合Schr\"{o}dinger-Boussinesq方程组的多种双行波解,包括双双曲正切函数解,双正切函数解,双有理函数解以及它们的混合解. 复合$(\frac{G''}{G})$展开方法不但直接有效地求出了两类非线性偏微分方程的双行波解,而且扩大了解的范围.这种新方法对于研究非线性偏微分方程具有广泛的应用意义.     

四阶具强阻尼非线性波动方程解的整体存在性与不存在性

研究四阶具强阻尼非线性波动方程utt-△u+△~2u-α△ut=f(u)的初边值问题.利用位势井方法,得到了问题整体解存在性与不存在性的最佳条件(门槛结果).     

Hausdorff measure of critical sets of solutions to magnetic schr?dinger equations

In this paper, we study the critical set of a complex-valued solution to a Schr?dinger equation involving the magnetic field and with a nonlinear term, where the critical set is ${\{x\in\Omega:~\psi(x)=0, ~\nabla\psi(x)=0\}}$ . We consider this equation in a bounded domain of ${\mathbb{R}^3}$ with the boundary condition: ${\nabla _{\mathbf{A}}\psi\cdot \nu=0}$ , and we establish a global 1-dimensional Hausdorff measure estimate for the critical sets. From the proof of global estimates, we find that our methods work as well for more general equations with a magnetic potential.     

The Cauchy Problem for the Average Field Equation Describing a Model of a Magnetic Solid

We study a model of a magnetic solid treated as a system of particles with mechanical moment $\vec s,\vec s \in S^2$ , and magnetic moment $\vec \mu ,\vec \mu = \vec s$ , interacting with one another via the magnetic field, which determines variations in the mechanical moment of each particle. We study the system of integro-differential equations describing the evolution of the one-particle distribution function for this system of particles. We prove existence and uniqueness theorems for the generalized and the classical solution of the Cauchy problem for this system of equations. We also prove that the generalized solution continuously depends on the initial conditions.     

一类带有铁磁材料参数的非线性涡流问题的A-φ有限元法

针对带有铁磁材料的非线性涡流问题,其非线性性通常体现在磁场强度和磁感应强度的关系上.本文提出了一种全离散的有限元A-φ格式,分别在时间和空间上采用向后欧拉公式以及节点有限元进行离散.首先,在合适的函数空间里给出时间上的半离散格式,通过考察其弱形式建立相应的适定性理论,并证明近似解收敛于弱解.其次,给出全离散格式并讨论其误差估计.最后,给出两个数值算例以验证理论结果.     

On an iterative method for a class of 2 point \& 3 point nonlinear SBVPs

In this article, we propose a novel modification to Quasi-Newton method, which is now a days popularly known as variation iteration method (VIM) and use it to solve the following class of nonlinear singular differential equations which arises in chemistry $-y''(x)-\frac{\alpha}{x}y''(x)=f(x,y),~x\in(0,1),$ where $\alpha\geq1$, subject to certain two point and three point boundary conditions. We compute the relaxation parameter as a function of Bessel and the modified Bessel functions. Since rate of convergence of solutions to the iterative scheme depends on the relaxation parameter, thus we can have faster convergence. We validate our results for two point and three point boundary conditions. We allow $\partial f/\partial y$ to take both positive and negative values.     

Biomagnetic viscoelastic fluid flow over a stretching sheet

Of concern in this paper is an investigation of biomagnetic flow of a non-Newtonian viscoelastic fluid over a stretching sheet under the influence of an applied magnetic field generated owing to the presence of a magnetic dipole. The viscoelasticity of the fluid is characterised by Walter’s B fluid model. The applied magnetic field has been considered to be sufficiently strong to saturate the ferrofluid. The magnetization of the fluid is considered to vary linearly with temperature as well as the magnetic field intensity. The theoretical treatment of the physical problem consists of reducing it to solving a system of non-linear coupled differential equations that involve six parameters, which are solved by developing a finite difference technique. The velocity profile, the skin-friction, the wall pressure and the rate of heat transfer at the sheet are computed for a specific situation. The study shows that the fluid velocity increases as the rate of heat transfer decreases, while the local skin-friction and the wall pressure increase as the magnetic field strength is increased. It is also revealed that fluid viscoelasticity has an enhancing effect on the local skin-friction. The study will have an important bearing on magnetic drug targeting and separation of red cells as well as on the control of blood flow during surgery.     

Inverse problems in magneto-electroscanning (in encephalography, for magnetic microscopes, etc.)

Contrary to the prevailing opinion about the incorrectness of the inverse MEEG-problem, we prove its unique solvability in the framework of the system of Maxwell''s equations [3]. The solution of this problem is the distribution of ${\bf y} \mapsto {\bf q}({\bf y})$ current dipoles of brain neurons that occupies the region $Y \subset \mathbb{R}^3 $. It is uniquely determined by the non-invasive measurements of the electric and magnetic fields induced by the current dipoles of neurons on the patient"s head. The solution can be represented in the form ${\bf q}={\bf q}_0+{\bf p}_0\delta\Big|_{\partial Y}$, where ${\bf q}_0$ is the usual function defined in $Y,$ and ${\bf p}_0\delta\Big|_{\partial Y} $ is a $\delta$-function on the boundary of the domain $Y$ with a certain density ${\bf p}_0$. It is essential that ${\bf p}_0$ and ${\bf q}_0$ are interrelated. This ensures the correctness of the inverse MEEG-problem. However, the components of the required 3-dimensional distribution $ {\bf q} $ must turn out to be linearly dependent if only the magnetic field ${\bf B}$ is taken into account. This question is considered in detail in a flat model of the situation.     

奇异非线性$p-$调和方程的一类正整体解

设p>1,β≥0是常数, n是自然数, 是一个连续函数.本文研究形如的奇异非线性p-调和方程的正整体解,给出了该类方程具有无穷多个其渐近阶刚好为|x|(2n-2)(当|x|→∞时)的径向对称的正整体解的若干充分条件.     

Optimal error estimate of the penalty finite element method for the time-dependent Navier-Stokes equations

A fully discrete penalty finite element method is presented for the two-dimensional time-dependent Navier-Stokes equations. The time discretization of the penalty Navier-Stokes equations is based on the backward Euler scheme; the spatial discretization of the time discretized penalty Navier-Stokes equations is based on a finite element space pair which satisfies some approximate assumption. An optimal error estimate of the numerical velocity and pressure is provided for the fully discrete penalty finite element method when the parameters and are sufficiently small.

     

求解非线性方程组的一种新的全局收敛的Levenberg-Marquardt算法

本文提出了求解非线性方程组的一种新的全局收敛的Levenberg-Marquardt算法,即μk=ακ(θ||F_k|| (1-θ)||J_k~TF_k||),θ∈[0,1],其中ακ利用信赖域技巧来修正.在不必假设雅可比矩阵非奇异的局部误差界条件下,证明了该算法是全局收敛和局部二次收敛的.数值试验表明该算法能有效地求解奇异非线性方程组问题.     

Fully discrete linear approximation scheme for electric field diffusion in type-II superconductors

Superconductors are attracting physicists thanks to their ability to conduct electric current with virtually zero resistance. Their nonlinear behaviour opens, on the other hand, challenging problems for mathematicians. Our model of the diffusion of electric field in superconductors is based on three pillars: the eddy-current version of Maxwell’s equations, power law model of type-II superconductivity and linear dependence of magnetic induction on magnetic field. This leads to a time-dependent nonlinear degenerate partial differential equation. We propose a linear fully discrete approximation scheme to solve it. We have proven the convergence of the method and derived the error estimates describing the dependence of the error on the discretization parameters. These theoretical results were successfully confronted with numerical experiments.     

${\mbox{\boldmath $R$}}^N$上奇异非线性多调和方程的正整体解

本文研究形如△((△nu)(p-1) )=f(|x|,u,|(?)u|)u-β,x∈RN的奇异非线性多调和方程在RN上的正整体解,此处P>1,β≥0是常数,n是自然数,f:R × R ×R →R 是一个连续函数, ξδ*:=sign(ξ)·|ξ|δ,,ξ∈R,δ>0,给出了该类方程具有无穷多个其渐进阶刚好为|x|2n的正整体解的充分条件与必要条件.这些结论可以推广到更一般的方程.     

A nonconforming Morley finite element method for the fully nonlinear Monge-Ampère equation

In this paper, we study finite element approximations of the viscosity solution of the fully nonlinear Monge-Ampère equation, det(D ² u) = f (> 0) using the well-known nonconforming Morley element. Our approach is based on the vanishing moment method, which was recently proposed as a constructive way to approximate fully nonlinear second order equations by the author and Feng (J Sci Comput 38(1):74–98, 2009). The vanishing moment method approximates the Monge-Ampère equation by the fourth order quasilinear equation -eD²u^e + det(D²u^e) = f{-\epsilon\Delta^2u^\epsilon + {\rm det}(D^2u^\epsilon) = f} with appropriate boundary conditions. We develop a finite element scheme using the n-dimensional Morley element introduced in Wang and Xu (Numer Math 103:155–169, 2006) to approximate the regularized fourth order problem in two and three dimensions, and then derive optimal order error estimates.