A tentative approach to the direct simulation of drag reduction by polymers

An efficient technique for drag reduction uses dilute solutions of a few p.p.m. of polymers. A possible reduction in drag of up to 80% is achieved. Several experimental observations have been made which tend to indicate that the polymers modify the turbulence structures within the buffer layer. Flow visualisations have shown that the changes consist of a weakening of the strength of the streamwise vortices. Existing literature reveals no attempts of numerical simulation of this phenomenon. In this paper an approach is pursued by using a constitutive equation which relates the elongation viscosity to the local properties of the flow. According to this model this viscosity is large in zones where the amount of strain rate is greater than the amount of vorticity, and is zero when the vorticity exceeds the strain rate. Simulations have been performed in a “minimal channel” to give good resolution with a limited number of grid points. The accuracy of the method is tested by comparison with the results of other techniques. For simulations with polymers, quantitative comparisons cannot be made, but the results reproduce the qualitative outputs of the experiments. The mean streamwise velocity is modified in the buffer layer; the peak of the streamwise turbulent intensity, in wall units, increases and its maximum moves far from the wall; and the vertical turbulent intensity is largely reduced in the wall region. An interesting outcome from both the simulation and the experiments is that the strength of the longitudinal vortices is reduced when the polymers are present.     

Meshless simulation of equilibrium swelling/deswelling of pH‐sensitive hydrogels

Hydrogels have been widely used in microelectromechanical systems (MEMS) and Bio‐MEMS devices. In this article, the equilibrium swelling/deswelling of the pH‐stimulus cylindrical hydrogel in the microchannel is studied and simulated by the meshless method. The multi‐field coupling model, called multi‐effect‐coupling pH‐stimulus (MECpH) model, is presented and used to describe the chemical field, electric field, and the mechanical field involved in the problem. The partial differential equations (PDEs) describing these three fields are either nonlinear or coupled together. This multi‐field coupling and high nonlinear characteristics produce difficulties for the conventional numerical methods (e.g., the finite element method or the finite difference method), so an alternative—meshless method is developed to discretize the PDEs, and the efficient iteration technique is adopted to solve the nonlinear problem. The computational results for the swelling/deswelling diameter of the hydrogel under the different pH values are firstly compared with experimental results, and they have a good agreement. The influences of other parameters on the mechanical properties of the hydrogel are also investigated in detail. It is shown that the multi‐field coupling model and the developed meshless method are efficient, stable, and accurate for simulation of the properties of the stimuli‐sensitive hydrogel. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 326–337, 2006     

湍流燃烧二阶矩模型的大涡模拟验证

对美国Sandia国家实验室测量的甲烷-空气湍流射流火焰进行了三维大涡模拟(LES),其统计的时均量以及脉动量符合实验测量值,且和RANS代数二阶矩湍流燃烧模型模拟结果大致符合。LES统计得到的化学反应率系数脉动和甲烷浓度脉动的互关联与对应的时均量梯度的乘积分布趋势相同,大小幅值的位置也相对应,与二阶矩湍流燃烧模型的封闭假设相符,验证了二阶矩代数燃烧模型的合理性。     

Planned positron experiments at FRM-II

The new research reactor FRM-II near Munich has a strong positron source, which delivers an intense, nearly monoenergetic positron beam. Our positron systems, the pulsed low energy positron source (PLEPS) and the scanning positron microscope (SPM) will be operated at this beam. Some aspects of matching these systems to the new positron source will be discussed.Considerable improvements are expected, e.g. more than 10⁵ s⁻¹ recorded events at PLEPS and sub-micrometre resolution at SPM. They will enable investigations in so far inaccessible problems like the evaluation of annihilation characteristics and trapping constants of individual defects or studies of fast dynamical processes. In applied materials science complex defect structures will be studied which demand a resolution into many differing lifetimes, e.g. fractured specimens, wear, corrosion, etc. Also large series of measurements at small systematic modifications are planned. There is also the opportunity to analyse in addition the chemical microstructure of the specimens by means of a hydrogen microprobe and other ion beam techniques available close to FRM-II at the Technical University of Munich.     

Recognition of ring pairs in the data analysis of Cherenkov detectors

We present a new algorithm for the identification of ring pairs in Cherenkov detectors. The approach combines linear regression methods for the construction of an appropriate reference quantity, methods from statistics for the evaluation of hypotheses and Monte Carlo simulation runs for the tuning of threshold levels.     

