


Том 57, № 6 (2017)
- Год: 2017
- Статей: 13
- URL: https://journals.rcsi.science/0965-5425/issue/view/11143
Article
Computation of zeros of the alpha exponential function
Аннотация
This paper deals with the function F(α; z) of complex variable z defined by the expansion \(F\left( {\alpha ;z} \right) = \sum\nolimits_{k = 0}^\infty {\frac{{{z^k}}}{{{{\left( {k!} \right)}^\alpha }}}} \) which is a natural generalization of the exponential function (hence the name). Primary attention is given to finding relations concerning the locations of its zeros for α ∈ (0,1). Note that the function F(α; z) arises in a number of modern problems in quantum mechanics and optics. For α = 1/2, 1/3,..., approximations of F(α; z) are constructed using combinations of degenerate hypergeometric functions 1F1(a; c; z) and their asymptotic expansions as z → ∞. These approximations to F(α; z) are used to approximate the countable set of complex zeros of this function in explicit form, and the resulting approximations are improved by applying Newton’s high-order accurate iterative method. A detailed numerical study reveals that the trajectories of the zeros under a varying parameter α ∈ (0,1] have a complex structure. For α = 1/2 and 1/3, the first 30 complex zeros of the function are calculated to high accuracy.



Synthesis of a rational filter in the presence of complete alternance
Аннотация
The construction of a rational function that is nonnegative on two intervals of which one is infinite is considered. It is assumed that the maximum deviation of the function from zero on the infinite interval takes the minimum possible value under the condition that the values of the function on the finite interval are within the given bounds. It is assumed that the rational function (fraction) has the complete alternance. In this case, the original problem is reduced to solving a system of nonlinear equations. For solving this system, a two-stage method is proposed. At the first stage, a subsystem is selected and used to find a good approximation for the complete system. At the second stage, the complete system of nonlinear equations is solved. The solution is explained in detail for the case when the order of the fraction is between one and four. Numerical results for a fraction of order ten are presented.



Parametrization of the Cauchy problem for systems of ordinary differential equations with limiting singular points
Аннотация
The method of solution continuation with respect to a parameter is used to solve an initial value problem for a system of ordinary differential equations with several limiting singular points. The solution is continued using an argument (called the best) measured along the integral curve of the problem. Additionally, a modified argument is introduced that is locally equivalent to the best one in the considered domain. Theoretical results are obtained concerning the conditioning of the Cauchy problem parametrized by the modified argument in a neighborhood of each point of its integral curve.






Parabolic equations with unknown time-dependent coefficients
Аннотация
The solvability of inverse problems of finding the coefficients of a parabolic equation together with solving this equation is studied. In these problems, certain additional conditions on the boundary are used as overdetermination conditions. Existence and uniqueness theorems for regular solutions of such problems are proven.



On power series representing solutions of the one-dimensional time-independent Schrödinger equation
Аннотация
For the equation χ″(x) = u(x)χ(x) with infinitely smooth u(x), the general solution χ(x) is found in the form of a power series. The coefficients of the series are expressed via all derivatives u(m)(y) of the function u(x) at a fixed point y. Examples of solutions for particular functions u(x) are considered.



Well-posedness analysis and numerical implementation of a linearized two-dimensional bottom sediment transport problem
Аннотация
A two-dimensional linearized model of coastal sediment transport due to the action of waves is studied. Up till now, one-dimensional sediment transport models have been used. The model under study makes allowance for complicated bottom relief, the porosity of the bottom sediment, the size and density of sediment particles, gravity, wave-generated shear stress, and other factors. For the corresponding initial–boundary value problem the uniqueness of a solution is proved, and an a priori estimate for the solution norm is obtained depending on integral estimates of the right-hand side, boundary conditions, and the norm of the initial condition. A conservative difference scheme with weights is constructed that approximates the continuous initial–boundary value problem. Sufficient conditions for the stability of the scheme, which impose constraints on its time step, are given. Numerical experiments for test problems of bottom sediment transport and bottom relief transformation are performed. The numerical results agree with actual physical experiments.



