Vibrational control of crystal growth from liquid phase

Modeling and numerical simulations of the convective flows induced by the vibration of the monocrystal during crystal growth have been performed for two configurations simulating the Cz and FZ methods. This permitted to emphasize the role of different vibrational mechanisms in the formation of the average flows. It is shown that an appropriate combination of these mechanisms can be used to counteract the usual convective flows (buoyancy- and/or thermocapillary-driven) inherent to crystal growth processes from the liquid phase. While vibrational convection is rather complex due to these identified mechanisms, the new modeling used in the present paper opens up very promising perspectives to efficiently control heat and mass transfer during real industrial applications of crystal growth from the liquid phase.     

Solution of the hyperbolic mild‐slope equation using the finite volume method

A finite volume solver for the 2D depth‐integrated harmonic hyperbolic formulation of the mild‐slope equation for wave propagation is presented and discussed. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov‐type second‐order finite volume scheme, whereby the numerical fluxes are computed using Roe's flux function. The eigensystem of the mild‐slope equations is derived and used for the construction of Roe's matrix. A formulation that updates the unknown variables in time implicitly is presented, which produces a more accurate and reliable scheme than hitherto available. Boundary conditions for different types of boundaries are also derived. The agreement of the computed results with analytical results for a range of wave propagation/transformation problems is very good, and the model is found to be virtually paraxiality‐free. Copyright © 2003 John Wiley & Sons, Ltd.     

Evaluation of a complete HPLC solvent optimization system involving piece-wise quadratic modelling

Summary Interpretive methods are accepted to give the best possible results for selectivity optimization in HPLC. However the methods are very complex, and most work so far has been detailed academic studies. This paper describes an evaluation of a complete integrated system incorporating peak labelling, modelling of retention behaviour and calculation of response surfaces, with particular emphasis on the retention modelling. The peak labelling section has been discussed previously.A piece-wise quadratic function is investigated for the modelling of retention times across an isoeluotropic plane to effect selectivity optimization in HPLC. This requires 10 data-points on the isoeluotropic plane. The predicted global optimum and local optima are evaluated by comparison of calculated and experimental retention data, for a nine component sample. Seven interstitial points, distributed across the whole plane between the data-points, are similarly evaluated for a related sample. The typical error (in retention time) is less than 2%, often 1%, and the maximum error is 4.2%. At the global optimum the error was found be less than 1.3% for all 9 peaks.     

Pervaporative transport modelling in a ternary system: ethyltertiarybutylether/ethanol/polyurethaneimide

The study deals with the analysis of diffusion and mass transfer modelling during pervaporation in a true ternary system involving a polar liquid mixture (ETBE/EtOH) and a polar block copolymer (polyurethaneimide or PUI). A survey of methods of pervaporative transfer modelling in ternary systems is first developed. From differential permeation experiments carried out with both pure liquids, it appears that both permeants obey a Fickian law. Moreover, the diffusional behaviour is consistent with Long's model, which has thus been assumed for the related ternary system. An extension of the Brun's model is then derived, which takes into account the diffusion coupling as well the significant deviation from sorption ideality. From a practical point of view, the calculated values of fluxes show generally good agreement with the experimental results, although a small deviation occurs for mixtures of low ethanol content. Diffusion coefficients of both pure solvents corresponding to transient or steady state are compared. A very good agreement is found for the aprotic permeant (ETBE). whereas the diffusion coefficient of ethanol in transient state is only the quarter of the value corresponding to steady state. The results are discussed in comparison with related investigations in the literature, involving specific liquid-polymer interactions.     

Kinetic modelling of cadmium removal from phosphoric acid by non-dispersive solvent extraction

