Viscous detonation in H2‒O2‒Ar using intrinsic low-dimensional manifolds and wavelet adaptive multilevel representation

A standard ignition delay problem for a mixture of hydrogen-oxygen-argon in a shock tube is extended to the viscous regime and solved using the method of intrinsic low-dimensional manifolds (ILDM) coupled with a wavelet adaptive multilevel representation (WAMR) spatial discretization technique. An operator-splitting method is used to describe the reactions as a system of ordinary differential equations at each spatial point. The ILDM method is used to eliminate the stiffness associated with the chemistry by decoupling processes which evolve on fast and slow time scales. The fast time scale processes are systematically equilibrated, thereby reducing the dimension of the phase space required to describe the reactive system. The WAMR technique captures the detailed spatial structures automatically with a small number of basis functions thereby further reducing the number of variables required to describe the system. A maximum of only 300 collocation points and 15 scale levels yields results with striking resolution of fine-scale viscous and induction zones. Additionally, the resolution of physical diffusion processes minimizes the effects of potentially reaction-inducing artificial entropy layers associated with numerical diffusion.     

The MAST FV/FE scheme for the simulation of two-dimensional thermohaline processes in variable-density saturated porous media

A novel methodology for the simulation of 2D thermohaline double diffusive processes, driven by heterogeneous temperature and concentration fields in variable-density saturated porous media, is presented. The stream function is used to describe the flow field and it is defined in terms of mass flux. The partial differential equations governing system is given by the mass conservation equation of the fluid phase written in terms of the mass-based stream function, as well as by the advection–diffusion transport equations of the contaminant concentration and of the heat. The unknown variables are the stream function, the contaminant concentration and the temperature. The governing equations system is solved using a fractional time step procedure, splitting the convective components from the diffusive ones. In the case of existing scalar potential of the flow field, the convective components are solved using a finite volume marching in space and time (MAST) procedure; this solves a sequence of small systems of ordinary differential equations, one for each computational cell, according to the decreasing value of the scalar potential. In the case of variable-density groundwater transport problem, where a scalar potential of the flow field does not exist, a second MAST procedure has to be applied to solve again the ODEs according to the increasing value of a new function, called approximated potential. The diffusive components are solved using a standard Galerkin finite element method. The numerical scheme is validated using literature tests.     

Passive ranging and a three-dimensional imaging system through wavefront coding

A single-image passive ranging and three-dimensional(3 D) imaging system with chiral phase encoding was proposed in 2011 [Opt. Lett. 36, 115(2011)]. A new theoretical analysis of the system in space domain is presented in this paper. We deduce the analytic relationships between the object distance and the point spread function, and between the object distance and the encoded image, respectively. Both the point spread function and the processed spectrum of the encoded image have two spots, which will rotate with the variation of the object distance. Then the depth map is extracted from the encoded image and it can be used to set up 3 D images. The theoretical analysis is verified by a wavefront coding system with a chiral phase which is generated by a phase-only liquid-crystal spatial light modulator. The phase generated by the liquid-crystal spatial light modulator is more flexible than the fixed phase mask and can be adjusted in real time. It is especially suitable for observing the object with a large depth of field.     

双曲正割光束在弱非局域非线性体材料中的传输特性研究

利用变分法研究了(2+1)维圆对称双曲正割光束在弱非局域非线性介质中的传输,得到了描述光束束宽、相位、波前曲率、振幅演化的一组微分方程,并得到了光束做孤子传输的临界功率;通过稳定性分析给出了弱非局域情形非局域效应对光束传输的稳定作用的定量描述,从而自洽地阐述了由不稳定的(2+1)维克尔孤子到稳定的(2+1)维弱非局域孤子的过渡情形. 数值模拟的结果验证了变分计算结果的正确性,并说明圆对称的双曲正割函数是(2+1)维弱非局域空间孤子的很好的近似. 关键词： 双曲正割光束 弱非局域非线性介质 空间光孤子     

Selkov模型不同扩散和流速下非Turing不稳定化学反应机制

建立了Selkov模型中间反应物具有不同扩散和不同流速条件下的反应 扩散 流动方程 ,理论分析了非Turing不稳定形成的条件 ,求得其参数区间 ,对Andresen的结论作了拓展 .研究还发现 ,在振荡Hopf区域之外 ,静止波动 (空间周期结构FDS)仍然可以存在 .因而 ,此结构存在的参数空间大于Andresen的结果 .同时 ,还将此种不稳定参数区间与Turing不稳定和差速流动引起不稳定 (DIFI)的结果进行了比较 ,结果发现静态FDS值总是处于DIFI临界曲线相应的最小值之上 ,这表明动力学机制是由DIFI不稳定造成的 ,DIFI不稳定区是产生静止波FDS不稳定结构的必要条件     

