Nonmonotone stabilization methods for nonlinear equations

We are concerned with defining new globalization criteria for solution methods of nonlinear equations. The current criteria used in these methods require a sufficient decrease of a particular merit function at each iteration of the algorithm. As was observed in the field of smooth unconstrained optimization, this descent requirement can considerably slow the rate of convergence of the sequence of points produced and, in some cases, can heavily deteriorate the performance of algorithms. The aim of this paper is to show that the global convergence of most methods proposed in the literature for solving systems of nonlinear equations can be obtained using less restrictive criteria that do not enforce a monotonic decrease of the chosen merit function. In particular, we show that a general stabilization scheme, recently proposed for the unconstrained minimization of continuously differentiable functions, can be extended to methods for the solution of nonlinear (nonsmooth) equations. This scheme includes different kinds of relaxation of the descent requirement and opens up the possibility of describing new classes of algorithms where the old monotone linesearch techniques are replace with more flexible nonmonotone stabilization procedures. As in the case of smooth unconstrained optimization, this should be the basis for defining more efficient algorithms with very good practical rates of convergence.This material is partially based on research supported by the Air Force Office of Scientific Research Grant AFOSR-89-0410, National Science Foundation Grant CCR-91-57632, and Istituto di Analisi dei Sistemi ed Informatica del CNR.     

Existence and asymptotic results for a system of integro-partial differential equations

A non-Fourier phase field model is considered. A global existence result for a Dirichlet, or generalized Neumann, initial-boundary value problem is obtained, followed by a discussion of the regularity and asymptotic properties of solutions ast.This research was supported in part by the National Science Foundation under Grant DMS 91-11794 and in part by the Italian M.U.R.S.T. project Problemi non lineari...Part of this author's work was done while visiting Ohio University.     

Interior-point algorithm for quadratically constrained entropy minimization problems

In this paper, we develop an interior point algorithm for quadratically constrained entropy problems. The algorithm uses a variation of Newton's method to follow a central path trajectory in the interior of the feasible set. The primal-dual gap is made less than a given in at most steps, wheren is the dimension of the problem andm is the number of quadratic inequality constraints.     

A generalized Dantzig—Wolfe decomposition principle for a class of nonconvex programming problems

Since Dantzig—Wolfe's pioneering contribution, the decomposition approach using a pricing mechanism has been developed for a wide class of mathematical programs. For convex programs a linear space of Lagrangean multipliers is enough to define price functions. For general mathematical programs the price functions could be defined by using a subclass of nondecreasing functions. However the space of nondecreasing functions is no longer finite dimensional. In this paper we consider a specific nonconvex optimization problem min {f(x):h _j (x)g(x),j=1, ,m, xX}, wheref(·),h _j (·) andg(·) are finite convex functions andX is a closed convex set. We generalize optimal price functions for this problem in such a way that the parameters of generalized price functions are defined in a finite dimensional space. Combining convex duality and a nonconvex duality we can develop a decomposition method to find a globally optimal solution.This paper is dedicated to Phil Wolfe on the occasion of his 65th birthday.     

On the global optimization properties of finite-difference local descent algorithms

The paper is devoted to the convergence properties of finite-difference local descent algorithms in global optimization problems with a special -convex structure. It is assumed that the objective function can be closely approximated by some smooth convex function. Stability properties of the perturbed gradient descent and coordinate descent methods are investigated. Basing on this results some global optimization properties of finite-difference local descent algorithms, in particular, coordinate descent method, are discovered. These properties are not inherent in methods using exact gradients.The paper was presented at the II. IIASA-Workshop on Global Optimization, Sopron (Hungary), December 9–14, 1990.     

Synthesis of substituted polymethylenes by radical polymerization of N,N,N′,N′-tetraalkylfumaramides and their characterization

