Study of shape isomer yields of ~（240）Am in the framework of a dynamical-statistical model

A dynamical statistical model is used to analyze the experimental shape isomer yields data in the reaction d＋240Pu at E=20–29 MeV.The possibility of determining the nuclear dissipation is discussed.Comparison of the experimental data with the calculations leads to a value of the reduced dissipation coeffcient β=0.45×1021 s-1 for the Am isotopes.     

Autocorrelation function formulations and the turbulence dissipation rate: Application to dispersion models

The classical statistical diffusion theory and the binomial autocorrelation function are used to obtain a new formulation for the turbulence dissipation rate ε. The approach employs the Maclaurin series expansion of a logarithm function contained in the dispersion parameter formulation. The numerical coefficient of this new relation for ε is 100% larger than the numerical coefficient of the classical relation derived from the exponential autocorrelation function. A similar approach shows that the dispersion parameter obtained from the even exponential autocorrelation function does not result in a relation for ε and, therefore, is not suitable for application in dispersion models. In addition, a statistical comparison to experimental ground-level concentration data demonstrates that this newly derived relation for ε as well as other formulations for the turbulence dissipation rate are suitable for application in Lagrangian stochastic dispersion models. Therefore, the analysis shows that there is an uncertainty regarding the turbulence dissipation rate function form and the autocorrelation function form.     

A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media

We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic–acoustic media. A velocity–strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic–acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic–acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.     

Effect of instrumental particle sizing resolution on the modelling of aerosol radiative parameters

A more realistic estimation of the scattering and hemispheric backscattering coefficients, σ_sp and σ_bsp, and their respective optical cross section, C_sca and C_bk, of aerosol particles is presented on the basis of the exact resolution of the width of the size bins of the particle counter instruments when size distribution measurements are used, and, with the exact optical detector instruments ability. The scattering and hemispheric backscattering cross sections, C_sca and C_bk, of the particles are averaged over the full size bins of the particle counter instrument, while these quantities are usually estimated only on the value of the mean geometric diameter of each size bin. Six instruments, the APS, ASASP-X, DMPS, FSSP-100, ELPI, and SMPS frequently used in particle size distribution measurements are reviewed, for spherical sea-salt particles at a wavelength λ=0.55 μm. The comparison using the conventional geometric mean diameter versus the use of the full size bin leads to large amount of errors for the optical cross section with non-negligible effects on their respective optical coefficients. The maximal accuracy expected for these optical quantities depend on the particle diameter as well as on the channel width of the instruments, and are also function of the angular detector probe used to measure them.     

A fractional order HIV-TB coinfection model with nonsingular Mittag-Leffler Law

The biological models for the study of human immunodeficiency virus (HIV) and its advanced stage acquired immune deficiency syndrome (AIDS) have been widely studied in last two decades. HIV virus can be transmitted by different means including blood, semen, preseminal fluid, rectal fluid, breast milk, and many more. Therefore, initiating HIV treatment with the TB treatment development has some advantages including less HIV-related losses and an inferior risk of HIV spread also having difficulties including incidence of immune reconstitution inflammatory syndrome (IRIS) because of a large pill encumbrance. It has been analyzed that patients with HIV have more weaker immune system and are susceptible to infections, for example, tuberculosis (TB). Keeping the importance of the HIV models, we are interested to consider an analysis of HIV-TB coinfected model in the Atangana-Baleanu fractional differential form. The model is studied for the existence, uniqueness of solution, Hyers-Ulam (HU) stability and numerical simulations with assumption of specific parameters.     

An efficient nonpolynomial spline method for distributed order fractional subdiffusion equations

In this paper, we propose an efficient numerical method for a distributed order fractional subdiffusion problem using nonpolynomial spline approach. The solvability, stability, and convergence of the scheme are established rigorously, and it is shown that the spatial convergence order improves some previous work done. Simulation is then conducted to verify the accuracy of the proposed scheme as well as to compare with earlier work.     

水凝胶类软物质材料理论中的数学问题

 水凝胶是一种具有广泛应用前景的软物质材料,一直是材料学家、物理化学家们关心的热点,目前有很多实验和专利产品.其机理和模型的研究处于初期阶段,本文综述该领域的部分进展,包括形变、微结构和宏观性质等的数学模型和相关计算方法,一并列出一些亟待解决的问题.     

Mathematical and numerical analysis of a three-species predator-prey model with herd behavior and time fractional-order derivative

In this paper, we investigate with a time fractional-order derivative in a three-species predator-prey model with the presence of prey social behavior. A new approximation for predator-prey interaction in the presence of prey social behavior has been considered. For the model analysis, the study has been divided into two principal parts. First of all, we study the local stability of the equilibria and the existence of Hopf bifurcation. Then, for the numerical analysis, the Caputo fractional derivative operator is utilized to approximate the numerical solution of the model. An excellent agreement is seen between the numerical results and the theoretical predictions.     

A comparison of nine PLS1 algorithms

Nine PLS1 algorithms were evaluated, primarily in terms of their numerical stability, and secondarily their speed. There were six existing algorithms: (a) NIPALS by Wold; (b) the non‐orthogonalized scores algorithm by Martens; (c) Bidiag2 by Golub and Kahan; (d) SIMPLS by de Jong; (e) improved kernel PLS by Dayal; and (f) PLSF by Manne. Three new algorithms were created: (g) direct‐scores PLS1 based on a new recurrent formula for the calculation of basis vectors yielding scores directly from X and y; (h) Krylov PLS1 with its regression vector defined explicitly, using only the original X and y; (i) PLSPLS1 with its regression vector recursively defined from X and the regression vectors of its previous recursions. Data from IR and NIR spectrometers applied to food, agricultural, and pharmaceutical products were used to demonstrate the numerical stability. It was found that three methods (c, f, h) create regression vectors that do not well resemble the corresponding precise PLS1 regression vectors. Because of this, their loading and score vectors were also concluded to be deviating, and their models of X and the corresponding residuals could be shown to be numerically suboptimal in a least squares sense. Methods (a, b, e, g) were the most stable. Two of them (e, g) were not only numerically stable but also much faster than methods (a, b). The fast method (d) and the moderately fast method (i) showed a tendency to become unstable at high numbers of PLS factors. Copyright © 2009 John Wiley & Sons, Ltd.