The solution of linear and nonlinear systems of Volterra functional equations using Adomian–Pade technique

The purpose of this study is to implement Adomian–Pade (Modified Adomian–Pade) technique, which is a combination of Adomian decomposition method (Modified Adomian decomposition method) and Pade approximation, for solving linear and nonlinear systems of Volterra functional equations. The results obtained by using Adomian–Pade (Modified Adomian–Pade) technique, are compared to those obtained by using Adomian decomposition method (Modified Adomian decomposition method) alone. The numerical results, demonstrate that ADM–PADE (MADM–PADE) technique, gives the approximate solution with faster convergence rate and higher accuracy than using the standard ADM (MADM).     

A Class of Nonmonotone Conjugate Gradient Methods for Unconstrained Optimization

In this paper, we introduce a class of nonmonotone conjugate gradient methods, which include the well-known Polak–Ribière method and Hestenes–Stiefel method as special cases. This class of nonmonotone conjugate gradient methods is proved to be globally convergent when it is applied to solve unconstrained optimization problems with convex objective functions. Numerical experiments show that the nonmonotone Polak–Ribière method and Hestenes–Stiefel method in this nonmonotone conjugate gradient class are competitive vis-à-vis their monotone counterparts.     

On the solution of Dirichlet’s problem of the complex Monge–Ampère equation for the Cartan–Hartogs domain of the first type

The complex Monge–Ampère equation is a nonlinear equation with high degree; therefore getting its solution is very difficult. In the present paper how to get the solution of Dirichlet’s problem of the complex Monge–Ampère equation on the Cartan–Hartogs domain of the first type is discussed, using an analytic method. Firstly, the complex Monge–Ampère equation is reduced to a nonlinear ordinary differential equation, then the solution of Dirichlet’s problem of the complex Monge–Ampère equation is reduced to the solution of a two-point boundary value problem for a nonlinear second-order ordinary differential equation. Secondly, the solution of Dirichlet’s problem is given as a semi-explicit formula, and in a special case the exact solution is obtained. These results may be helpful for a numerical method approach to Dirichlet’s problem of the complex Monge–Ampère equation on the Cartan–Hartogs domain of the first type.     

Global Convergence of Conjugate Gradient Methods without Line Search

Global convergence results are derived for well-known conjugate gradient methods in which the line search step is replaced by a step whose length is determined by a formula. The results include the following cases: (1) The Fletcher–Reeves method, the Hestenes–Stiefel method, and the Dai–Yuan method applied to a strongly convex LC ¹ objective function; (2) The Polak–Ribière method and the Conjugate Descent method applied to a general, not necessarily convex, LC ¹ objective function.     

Chebyshev pseudospectral method for computing numerical solution of convection–diffusion equation

A method for computing highly accurate numerical solutions of 1D convection–diffusion equations is proposed. In this method, the equation is first discretized with respect to the spatial variable, transforming the original problem into a set of ordinary differential equations, and then the resulting system is integrated in time by the fourth-order Runge–Kutta method. Spatial discretization is done by using the Chebyshev pseudospectral collocation method. Before describing the method, we review a finite difference-based method by Salkuyeh [D. Khojasteh Salkuyeh, On the finite difference approximation to the convection–diffusion equation, Appl. Math. Comput. 179 (2006) 79–86], and, contrary to the proposal of the author, we show that this method is not suitable for problems involving time dependent boundary conditions, which calls for revision. Stability analysis based on pseudoeigenvalues to determine the maximum time step for the proposed method is also carried out. Superiority of the proposed method over a revised version of Salkuyeh’s method is verified by numerical examples.     

Ultrasonic method of determining the elastic characteristics of anisotropic materials

An ultrasonic method of determining the elastic characteristics, developed for isotropic materials and based on the measurement of the angles of total internal reflection of ultrasound at a liquid-solid interface, is applied to an anisotropic Compreg of the DSPA type. The results obtained for the shear modulus using the ultrasonic method at 2.5 MHz and an independent method based on measuring the frequency of the damped torsional vibrations of a prismatic bar at 10–30 Hz are compared. It is shown that the shear modulus varies only slightly (10–12%) on the frequency interval investigated.Leningrad Kirov Forest Engineering Academy. Translated from Mekhanika Polimerov, No. 4, pp. 728–731, July–August, 1969.     