Numerical evaluation of linearized image reconstruction based on finite element method for biomedical photoacoustic imaging

An image reconstruction algorithm for biomedical photoacoustic imaging is discussed. The algorithm solves the inverse problem of the photoacoustic phenomenon in biological media and images the distribution of large optical absorption coefficients, which can indicate diseased tissues such as cancers with angiogenesis and the tissues labeled by exogenous photon absorbers. The linearized forward problem, which relates the absorption coefficients to the detected photoacoustic signals, is formulated by using photon diffusion and photoacoustic wave equations. Both partial differential equations are solved by a finite element method. The inverse problem is solved by truncated singular value decomposition, which reduces the effects of the measurement noise and the errors between forward modeling and actual measurement systems. The spatial resolution and the robustness to various factors affecting the image reconstruction are evaluated by numerical experiments with 2D geometry.     

Structure and dynamic response of radiative diffusion flames

Nitrogen-diluted hydrogen burning in air is modeled numerically using a constant density and one-step reaction model in a plane two-dimensional counterflow configuration. An optically thin assumption is used to investigate the effects of radiation on the dynamics, structure, and extinction of diffusion flames. While there exist dual steady-state extinction limits for the 1D radiative flame response, it is found that as the 1D radiative extinction point is approached the 1D low-stretch diffusion flame exhibits oscillatory response, even with sub-unity Lewis number fuel. These radiation-induced limit cycle oscillations are found to have increasing amplitude and decreasing frequency as the stretch rate is reduced. Flame oscillation eventually leads to permanent extinction at the stretch rate which is larger than the steady-state radiative extinction value. Along the 1D radiative response curve, the transition from 1D flame to 2D structure and the differences in the resulting 2D flame patterns are also examined using a variety of initial profiles, with special emphasis on the comparison of using the initial profiles with and without a flame edge. Similar to the previous studies on the high-stretch adiabatic edge flames using the same configuration, the high-stretch radiative flames are found to resist 1D blow-off quenching through various 2D structures, including propagating front and steady cellular flames for initial profiles with and without flame edges. For all initial profiles studied, the low-stretch radiative flames are also found to exhibit different 2D flame phenomena near the 1D radiative extinction limit, such as transient cellular structures, steady cellular structures, and pulsating ignition fronts. Although the results demonstrate the presence of low-stretch and high-stretch 2D bifurcation branches close to the corresponding 1D extinction limits irrespective of the initial profile used, particular 2D flame structures in certain stretch rate range are initial profile dependent. The existence of two-dimensional flame structures beyond the 1D steady-state radiative extinction limit suggests that the flammable range is expanded as compared to that predicted by the 1D model. Hence, multi-dimensional flame patterns need to be accounted for when determining the flammability limits for a given system.     

Nonlinear diffusion under a time dependent external force: q-maximum entropy solutions

Nonlinear diffusion equations provide useful models for a number of interesting phenomena, such as diffusion processes in porous media. We study here a family of nonlinear Fokker-Planck equations endowed both with a power-law nonlinear diffusion term and a drift term with a time dependent force linear in the spatial variable. We show that these partial differential equations exhibit exact time dependent particular solutions of the Tsallis maximum entropy (q-MaxEnt) form. These results constitute generalizations of previous ones recently discussed in the literature [C. Tsallis, D.J. Bukman, Phys. Rev. E 54, R2197 (1996)], concerning q-MaxEnt solutions to nonlinear Fokker-Planck equations with linear, time independent drift forces. We also show that the present formalism can be used to generate approximate q-MaxEnt solutions for nonlinear Fokker-Planck equations with time independent drift forces characterized by a general spatial dependence. Received 25 April 2001 and Received in final form 6 June 2001     

Quantization of the O（N） Nonlinear Sigma Model

The Hamilton-Jacobi method of quantizing singular systems is discussed.The equations of motion are obtained as total differential equations in many variables.It is shown that if the system is integrable,one can obtain the canonical phase space coordinates and set of canonical Hamilton-Jacobi partial differential equations without any need to introduce unphysical auxiliary fields.As an example we quantize the O(2) nonlinear sigma model using two different approaches:the functional Schrodinger method to obtain the wave functionals for the ground and the exited state and then we quantize the same model using the canonical path integral quantization as an integration over the canonical phase-space coordinates.     