In this work it is reported that the kinetic modelling of the separation of cadmium from phosphoric acid by non-dispersive solvent extraction. Using Cyanex 302 as selective extractant, the extraction step was carried out in a hollow fibre module containing polypropylene fibres, whereas the concentration step required a ceramic module with tubular channels due to the high acidity of the backextraction agent. Application of the methodology previously reported by the authors led to the development of a kinetic model with three design parameters, i.e., equilibrium constant of the extraction reaction (K'_e = 6 × 10³ mol⁻²/l⁻²), membrane mass transport coefficient of the extraction module K_me=8.33×10⁻⁸ m/s) and of the backextraction module (K_ms=3.33×10⁻⁸ m/s), that described satisfactorily the behaviour of the separation-concentration system. Thus, in this work a new application of the non-dispersive solvent extraction technology is presented, characterising at the same time the behaviour and parameters of a new type of contactor, i.e., a tubular ceramic module.     

Computer simulation of the steady state currents at enzyme doped carbon paste electrodes

The plate-gap model of porous enzyme doped electrode has been proposed and analyzed. It was suggested that reaction diffusion conditions in pores of bulk electrode resemble particular conditions in thin gap between parallel conducting plates. The model is based on the diffusion equations containing a nonlinear term related to the Michaelis–Menten kinetic of the enzymatic reaction inside gap. Steady state current was calculated for the wide range of given parameters and substrate concentrations. All dependences of current on substrate concentration were approximated by hyperbolas in order to obtain “apparent” parameters (maximal currents and apparent Michaelis constants) of modelled biosensors. Simple approximate relationships between given and apparent parameters were derived. The applicability of theoretical plate-gap model was tested for the case of carbon paste electrodes which were doped with PQQ – dependent glucose dehydrogenase. It was found, that soluble glucose dehydrogenase based biosensors exhibit characteristic features of the theoretical plate-gap biosensors.     

Functionality map analysis of the active site cleft of human thrombin

Summary The Multiple Copy Simultaneous Search methodology has been used to construct functionality maps for an extended region of human thrombin, including the active site. This method allows the determination of energetically favorable positions and orientations for functional groups defined by the user on the three-dimensional surface of a protein. The positions of 10 functional group sites are compared with those of corresponding groups of four thrombin-inhibitor complexes. Many, but not all features, of known thrombin inhibitors are reproduced by the method. The results indicate that certain aspects of the binding modes of these inhibitors are not optimal. In addition, suggestions are made for improving binding by interaction with functional group sites on the thrombin surface that are not used by the thrombin inhibitors. Abbreviations: MCSS, multiple copy simultaneous search; PPACK, d-phenylalanyl-l-propyl-l-arginine chloromethane; NAPAP, N -(2-naphthylsulfonylglycyl)-d-para-amidinophenylalanylpiperidine; argatroban, (2R,4R)-4-methyl-1-[N -(3-methyl-1,2,3,4-tetrahydro-8-quinolinylsulfonyl)-l-arginyl]-2-piperidine carboxylic acid; rms, root mean square. The thrombin residues are numbered according to the chymotrypsin-based numbering by Bode et al. [8]. P1, P2, P3, etc., denote the peptide inhibitor residues on the amino-terminal side of the scissile peptide bond, and S1, S2, S3, etc., the corresponding subsites of thrombin     

A molecular mechanics approach for analyzing tensile nonlinear deformation behavior of single-walled carbon nanotubes

In this paper, by capturing the atomic information and reflecting the behaviour governed by the nonlinear potential function, an analytical molecular mechanics approach is proposed. A constitutive relation for single-walled carbon nanotubes (SWCNT’s) is established to describe the nonlinear stress-strain curve of SWCNT’s and to predict both the elastic properties and breaking strain of SWCNT’s during tensile deformation. An analysis based on the virtual internal bond (VIB) model proposed by P. Zhang et al. is also presented for comparison. The results indicate that the proposed molecular mechanics approach is indeed an acceptable analytical method for analyzing the mechanical behavior of SWCNT’s. The project supported by the National Natural Science Foundation of China (10121202, 90305015 and 10328203), the Key Grant Project of Chinese Ministry of Education (0306) and the Research Grants Council of the Hong Kong Special Administrative Region, China (HKU 7195/04E).     

Nonlinear Modeling Study of Aerodynamic Characteristics of an X38-like Vehicle at Strong Viscous Interaction Regions

Strong viscous interaction and multiple flow regimes exist when vehicles fly at high altitude and high Mach number conditions. The Navier–Stokes(NS) solver is no longer applicable in the above situation. Instead, the direct simulation Monte Carlo (DSMC) method or Boltzmann model equation solvers are usually needed. However, they are computationally more expensive than the NS solver. Therefore, it is of great engineering value to establish the aerodynamic prediction model of vehicles at high altitude and high Mach number conditions. In this paper, the hypersonic aerodynamic characteristics of an X38-like vehicle in typical conditions from 70 km to 110 km are simulated using the unified gas kinetic scheme (UGKS), which is applicable for all flow regimes. The contributions of pressure and viscous stress on the force coefficients are analyzed. The viscous interaction parameters, Mach number, and angle of attack are used as independent variables, and the difference between the force coefficients calculated by UGKS and the Euler solver is used as a dependent variable to establish a nonlinear viscous interaction model between them in the range of 70–110 km. The evaluation of the model is completed using the correlation coefficient and the relative orthogonal distance. The conventional viscous interaction effect and rarefied effect are both taken into account in the model. The model can be used to quickly obtain the hypersonic aerodynamic characteristics of X38-like vehicle in a wide range, which is meaningful for engineering design.