Stability of the Kolmogorov flow and its modifications
Аннотация
Recurrence formulas are obtained for the kth term of the long wavelength asymptotics in the stability problem for general two-dimensional viscous incompressible shear flows. It is shown that the eigenvalues of the linear eigenvalue problem are odd functions of the wave number, while the critical values of viscosity are even functions. If the velocity averaged over the long period is nonzero, then the loss of stability is oscillatory. If the averaged velocity is zero, then the loss of stability can be monotone or oscillatory. If the deviation of the velocity from its period-average value is an odd function of spatial variable about some x0, then the expansion coefficients of the velocity perturbations are even functions about x0 for even powers of the wave number and odd functions about for x0 odd powers of the wave number, while the expansion coefficients of the pressure perturbations have an opposite property. In this case, the eigenvalues can be found precisely. As a result, the monotone loss of stability in the Kolmogorov flow can be substantiated by a method other than those available in the literature.



Special discontinuities in nonlinearly elastic media
Аннотация
Solutions of a nonlinear hyperbolic system of equations describing weakly nonlinear quasitransverse waves in a weakly anisotropic elastic medium are studied. The influence of small-scale processes of dissipation and dispersion is investigated. The small-scale processes determine the structure of discontinuities (shocks) and a set of discontinuities with a stationary structure. Among the discontinuities with a stationary structure, there are special ones that, in addition to relations following from conservation laws, satisfy additional relations required for the existence of their structure. In the phase plane, the structure of such discontinuities is represented by an integral curve joining two saddles. Special discontinuities lead to nonunique self-similar solutions of the Riemann problem. Asymptotics of non-self-similar problems for equations with dissipation and dispersion are found numerically. These asymptotics correspond to self-similar solutions of the problems.



Lebesgue averaging method in serial computations of atmospheric radiation
Аннотация
The Lebesgue averaging method was applied to the numerical simulation of the radiative transfer equation. It was found that the method ensures good accuracy, while the amount of computations with respect to the energy variable is reduced by more than three orders of magnitude. “Fast” simplified techniques for the Lebesgue processing of photon absorption cross sections in serial computations of atmospheric radiation were examined. Attention was given to the convenience of using the techniques, including by experienced users.



Application of numerical schemes with singling out the boundary layer for the computation of turbulent flows using eddy-resolving approaches on unstructured grids
Аннотация
The use of eddy-resolving approaches to solving problems on arbitrary unstructured grids is investigated. The applications of such approaches requires the use of low dissipation numerical schemes, which can lead to numerical oscillations of the solution on unstructured grids. Numerical oscillations typically occur in domains with large gradients of velocities, in particular, in the near-wall layer. It is proposed to single out the boundary layer and use a numerical scheme with increased numerical dissipation in it. The algorithm for singling out the boundary layer uses a switching function to change the parameters of the numerical scheme. This algorithm is formulated based on the BCD scheme from the family NVD. Its validity and advantages are investigated using the zonal RANS–LES approach for solving some problems of turbulent flow of incompressible fluids.



Solution of the Wang Chang–Uhlenbeck equation for molecular hydrogen
Аннотация
Molecular hydrogen is modeled by numerically solving the Wang Chang–Uhlenbeck equation. The differential scattering cross sections of molecules are calculated using the quantum mechanical scattering theory of rigid rotors. The collision integral is computed by applying a fully conservative projection method. Numerical results for relaxation, heat conduction, and a one-dimensional shock wave are presented.



Convective heat and mass transfer in three-dimensional mixed convection flow of viscoelastic fluid in presence of chemical reaction and heat source/sink
Аннотация
Heat and mass transfer effects in the three-dimensional mixed convection flow of a viscoelastic fluid with internal heat source/sink and chemical reaction have been investigated in the present work. The flow generation is because of an exponentially stretching surface. Magnetic field normal to the direction of flow is considered. Convective conditions at the surface are also encountered. Appropriate similarity transformations are utilized to reduce the boundary layer partial differential equations into the ordinary differential equations. The homotopy analysis method is used to develop the solution expressions. Impacts of different controlling parameters such as ratio parameter, Hartman number, internal heat source/sink, chemical reaction, mixed convection, concentration buoyancy parameter and Biot numbers on the velocity, temperature and concentration profiles are analyzed. The local Nusselt and Sherwood numbers are sketched and examined.


