Hydrodynamic interaction of two unsteady model microorganisms.

The study of pair-wise interactions between swimming microorganisms is fundamental to the understanding of the rheological and transport properties of semi-dilute suspensions. In this paper, the hydrodynamic interaction of two ciliated microorganisms is investigated numerically using a boundary-element method, and the microorganisms are modeled as spherical squirmers that swim by time-dependent surface deformations. The results show that the inclusion of the unsteady terms in the ciliary propulsion model has a large impact on the trajectories of the interacting cells, and causes a significant change in scattering angles with potential important consequences on the diffusion properties of semi-dilute suspensions. Furthermore, the analysis of the shear stress acting on the surface of the microorganisms revealed that the duration and the intensity of the near-field interaction are significantly modified by the presence of unsteadiness. This observation may account for the hydrodynamic nature of randomness in some biological reactions, and supersedes the distinction between intrinsic randomness and hydrodynamic interactions, adding a further element to the understanding and modeling of interacting microorganisms.

Estimating material elasticity by spherical indentation load-relaxation tests on viscoelastic samples of finite thickness

A two-step viscoelastic spherical indentation method is proposed to compensate for 1) material relaxation and 2) sample thickness. In the first step, the indenter is moved at a constant speed and the reaction force is measured. In the second step, the indenter is held at a constant position and the relaxation response of the material is measured. Then the relaxation response is fit with a multi-exponential function which corresponds to a three-branch general Maxwell model. The relaxation modulus is derived by correcting the finite ramp time introduced in the first step. The proposed model takes into account the sample thickness, which is important for applications in which the sample thickness is less than ten times the indenter radius. The model is validated numerically by finite element simulations. Experiments are carried out on a 10% gelatin phantom and a chicken breast sample with the proposed method. The results for both the gelatin phantom and the chicken breast sample agree with the results obtained from a surface wave method. Both the finite element simulations and experimental results show improved elasticity estimations by incorporating the sample thickness into the model. The measured shear elasticities of the 10% gelatin sample are 6.79 and 6.93 kPa by the proposed finite indentation method at sample thickness of 40 and 20 mm, respectively. The elasticity of the same sample is estimated to be 6.53 kPa by the surface wave method. For the chicken breast sample, the shear elasticity is measured to be 4.51 and 5.17 kPa by the proposed indentation method at sample thickness of 40 and 20 mm, respectively. Its elasticity is measured by the surface wave method to be 4.14 kPa. © 2011 IEEE.

Simulation of spherical turbulent premixed flames

Outwardly-propagating spherical turbulent premixed flames are studied using Reynolds-Averaged Navier-Stokes methodology. The reaction rate closure is based on recently developed strained flamelets, where the flamelets are parametrised using the scalar dissipation rate. It is shown that the leading edge speed for the spherical flame is higher than for the planar case, for a given turbulence and thermo-chemical conditions. In addition, it is shown that including the mean curvature effects in the reaction rate closure do not influence the spherical flame speeds when compared with a model that excludes these effects.

Eccentricity effects on acoustic radiation from a spherical source suspended within a thermoviscous fluid sphere.

Acoustic radiation from a spherical source undergoing angularly periodic axisymmetric harmonic surface vibrations while eccentrically suspended within a thermoviscous fluid sphere, which is immersed in a viscous thermally conducting unbounded fluid medium, is analyzed in an exact fashion. The formulation uses the appropriate wave-harmonic field expansions along with the translational addition theorem for spherical wave functions and the relevant boundary conditions to develop a closed-form solution in form of infinite series. The analytical results are illustrated with a numerical example in which the vibrating source is eccentrically positioned within a chemical fluid sphere submerged in water. The modal acoustic radiation impedance load on the source and the radiated far-field pressure are evaluated and discussed for representative values of the parameters characterizing the system. The proposed model can lead to a better understanding of dynamic response of an underwater acoustic lens. It is equally applicable in miniature transducer analysis and design with applications in medical ultrasonics.

The modes of a tethered, spherical balloon

Observations of a tethered meteorological balloon show that a strong vibration coupling exists between axial forcing of the tether and ovalling deformations of the balloon. Such coupling may lead to system instabilities and fatigue failure in a tethered-balloon system. This is particularly relevant in the case of a balloon launched from a moving vessel, as is proposed as part of the SPICE geoengineering project. This paper investigates the vibration characteristics of a tethered, spherical balloon using a simple analytical model: a tensioned, spherical membrane attached to a spring. The analytical solution for the natural frequencies and modeshapes of this system is compared to transfer functions obtained by laser vibrometry. These results are then used to determine the most suitable method of modelling the dynamic response of a tethered balloon.

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

Theoretical study on spherical proton emission

