The Rayleigh–Stokes problem for an edge in a generalized Oldroyd-B fluid

The velocity field corresponding to the Rayleigh–Stokes problem for an edge, in an incompressible generalized Oldroyd-B fluid has been established by means of the double Fourier sine and Laplace transforms. The fractional calculus approach is used in the constitutive relationship of the fluid model. The obtained solution, written in terms of the generalized G-functions, is presented as a sum of the Newtonian solution and the corresponding non-Newtonian contribution. The solution for generalized Maxwell fluids, as well as those for ordinary Maxwell and Oldroyd-B fluids, performing the same motion, is obtained as a limiting case of the present solution. This solution can be also specialized to give the similar solution for generalized second grade fluids. However, for simplicity, a new and simpler exact solution is established for these fluids. For β → 1, this last solution reduces to a previous solution obtained by a different technique.     

The Rayleigh–Stokes problem for an edge in a generalized Oldroyd-B fluid

The velocity field corresponding to the Rayleigh–Stokes problem for an edge, in an incompressible generalized Oldroyd-B fluid has been established by means of the double Fourier sine and Laplace transforms. The fractional calculus approach is used in the constitutive relationship of the fluid model. The obtained solution, written in terms of the generalized G-functions, is presented as a sum of the Newtonian solution and the corresponding non-Newtonian contribution. The solution for generalized Maxwell fluids, as well as those for ordinary Maxwell and Oldroyd-B fluids, performing the same motion, is obtained as a limiting case of the present solution. This solution can be also specialized to give the similar solution for generalized second grade fluids. However, for simplicity, a new and simpler exact solution is established for these fluids. For β → 1, this last solution reduces to a previous solution obtained by a different technique.      

On Torsion of Prismatic Bodies

The method of decaying residual solution is applied to obtain an approximate interior solution for the torsion of slender prismatic elastic bodies under different end conditions. The approximate solution is generally accurate up to terms that are exponentially small in the length-to-cross-sectional-width ratio. For stress end conditions, the result is identical to the classical Saint-Venant torsion solution. Similar types of simple solutions, not known previously, are obtained for different types of mixed end conditions. For displacement conditions at both ends, the corresponding Saint-Venant type result requires an accurate solution of a canonical problem for a semi-infinite prismatic body that is to be obtained once and for all. The solution of the canonical problem is elementary for a circular cross section. The approximate interior solution in that case is identical to the known exact interior solution.     

Effect of the dynamic pressure on the similarity solution of cylindrical shock waves in a rarefied polyatomic gas

The similarity solution for a strong cylindrical shock wave in a rarefied polyatomic gas is analyzed on the basis of Rational Extended Thermodynamics with six independent fields; the mass density, the velocity, the pressure and the dynamic pressure. A new ODE system for the similarity solution is derived in a systematic way by using the method based on the Lie group theory proposed in the context of the spherical shock wave in a rarefied monoatomic gas in Donato and Ruggeri (J Math Anal Appl 251:395, 2000). The boundary conditions are also specified from the Rankine–Hugoniot conditions for the sub-shock. The derived similarity solution is characterized by only one dimensionless parameter \(\alpha \) related to the relaxation time for the dynamic pressure. The numerical analysis of the similarity solution is also performed. The solution agrees with the well-known Sedov–von Neumann–Taylor (SNT) solution when \(\alpha \) is small. When \(\alpha \) is larger, due to the presence of the dynamic pressure, the deviation from the SNT solution is evident; the strength of a peak near the shock front becomes smaller and the profile becomes broader.

     

ANALYSIS OF THE LIMIT CYCLE OF A GENERALISED VAN DER POL EQUATION BY A TIME TRANSFORMATION METHOD

Abstract

The properties of the limit cycle of a generalised van der Pol equation of the form ü + u = ε (1—u²ⁿ)u, where ε is small and n is any positive integer, are investigated by applying a time transformation perturbation method due to Burton. It is found that as n increases the amplitude of the limit cycle oscillation decreases and its period increases. The time transformation solution is compared with the solution derived using the method of multiple scales and with a numerical solution. It is found that, to first order in ε, the time transformation solution for the limit cycle agrees better with the numerical solution than the multiple scales solution. Both perturbation solutions give the same result for the period of the limit cycle to second order in ε. The accuracy of the time transformation solution decreases as n increases.     

Numerical analysis for a conservative difference scheme to solve the Schrödinger-Boussinesq equation

In this paper, we present a finite difference scheme for the solution of an initial-boundary value problem of the Schrödinger-Boussinesq equation. The scheme is fully implicit and conserves two invariable quantities of the system. We investigate the existence of the solution for the scheme, give computational process for the numerical solution and prove convergence of iteration method by which a nonlinear algebra system for unknown Vⁿ⁺¹ is solved. On the basis of a priori estimates for a numerical solution, the uniqueness, convergence and stability for the difference solution is discussed. Numerical experiments verify the accuracy of our method.     

