Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

Microstructure and mechanical properties at TiCp/Ni-alloy interfaces in laser-synthesized coatings

Titanium carbide particle (TiCp) reinforced Ni alloy composite coatings were synthesized by laser cladding using a cw 3 kW CO2 laser. Two kinds of coatings were present in terms of TiCp origins, i.e. undissolved and in situ reacted TiCp, respectively. The former came from the TiCp pre-coated on the sample, whereas the latter from in situ reaction between titanium and graphite in the molten pool during laser irradiation. Conventional and high-resolution transmission electron microscope observations showed the epitaxial growth of TiC, the precipitation of CrB, and the chemical reaction between Ti and B elements around phase interfaces of undissolved TiCp. The hardness, H, and elastic modulus, E, were measured by nanoindentation of the matrix near the TiCp interface. For undissolved TiCp, the loading curve revealed pop-in phenomena caused by the plastic deformation of the crack formation or debounding of TiCp from the matrix. As for in situ generated TiCp, no pop-in mark appears. On the other hand, in situ reacted TiCp led to much higher hardness and modulus than that in the case of undissolved TiCp. The coating reinforced by in situ generated TiCp displayed the highest impact wear resistance at both low and high impact conditions, as compared with coatings with undissolved TiCp and without TiCp. The impact wear resistance of the coating reinforced by undissolved TiCp increases at a low impact work but decreases at a high impact work, as compared with the single Ni alloy coating. The degree of wear for the composite coating depends primarily on the debonding removal of TiCp.

Enhanced transduction methods for electrostatically driven MEMS resonators

Electrically addressed silicon bulk acoustic wave microresonators offer high Q solutions for applications in sensing and signal processing. However, the electrically transduced motional signal is often swamped by parasitic feedthrough in hybrid technologies. With the aim of enhancing the ratio of the motional to feedthrough current at nominal operating voltages, this paper benchmarks a variety of drive and detection principles for electrostatically driven square-extensional mode resonators operating in air and in a foundry MEMS process utilizing 2μm gaps. A new detection technique, combining second harmonic capacitive actuation and piezoresistive detection, outperforms previously reported methods utilizing voltages as low as ± 3V in air providing a promising solution for low voltage CMOS-MEMS integration. ©2009 IEEE.

Enhancing parametric sensitivity using mode localization in electrically coupled mems resonators

We use vibration localization as a sensitive means of detecting small perturbations in stiffness in a pair of weakly coupled micromechanical resonators. For the first time, the variation in the eigenstates is studied by electrostatically coupling nearly identical resonators to allow for stronger localization of vibrational energy due to perturbations in stiffness. Eigenstate variations that are orders of magnitude greater than corresponding shifts in resonant frequency for an induced stiffness perturbation are experimentally demonstrated. Such high, voltagetunable parametric sensitivities together with the added advantage of intrinsic common mode rejection pave the way to a new paradigm of mechanical sensing. ©2009 IEEE.

Studies of the Effect of a Coal Concentrator on NO Formation in Swirling Coal Combustion

To develop low-pollution burners, the effect of a coal concentrator on NO formation in swirling coal combustion is studied using both numerical simulation and experiments. The isothermal gas-particle two-phase velocities and particle concentration in a cold model of swirl burners with and without coal concentrators were measured using the phase Doppler particle anemometer (PDPA). A full two-fluid model of reacting gas-particle flows and coal combustion with an algebraic unified second-order moment (AUSM) turbulence-chemistry model for the turbulent reaction rate of NO formation are used to simulate swirling coal combustion and NO formation with different coal concentrators. The results give the turbulent kinetic energy, particle concentration, temperature and NO concentration in cases of with and without coal concentrators. The predicted results for cold two-phase flows are in good agreement with the PDPA measurement results, showing that the coal concentrator increases the turbulence and particle concentration in the recirculation zone. The combustion modeling results indicate that although the coal concentrator increases the turbulence and combustion temperature, but still can remarkably reduce the NO formation due to creating high coal concentration in the recirculation zone.

Numerical Investigation On Detonation Cell Evolution In A Channel With Area-Changing Cross Section

The two-dimensional cellular detonation propagating in a channel with area-changing cross section was numerically simulated with the dispersion-controlled dissipative scheme and a detailed chemical reaction model. Effects of the flow expansion and compression on the cellular detonation cell were investigated to illustrate the mechanism of the transverse wave development and the cellular detonation cell evolution. By examining gas composition variations behind the leading shock, the chemical reaction rate, the reaction zone length, and thermodynamic parameters, two kinds of the abnormal detonation waves were identified. To explore their development mechanism, chemical reactions, reflected shocks and rarefaction waves were discussed, which interact with each other and affect the cellular detonation in different ways.