Radical polymerization of N,N,N′,N′-tetraalkylfumaramides (TRFAm) bearing methyl, ethyl, n-propyl, isopropyl, and isobutyl groups as N-substituents (TMFAm, TEFAm, TnPFAm, TIPFAm, and TIBFAm, respectively) was investigated. In the polymerization of TEFAm initiated with 1,1′-azobiscyclohexane-1-carbonitrile (ACN) in benzene, the polymerization rate (R_p) was expressed as follows: R_p = k [ACN]^0.28 [TEFAm]^1.26, and the overall activation energy was 102.1 kJ/mol. The introduction of a bulky alkyl group into N-substituent of TRFAm decreased the R_p in the following order: TMFAm > TEFAm > TnPFAm > TIBFAm > TIPFAm ～ 0. The relative reactivities of these monomers were also investigated in radical copolymerization with styrene (St) and methyl methacrylate (MMA). In copolymerization of TRFAm (M₂) with St (M₁), monomer reactivity ratios were determined to be r₁ = 1.07 and r₂ = 0.20 for St–TMFAm, and r₁ = 1.88 and r₂ = 0.11 for St–TEFAm, from which Q₂ and e₂ values were estimated to be 0.35 and 0.44 for TMFAm, and 0.19 and 0.47 for TEFAm, respectively. The other TRFAm were also copolymerized with St, but copolymerization with MMA gave polymers containing a small amount of TRFAm units. The polymer from TRFAm consists of a less-flexible poly(N,N-dialkylaminocarbonylmethylene) structure. The solubility and thermal property of the polymers were also investigated.     

The relationship between the linear (oscillatory) and nonlinear (steady-state) flow properties of a series of polymer and colloidal systems

The Cox-Merz empirical relationship between the linear (oscillatory) and nonlinear (steady-state) viscosities has been shown to be valid for many polymeric systems. Here, we present an equivalent expression to relate the linear (G) and nonlinear (N ₁) elastic properties of viscoelastic systems. Like the Cox-Merz relationship, it uses a combination of elastic and viscous parameters. The modified form of the storage modulus is then equivalent to the Cox-Merz complex viscosity. It can be used to correlate with (half) the normal force at numerically equal circular frequency and shear rate, respectively.This new expression and the Cox-Merz rule are tested for a range of polymeric and colloidal systems. It is found that both expression work for the polymeric systems considered, but fail for the colloidal systems. In the latter, the steady state values of viscosity and elasticity are consistently low, and replacing them by the complex viscosity and our new elastic expression only makes matters worse.For polymer systems, we suggest this is a general but not universal observation, since we are aware of exceptions to the rule that polymeric systems obey the Cox-Merz rule for viscosity and our rule for elasticity. For colloidal systems we find that neither rule is obeyed for any of our systems.     

Nitroxide spin probe study of probe size,hydrogen bonding and polymer matrix rigidity effects on poly(acrylic acid)/poly(ethylene oxide) complexes

An electron spin resonance (ESR) spin probe study was performed on 1 : 1 by weight poly(acrylic acid) (PAA)/poly(ethylene oxide) (PEO) complex over the 100–450 K temperature range with a series of tetramethylpiperidyloxy‐based spin probes. Measurements of the parameters T_5mT, Ta and Td demonstrated the effects of probe size and the strength of hydrogen bonding. The probes in the series Tempone, Tempo, Tempol and Tamine (respectively 4‐oxo‐, unsubstituted, 4‐hydroxy‐ and 4‐amino‐2,2,6,6,‐tetramethylpiperidine ‐1‐oxyl) displayed noticeable increases in the hydrogen‐bonding effect, as indicated by Ta and Td. These increases correlated with increasing hydrogen bond acceptor strength. On the other hand, as the probe size became larger, T_5mT gradually increased due to the free volume decrease. These effects were analyzed using the established theoretical relationship of T_5mT to probe volume expressed by f. Meanwhile, in order to investigate the effect of polymer matrix rigidity, a similar study was performed with a nitroxide spin probe, 2,2,6,6‐tetramethyl‐1‐piperidine‐1‐oxyl (Tempo), on PAA/PEO complexes of different weight compositions. The quantitative fast motion fraction in the composite ESR spectrum was calculated. The influence of changes in the composition of PAA on the molecular mobility was characterized by changes of the spectral parameters and τ_c. The molecular mobility was shown to diminish with increasing content of PAA in PAA/PEO blends duo to the restriction of the polymer matrix rigidity increase. Copyright © 2003 John Wiley & Sons, Ltd.     