Globally Convergent Broyden-Like Methods for Semismooth Equations and Applications to VIP, NCP and MCP

In this paper, we propose a general smoothing Broyden-like quasi-Newton method for solving a class of nonsmooth equations. Under appropriate conditions, the proposed method converges to a solution of the equation globally and superlinearly. In particular, the proposed method provides the possibility of developing a quasi-Newton method that enjoys superlinear convergence even if strict complementarity fails to hold. We pay particular attention to semismooth equations arising from nonlinear complementarity problems, mixed complementarity problems and variational inequality problems. We show that under certain conditions, the related methods based on the perturbed Fischer–Burmeister function, Chen–Harker–Kanzow–Smale smoothing function and the Gabriel–Moré class of smoothing functions converge globally and superlinearly.     

A higher-order moment method of the lattice Boltzmann model for the conservation law equation

In this paper, we proposed a higher-order moment method in the lattice Boltzmann model for the conservation law equation. In contrast to the lattice Bhatnagar–Gross–Krook (BGK) model, the higher-order moment method has a wide flexibility to select equilibrium distribution function. This method is based on so-called a series of partial differential equations obtained by using multi-scale technique and Chapman–Enskog expansion. According to Hirt’s heuristic stability theory, the stability of the scheme can be controlled by modulating some special moments to design the third-order dispersion term and the fourth-order dissipation term. As results, the conservation law equation is recovered with higher-order truncation error. The numerical examples show the higher-order moment method can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the conservation law equation.     

Primal–Dual–Infeasible Newton Approach for the Analytic Center Deep-Cutting Plane Method

The convergence and complexity of a primal–dual column generation and cutting plane algorithm from approximate analytic centers for solving convex feasibility problems defined by a deep cut separation oracle is studied. The primal–dual–infeasible Newton method is used to generate a primal–dual updating direction. The number of recentering steps is O(1) for cuts as deep as half way to the deepest cut, where the deepest cut is tangent to the primal–dual variant of Dikin's ellipsoid.     

On a Commutative Class of Search Directions for Linear Programming over Symmetric Cones

The commutative class of search directions for semidefinite programming was first proposed by Monteiro and Zhang (Ref. 1). In this paper, we investigate the corresponding class of search directions for linear programming over symmetric cones, which is a class of convex optimization problems including linear programming, second-order cone programming, and semidefinite programming as special cases. Complexity results are established for short-step, semilong-step, and long-step algorithms. Then, we propose a subclass of the commutative class for which we can prove polynomial complexities of the interior-point method using semilong steps and long steps. This subclass still contains the Nesterov–Todd direction and the Helmberg–Rendl–Vanderbei–Wolkowicz/Kojima–Shindoh–Hara/Monteiro direction. An explicit formula to calculate any member of the class is also given.     

Practical enhancements to the Magnanti–Wong method

The Magnanti–Wong method–accelerating Benders decomposition–is shown to exhibit difficulties due to its dependence on the subproblem; an independent version is therefore introduced. The method additionally requires a–sometimes intractable–master problem core point; for several applications it is proved and experimentally verified that alternative points may be used.     

Interior penalty discontinuous approximations of convection–diffusion problems with parabolic layers

Summary A nonsymmetric discontinuous Galerkin finite element method with interior penalties is considered for two–dimensional convection–diffusion problems with regular and parabolic layers. On an anisotropic Shishkin–type mesh with bilinear elements we prove error estimates (uniformly in the perturbation parameter) in an integral norm associated with this method. On different types of interelement edges we derive the values of discontinuity–penalization parameters. Numerical experiments complement the theoretical results.     

A nonsmooth Levenberg–Marquardt method for solving semi-infinite programming problems

In this paper, we first transform the semi-infinite programming problem into the KKT system by the techniques in [D.H. Li, L. Qi, J. Tam, S.Y. Wu, A smoothing Newton method for semi-infinite programming, J. Global. Optim. 30 (2004) 169–194; L. Qi, S.Y. Wu, G.L. Zhou, Semismooth Newton methods for solving semi-infinite programming problems, J. Global. Optim. 27 (2003) 215–232]. Then a nonsmooth and inexact Levenberg–Marquardt method is proposed for solving this KKT system based on [H. Dan, N. Yamashita, M. Fukushima, Convergence properties of the inexact Levenberg–Marquardt method under local error bound conditions, Optimim. Methods Softw., 11 (2002) 605–626]. This method is globally and superlinearly (even quadratically) convergent. Finally, some numerical results are given.     

