Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

Absorbing boundary conditions for reaction-diffusion equations

We treat here of the question of absorbing boundary conditionsfor nonlinear diffusion equations. We use the conditions designedfor the linear equation, we prove them to be well posed forthe nonlinear problem, and through numerical experiments thatthey are well suited for reaction–diffusion equations.     

Statistical thermodynamics evaluation of polymer–polymer miscibility

The Simha and Somcynsky (S–S) statistical thermodynamics theory was used to compute the solubility parameters as a function of temperature and pressure [δ = δ(T, P)], for a series of polymer melts. The characteristic scaling parameters required for this task, P*, T*, and V*, were extracted from the pressure–temperature–volume (PVT) data. To determine the potential polymer–polymer miscibility, the dependence of δ versus T (at ambient pressure) was computed for 17 polymers. Close proximity of the δ versus T curves for four miscible polymer pairs: PPE/PS, PS/PVME, and PC/PMMA signaled the usefulness of this approach. It is noteworthy, that the tabulated solubility parameters (derived from the solution data under ambient conditions) propounded the immiscibility of the PVC/PVAc pair. The computed values of δ also suggested miscibility for polymer pairs of unknown miscibility, namely PPE/PVC, PPE/PVAc, and PET/PSF. In recognizing the limitations of the solubility parameter approach (the omission of several thermodynamic contributions), these preliminary results are auspicious because they indicate a new route for estimating the miscibility of any polymeric material at a given temperature and pressure. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2909–2915, 2004     

Selective quenching of tryptophanyl fluorescence in bovine serum albumin by the iodide ion

Modified Stern-Volmer equation is obeyed by bovine serum albumin (BSA)-iodide system showing selective quenching of tryptophanyl fluorescence of BSA. The fraction of accessible protein fluorescence is 0.56 and the effective Stern-Volmer constant is 290 M^-1 at pH 7.4 in 0.005 M phosphate buffer at 25°C. Collisional quenching is operative both in the BSA -I⁻¹ system and the model system, tryptophan-I⁻¹. It is supported by the observed relationship between the ratio of quenching rate constants (k _q ) and diffusion coefficients and alsok _q with bulk viscosity.     

Needle Variations and Almost Lower Semicontinuous Differential Inclusions

It is shown that the construction of needle variations can be carried out for almost lower semicontinuous differential inclusions rather than for the case of ordinary single-valued continuously differentiable vector fields usually considered in the literature. The construction leads to needle variations whose flows are in general set-valued but still have good differentiability properties. The variations are constructed by using single-valued selections that are not necessarily continuous with respect to the state variable, but have instead a much weaker 'integral continuity' property, somewhat more general that the 'directional continuity' considered in previous work by A. Cambini and S. Querci, A. Pucci, and A. Bressan. The existence of many such selections is proved by slightly adapting an argument due to Bressan, extending it from the lower semicontinuous to the almost lower semicontinuous case, and strengthening it to yield not only directional continuity at all points but also full continuity at a specified point.     

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.     

Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

Observation of hydrogen induced intermediate borides in PrFeB based alloys by Mössbauer effect spectroscopy

In the present work, a quantitative analysis of the phase compositions by Mössbauer effect spectroscopy of solid and conventional hydrogen disproportionated Pr_13.7Fe_80.3B_6.0 and Pr_13.7Fe_63.5Co_16.7Zr_0.1B_6.0 alloys was carried out. Significant amounts of intermediate borides t-Fe₃B and Pr(Fe, Co)₁₂B₆ were detected after solid hydrogen disproportionation treatment in Pr_13.7Fe_80.3B_6.0 and Pr_13.7Fe_63.5Co_16.7Zr_0.1B_6.0 alloys, respectively. After conventional hydrogenation–disproportionation–desorption–recombination treatment these phases were not detected and in no case residual Pr₂Fe₁₄B-phase was found. It was observed that the amount of intermediate borides after disproportionation can be correlated with the degree of texture after recombination at various temperatures.     

Properties of polystyrene-poly(methylmethacrylate) random and diblock copolymers in dilute and semidilute solutions

Solution properties for random and diblock copolymers of polystyrene (PS) and poly(methyl methacrylate) (PMMA) have been measured by dynamic and total intensity light scattering in solvents of differing quality. The results are compared with the corresponding properties for PS and PMMA homopolymers of similar molecular weight, in order to determine if interactions between unlike monomers are significant. The hydrodynamic radius (R_h) and diffusion second virial coefficient (k_d) for the random copolymer are found to be larger than the corresponding values for the homopolymers in a solvent which is near-theta for the two homopolymers, whereas no such effect is observed for the block copolymer. This suggests that most intrachain interactions occur a relatively short distance along the chain backbone. In a mutual good solvent R_h and k_d of the random copolymer are comparable to the average of the values for the homopolymers, indicating that in a good solvent monomer/solvent interactions dominate over monomer/monomer interactions. For an isolated diblock copolymer in a mutual good solvent, there is no evidence that interactions between unlike monomers lead to additional expansion of the entire molecule, as measured by R_h, nor expansion of the individual blocks as probed by light scattering with one block optically masked. However, at low but finite concentration there is evidence (the coefficients of the binary interaction terms in the viscosity and the mutual diffusion coefficient, and the second and third virial coefficients) that a weak ordering effect may exist in block copolymer solutions, far from the conditions where microphase separation occurs. Finally, measurements of ternary polymer-polymer-solvent solutions show no dependence on monomer composition or monomer distribution for the tracer diffusion of probe PS-PMMA copolymers in a PMMA/toluene matrix. This indicate that the frictional interaction is largely unaffected by interactions between unlike monomers. However, there is evidence that the thermodynamic interaction is more unfavorable between a random copolymer and the homopolymer matrix than between a diblock and the matrix. © 1994 John Wiley & Sons, Inc.     

Influence of concentration on the activation energy for diffusion in polymer-solvent systems

The self-diffusion of benzene, toluene, and ethylbenzene in polystyrene have been analyzed using the Vrentas/Duda free-volume diffusion model. Diffusion coefficient predictions suggest an exponential concentration dependence of the activation energy required to overcome attractive forces, E. Without the use of any diffusion data approximating E as zero over the entire concentration range yields self-diffusion coefficient predictions which are in good agreement with experimental data. © 1992 John Wiley & Sons, Inc.