The proton radioactivity half-lives of spherical proton emitters are investigated within a generalized liquid drop model (GLDM), including the proximity effects between nuclei in a neck and the mass and charge asymmetry. The penetrability is calculated in the WKB approximation and the assault frequency is estimated by the quantum mechanism method considering the structure of the parent nucleus. The spectroscopic factor is taken into account in half-life calculation, which is obtained by employing the relativistic mean field (RMF) theory. The half-lives within the GLDM are compared with the experimental data and other theoretical values. The results show that the GLDM works quite well for spherical proton emitters when the assault frequency is estimated by the quantum mechanical method and the spectroscopic factor is considered.

Shell model study of single-particle and collective structure in neutron-rich Cr isotopes

The structure of neutron-rich Cr isotopes is systematically investigated by using the spherical shell model. The calculations reproduce well the known energy levels for the even-even Cr52-62 and odd-mass Cr53-59 nuclei, and predict a lowering of excitation energies around neutron number N = 40. The calculated B(E2; 2(1)(+) -> 0(1)(+)) systematics shows a pronounced collectivity around N = 40; a similar characteristic behavior has been suggested for Zn and Ge isotopes. Causes for the sudden drop of the 9/2(1)(+) energy in Cr-59 and the appearance of very low 0(2)(+) states around N = 40 are discussed. We also predict a new band with strong collectivity built on the 0(2)(+) state in the N = 40 isotope Cr-64.

Unified fission model for proton emission

The unified fission model (UFM) combining with the phenomenological assault frequency has been carried out to investigate the proton-radioactivity half-lives of spherical proton emitters. The results are in good agreement with the experimental data and other theoretical values, and newly observed spherical proton emitters have been analyzed. Finally, the effect of angular momentum transfer on half-life of proton emission has been discussed in detail and a formula can be used to describe this relationship.

Fabrication, characterization of spherical CaWO4:Ln@MCM-41(Ln=Eu3+, Dy3+, Sm3+, Er3+) composites and their applications as drug release systems

Spherical MCM-41 particles with a diameter of about 150 nm have been successfully coated with CaWO4:Ln (Ln = Eu3+, Dy3+, Sm3+, Er3+) phosphor layers through a simple Pechini sol-gel process. The obtained CaWO4:Ln@MCM-41 composites, which kept the mesoporous structure of MCM-41 and the luminescent properties of phosphors, were investigated as a drug delivery system using aspirin (ASPL) as a model drug.

Controlling electric field distribution by graded spherical core-shell metamaterials

This paper investigates analytically the electric field distribution of graded spherical core-shell metamaterials, whose permittivity is given by the graded Drude model. Under the illumination of a uniform incident optical field, the obtained results show that the electrical field distribution in the shell region is controllable and the electric field peak's position inside the spherical shell can be confined in a desired position by varying the frequency of the optical field as well as the parameters of the graded dielectric profiles. It has also offered an intuitive explanation for controlling the local electric field by graded metamaterials.

Thermal conductivity of graded spherical composites with contact resistance

The effective thermal conductivity of graded composites with contact resistance on the inclusion surface is investigated. As an example, we have considered the graded composite media with a spherical particle embedded in a homogeneous matrix, where the thermal conductivity of spherical inclusion is an exponential function k(i) = c exp(betar) (where r is the inside distance of a point in particle from the center of the spherical particle in a spherical coordinate). For both heat contact resistance and perfect contact cases, we have given a reasonable effective medium approximation to calculate the effective conductivity. (C) 2003 Elsevier B.V. All rights reserved.

A Perturbation Scheme for Spherical Arrangements with Application to Molecular Modeling

We describe a software package for computing and manipulating the subdivision of a sphere by a collection of (not necessarily great) circles and for computing the boundary surface of the union of spheres. We present problems that arise in the implementation of the software and the solutions that we have found for them. At the core of the paper is a novel perturbation scheme to overcome degeneracies and precision problems in computing spherical arrangements while using floating point arithmetic. The scheme is relatively simple, it balances between the efficiency of computation and the magnitude of the perturbation, and it performs well in practice. In one O(n) time pass through the data, it perturbs the inputs necessary to insure no potential degeneracies and then passes the perturbed inputs on to the geometric algorithm. We report and discuss experimental results. Our package is a major component in a larger package aimed to support geometric queries on molecular models; it is currently employed by chemists working in "rational drug design." The spherical subdivisions are used to construct a geometric model of a molecule where each sphere represents an atom. We also give an overview of the molecular modeling package and detail additional features and implementation issues.

A Biomechanical Model of Human Oculomotor Plant Kinematics Based Upon Geometric Algebra

A biomechanical model of the human oculomotor plant kinematics in 3-D as a function of muscle length changes is presented. It can represent a range of alternative interpretations of the data as a function of one parameter. The model is free from such deficits as singularities and the nesting of axes found in alternative formulations such as the spherical wrist (Paul, l98l). The equations of motion are defined on a quaternion based representation of eye rotations and are compact atnd computationally efficient.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.