Solving Linear Systems Whose Elements Are Nonlinear Functions of Intervals

The paper addresses the problem of solving linear algebraic systems the elements of which are, in the general case, nonlinear functions of a given set of independent parameters taking on their values within prescribed intervals. Three kinds of solutions are considered: (i) outer solution, (ii) interval hull solution, and (iii) inner solution. A simple direct method for computing a tight outer solution to such systems is suggested. It reduces, essentially, to inverting a real matrix and solving a system of real linear equations whose size n is the size of the original system. The interval hull solution (which is a NP-hard problem) can be easily determined if certain monotonicity conditions are fulfilled. The resulting method involves solving n+1 interval outer solution problems as well as 2n real linear systems of size n. A simple iterative method for computing an inner solution is also given. A numerical example illustrating the applicability of the methods suggested is solved.     

On the Convergence of Coderivative of SAA Solution Mapping for a Parametric Stochastic Generalized Equation

The aim of this paper is to investigate the convergence properties for Mordukhovich’s coderivative of the solution map of the sample average approximation (SAA) problem for a parametric stochastic generalized equation. It is demonstrated that, under suitable conditions, both the cosmic deviation and the ρ-deviation between the coderivative of the solution mapping to SAA problem and that of the solution mapping to the parametric stochastic generalized equation converge almost surely to zero as the sample size tends to infinity. Moreover, the exponential convergence rate of coderivatives of the solution maps to the SAA parametric generalized equations is established. The results are used to develop sufficient conditions for the consistency of the Lipschitz-like property of the solution map of SAA problem and the consistency of stationary points of the SAA estimator for a stochastic mathematical program with complementarity constraints.     

SOME RECENT PROGRESS ON STOCHASTIC HEAT EQUATIONS

This article attempts to give a short survey of recent progress on a class of elementary stochastic partial differential equations(for example, stochastic heat equations)driven by Gaussian noise of various covariance structures. The focus is on the existence and uniqueness of the classical(square integrable) solution(mild solution, weak solution). It is also concerned with the Feynman-Kac formula for the solution; Feynman-Kac formula for the moments of the solution; and their applications to the asymptotic moment bounds of the solution. It also briefly touches the exact asymptotics of the moments of the solution.     

An approximation of stochastic hyperbolic equations: case with Wiener process

In the present paper, the two‐step difference scheme for the Cauchy problem for the stochastic hyperbolic equation is presented. The convergence estimate for the solution of the difference scheme is established. In applications, the convergence estimates for the solution of difference schemes for the numerical solution of four problems for hyperbolic equations are obtained. The theoretical statements for the solution of this difference scheme are supported by the results of the numerical experiment. Copyright © 2012 John Wiley & Sons, Ltd.     

Exact Closed Form Solution for Two Coupled Inhomogeneous First Order Difference Equations with Variable Coefficients

We present an exact closed form solution for two coupled, homogeneous as well as inhomogeneous, first order difference equations with variable coefficients. The solution is obtained by using the graph theoretic, discrete path formalism. The central parameters in the solution are the crossing index and the crossing number. The transition from an enumerative graph theoretic solution to a closed form combinatorial solution is made possible by an isomorphism in-between paths on the signal flow graph, and n -tuplets of binary numbers.     

Two-dimensional wave equations with fractal boundaries

This paper focuses on two cases of two-dimensional wave equations with fractal boundaries. The first case is the equation with classical derivative. The formal solution is obtained. And a definition of the solution is given. Then we prove that under certain conditions, the solution is a kind of fractal function, which is continuous, differentiable nowhere in its domain. Next, for specific given initial position and 3 different initial velocities, the graphs of solutions are sketched. By computing the box dimensions of boundaries of cross-sections for solution surfaces, we evaluate the range of box dimension of the vibrating membrane. The second case is the equation with p-type derivative. The corresponding solution is shown and numerical example is given.     

Embedding evolutionary strategy in ordinal optimization for hard optimization problems

This work proposes a method for embedding evolutionary strategy (ES) in ordinal optimization (OO), abbreviated as ESOO, for solving real-time hard optimization problems with time-consuming evaluation of the objective function and a huge discrete solution space. Firstly, an approximate model that is based on a radial basis function (RBF) network is utilized to evaluate approximately the objective value of a solution. Secondly, ES associated with the approximate model is applied to generate a representative subset from a huge discrete solution space. Finally, the optimal computing budget allocation (OCBA) technique is adopted to select the best solution in the representative subset as the obtained “good enough” solution. The proposed method is applied to a hotel booking limits (HBL) problem, which is formulated as a stochastic combinatorial optimization problem with a huge discrete solution space. The good enough booking limits, obtained by the proposed method, have promising solution quality, and the computational efficiency of the method makes it suitable for real-time applications. To demonstrate the computational efficiency of the proposed method and the quality of the obtained solution, it is compared with two competing methods – the canonical ES and the genetic algorithm (GA). Test results demonstrate that the proposed approach greatly outperforms the canonical ES and GA.     

