Analytical solution and numerical model for the interface in a stratified long wave system

For a two layered long wave propagation, linearized governing equations, which were derived earlier from the Euler equations of mass and momentum assuming negligible friction and interfacial mixing are solved analytically using Fourier transform. For the solution, variations of upper layer water level is assumed to be sinosoidal having known amplitude and variations of interface level is solved. As the governing equations are too complex to solve it analytically, density of upper layer fluid is assumed as very close to the density of lower layer fluid to simplify the lower layer equation. A numerical model is developed using the staggered leap-forg scheme for computation of water level and discharge in one dimensional propagation having known amplitude for the variations of upper layer water level and interface level to be solved. For the numerical model, water levels (upper layer and interface) at both the boundaries are assumed to be known from analytical solution. Results of numerical model are verified by comparing with the analytical solutions for different time period. Good agreements between analytical solution and numerical model are found for the stated boundary condition. The reliability of the developed numerical model is discussed, using it for different a (ratio of density of fluid in the upper layer to that in the lower layer) and p (ratio of water depth in the lower layer to that in the upper layer) values. It is found that as ‘CX’ increases amplification of interface also increases for same upper layer amplitude. Again for a constant lower layer depth, as ‘p’ increases amplification of interface. also increases for same upper layer amplitude.

Separation of mercury (II), lead (II), and silver (I) from aqueous solution using an aminopolycarboxylate silica

A quartz crystal microbalance modified by the attachment of silica particles derivatized with the aminopolycarboxylate ligand N-[(3-trimethoxysilyl)propyl]ethylenediamine-N,N',N'-triacetic acid has been employed to assess conditions under which mercury (II), lead (II), and silver (I) nitrates may be separated in aqueous solution. The separation protocol, which involved removal of Hg(II), as [HgI4](2-), and Pb(II) with H+ was successfully applied to a batchwise separation of the 3 metal ions.

A new wavelet-based method for the solution of the population balance equation

A new wavelet-based method for solving population balance equations with simultaneous nucleation, growth and agglomeration is proposed, which uses wavelets to express the functions. The technique is very general, powerful and overcomes the crucial problems of numerical diffusion and stability that often characterize previous techniques in this area. It is also applicable to an arbitrary grid to control resolution and computational efficiency. The proposed technique has been tested for pure agglomeration, simultaneous nucleation and growth, and simultaneous growth and agglomeration. In all cases, the predicted and analytical particle size distributions are in excellent agreement. The presence of moving sharp fronts can be addressed without the prior investigation of the characteristics of the processes. (C) 2001 Published by Elsevier Science Ltd.

Diverse solid-state and solution structures within a series of hexaamine dicopper(II) complexes

Mono- and dicopper(II) complexes of a series of potentially bridging hexaamine ligands have been prepared and characterized in the solid state by X-ray crystallography. The crystal structures of the following Cu-II complexes are reported: [Cu(HL3)](ClO4)(3), C11H31Cl3CuN6O12, monoclinic, P2(1)/n, a = 8.294(2) Angstrom, b = 18.364(3) Angstrom, c = 15.674(3) Angstrom, beta = 94.73(2)degrees, Z = 4; {[Cu-2(L-4)(CO3)](2)}(ClO4)(4). 4H(2)O, C40H100Cl4Cu4N12O26, triclinic, P (1) over bar, a = 9.4888(8) Angstrom, b=13.353(1) Angstrom,. c = 15.329(1) Angstrom, alpha = 111.250(7)degrees, beta = 90.068(8)degrees, gamma = 105.081(8)degrees, Z=1; [Cu-2(L-5)(OH2)(2)](ClO4)(4), C(13)H(36)Cl(4)Cu(2)Z(6)O(18), monoclinic, P2(1)/c, a = 7.225(2) Angstrom. b = 8.5555(5) Angstrom, c = 23.134(8) Angstrom, beta = 92.37(1)degrees, Z = 2; [Cu-2(L-6)(OH2)(2)](ClO4)(4). 3H(2)O, C14H44Cl4Cu2N6O21, monoclinic, P2(1)/a, a = 15.204(5) Angstrom, b = 7.6810(7) Angstrom, c = 29.370(1) Angstrom, beta = 100.42(2)degrees, Z = 4. Solution spectroscopic properties of the bimetallic complexes indicate that significant conformational changes occur upon dissolution, and this has been probed with EPR spectroscopy and molecular mechanics calculations.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Generalized quasi mode theory of macroscopic canonical quantization in cavity quantum electrodynamics and quantum optics I. Theory