Concentration wave description of an early stage kinetics of clustering and ordering in ordered phase

Kinetic theory of early stages of coherent clustering and continuous ordering in ordered phase is developed using the concept of concentration waves on the assumption of mean field energy. Linear differential equations of atom diffusion in ordered phase are solved by diagonalizing the coefficient matrix of the equations. It is shown that at the early stages of the reaction the concentration (fluctuation) waves evolve exponentially with respect to time on the sublattices of the ordered phase.The theory is applied to the Fe-10.9 at. %Si alloy which exhibits a sequence of transitions on aging at 775K. It is proved that the first transition from the disordered to the B2 phase occurs purely due to the kinetic reason that the amplification factors of the fluctuation waves have a positive maximum value only at the superlattice point of the B2 structure in the first Brillouin zone of the disordered phase. The second transition, which is from the metastable B2 phase to the stable phase, i.e. the spinodally decomposed D03 phase, is also explained consistently by the present theory that the amplifications calculated in the B2 phase exhibit two positive maxima in the first Brillouin zone. The maximum at point $\text{(q) = (}\text{1}\text{2}\text{1}\text{2}\text{1}\text{2}\text{)}$ represents the amplification factor of the D03 ordering-wave in the B2 phase, while the other close to the Brillouin zone center represents the amplification factor of the fluctuation wave of the spinodal decomposition.     

Analysis of a spray diffusion flame and its extinction in a shear layer

A new semi-analytical solution for a laminar spray diffusion flame in the shear layer between fuel and oxidant streams is developed. The Stokes number is identified as a small spray droplet-related parameter to be used in a perturbation analysis of the liquid phase governing equations. Appropriate specification of an additional parameter ensures that similarity is achieved so that the concentration field of the liquid in the spray can be readily found. The coupled gas-phase equations are treated using the usual inverse of the large Zeldovitch number for the asymptotic analysis. Numerical results demonstrate how the distribution of the liquid phase in the developing shear layer between two unidirectional gas streams flowing over one another with (the possibility of) dissimilar velocities in their respective free-streams influences the flame shape, location, fuel vapour and temperature fields. An extinction analysis enables a parametric mapping of conditions for extinguishment of the spray diffusion flames to be drawn.     

On the down-sampling of diffusion MRI data along the angular dimension

BackgroundIt has been established that the diffusion gradient directions in diffusion MRI should be uniformly distributed in 3D spherical space, so that orientation-dependent diffusion properties (e.g., fractional anisotropy or FA) can be properly quantified. Sometimes the acquired data need to be down-sampled along the angular dimension before computing diffusion properties (e.g., to exclude data points corrupted by motion artifact; to harmonize data obtained with different protocols). It is important to quantitatively assess the impact of data down-sampling on measurement of diffusion properties.Materials and methodsHere we report 1) a numerical procedure for down-sampling diffusion MRI (e.g., for data harmonization), and 2) a spatial uniformity index of diffusion directions, aiming to predict the quality of the chosen down-sampling schemes (e.g., from data harmonization; or rejection of motion corrupted data points). We quantitatively evaluated human diffusion MRI data, which were down-sampled from 64 or 60 diffusion gradient directions to 30 directions, in terms of their 1) FA value accuracy (using fully-sampled data as the ground truth), 2) FA fitting residuals, and 3) spatial uniformity indices.ResultsOur experimental data show that the proposed spatial uniformity index is correlated with errors in FA obtained from down-sampled diffusion MRI data. The FA fitting residuals that are typically used to assess diffusion MRI quality are not correlated with either FA errors or spatial uniformity index.ConclusionsThese results suggest that the spatial uniformity index could be more valuable in assessing quality of down-sampled diffusion MRI data, as compared with FA fitting residual measures. We expect that our implemented software procedure should prove valuable for 1) guiding data harmonization for multi-site diffusion MRI studies, and 2) assessing the impact of rejecting motion corrupted data points on the accuracy of diffusion measures.     

Numerical simulation for two-dimensional Riesz space fractional diffusion equations with a nonlinear reaction term

