Finite Element Methods for a Model for Full Waveform Acoustic Logging

A model is defined to simulate the propagation of waves in aradially symmetric, isotropic, composite system consisting ofa fluid-filled well bore ^f through a fluid-saturated poroussolid ^p. Biot's equations of motion are chosen to describe thepropagation of waves in ^p, while the standard equation of motionfor compressible inviscid fluids is used for ^f, with appropriateboundary conditions at the contact surface between ^f and ^p.Also, absorbing boundary conditions for the artificial boundariesof ^p are derived for the model, their effect being to make themtransparent for waves arriving normally First, results on the existence and uniqueness of the solutionof the differential problem are given and then a discrete-time,explicit finite element procedure is defined and analysed, withfinite element spaces suited for radially symmetric problemsbeing used for the spatial discretisation.     

Modelling,simulation, animation,and real-time control (MoSART) for a class of electromechanical systems: a system-theoretic approach

This paper describes an interactive modelling, simulation, animation, and real-time control (MoSART) environment for a class of ‘cart-pendulum’ electromechanical systems that may be used to enhance learning within differential equations and linear algebra classes. The environment is useful for conveying fundamental mathematical/systems concepts through computer-aided analysis, design, graphical visualization, and 3D animation. Referred to as Cart-Pendulum Control3D-Lab1Control3D-Lab-short for control 3D animation laboratory. The author to whom correspondence should be addressed., the environment is based on Microsoft Windows/Visual C++/Direct-3D and MATLAB/Simulink2MATLAB and Simulant are trademarks of The Masterworks, Inc.. Pull-down menus provide access to systems/models/control laws, exogenous signals, parameters, animation models and visual indicators, a suite of (easy-to-modify) Simulant diagrams containing models and control laws, MATLAB m-files for detailed analysis and design, and detailed documentation for each of the above. Three blocks enable animation, joystick inputs, and real-time animation within Simulant. Examples are presented to illustrate the utility of the environment as a valuable tool for analysing/visualizing the above class of electromechanical systems and for enhancing mathematics instruction.     

Calculation of surface tension via area sampling

We examine the performance of several molecular simulation techniques aimed at evaluation of the surface tension through its thermodynamic definition. For all methods explored, the surface tension is calculated by approximating the change in Helmholtz free energy associated with a change in interfacial area through simulation of a liquid slab at constant particle number, volume, and temperature. The methods explored fall within three general classes: free-energy perturbation, the Bennett acceptance-ratio scheme, and the expanded ensemble technique. Calculations are performed for both the truncated Lennard-Jones and square-well fluids at select temperatures spaced along their respective liquid-vapor saturation lines. Overall, we find that Bennett and expanded ensemble approaches provide the best combination of accuracy and precision. All of the methods, when applied using sufficiently small area perturbation, generate equivalent results for the Lennard-Jones fluid. However, single-stage free-energy-perturbation methods and the closely related test-area technique recently introduced by Gloor et al. [J. Chem. Phys. 123, 134703 (2005)] generate surface tension values for the square-well fluid that are not consistent with those obtained from the more robust expanded ensemble and Bennett approaches, regardless of the size of the area perturbation. Single-stage perturbation methods fail also for the Lennard-Jones system when applied using large area perturbations. Here an analysis of phase-space overlap produces a quantitative explanation of the observed inaccuracy and shows that the satisfactory results obtained in these cases from the test-area method arise from a cancellation of errors that cannot be expected in general. We also briefly analyze the variation in method performance with respect to the adjustable parameters inherent to the techniques.     

The Saffman-Taylor problem for an extremely shear-thinning fluid

We consider a steady flow driven by pushing a finger of gasinto a highly shear-thinning power-law fluid, with exponentn, in a Hele-Shaw channel. We formulate the problem in termsof the streamfunction , which satisfies the p-Laplacian equation (with ), and investigate travelling wave solutions in the large-n (extreme shear-thinning) limit.We take a Legendre transform of the free-boundary problem for, which reduces it to a linear problem on a fixed domain. The solution to this problem is foundby using matched asymptotic expansions and the resulting shapeof the finger deduced (being, to leading order, a semi-infinitestrip). The nonlinear problem for the streamfunction is alsotreated using matched asymptotic expansion in the physical plane.The finger-width selection problem is briefly discussed in termsof our results.     

Supramolecular circular helicates formed by destabilisation of supramolecular dimers

