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.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

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.

A new equation for macroscopic description of capillary rise in porous media

Capillary rise in porous media is frequently modeled using the Washburn equation. Recent accurate measurements of advancing fronts clearly illustrate its failure to describe the phenomenon in the long term. The observed underprediction of the position of the front is due to the neglect of dynamic saturation gradients implicit in the formulation of the Washburn equation. We consider an approximate solution of the governing macroscopic equation, which retains these gradients, and derive new analytical formulae for the position of the advancing front, its speed of propagation, and the cumulative uptake. The new solution properly describes the capillary rise in the long term, while the Washburn equation may be recovered as a special case. (C) 2004 Elsevier Inc. All rights reserved.

Entropy generation for forced convection in a porous saturated circular tube with uniform wall temperature

A numerical study is reported to investigate both the First and the Second Law of Thermodynamics for thermally developing forced convection in a circular tube filled by a saturated porous medium, with uniform wall temperature, and with the effects of viscous dissipation included. A theoretical analysis is also presented to study the problem for the asymptotic region applying the perturbation solution of the Brinkman momentum equation reported by Hooman and Kani [1]. Expressions are reported for the temperature profile, the Nusselt number, the Bejan number, and the dimensionless entropy generation rate in the asymptotic region. Numerical results are found to be in good agreement with theoretical counterparts.

On a semilinear Schrodinger equation with critical Sobolev exponent

We consider the semilinear Schrodinger equation -Deltau+V(x)u= K(x) \u \ (2*-2 u) + g(x; u), u is an element of W-1,W-2 (R-N), where N greater than or equal to4, V, K, g are periodic in x(j) for 1 less than or equal toj less than or equal toN, K>0, g is of subcritical growth and 0 is in a gap of the spectrum of -Delta +V. We show that under suitable hypotheses this equation has a solution u not equal 0. In particular, such a solution exists if K equivalent to 1 and g equivalent to 0.

Analysis of heat transfer and entropy generation for a thermally developing Brinkman–Brinkman forced convection problem in a rectangular duct with isoflux walls

Heat transfer and entropy generation analysis of the thermally developing forced convection in a porous-saturated duct of rectangular cross-section, with walls maintained at a constant and uniform heat flux, is investigated based on the Brinkman flow model. The classical Galerkin method is used to obtain the fully developed velocity distribution. To solve the thermal energy equation, with the effects of viscous dissipation being included, the Extended Weighted Residuals Method (EWRM) is applied. The local (three dimensional) temperature field is solved by utilizing the Green’s function solution based on the EWRM where symbolic algebra is being used for convenience in presentation. Following the computation of the temperature field, expressions are presented for the local Nusselt number and the bulk temperature as a function of the dimensionless longitudinal coordinate, the aspect ratio, the Darcy number, the viscosity ratio, and the Brinkman number. With the velocity and temperature field being determined, the Second Law (of Thermodynamics) aspect of the problem is also investigated. Approximate closed form solutions are also presented for two limiting cases of MDa values. It is observed that decreasing the aspect ratio and MDa values increases the entropy generation rate.

Refining the risk-benefit equation for thrombolysis: How to identity the low risk patient before administering thrombolytic therapy

In view of the relative risk of intracranial haemorrhage and major bleeding with thrombolytic therapy, it is important ro identify as early as possible the low risk patient who may not have a net clinical benefit from thrombolysis in the setting of acute myocardial infarction. An analysis of 5434 hospital-treated patients with myocardial infarction in the Perth MONICA study showed that age below 60 and absence of previous infarction or diabetes, shock, pulmonary oedema, cardiac arrest and Q-wave or left bundle branch block on the initial ECG identified a large group of patients with a 28 day mortality of only 1%, and one year mortality of only 2%. Identification of baseline risk in this way helps refine the risk-benefit equation for thrombolytic therapy, and may help avoid unnecessary use of thrombolysis in those unlikely to benefit.

A pore-size-dependent equation of state for multilayer adsorption in cylindrical mesopores

The chemical potential of adsorbed film inside cylindrical mesopores is dependent on the attractive interactions between the adsorbed molecules and adsorbent, the curvature of gas/adsorbed phase interface, and surface tension. A state equation of the adsorbed film is proposed to take into account the above factors. Nitrogen adsorption on model adsorbents, MCM-41, which exhibit uniform cylindrical channels, are used to verify the theoretical analysis. The proposed theory is capable of describing the important features of adsorption processes in cylindrical mesopores. According to this theory, at a given relative pressure, the smaller the pore radius is, the thicker the adsorbed film will be. The thickening of adsorbed films in the pores as the vapor pressure increases inevitably causes an increase in the interface curvature, which consequently leads to capillary condensation. Besides, this study confirmed that the interface tension depends substantially on the interface curvature in small mesopores. A quantitative relationship between the condensation pressure and the pore radius can be derived from the state equation and used to predict the pore radius from a condensation pressure, or vice versa.

Groundwater waves in a coastal aquifer: A new governing equation including vertical effects and capillarity

Groundwater waves, that is, water table fluctuations, are a natural phenomenon in coastal aquifers. They represent an important part of the interaction between the ocean and aquifer and affect the mass exchange between them. This paper presents a new groundwater wave equation. Because it includes the effects of vertical flows and capillarity, the new equation is applicable to both intermediate-depth aquifers and high-frequency waves. Compared with the wave equation derived by Nielsen ed al. [1997], the present equation provides a closer representation of groundwater waves. In particular, it predicts high-frequency water table fluctuations as observed in the field. A validation of the new equation has been carried out by comparing the analytical solutions to it with predictions from direct simulations using the numerical model SUTRA. The effects of various physical parameters and their relative importance are also discussed.

Asymptotic solutions to the Gross-Pitaevskii gain equation: Growth of a Bose-Einstein condensate

We give an asymptotic analytic solution for the generic atom-laser system with gain in a D-dimensional trap, and show that this has a non-Thomas-Fermi behavior. The effect is due to Bose-enhanced condensate growth, which creates a local-density maximum and a corresponding outward momentum component. In addition, the solution predicts amplified center-of-mass oscillations, leading to enhanced center-of-mass temperature.

Transfer matrix eigenvalues of the anisotropic multiparametric U model

A multiparametric extension of the anisotropic U model is discussed which maintains integrability. The R-matrix solving the Yang-Baxter equation is obtained through a twisting construction applied to the underlying U-q(sl (2/1)) superalgebraic structure which introduces the additional free parameters that arise in the model. Three forms of Bethe ansatz solution for the transfer matrix eigenvalues are given which we show to be equivalent.

Similitude applied to centrifugal scaling of unsaturated flow

Centrifuge experiments modeling single-phase flow in prototype porous media typically use the same porous medium and permeant. Then, well-known scaling laws are used to transfer the results to the prototype. More general scaling laws that relax these restrictions are presented. For permeants that are immiscible with an accompanying gas phase, model-prototype (i.e., centrifuge model experiment-target system) scaling is demonstrated. Scaling is shown to be feasible for Miller-similar (or geometrically similar) media. Scalings are presented for a more, general class, Lisle-similar media, based on the equivalence mapping of Richards' equation onto itself. Whereas model-prototype scaling of Miller-similar media can be realized easily for arbitrary boundary conditions, Lisle-similarity in a finite length medium generally, but not always, involves a mapping to a moving boundary problem. An exception occurs for redistribution in Lisle-similar porous media, which is shown to map to spatially fixed boundary conditions. Complete model-prototype scalings for this example are derived.