Improvement of eigenvalues of Sturm–Liouville problem with -periodic boundary conditions

The computation of eigenvalues of Sturm–Liouville problem with t-periodic boundary conditions is considered. Using the Richardson extrapolation based on the finite difference, the accuracy of the eigenvalues are improved. Numerical results demonstrate the usefulness of the correction.     

On the Global Stability of Solutions of Moment Systems in Nonequilibrium Thermodynamics

In this paper, we study the linearization of the Cauchy problem and the mixed problem for the system of Grad--Hermite moments in nonequilibrium thermodynamics in the neighborhood of the equilibrium state. Stability conditions for solutions of the Cauchy problem are proved as a generalization of the classical Hermite--Biller theorem on stable polynomials. For the mixed problem, we prove an analog of the Vishik--Lyusternik theorem on small singular perturbations of general elliptic problems. The last observation allows us to introduce the Shapiro--Lopatinskii condition, which implies the well-posedness of the mixed problem.     

On the multidimensional distribution of the subset sum generator of pseudorandom numbers

We show that for a random choice of the parameters, the subset sum pseudorandom number generator produces a sequence of uniformly and independently distributed pseudorandom numbers. The result can be useful for both cryptographic and quasi-Monte Carlo applications and relies on bounds of exponential sums.

     

Monte Carlo ray-tracing simulation for radiative heat transfer in turbulent fluctuating media under the optically thin fluctuation approximation

The Monte Carlo ray-tracing method (MCRT) based on the concept of radiation distribution factor is extended to solve radiative heat transfer problem in turbulent fluctuating media under the optically thin fluctuation approximation. A one-dimensional non-scattering turbulent fluctuating media is considered, in which the mean temperature and absorption coefficient distribution are assumed and the shape of probability density function is given. The distribution of the time-averaged volume radiation heat source is solved by MCRT and direct integration method. It is shown that the results of MCRT based on the concept of radiation distribution factor agree with these of integration solution very well, but results of MCRT based on the concept of radiative transfer coefficient do not agree with these of integration solution. The solution of time-averaged radiative transfer equation by the concept of radiative transfer coefficient should be treated with caution.     

Lorentz Invariant Decompositions of the State Vector Spaces and the Basis Problem

We consider a representation of the state reduction which depends neither on its reality nor on the details of when and how it emerges. Then by means of the representation we find necessary conditions, even if not the sufficient ones, for a decomposition of the state vector space to be a solution to the basis problem. The conditions are that the decomposition should be Lorentz invariant and orthogonal and that the associated projections should be continuous. They are shown to be able to determine a decomposition in each of a few examples considered if the other circumstances are taken into account together.     

Capacity Constrained Transit Assignment with Common Lines

This paper proposes the use of absorbing Markov chains to solve the capacity constrained transit network loading problem taking common lines into account. The approach handles congested transit networks, where some passengers will not be able to board because of the absence of sufficient space. The model also handles the common lines problem, where choice of route depends on frequency of arrivals. The mathematical formulation of the problem is presented together with a numerical example. This revised version was published online in July 2006 with corrections to the Cover Date.     

Comment on “A Hierarchy of Integrable Nonlinear Lattice Equations and New Integrable Symplectic Map“

Comment on a recent paper on Commun. Theor. Phys. (Beijing, China) 38 (2002) pp. 523-528.     

THE COUPLING OF NATURAL BOUNDARY ELEMENT AND FINITE ELEMENT METHOD FOR 2D HYPERBOLIC EQUATIONS

In this paper, we investigate the coupling of natural boundary element and finite element methods of exterior initial boundary value problems for hyperbolic equations. The governing equation is first discretized in time, leading to a time-step scheme, where an exterior elliptic problem has to be solved in each time step. Second, a circular artificial boundary FR consisting of a circle of radius R is introduced, the original problem in an unbounded domain is transformed into the nonlocal boundary value problem in abounded subdomain. And the natural integral equation and the Poisson integral formula are obtained in the infinite domainΩ2 outside circle of radius R. The coupled variational formulation is given. Only the function itself, not its normal derivative at artificial boundary ΓR, appears in the variational equation, so that the unknown numbers are reducedand the boundary element stiffness matrix has a few different elements. Such a coupled method is superior to the one based on direct boundary element method. This paper discusses finite element discretization for variational problem and its corresponding numerical technique, and the convergence for the numerical solutions. Finally, the numerical example is presented to illustrate feasibility and efficiency of this method.     

风冷热泵冷热水机组动态过程仿真

本文给出了风冷热泵冷热水机组动态过程模型。系统的高、低压段被分成制冷剂侧、管壁、以及水或空气侧三个控制容积。针对每个控制容积建立了质量和能量平衡的微分方程。由于压缩机、热力膨胀阀阀体、以及四通换向阀具有较小的热惯性,因而采用稳态模型描述压缩过程、膨胀过程、以及内泄露过程。膨胀阀感温包中的制冷剂温度对于蒸发器出口温度的延迟用一阶惯性环节来描述。通过“预测-校正”和“自适应步长”方法实现系统仿真。仿真结果与测试数据吻合良好。     

Analysis of Two-Dimensional Optical Fields Using an Additional Shift

A solution to the phase problem in optics is considered within the context of registration and analysis of two-dimensional stationary optical fields transformed by the object under study or fields forming an image. To obtain information on amplitude and phase distributions of the light field analyzed, a method of registration of two intensity distributions is used. The first distribution corresponds directly to the amplitude distribution. The other is formed for the sum of the initial field and the field shifted along a certain direction. The intensity distributions obtained allow one to calculate the two-dimensional structure of the field under study. It is noteworthy that the method requires no iteration procedures in solving the problem. This leads to speeding up of the processing and analysis of the information. Two variants of optical schemes for the analysis of light fields are considered. The first one corresponds to registration of the image of the analyzed plane and the second to registration of the spectrum of the spatial frequencies.