通过脱碳控制全球平均温度

Establishing a reliable method to predict the global mean temperature (T_e) is of great importance because CO₂ reduction activities require political and global cooperation and significant financial resources. The current climate models all seem to predict that the earth's temperature will continue to increase, mainly based on the assumption that CO₂ emissions cannot be lowered significantly in the foreseeable future. Given the earth's multifactor climate system, attributing atmospheric CO₂ as the only cause for the observed temperature anomaly is most likely an oversimplification; the presence of water (H₂O) in the atmosphere should at least be considered. As such, T_e is determined by atmospheric water content controlled by solar activity, along with anthropogenic CO₂ activities. It is possible that the anthropogenic CO₂ activities can be reduced in the future. Based on temperature measurements and thermodynamic data, a new model for predicting T_e has been developed. Using this model, past, current, and future CO₂ and H₂O data can be analyzed and the associated T_e calculated. This new, esoteric approach is more accurate than various other models, but has not been reported in the open literature. According to this model, by 2050, T_e may increase to 15.5 ℃ under "business-as-usual" emissions. By applying a reasonable green technology activity scenario, T_e may be reduced to approximately 14.2 ℃. To achieve CO₂ reductions, the scenario described herein predicts a CO₂ reduction potential of 513 gigatons in 30 years. This proposed scenario includes various CO₂ reduction activities, carbon capturing technology, mineralization, and bio-char production; the most important CO₂ reductions by 2050 are expected to be achieved mainly in the electricity, agriculture, and transportation sectors. Other more aggressive and plausible drawdown scenarios have been analyzed as well, yielding CO₂ reduction potentials of 1051 and 1747 gigatons, respectively, in 30 years, but they may reduce global food production. It is emphasized that the causes and predictions of the global warming trend should be regarded as open scientific questions because several details concerning the physical processes associated with global warming remain uncertain. For example, the role of solar activities coupled with Milankovitch cycles are not yet fully understood. In addition, other factors, such as ocean CO₂ uptake and volcanic activity, may not be negligible.     

Efficient algorithms for solving rank two and rank three bilinear programming problems

Finitely convergent algorithms for solving rank two and three bilinear programming problems are proposed. A rank k bilinear programming problem is a nonconvex quadratic programming problem with the following structure: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaaieaacaWFTbGaa8% xAaiaa-5gacaWFPbGaa8xBaiaa-LgacaWF6bGaa8xzaiaa-bcacaWF% 7bacbiGaa43yamaaDaaaleaacaGFWaaabaGaa4hDaaaakiaa+Hhaca% GFRaGaa4hzamaaDaaaleaacaGFWaaabaGaa4hDaaaakiaa+LhacaGF% RaWaaabuaeaacaGFJbWaa0baaSqaaiaa+PgaaeaacaGF0baaaOGaam% iEaiabl+y6NjaadsgadaqhaaWcbaGaamOAaaqaaiaadshaaaGccaWG% 5bGaaiiFaaWcbaGaa8NAaiaa-1dacaWFXaaabeqdcqGHris5aOGaa4% hEaiabgIGiolaa+HfacaGFSaGaa4xEaiabgIGiolaa+LfacaWF9bGa% a8hlaaaa!5D2E!\[minimize \{ c_0^t x + d_0^t y + \sum\limits_{j = 1} {c_j^t xd_j^t y|} x \in X,y \in Y\} ,\]where X Rⁿ¹ and Y R ⁿ² are non-empty and bounded polytopes. We show that a variant of parametric simplex algorithm can solve large scale rank two bilinear programming problems efficiently. Also, we show that a cutting-cake algorithm, a more elaborate variant of parametric simplex algorithm can solve medium scale rank three problems.This research was supported in part by Grant-in-Aid for Scientific Research of the Ministry of Education, Science and Culture, Grant No. 63490010.