Confined diffusion

Summary An equation for the unidimensional confined diffusion is proposed. The equation coincides with the well-known homogeneous equation except the presence of a source term. This term which has the form of a dipole distribution is located on a moving front which sharply separates two distinct regions. In the first region (from the boundary up to the front) the confined solution coincides with a suitable solution of the homogeneous equation; in the second region (besides the front) it vanishes. The source term, moreover, switches off the diffusing flux at the front. The sharp confinement allows to relax the original boundary conditions of the homogeneous equation. Precisely, to the function depending on the time at the boundary, another arbitrary function depending on the space at the initial time is added. This new function (provided not vanishing) allows to obtain in general an acceptable evolution of the front and does not prevent the validity of the conservation law: flux at the boundary is equal to the time variation of the diffusing quantity contained between the boundary and the front. By a suitable choice of this new function, so that it results to be connected to the other boundary condition (that depending on time) it is possible to arrive at an evolution of the front such as: , where λ,K, corresponding, respectively, to a dimensionless parameter and diffusivity, depend on the medium. Under such simplifying assumption, it is possible to obtain an analytical expression for the confined solution. This solution, evaluated in a point of the space, arrives asymptotically at the same value reached by the solution of the homogeneous equation.     

Zero Viscosity Limit for Analytic Solutions of the Navier-Stokes Equation on a Half-Space.¶ II. Construction of the Navier-Stokes Solution

This is the second of two papers on the zero-viscosity limit for the incompressible Navier-Stokes equations in a half-space in either 2D or 3D. Under the assumption of analytic initial data, we construct solutions of Navier-Stokes for a short time which is independent of the viscosity. The Navier-Stokes solution is constructed through a composite asymptotic expansion involving the solutions of the Euler and Prandtl equations, which were constructed in the first paper, plus an error term. This shows that the Navier-Stokes solution goes to an Euler solution outside a boundary layer and to a solution of the Prandtl equations within the boundary layer. The error term is written as a sum of first order Euler and Prandtl corrections plus a further error term. The equation for the error term is weakly nonlinear; its linear part is the time dependent Stokes equation. This error equation is solved by inversion of the Stokes equation, through expressing the solution as a regular (Euler-like) part plus a boundary layer (Prandtl-like) part. The main technical tool in this analysis is the Abstract Cauchy-Kowalewski Theorem. Received: 5 September 1996 / Accepted: 14 July 1997     

Conservative method for simulation of a high-order nonlinear Schrdinger equation with a trapped term

We propose a new scheme for simulation of a high-order nonlinear Schrodinger equation with a trapped term by using the mid-point rule and Fourier pseudospectral method to approximate time and space derivatives, respectively. The method is proved to be both charge- and energy-conserved. Various numerical experiments for the equation in different cases are conducted. From the numerical evidence, we see the present method provides an accurate solution and conserves the discrete charge and energy invariants to machine accuracy which are consistent with the theoretical analysis.     

Interacting New Agegraphic Viscous Dark Energy with Varying <Emphasis Type="Italic">G</Emphasis>

We consider the new agegraphic model of dark energy with a varying gravitational constant, G, in a non-flat universe. We obtain the equation of state and the deceleration parameters for both interacting and noninteracting new agegraphic dark energy. We also present the equation of motion determining the evolution behavior of the dark energy density with a time variable gravitational constant. Finally, we generalize our study to the case of viscous new agegraphic dark energy in the presence of an interaction term between both dark components.     

Energy-transport by radiative diffusion

A new derivation of the general-relativistic Fourier equation is given for radiation transport by using the principle of conservation of momentum plus some rather simple assumptions. The Fourier equation at which I arrive is not the usual one but has an additional term. For this reason it leads to a hyperbolic equation for heat conduction, thus avoiding the paradox of infinite velocity of heat propagation, which is a consequence of the usual Fourier equation, as the latter one leads to a parabolic equation for heat conduction. The new Fourier equation is compared with the one that was given by Kranys by using ad hoc assumptions.     

The evolution equation of the scalar field in the new inflationary universe

We derive an effective evolution equation for the scalar field driving inflation in the new inflationary universe. We use a perturbative calculation scheme proposed recently by Morikawa and Sasaki. The relevant initial conditions and dynamical constraints for the Coleman-Weinberg effective potential to appear in the evolution equation are discussed as well as the form of the particle production damping term. The validity of these conditions in the new inflationary universe model is discussed.     

正压大气环流中的曲面周期波和孤波