Variation of Curvatures of Subanalytic Spaces and Schläfli-Type Formulas

We prove a strong variational formula for Lipschitz–Killing curvaturesof subanalytic sets. As corollaries, we reprove the Chern–Gauss–Bonnettheorem and higher Schläfli formulas. The proof of the variationalformula uses normal cycles of subanalytic sets and a new method allowinga reduction from the difficult singular geometry to computations withdifferential forms.     

Extended Gauss–Markov Theorem for Nonparametric Mixed-Effects Models

The Gauss–Markov theorem provides a golden standard for constructing the best linear unbiased estimation for linear models. The main purpose of this article is to extend the Gauss–Markov theorem to include nonparametric mixed-effects models. The extended Gauss–Markov estimation (or prediction) is shown to be equivalent to a regularization method and its minimaxity is addressed. The resulting Gauss–Markov estimation serves as an oracle to guide the exploration for effective nonlinear estimators adaptively. Various examples are discussed. Particularly, the wavelet nonparametric regression example and its connection with a Sobolev regularization is presented.     

A new method of measuring the viscoelastic properties of polymer systems in fluid and highly elastic states using forced vibrations in the ultralow-frequency region

A method is considered for measuring the dynamic mechanical characteristics of polyer systems using forced vibrations in the ultralow-frequency region down to 10^–6 Hz. The essence of this method lies in the fact that during the measurement process the motion of the polymer sample is controlled with a set amplitude and frequency, and the dynamic characteristics of the polymer are determined from the controlling mechanical stress. The method is illustrated using the results of measurements of the components of the complex modulus for polybutadiene over a wide range of temperatures and frequencies without making use of superimposition.The Central Constructional Bureau for Prototype Instrument Construction, Academy of Sciences of the USSR, Moscow. Translated from Mekhanika Polimerov, No. 2, pp. 344–347, March–April, 1976.     

Procedure for measuring the thermophysical properties of polymers at high hydrostatic pressures

A procedure is suggested for determining the thermophysical properties of solids under hydrostatic pressure conditions by the instantaneous cylindrical heat source method without allowance for heat exchange at the sample-liquid boundary. A calculation performed showed that the temperature field of a cylinder of finite radius is identical to the temperature field of an infinite body up to Fo=0.05–0.07 with an accuracy far better than 0.5%. It is shown that the sample dimensions can be small thereupon at measurement times sufficient for performing the measurements on a standard apparatus of the order of 40–70 sec.Moscow Chemical Engineering Institute. Translated from Mekhanika Polimerov, Vol. 9, No. 3, pp. 568–570, May–June, 1973.     

Resistance of composite flywheel rim to radial tensile stresses from centrifugal forces based on results of pure bending of a rim segment

A method is proposed for evaluating the resistance of a flywheel rim to radial stresses in free rotation. The method is based on loading a rim segment in pure bending and calculating the limiting moment and the corresponding limiting angular velocity. Applicability of the method is substantiated theoretically by investigating the similarity of the radial stress diagrams in rotation and pure bending. The method is verified experimentally for the strained state of a rim segment in pure bending.Translated from Mekhanika Kompozitnykh Materialov, Vol. 29, No. 4, pp. 521–526, July–August, 1993.     

Creep of polymer materials under periodically varying loads

A method is proposed for constructing the creep curves of a material whose nonlinear memory properties are described by Rozovskii's nonlinear integral equation [2] (with allowance for the stress dependence of the relaxation time) under given periodic loading from known creep curves recorded at constant stress. In deriving the theoretical relation certain simplifying assumptions are made (the creep strain accumulated in 1–2 cycles is small, no vibration [4–6]). An experimental check shows that the proposed method can be used to predict the behavior of a material under periodic loading with an accuracy sufficient for practical purposes.Mekhanika Polimerov, Vol. 2, No. 3, pp. 330–336, 1966