The quasi mode theory of macroscopic quantization in quantum optics and cavity QED developed by Dalton, Barnett and Knight is generalized. This generalization allows for cases in which two or more quasi permittivities, along with their associated mode functions, are needed to describe the classical optics device. It brings problems such as reflection and refraction at a dielectric boundary, the linear coupler, and the coupling of two optical cavities within the scope of the theory. For the most part, the results that are obtained here are simple generalizations of those obtained in previous work. However the coupling constants, which are of great importance in applications of the theory, are shown to contain significant additional terms which cannot be 'guessed' from the simpler forms. The expressions for the coupling constants suggest that the critical factor in determining the strength of coupling between a pair of quasi modes is their degree of spatial overlap. In an accompanying paper a fully quantum theoretic derivation of the laws of reflection and refraction at a boundary is given as an illustration of the generalized theory. The quasi mode picture of this process involves the annihilation of a photon travelling in the incident region quasi mode, and the subsequent creation of a photon in either the incident region or transmitted region quasi modes.

Generalized quasi mode theory of macroscopic canonical quantization in cavity quantum electrodynamics and quantum optics II. Application to reflection and refraction at a dielectric interface

The generalization of the quasi mode theory of macroscopic quantization in quantum optics and cavity QED presented in the previous paper, is applied to provide a fully quantum theoretic derivation of the laws of reflection and refraction at a boundary. The quasi mode picture of this process involves the annihilation of a photon travelling in the incident region quasi mode, and the subsequent creation of a photon in either the incident region or transmitted region quasi modes. The derivation of the laws of reflection and refraction is achieved through the dual application of the quasi mode theory and a quantum scattering theory based on the Heisenberg picture. Formal expressions from scattering theory are given for the reflection and transmission coefficients. The behaviour of the intensity for a localized one photon wave packet coming in at time minus infinity from the incident direction is examined and it is shown that at time plus infinity, the light intensity is only significant where the classical laws of reflection and refraction predict. The occurrence of both refraction and reflection is dependent upon the quasi mode theory coupling constants between incident and transmitted region quasi modes being nonzero, and it is seen that the contributions to such coupling constants come from the overlap of the mode functions in the boundary layer region, as might be expected from a microscopic theory.

Solution structures by H-1 NMR of the novel cyclic trypsin inhibitor SFTI-1 from sunflower seeds and an acyclic permutant

SFTI-1 is a recently discovered cyclic peptide trypsin inhibitor from sunflower seeds comprising 14 amino acid residues. It is the most potent known Bowman-Birk inhibitor and the only naturally occurring cyclic one. The solution structure of SFTI-1 has been determined by H-1-NMR spectroscopy and compared with a synthetic acyclic permutant. The solution structures of both are remarkably similar. The lowest energy structures from each family of 20 structures of cyclic and acyclic SFTI-1 have an rmsd over the backbone and heavy atoms of 0.29 Angstrom and 0.66 Angstrom, respectively. The structures consist of two short antiparallel beta -strands joined by an extended loop containing the active site at one end. Cyclic SFTI-1 also has a hairpin turn completing the cycle. Both molecules contain particularly stable arrangements of cross-linking hydrogen bonds between the beta -strands and a single disulfide bridge, making them rigid and well defined in solution. These stable arrangements allow both the cyclic and acyclic variants of SFTI-1 to inhibit trypsin with very high potencies (0.5 nM and 12.1 nM, respectively). The cyclic nature of SFTI-1 appears to have evolved to provide higher trypsin inhibition as well as higher stability. The solution structures are similar to the crystal structure of the cyclic inhibitor in complex with trypsin. The lack of a major conformational change upon binding suggests that the structure of SFTI-1 is rigid and already pre-organized for maximal binding due to minimization of entropic losses compared to a more flexible ligand. These properties make SFTI-1 an ideal platform for the design of small peptidic pharmaceuticals or pesticides. (C) 2001 Academic Press.

A critical comparison of the external and internal boron requirements for contrasting species in boron-buffered solution culture

Despite reports that boron (B) requirements differ among plant species there is a shortage of critical evidence to demonstrate unequivocally whether species differ in internal or external B requirements or both. The present research was conducted to establish the external and internal B requirements of three contrasting species, a woody dicot (marri), an herbaceous dicot (sunflower) and a monocot (wheat) using B-buffered solution culture. Boron-buffered solution culture provided satisfactory control of external B concentrations ranging from 0.04 to 30 muM throughout the 20- (sunflower and wheat) or 40-day (marri) growth period. At low external B concentrations (less than or equal to 0.13 muM), the growth of marri and sunflower was severely depressed but by contrast the vegetative growth of wheat plants was satisfactory and free of B deficiency symptoms. Marri and sunflower plants achieved total maximum shoot growth at greater than or equal to1.2 muM B in solutions while wheat plants did so at greater than or equal to 0.6 muM B. The critical B concentrations (mg kg(-1) dry matter) in the youngest open leaf blades of marri, sunflower and wheat plants were 17.9, 19.7 and 1.2 on 20, 10 and 10 days after transplanting (DAT), respectively. Lower internal and external B requirements of wheat were matched by a lower uptake rate of B compared to marri and sunflower.

Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.