大尺度正压大气环流的波动特征对理解气候变化具有重要的意义,而非线性浅水波方程组是描述大尺度正压大气环流的原始控制方程.本文对线性方程的复变函数解,通过二次适当的移植,求得浅水波方程组的发展方程的扰动位势的实变函数解,该实变函数解析解由基流项和波动项两部分组成.其中基流由波数、波速、β效应、变形半径和时间的任意函数共同决定;波动项与β效应有关.分析表明,在大尺度正压大气环流中扰动位势存在曲面的周期波和孤波的现象,周期波与孤波相互调制而呈现不稳定性;当多个周期孤波同时出现时,则彼此独立传播;扰动位势波动项中的时间任意函数对曲面周期孤波的波幅有调制作用,可控制波的产生、发展和消失.所得结果对研究大气波动现象和气候变化具有一定的理论参考价值.     

Numerical solution of under-resolved detonations

A new fractional-step method is proposed for the numerical solution of high speed reacting flows, where the chemical time scales are often much smaller than the fluid dynamical time scales. When the problem is stiff, because of insufficient spatial/temporal resolution, a well-known spurious numerical phenomenon occurs in standard finite volume schemes: the incorrect calculation of the speed of propagation of discontinuities. The new method is first illustrated considering a one-dimensional scalar hyperbolic advection/reaction equation with stiff source term, which may be considered as a model problem to under-resolved detonations. During the reaction step, the proposed scheme replaces the cell average representation with a two-value reconstruction, which allows us to locate the discontinuity position inside the cell during the computation of the source term. This results in the correct propagation of discontinuities even in the stiff case. The method is proved to be second-order accurate for smooth solutions of scalar equations and is applied successfully to the solution of the one-dimensional reactive Euler equations for Chapman–Jouguet detonations.     

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.     

On time-dependent radiative transfer

An integral equation is developed for application to time-dependent laboratory experiments in which partial redistribution effects are important. The equation of transport with the Heasley-Kneer emission coefficient and the equation of statistical equilibrium lead to a time-dependent redistribution function containing an absorption—re-emission term which decays exponentially in time and a scattering term which is instantaneous. This integral equation does not agree with an equation written by Payne et al. [Phys. Rev. A 9, 1050 (1974)] that has been used to compare theory with experiments. The difference between the Payne equation and the equation developed here needs to be examined in detail, since it might under some circumstances be on the same order as the difference between partial and complete redistribution.     

广义测不准关系与三维BTZ黑洞熵

通过应用在量子引力中,由广义测不准关系得出的新的态密度方程,研究三维BTZ背景下黑洞的熵.当取广义测不准关系中引入的,具有Planck量级与空间维数有关的常数λ为特定值时,得到BTZ黑洞Bekenstein-Hawking 熵和修正项.由于利用新的态密度方程,在计算中不存在用brick-wall模型计算黑洞熵时出现的发散项和小质量近似.所得结论,从量子统计力学角度给出了黑洞Bekenstein-Hawking 熵的修正值,使人们对黑洞熵的修正值有更深入的认识. 关键词： 广义测不准关系 量子统计 BTZ黑洞熵     

FRW-Cosmological Model for Conharmonically Flat Space Time

This paper deals with FRW-Cosmological model of the universe for conharmonically flat space time. Einstein field equations with variable cosmological term are solved by using a law of variation for Hubble’s parameter, which is related to average scale factor. A new class of exact solution of the field equation has been obtained in which cosmological-term decreases with cosmic time. A detailed study of physical and kinematical properties of the model is also carried out.     

Torque equations for spin and orbital angular momenta of radiation fields and Faraday effect

The spin angular momentum S of light has never been linked to the Faraday rotation of light traveling in an optically active medium possessing a rotational invariance of a crystal, because there was no helicity term associated with the phase shift in the previous torque equation for S. In order to relate the change in S with time to the Faraday rotation, therefore, we derived an exact torque equation for S. As a result, a magnetic helicity term appeared in a new torque equation for S, so that one-half of the phase shift derived from the helicity term was equivalent to the Faraday rotation angle. However, the orbital angular momentum L had no relation to the Faraday rotation. It was thus clarified that the change in S with time is related to the Faraday rotation angle of light traveling in an optically active medium, owing to the appearance of the helicity term without a rotational invariance around the optical axis. It was also demonstrated theoretically that the Faraday rotation is accompanied by a torque acting on the crystal so that the total angular momentum of light and matter is conserved.     