时滞均值回复θ过程及其数值解的收敛性

时滞均值回复θ过程用于描述受时间延迟影响的利率、波动率等金融特征,本文利用随机时滞微分方程理论证明了过程在1/2≤θ<1情况时解的存在唯一性和非负性.由于表示该过程的随机时滞微分方程没有显示解,所以数值近似解是研究过程的重要的方法,本文证明了时滞均值回复θ过程Euler-Maruyama数值解的p(p≥2)阶矩意义上的...     

Unique Solvability of the Inverse Problem with a Time-Dependent Boundary Condition for a Sorption Model

A mixed initial boundary-value problem is considered for nonequilibrium sorption dynamics with inner-diffusion kinetics. The problem allows for convection and longitudinal diffusion and has a time-dependent boundary condition. This condition contains the time derivative of a solution component and constitutes the balance equation for the absorbed mixture near the output boundary of the sorption region—inside the diffusion barrier. Bounds on the solution of the direct problem are obtained: nonnegativity of the solution and its first time derivatives, and boundedness of the solution by known functions. The inverse problem of estimating the nonlinear system parameter—the sorption isotherm—is considered and a uniqueness theorem is proved.     

Computational solution of a fractional generalization of the Schrödinger equation occurring in quantum mechanics

The object of this article is to present the solution of a fractional generalization of the Schrödinger equation of quantum mechanics in one dimension. The method followed in deriving the solution is that of joint Laplace and Fourier transforms. The solution is obtained in a compact and computational form in terms of the H-function. A result given earlier by Debnath for the solution of a generalized Schrödinger equation is obtained in an explicit form in terms of the H-function, as a special case of our findings.     

Steady solution and its stability for Navier–Stokes equations with general Navier slip boundary condition

The steady solution and the asymptotic behavior of the corresponding nonsteady solution are studied for Navier–Stokes equations under the general Navier slip boundary condition. The existence of a unique stationary solution is established. It is also proved that this solution is asymptotically stable under some restrictions on the data. Bibliography: 16 titles. Dedicated to Vsevolod Alekseevich Solonnikov on the occasion of his jubilee Published in Zapiski Nauchnykh Seminarov POMI, Vol. 362, 2008, pp. 153–175.     

Initial-boundary-value problems for the generalized multi-term time-fractional diffusion equation

In this paper, the initial-boundary-value problems for the generalized multi-term time-fractional diffusion equation over an open bounded domain G×(0,T), G∈Rn are considered. Based on an appropriate maximum principle that is formulated and proved in the paper, too, some a priory estimates for the solution and then its uniqueness are established. To show the existence of the solution, first a formal solution is constructed using the Fourier method of the separation of the variables. The time-dependent components of the solution are given in terms of the multinomial Mittag-Leffler function. Under certain conditions, the formal solution is shown to be a generalized solution of the initial-boundary-value problem for the generalized time-fractional multi-term diffusion equation that turns out to be a classical solution under some additional conditions. Another important consequence from the maximum principle is a continuously dependence of the solution on the problem data (initial and boundary conditions and a source function) that - together with the uniqueness and existence results - makes the problem under consideration to a well-posed problem in the Hadamard sense.     

Mathematical studies of the solution of Burgers' equations by Adomian decomposition method

In this paper, a novel Adomian decomposition method (ADM) is developed for the solution of Burgers' equation. While high level of this method for differential equations are found in the literature, this work covers most of the necessary details required to apply ADM for partial differential equations. The present ADM has the capability to produce three different types of solutions, namely, explicit exact solution, analytic solution, and semi-analytic solution. In the best cases, when a closed-form solution exists, ADM is able to capture this exact solution, while most of the numerical methods can only provide an approximation solution. The proposed ADM is validated using different test cases dealing with inviscid and viscous Burgers' equations. Satisfactory results are obtained for all test cases, and, particularly, results reported in this paper agree well with those reported by other researchers.     

First boundary-value problem in the half-strip for a parabolic-type equation with bessel operator and Riemann–Liouville derivative

The first boundary-value problem in the half-strip for a parabolic-type equation with Bessel operator and Riemann–Liouville derivative is studied. In the case of the zero initial condition, the representation of the solution in terms of the Fox H-function is obtained. The uniqueness of the solution for a class of functions vanishing at infinity is proved. It is shown that when the equation under consideration coincides with the Fourier equation, the obtained representation of the solution becomes the known representation of the solution of the corresponding problem.