A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

Stability and slightly supercritical oscillatory regimes of natural convection in a 8:1 cavity: solution of the benchmark problem by a global Galerkin method

The global Galerkin method is applied to the benchmark problem that considers an oscillatory regime of convection of air in a tall two‐dimensional rectangular cavity. The three most unstable modes of the linearized system of the Boussinesq equations are studied. The converged values of the critical Rayleigh numbers together with the corresponding oscillation frequencies are calculated for each mode. The oscillatory flow regimes corresponding to each of the three modes are approximated asymptotically. No direct time integration is applied. Good agreement with the previously published results obtained by solution of the time‐dependent Boussinesq equations is reported. Copyright © 2004 John Wiley & Sons, Ltd.     

High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Free boundary viscous flows at micro and mesolevel during liquid composites moulding process

Numerical simulation aspects, related to low Reynolds number free boundary viscous flows at micro and mesolevel during the resin impregnation stage of the liquid composite moulding process (LCM), are presented in this article. A free boundary program (FBP), developed by the authors, is used to track the movement of the resin front accurately by accounting for the surface tension effects at the boundary. Issues related to the global and local mass conservation (GMC and LMC) are identified and discussed. Unsuitable conditions for LMC and consequently GMC are uncovered at low capillary numbers, and hence a strategy for the numerical simulation of such flows is suggested. FBP encompasses a set of subroutines that are linked to modules in ANSYS. FBP can capture the void formation dynamics based on the analysis developed. We present resin impregnation dynamics in two dimensions. Extension to three dimensions is a subject for further research. Several examples are shown and efficiency of different stabilization techniques are compared. Copyright © 2004 John Wiley & Sons, Ltd.     

Composite Silica Particles Using a Mixture of Organosilane Monomers

Composite silica particles were synthesized by a two-step (acid-base) process in an aqueous solution with a mixture of organoalkoxysilane monomers. The two-step process separates the hydrolysis and condensation procedures to easily control condensation rate. In this study, the silane monomers used were phenyltrimethoxysilane (PTMS), vinyltrimethoxysilane (VTMS), methyltrimethoxysilane (MTMS), and tetraethyl-orthosilicate (TEOS). The physical properties of the resultant composite particles were investigated with the change in the molar ratio of monomers. The size of the particles increased with increasing the molar ratio of R_aSi(OR)₃/R_bSi(OR)₃ or R_aSi(OR)₃/TEOS (R_a: phenyl; R_b: vinyl, methyl).     

Two-Layered Poly-o-phenylenediamine/Polyaniline Composite: Electrosynthesis and Spectroelectrochemical Properties

Spectroelectrochemical properties of poly-o-phenylenediamine (PPD) and the synthesized composite of PPD and polyaniline—two chemically related polymers containing an amino-substituted benzene ring but having different conduction nature—are studied. The polyaniline synthesis on PPD-modified electrodes involves stages of the reaction initiation, the copolymer formation, and the formation of a polyaniline layer at the copolymer/solution interface.     

On General Mixed Quasivariational Inequalities

In this paper, we suggest and analyze several iterative methods for solving general mixed quasivariational inequalities by using the technique of updating the solution and the auxiliary principle. It is shown that the convergence of these methods requires either the pseudomonotonicity or the partially relaxed strong monotonicity of the operator. Proofs of convergence is very simple. Our new methods differ from the existing methods for solving various classes of variational inequalities and related optimization problems. Various special cases are also discussed.     

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.     

A relaxed model and its homogenization for nematic liquid crystals in composite materials

We analyse a model for equilibrium configurations of composite systems of nematic liquid crystal with polymer inclusions, in the presence of an external magnetic field. We assume that the system has a periodic structure, and consider the relaxed problem on the unit length constraint of the nematic director field. The relaxation of the Oseen–Frank energy functional is carried out by including bulk as well as surface energy penalty terms, rendering the problem fully non‐linear. We employ two‐scale convergence methods to obtain effective configurations of the system, as the size of the polymeric inclusions tends to zero. We discuss the minimizers of the effective energies for, both, the constrained as well as the unconstrained models. Copyright © 2004 John Wiley & Sons, Ltd.     

Diffusion in poro‐plastic media

A model is developed for the flow of a slightly compressible fluid through a saturated inelastic porous medium. The initial‐boundary‐value problem is a system that consists of the diffusion equation for the fluid coupled to the momentum equation for the porous solid together with a constitutive law which includes a possibly hysteretic relation of elasto‐visco‐plastic type. The variational form of this problem in Hilbert space is a non‐linear evolution equation for which the existence and uniqueness of a global strong solution is proved by means of monotonicity methods. Various degenerate situations are permitted, such as incompressible fluid, negligible porosity, or a quasi‐static momentum equation. The essential sufficient conditions for the well‐posedness of the system consist of an ellipticity condition on the term for diffusion of fluid and either a viscous or a hardening assumption in the constitutive relation for the porous solid. Copyright © 2004 John Wiley & Sons, Ltd.