The effect of changes in the angles at the connection points of linear/circular helicates is explored as a route to control the nuclearity and architecture of metallo-supramolecular arrays. This effect is probed by changing the geometry of the metal centre used to assemble bis-pyridylimine ligands that contain a 1,3-bis(aminomethyl) benzene spacer group. Tetrahedral metal ions favour linear dimers, whereas octahedral nickel(II) predominantly gives a triangular circular helicate. Five-coordinate copper(II) falls in the middle of these extremes and results in the formation of solvent-dependent mixtures of dimer and trimer. The trinuclear, triangular, circular helicate structures, which result from coordination to copper(II) and nickel(II), are structurally characterised by X-ray crystallography and reveal that the units can aggregate into hexagonal arrays that contain anion-filled tube-like channels in the solid state.     

Investigation of the phase behavior of an embedded charge protein model through molecular simulation

The phase behavior of an embedded-charge model for lysozyme developed by Carlsson and co-workers (J. Phys. Chem. B 2001, 105, 9040) is investigated using grand canonical transition matrix Monte Carlo simulation. Within this model, protein-protein interactions are approximated through a combination of hard-sphere repulsion, isotropic hydrophobic attraction, and screened electrostatic interactions through a series of embedded point charges located at the positions of charged amino acid groups within lysozyme. Liquid-liquid phase diagrams are constructed for a wide range of solution conditions and compared with experimental data. Our results indicate that the model is generally capable of describing qualitative trends in the evolution of protein phase behavior with variation of pH and ionic strength. From a quantitative perspective, model estimates for both the change in critical temperature with variation of the solution conditions and the critical concentration do not agree with experimental results. We find the width of model coexistence curves to be independent of solution conditions and narrow relative to experimentally obtained phase envelopes. Connections between the value of the second virial coefficient evaluated at the critical temperature and the location of the liquid-liquid phase envelope are also examined.     

HUMIC MACROMOLECULES IN NATURAL WATERS

Abstract

Humic substances are the major organic constituents of soils and sediments. They also occur in small concentrations in natural surface waters and groundwaters. They form through the breakdown of plant and animal tissues by chemical and biological processes that tend to produce complex chemical structures that are more stable than the original material from which they were derived. One of the more important characteristics of humic substances is their ability to form water-soluble and water-insoluble complexes with metal ions and hydrous oxides and to interact with clay minerals and various organic compounds such as alkanes, fatty acids, and toxic organic substances such as pesticides.     

Rational addition of capping groups to the phosphomolybdate Keggin anion [PMo12O40](3-) by mild, non-aqueous reductive aggregation

Controlled reductive assembly of capped Keggin anions [PMo(12)O(40)(ML(m))(n)](3-) has been achieved by reduction of [PMo(12)O(40)](3-) with sodium-mercury amalgam in the presence of metal halides, as exemplified by the rational syntheses of mono-capped [PMo(12)O(40){Co(MeCN)(2)}](3-) and bi-capped [PMo(12)O(40)(VO)(2)](3-) and [PMo(12)O(40)Sb(2)](3-).     

Quantitative link between single-particle dynamics and static structure of supercooled liquids

We present evidence via molecular simulation that the supercooled fluid states of SPC/E water as well as the "repulsive" and "attractive" supercooled fluid states of a recently introduced model for colloids with short-ranged attractions are characterized by the same functional relationship between self-diffusivity and the pair correlation function. We discuss how this simple relationship connects to an earlier finding that the temperature dependency of a supercooled fluid's single-particle dynamics tracks that of its excess entropy (relative to ideal gas). The generality of this observed structure-property relationship is supported by its ability to successfully describe the nontrivial behaviors of these very different types of model systems.     

Excess-entropy-based anomalies for a waterlike fluid

Many thermodynamic and dynamic properties of water display unusual behavior at low enough temperatures. In a recent study, Yan et al. [Phys. Rev. Lett. 95, 130604 (2005)] identified a spherically symmetric two-scale potential that displays many of the same anomalous properties as water. More specifically, for select parametrizations of the potential, one finds that the regions where isothermal compression anomalously (i) decreases the fluid's structural order, (ii) increases its translational self-diffusivity, and (iii) increases its entropy form nested domes in the temperature-density plane. These property relationships are similar to those found for more realistic models of water. In this work, the authors provide evidence that suggests that the anomalous regions specified above can all be linked through knowledge of the excess entropy. Specifically, the authors show how entropy scaling relationships developed by Rosenfeld [Phys. Rev. A 15, 2545 (1977)] can be used to describe the region of diffusivity anomalies and to predict the state conditions for which anomalous viscosity and thermal conductivity behavior might be found.