The search for exudates from Eurasian watermilfoil and hydrilla

Secondary metabolites are produced by aquatic plants, and in some instances, exudation of these metabolites into the surrounding water has been detected. To determine whether infestations of Eurasian watermilfoil or hydrilla produce such exudates, plant tissues and water samples were collected from laboratory cultures and pond populations and were analyzed using solid phase extraction, HPLC, and various methods of mass spectrometry including electrospray ionization, GC/MS, electron impact and chemical ionization. Previously reported compounds such as tellimagrandin II (from Eurasian watermilfoil) and a caffeic acid ester (from hvdrilla), along with a newly discovered flavonoid, cyanidin 3 dimalonyl glucoside (from hydrilla), were readily detected in plant tissues used in this research but were not detected in any of the water samples. If compounds are being released, as suggested by researchers using axenic cultures, we hypothesize that they may be rapidly degraded by bacteria and therefore undetectable.

Magnitude and extent of chemical contamination and toxicity in sediments of Biscayne Bay and vicinity

The toxicity of sediments in Biscayne Bay and many adjoining tributaries was determined as part of a bioeffects assessments program managed by NOAA’s National Status and Trends Program. The objectives of the survey were to determine: (1) the incidence and degree of toxicity of sediments throughout the study area; (2) the spatial patterns (or gradients) in chemical contamination and toxicity, if any, throughout the study area; (3) the spatial extent of chemical contamination and toxicity; and (4) the statistical relationships between measures of toxicity and concentrations of chemicals in the sediments. The survey was designed to characterize sediment quality throughout the greater Biscayne Bay area. Surficial sediment samples were collected during 1995 and 1996 from 226 randomly-chosen locations throughout nine major regions. Laboratory toxicity tests were performed as indicators of potential ecotoxicological effects in sediments. A battery of tests was performed to generate information from different phases (components) of the sediments. Tests were selected to represent a range in toxicological endpoints from acute to chronic sublethal responses. Toxicological tests were conducted to measure: reduced survival of adult amphipods exposed to solid-phase sediments; impaired fertilization success and abnormal morphological development in gametes and embryos, respectively, of sea urchins exposed to pore waters; reduced metabolic activity of a marine bioluminescent bacteria exposed to organic solvent extracts; induction of a cytochrome P-450 reporter gene system in exposures to solvent extracts; and reduced reproductive success in marine copepods exposed to solid-phase sediments. Contamination and toxicity were most severe in several peripheral canals and tributaries, including the lower Miami River, adjoining the main axis of the bay. In the open basins of the bay, chemical concentrations and toxicity generally were higher in areas north of the Rickenbacker Causeway than south of it. Sediments from the main basins of the bay generally were less toxic than those from the adjoining tributaries and canals. The different toxicity tests, however, indicated differences in severity, incidence, spatial patterns, and spatial extent in toxicity. The most sensitive test among those performed on all samples, a bioassay of normal morphological development of sea urchin embryos, indicated toxicity was pervasive throughout the entire study area. The least sensitive test, an acute bioassay performed with a benthic amphipod, indicated toxicity was restricted to a very small percentage of the area. Both the degree and spatial extent of chemical contamination and toxicity in this study area were similar to or less severe than those observed in many other areas in the U.S. The spatial extent of toxicity in all four tests performed throughout the bay were comparable to the “national averages” calculated by NOAA from previous surveys conducted in a similar manner. Several trace metals occurred in concentrations in excess of those expected in reference sediments. Mixtures of substances, including pesticides, petroleum constituents, trace metals, and ammonia, were associated statistically with the measures of toxicity. Substances most elevated in concentration relative to numerical guidelines and associated with toxicity included polychlorinated biphenyls, DDT pesticides, polynuclear aromatic hydrocarbons, hexachloro cyclohexanes, lead, and mercury. These (and other) substances occurred in concentrations greater than effects-based guidelines in the samples that were most toxic in one or more of the tests. (PDF contains 180 pages)

Mechanics of adhesive contact on a power-law graded elastic half-space

We consider adhesive contact between a rigid sphere of radius R and a graded elastic half-space with Young's modulus varying with depth according to a power law E = E-0(z/c(0))(k) (0 < k < 1) while Poisson's ratio v remaining a constant. Closed-form analytical solutions are established for the critical force, the critical radius of contact area and the critical interfacial stress at pull-off. We highlight that the pull-off force has a simple solution of P-cr= -(k+3)pi R Delta gamma/2 where Delta gamma is the work of adhesion and make further discussions with respect to three interesting limits: the classical JKR solution when k = 0, the Gibson solid when k --> 1 and v = 0.5, and the strength limit in which the interfacial stress reaches the theoretical strength of adhesion at pull-off. (C) 2009 Elsevier Ltd. All rights reserved.

Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.