Augmented Lagrangian Algorithms Based on the Spectral Projected Gradient Method for Solving Nonlinear Programming Problems

The spectral projected gradient method SPG is an algorithm for large-scale bound-constrained optimization introduced recently by Birgin, Martínez, and Raydan. It is based on the Raydan unconstrained generalization of the Barzilai-Borwein method for quadratics. The SPG algorithm turned out to be surprisingly effective for solving many large-scale minimization problems with box constraints. Therefore, it is natural to test its perfomance for solving the sub-problems that appear in nonlinear programming methods based on augmented Lagrangians. In this work, augmented Lagrangian methods which use SPG as the underlying convex-constraint solver are introduced (ALSPG) and the methods are tested in two sets of problems. First, a meaningful subset of large-scale nonlinearly constrained problems of the CUTE collection is solved and compared with the perfomance of LANCELOT. Second, a family of location problems in the minimax formulation is solved against the package FFSQP.     

Wavelets and adaptive grids for the discontinuous Galerkin method

In this paper, space adaptivity is introduced to control the error in the numerical solution of hyperbolic systems of conservation laws. The reference numerical scheme is a new version of the discontinuous Galerkin method, which uses an implicit diffusive term in the direction of the streamlines, for stability purposes. The decision whether to refine or to unrefine the grid in a certain location is taken according to the magnitude of wavelet coefficients, which are indicators of local smoothness of the numerical solution. Numerical solutions of the nonlinear Euler equations illustrate the efficiency of the method.     

Fast electrons generation by RF and random fields near the antenna I. One dimensional model and regular fields

Test particle motion and acceleration has been explored in strong radio frequency (RF) fields, for which quasilinear ponderomotive force approximation is not valid. By nonlinear acceleration in spatially varying wave amplitude of RF travelling wave, electrons may be accelerated to time averaged velocities significantly larger than the RF wave phase velocity, and than the boundary plasma thermal velocity, in RF fields of several Volts per centimeter at wave frequency of 7 MHz. It is also demonstrated that even weak spatial gradients, much weaker than those expected in experiments, of the RF wave field amplitude, have significant consequences for the particle motion. Estimates are presented of the total energy transferred from the near antenna RF field to the plasma due to the nonlinear electron acceleration effects.     

Transient binary nucleation from ethanol-hexanol vapour

We study condensation of ethanol-hexanol vapour by numerical solution of kinetic equations. The number of droplets formed in unit volume is computed within self-consistent classical model. It is shown that formation of ethanol-rich droplets prevails at the initial stage of nucleation process, but in the stationary state formation of droplets near the saddle point (on cluster formation energy surface) plays the dominant role. Presented at the 6th Joint Seminar “Development of Materials Science in Research and Education”, Karlštejn, Czech Republic, 17–19 September 1996. This work was supported by Grant No. A1010615 of the Grant Agency of the Academy of Sciences of the Czech Republic.     

Study of the gas-phase parameters affecting the silicon-oxide film deposition induced by an ArF laser

A study of the gas-phase parameters involved in ArF laser induced chemical vapour deposition of silicon-oxide thin films is presented. A complete set of experiments has been performed showing the influence of the concentration of the precursor gases, N₂O and SiH₄, and their influence on total and partial pressures on film growth and properties. In this paper we demonstrate the ability of this LCVD method to deposit silicon oxide films of different compositions and densities by appropriate control of gas composition and total pressure. Moreover, a material specific calibration plot comprising data obtained using different preparation techniques is presented, allowing determination of the stoichiometry of SiO _x films by using FTIR spectroscopy independently of the deposition method. For the range of processing conditions examined, the experimental results suggest that chemical processes governing deposition take place mainly in the gas phase.     

Molecular dynamics simulations of the oligonucleotide with the modified phosphate/phosphonate internucleotide linkage

Impact of the internucleotide linkage modification by inserting a methylene group to the P-O bond (—O—PO ₂ ^— —O— chain changed for —O—PO ₂ ^— —CH₂—O—), on the modified oligonucleotide binding ability to the natural DNA strand was studied by molecular dynamics simulation. Complex of (dT)₁₁ with a deoxyadenosine undecamer containing alternating modified and natural internucleotide linkage was studied as a model system. The Amber force field was completed by a set of new parameters needed to model the modified part of the nucleotide. The simulations confirmed existence of a double-helical complex the melting point of which is considerably higher than 300 K. While the thymidine (unmodified) strand possesses a B-type secondary structure, the conformation of the adenosine (modified) strand is not stable at 300 K. The -ggt conformation of the modified linkages is highly preferred, temporary jumps to the -g-gt and ggt conformations were, however, observed.     

Global Analysis of the Flows of Fluids with Pressure-Dependent Viscosities

 To describe the flows of fluids over a wide range of pressures, it is necessary to take into account the fact that the viscosity of the fluid depends on the pressure. That the viscosity depends on the pressure has been verified by numerous careful experiments. While the existence of solutions local-in-time to the equations governing the flows of such fluids are available for small, special data and rather unrealistic dependence of the viscosity on the pressure, no global existence results are in place. Our interest here is to establish the existence of weak solutions for spatially periodic three-dimensional flows that are global in time, for a large class of physically meaningful viscosity-pressure relationships. (Accepted May 1, 2002) Published online November 15, 2002 Communicated by S. S. ANTMAN     

A system of parabolic equations in nonequilibrium thermodynamics including thermal and electrical effects

The time-dependent equations for a charged gas or fluid consisting of several components, exposed to an electric field, are considered. These equations form a system of strongly coupled, quasilinear parabolic equations which in some situations can be derived from the Boltzmann equation. The model uses the duality between the thermodynamic fluxes and the thermodynamic forces. Physically motivated mixed Dirichlet-Neumann boundary conditions and initial conditions are prescribed.The existence of weak solutions is proven. The key of the proof is (i) a transformation of the problem by using the entropic variables, or electro-chemical potentials, which symmetrizes the equations, and (ii) a priori estimates obtained by using the entropy function. Finally, the entropy inequality is employed to show the convergence of the solutions to the thermal equilibrium state as the time tends to infinity.