Antimicorbial effects of newly synthesized organotin(IV) and organolead(IV) derivatives

Triorganotin(IV) and triorganolead(IV) derivatives of the types Me₃Sn(SCZ) and Ph₃Pb(SCZ) (where SCZ^? is the anion of a semicarbanzone ligand) have been synthesized by substitution reactions of trimethyltin chloride and triphenyl-lead chloride with semicarbazones derived from heterocyclic ketones. The resulting complexes have been characterized by elemental analyses, molecualr weight determinations and conductivity measurements. The mode of bonding has been established on the basis of infrared and ¹H, ¹³C and ¹¹⁹Sn NMR spectroscopic studies. Some respresentative complexes have also been evaluated for their antimicrobial effects on different species of pathogenic fungi and bacteria; the results of these investigations have been reported in the present paper.     

An evolutionary method for optimization of plate buckling resistance

Optimization of plate buckling resistance is very complicated, because the in-plane stress resultants in the prebuckled state of a plate are functions of thickness distribution. This paper discusses the problem of finding the optimum thickness distribution of isotropic plate structures, with a given volume and layout, that would maximise the buckling load. A simple numerical method using the finite-element analysis is presented to obtain the optimum thickness distribution. Optimum designs of compression-loaded rectangular plates with different boundary conditions and plate aspect ratios are obtained by using the proposed method. Optimum designs from earlier studies and the methods of buckling analysis used to attain these results are discussed and compared with the designs from the proposed method. This paper also examines the reliability of the optimality criterion generally used for plate buckling optimization, which is based on the uniform strain energy density.     

Surface properties of finite classical Coulomb systems: Debye-Hückel approximation and computer simulations

We report analytical and numerical studies of surface correlations in finite, homogeneously polarizable, classical Coulomb systems placed in an insulating or conducting environment. Their purpose is to understand the phenomenological, shape-dependent laws of electrostatics, from the point of view of statistical mechanics; we focus on the knowledge of the dielectric susceptibility of the system, a quantity proportional to the equilibrium fluctuation of the system's instantaneous polarization per unit volume. This goal has been achieved for a system in a conducting state. The picture is that the shape-dependent part of the susceptibilities results from the action of unbounded observables (the second moments of the instantaneous polarization of the system) on long-range surface correlations and that the relations of electrostatics are verified by means of shape-dependent thermodynamic limits. This picture is supported (i) by exact solutions and asymptotic analysis of the Debye-Hückel approximation of multicomponent plasmas in disks and spheres with insulating and conducting environment and also in ellipses in a vacuum, and (ii) by computer simulations of a one-component plasma in a disk with different environments, notably a conducting environment with permeable and impermeable wall. These observations have revealed for the first time the reason why the susceptibility of a conducting disk in a conductor with impermeable walls diverges linearly with the radius of the disk: this is due to the occurrence of long-range radial correlations in the conductor. These findings are quantitatively interpreted in terms of a novel canonical Debye-Huckel approximation as contrasted to the ordinary grand canonical version. Lastly a fresh look at the problem of the surface correlations of a conductor in a vacuum, which places the observer close to the surface of the conductor but in the vacuum, is presented and applied to the disk, the ellipse, the cylinder, the sphere, and the wedge.     

A method for selection and use of numerical integration points for molecules with cylindrical symmetry

Boys and Handy [1] have discussed the solution of the bivariational equations with restricted numerical integration. One of the weaknesses of the method was that in the numerical summations over points, some points arose with r _ij= 0 and non-zero weights. This makes the method quite impractical for the Schrodinger Hamiltonian (because of the singularity at r _ij= 0), and it cannot be advantageous for the transcorrelated Hamiltonian C^–1HC because there will be some discontinuous higher derivatives at r _ij=0. Here it is shown how the symmetry of cylindrically symmetric molecules can be used to eliminate such points, without losing any of the advantages of the overall method, such as the convergence of the eigensolutions. It is also shown how the primary numerical integration points (z _i, r_i) may be chosen in any calculation such that each is associated with an equal amount of one-electron density. The choice of the angular coordinates are governed by the removal of the r _ij=0 points and maintaining the natural orthogonality between orbitals of different symmetry types. The method has been programmed and found to be practical, although no new molecular calculations have yet been performed. It is to be hoped that these points will give a basis for new transcorrelated calculations on diatomic molecules.This paper was presented during the session on numerical integration methods for molecules of the 1970 Quantum Theory Conference in Nottingham. It has been revised in the light of the interesting discussion which followed.