Solution structure of BSTI: A new trypsin inhibitor from skin secretions of Bombina bombina

The three-dimensional solution structure of BSTI, a trypsin inhibitor from the European frog Bombina bombina, has been solved using H-1 NMR spectroscopy. The 60 amino acid protein contains five disulfide bonds, which were unambiguously determined to be Cvs (4-38), Cys (13-34), Cys (17-30), Cys (21-60), and Cys (40-54) by experimental restraints and subsequent structure calculations. The main elements of secondary structure are four beta -strands, arranged as two small antiparallel beta -sheets, The overall fold of BSTI is disk shaped and is characterized by the lack of a hydrophobic core. The presumed active site is located on a loop comprising residues 21-34, which is a relatively disordered region similar to that seen in many other protease inhibitors. However, the overall fold is different to other known protease inhibitors with the exception of a small family of inhibitors isolated from nematodes of the family Ascaris and recently also from the haemolymph of Apis mellifera. BSTI may thus be classified as a new member of this recently discovered family of protease inhibitors.

The solution structure of C1-T1, a two-domain proteinase inhibitor derived from a circular precursor protein from Nicotiana alata

A two-domain portion of the proteinase inhibitor precursor from Nicotiana alata (NaProPI) has been expressed and its structure determined by NMR spectroscopy. NaProPI contains six almost identical 53 amino acid repeats that fold into six highly similar domains; however, the sequence repeats do nut coincide with the structural domains. Five of the structural domains comprise the C-terminal portion of one repeat and the N-terminal portion of the next. The sixth domain contains the C-terminal portion of the sixth repeat and the N-terminal portion of the first repeat. Disulphide bonds link these C and N-terminal fragments to generate the clasped-bracelet fold of NaProPI. The three-dimensional structure of NaProPI is not known, but it is conceivable that adjacent domains in NaProPI interact to generate the circular bracelet with the N and C termini in close enough proximity to facilitate formation of the disulphide bonds that form the clasp The expressed protein, examined in the current study, comprises residues 25-135 of NaProPI and encompasses the first two contiguous structural domains, namely the chymotrypsin inhibitor C1 and the trypsin inhibitor T1, joined by a five-residue linker, and is referred to as C1-T1. The tertiary structure of each domain in C1-T1 is identical to that found in the isolated inhibitors. However, no nuclear Overhauser effect contacts are observed between the two domains and the five-residue linker adopts an extended conformation. The absence of interactions between the domains indicates that adjacent domains do not specifically interact to drive the circularisation of NaProPI. These results are in agreement with recent data which describe similar PI precursors from other members of the Solanaceae having two, three, or four repeats. The lack of strong interdomain association is likely to be important for the function of individual inhibitors by ensuring that there is no masking of reactive sites upon release from the precursor. (C) 2001 Academic Press.

Application of heterogeneous vacancy solution theory to characterization of microporous solids

The vacancy solution theory of adsorption is re-formulated here through the mass-action law, and placed in a convenient framework permitting the development of thermodynamic ally consistent isotherms. It is shown that both the multisite Langmuir model and the classical vacancy solution theory expression are special cases of the more general approach when the Flory-Huggins activity coefficient model is used, with the former being the thermodynamically consistent result. The improved vacancy solution theory approach is further extended here to heterogeneous adsorbents by considering the pore-width dependent potential along with a pore size distribution. However, application of the model to numerous hydrocarbons as well as other adsorptives on microporous activated carbons shows that the multisite model has difficulty in the presence of a pore size distribution, because pores of different sizes can have different numbers of adsorbed layers and therefore different site occupancies. On the other hand, use of the classical vacancy solution theory expression for the local isotherm leads to good simultaneous fit of the data, while yielding a site diameter of about 0.257 nm, consistent with that expected for the potential well in aromatic rings on carbon pore surfaces. It is argued that the classical approach is successful because the Flory-Huggins term effectively represents adsorbate interactions in disguise. When used together with the ideal adsorbed solution theory the heterogeneous vacancy solution theory successfully predicts binary adsorption equilibria, and is found to perform better than the multisite Langmuir as well as the heterogeneous Langmuir model. (C) 2001 Elsevier Science Ltd. All rights reserved.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.