Fractional differential equations have attracted considerable interest because of their ability to model anomalous transport phenomena. Space fractional diffusion equations with a nonlinear reaction term have been presented and used to model many problems of practical interest. In this paper, a two-dimensional Riesz space fractional diffusion equation with a nonlinear reaction term (2D-RSFDE-NRT) is considered. A novel alternating direction implicit method for the 2D-RSFDE-NRT with homogeneous Dirichlet boundary conditions is proposed. The stability and convergence of the alternating direction implicit method are discussed. These numerical techniques are used for simulating a two-dimensional Riesz space fractional Fitzhugh-Nagumo model. Finally, a numerical example of a two-dimensional Riesz space fractional diffusion equation with an exact solution is given. The numerical results demonstrate the effectiveness of the methods. These methods and techniques can be extended in a straightforward method to three spatial dimensions, which will be the topic of our future research.     

Semi-Lagrangian multistep exponential integrators for index 2 differential–algebraic systems

Implicit-explicit (IMEX) multistep methods are very useful for the time discretization of convection diffusion PDE problems such as the Burgers equations and the incompressible Navier–Stokes equations. In the latter as well as in PDE models of plasma physics and of electromechanical systems, semi-discretization in space gives rise to differential–algebraic (DAE) system of equations often of index higher than 1. In this paper we propose a new class of exponential integrators for index 2 DAEs arising from the semi-discretization of PDEs with a dominating and typically nonlinear convection term. This class of problems includes the incompressible Navier–Stokes equations. The integration methods are based on the backward differentiation formulae (BDF) and they can be applied without modifications in the semi-Lagrangian integration of convection diffusion problems. The approach gives improved performance at low viscosity regimes.     

Diffusion of Au atoms and (5×2) phase formation on the Si(111) (7×7) surface

In this article we investigate the complex 1D mesoscopic model of adatom diffusion and the evolution of an ordered phase on the substrate surface. The analysis of the theoretical model is compared with the experimental results of the spreading of Au adatoms on Si(111)-(7×7) surface. The steady state solutions and their stability conditions are determined within the concept of the traveling-wave solution. It is shown that the formation of the ordered phase (5×2) and the difference in the diffusion of Au on (7×7) and on (5×2) structure results in a sharp edge of diffusion front which corresponds to the coverage of a saturated (5×2) phase. This edge moves linearly in time and α can be determined by experiment. The system of model equations enables the damped waves solution or temporary evolution of two steps.     

Higher-order dynamics in lattice-based models using the chapman-enskog method

In this paper, we investigate the existence of higher-order dynamics in lattice-based models. We have identified two conditions that determine whether a model would allow some Burnett-like equations when the Chapman-Enskog expansion is used. These two conditions are the number of the conserved quantity as well as the space and time discretization. We shall demonstrate these conditions by discussing (1) pure diffusion equation and (2) hydrodynamic equations. While the fact that diffusion equation allows the higher-order dynamics can be shown easily, we will illustrate that care must be taken when deriving Burnett-like equations for lattice-based hydrodynamics models using the Chapman-Enskog method.     

Polarization changes and beam profile deformation of light during acousto-optic interaction in isotropic media

A spatial Fourier transform approach is proposed to investigate the effects of polarization changes and beam profile deformation of light during acousto-optic (AO) interaction in isotropic media. The behaviour of the total scattered optical fields inside the AO cell can be properly described by a vector wave equation of which the permittivity is perturbed by an acoustic wave propagating inside the medium. In the Bragg regime, using a spatial Fourier transform approach, two coupled differential equations can be derived from the wave equation to depict AO interaction in the spatial frequency domain. Analytic solutions, which comprise the effects of changing polarization, beam deformation and propagating diffraction, can be found from the coupled equations. Detailed numerical simulations, including Fourier transforming the incident light profile to calculate the spectra of the scattered light beams and, hence, their profiles in space using the inverse transform, are presented.     

The Canonical Path Integral Quantization of CP1 Model

The Hamilton-Jacobi method of quantizing singular systems is discussed. The equations of motion are obtained as total differential equations in many variables. It is shown that if the system is integrable, then one can obtain the canonical phase space coordinates and the set of the canonical Hamilton-Jacobi partial differential equations without any need to introduce unphysical auxiliary fields. As an example we quantize the CP¹ model using the canonical path integral quantization formalism to obtain the path integral as an integration over the canonical phase-space coordinates.     

New cases of quasiperiodic motions in reversible systems

A theorem on the existence of invariant D-dimensional tori in reversible mappings near surfaces foliated into invariant tori of dimension d is announced, where d相似文献