弹粘塑性双球壳塌缩热点反应模型

 基于Kim的弹粘塑性单球壳塌缩模型,考虑PBX炸药中的粘结剂效应,假设炸药和粘结剂均为弹粘塑性材料,建立了弹粘塑性双球壳塌缩热点反应模型,给出了炸药球壳在冲击压力作用下的速度、应变、温度和化学反应速率的时空分布,以及新的热点反应速率理论表达式。把新的热点反应项与Kim的低压下慢反应项和张震宇提出的高压反应速率方程相结合,得到了新的冲击起爆三项式细观反应速率模型。把该模型加入DYNA2D中,模拟了PBX-9501炸药的一维冲击起爆过程,结果表明：该模型除了可以解释炸药颗粒度和孔隙度的影响外,还可以较好地描述粘结剂强度和含量对PBX炸药冲击起爆感度的影响。     

New longitudinal waves in electron-positron-ion quantum plasmas

A general quantum dispersion equation for electron-positron(hole)-ion quantum plasmas is derived and studied for some interesting cases. In an electron-positron-ion degenerate Fermi gas, with or without the Madelung term, a new type of zero sound waves are found. Whereas in an electron-hole-ion plasmas a new longitudinal quantum waves are revealed, which have no analogies in quantum electron-ion plasmas. The excitation of these quantum waves by a low-density monoenergetic straight electron beam is examined. Furthermore, the Korteweg-de Vries (KdV) equation for novel quantum waves is derived and the contribution of the Madelung term in the formation of the KdV solitons is discussed.     

A robust numerical model for premixed flames with high density ratios based on new pressure correction and IMEX schemes

In this study we present a model for the interaction of premixed flames with obstacles in a channel flow. Although the flow equations are solved with Direct Numerical Simulation using a low Mach number approximation, the resolution used in the computation is limited (∼1 mm) hence the inner structure of the flame and the chemical scales are not solved. The species equations are substituted with a source term in the energy equation that simulates a one-step global reaction. A level set method is applied to track the position of the flame and its zero level is used to activate the source term in the energy equation only at the flame front. An immersed boundary method reproduces the geometry of the obstacles. The main contribution of the paper is represented by the proposed numerical approach: an IMEX (implicit–explicit) Runge–Kutta scheme is used for the time integration of the energy equation and a new pressure correction algorithm is introduced for the time integration of the momentum equations. The approach presented here allows to calculate flames which produce high density ratios between burnt and unburnt regions. The model is verified by simulating first simple solutions for one- and two-dimensional flames. At last, the experiments performed by Masri and Ibrahim with square and rectangular bodies are calculated.     

Nucleus-nucleus collisions: Intrinsic motion with external correlations

A consistent treatment of the relative and intrinsic motion which goes beyond the mean-field approach allows to include the fluctuations of the time-dependent mean field for the intrinsic as well as for the relative motion. Starting with the v. Neumann equation for the total density matrix, we derive a modified equation for the intrinsic many-body density matrix. This equation is used to obtain the quantum kinetic equations for the one-body density matrix and the two-body correlator. In the time-dependent single-particle basis, the occupation numbers change in time due to a collision term originating from residual two-body interactions which account for equilibration, and due to the fluctuations of the external mean field. The relations to TDHF with collision term are discussed. Special attention is paid to the conditions for a diabatic evolution of the single-particle states and to finite size effects which play an important role to make two-body collisions operative in finite nuclei.     

Results for a fractional diffusion equation with a nonlocal term in spherical symmetry

We obtain time dependent solutions for a fractional diffusion equation containing a nonlocal term by considering the spherical symmetry and using the Green function approach. The nonlocal term incorporated in the diffusion equation may also be related to the spatial and time fractional derivative and introduces different regimes of spreading of the solution with the time evolution. In addition, a rich class of anomalous diffusion processes may be described from the results obtained here.     

双曲函数型辅助方程构造具5次强非线性项的波方程的新精确孤波解

在辅助方程法的基础上利用两种函数变换和一种双曲函数型辅助方程，通过符号计算系统Mathematica构造了在力学当中一个重要的模型，有5次强非线性项的波方程的新三角函数型和双曲函数型精确孤波解.这种方法寻找其他具5次强非线性项的非线性发展方程的新精确解方面具有普遍意义. 关键词： 双曲函数型辅助方程 函数变换 具5次强非线性项的波方程 精确孤波解     

Effect of Torsion in Dirac Equation for Coulomb Potential in Robertson–Walker Space–Time

The Dirac equation with Coulomb-like potential and self-interaction term, that originates from torsion, is studied in the Robertson–Walker space–time. It is shown that the angular dependence of the equation can be separated also in presence of nonlinear terms. Under reasonable physical assumptions, the time dependence is also separated. An extended perturbative calculation can then be applied qualitatively. The conclusion is that the perturbation of the energy levels of the system, as consequence of the self-interacting term, is not relevant on physical grounds.