A new approach to potential flow around axisymmetric bodies

A new method is introduced to solve potential flow problems around axisymmetric bodies. The approach relies on expressing the infinite series expansion of the Laplace equation solution in terms of a finite sum which preserves the Laplace solution for the potential function under a Neumann-type boundary condition. Then the coefficients of the finite sum are calculated in a least squares approximation sense using the Gram-Schmidt orthonormalization method. Sample benchmark problems are presented and discussed in some detail. The solutions are accurate and converged faster when a rather small number of terms were used. The method is simple and can be easily programmed.     

A hybrid finite-element–volume-of-fluid method for simulating free surface flows and interfaces

A numerical technique is developed for the simulation of free surface flows and interfaces. This technique combines the strength on the finite element method (FEM) in calculating the field variables for a deforming boundary and the versatility of the volume-of-fluid (VOF) technique in advection of the fluid interfaces. The advantage of the VOF technique is that it allows the simulation of interfaces with large deformations, including surface merging and breaking. However, its disadantage is that is solving the flow equations, it cannot resolve interfaces smaller than the cell size, since information on the subgrid scale is lost. Therefore the accuracy of the interface reconstruction and the treatment of the boundary conditions (i.e. viscous stresses and surface tension forces) become grid-size-dependent. On the other hand, the FEM with deforming interface mesh allows accurate implementation of the boundary conditions, but it cannot handle large surface deformations occurring in breaking and merging of liquid regions. Combining the two methods into a hybrid FEM-VOF method eliminates the major shortcomings of both. The outcome is a technique which can handle large surface deformations with accurate treatment of the boundary conditions. For illustration, two computational examples are presented, namely the instability and break-up of a capillary jet and the coalescence collision of two liquid drops.     

Measurement and prediction of solubilities and diffusion coefficients of carbon dioxide in starch–water mixtures at elevated pressures

The solubility and diffusion coefficient of carbon dioxide in intermediate‐moisture starch–water mixtures were determined both experimentally and theoretically at elevated pressures up to 16 MPa at 50 °C. A high‐pressure decay sorption system was assembled to measure the equilibrium CO₂ mass uptake by the starch–water system. The experimentally measured solubilities accounted for the estimated swollen volume by Sanchez–Lacombe equation of state (S‐L EOS) were found to increase almost linearly with pressure, yielding 4.0 g CO₂/g starch–water system at 16 MPa. Moreover, CO₂ solubilities above 5 MPa displayed a solubility increase, which was not contributed by the water fraction in the starch–water mixture. The solubilities, however, showed no dependence on the degree of gelatinization (DG) of starch. The diffusion coefficient of CO₂ was found to increase with concentration of dissolved CO₂, which is pressure‐dependent, and decrease with increasing DG in the range of 50–100%. A free‐volume‐based diffusion model proposed by Areerat was employed to predict the CO₂ diffusivity in terms of pressure, temperature, and the concentration of dissolved CO₂. S‐L EOS was once more used to determine the specific free volume of the mixture system. The predicted diffusion coefficients showed to correlate well with the measured values for all starch–water mixtures. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 607–621, 2006     

Freezing temperatures and densities of solutions of some hydrocarbons and sodium oleate inN-methylacetamide. Evidence for a solvophobic interaction

The freezing temperatures and densities (at 31°C) of solutions of octane, nonane, decane, 3,3-diethylpentane, and sodium oleate inN-methylacetamide (NMA) have been measured. The molality of the freezing solution was calculated from the density. The solubilities of octane, nonane, and decane inN-methylacetamide are also reported. Apparent molal volumes calculated from the densities are close to the values in the pure hydrocarbons and are not strong functions of the concentration. This indicates the absence of any unusual packing effect. The calculated free energies of transfer of the hydrocarbons from pure hydrocarbon to NMA solution are much less negative than the corresponding values for water, showing that the bulk solvophobic interaction inN-methylacetamide is smaller than in water. This is consistent with the freezing temperatures of sodium oleate which show that micelles do not form below 0.1 mole-kg^–1. The osmotic coefficients of the hydrocarbons calculated from the freezing temperatures showed negative deviations from ideality that were larger for the hydrocarbons with the higher molecular weights. Two estimates of the pairwise solvophobic interaction inN-methylacetamide indicate that it is also smaller than in water. The solvophobic effect in this solvent does not include the large entropy and enthalpy effects found in aqueous solutions.     

Flow injection-spectrophotometric determination of cresol compounds in water by reaction withp-aminophenol

A spectrophotometric method has been developed for the simultaneous determination ofo-cresol andm-cresol in water by reaction withp-aminophenol (PAP). Three different methodologies have been assayed; (i) batch analysis, after reaction in an alkaline medium in the presence of dissolved molecular oxygen as oxidizing agent, (ii) a stopped-flow procedure, carried out in the presence of KIO₄ and (iii) a flow injection method based on the same approach. The batch procedure requires 22 min for the full development of colour witho-cresol and 12 min form-cresol. In the stopped-flow mode, using KIO₄ and a reaction time of 12 min, better sensitivity can be obtained for both compounds and limits of detection of 10 g 1^–1 foro-cresol and 30 g 1^–1 form-cresol were found. The flow injection method has a lower sensitivity but permits more than 80 injections per hour. Based on the different maximum absorbance wavelengths obtained for the reaction products of PAP witho-cresol (614 nm) andm-cresol (632 nm), both compounds can be simultaneously determined in water samples and recoveries of 90 to 115% were found in spiked water samples of different types.     

Generalized solutions of a stochastic partial differential equation

We discuss the Cauchy problem of a certain stochastic parabolic partial differential equation arising in the nonlinear filtering theory, where the initial data and the nonhomogeneous noise term of the equation are given by Schwartz distributions. The generalized (distributional) solution is represented by a partial (conditional) generalized expectation ofT(t)° _0,t ^–1 , whereT(t) is a stochastic process with values in distributions and _s,t is a stochastic flow generated by a certain stochastic differential equation. The representation is used for getting estimates of the solution with respect to Sobolev norms.Further, by applying the partial Malliavin calculus of Kusuoka-Stroock, we show that any generalized solution is aC -function under a condition similar to Hörmander's hypoellipticity condition.     

Approximate minimum-cost multicommodity flows in\tilde O(ɛ −2 KNM) timetime

We show that an ε-approximate solution of the cost-constrainedK-commodity flow problem on anN-nodeM-arc network,G can be computed by sequentially solving O(K(? ^?2+logGK) logGM log (G? ^?1 GK)) single-commodity minimum-cost flow problems on the same network. In particular, an approximate minimum-cost multicommodity flow can be computed in $\tilde O$ (G? ^?2 GKNM) running time, where the notation Õ(·) means “up to logarithmic factors”. This result improves the time bound mentioned by Grigoriadis and Khachiyan [4] by a factor ofM/N and that developed more recently by Karger and Plotkin [8] by a factor of? ^?1. We also provide a simple $\tilde O$ (NM)-time algorithm for single-commodity budget-constrained minimum-cost flows which is $\tilde O$ (? ^?3) times faster than the algorithm developed in the latter paper.     

THE PARAMETER KERNELS FOR THE LINEAR INVERSION OF THE FREE OSCILLATION OF THE EARTH (Ⅰ)

In the linear inversion of the radial variation of the parameters of the Earth by usingthe observed frequencies of various normal modes of free oscillation of the earth, it is neces-sary to know the values of the kernels of the parameters ρ, μ and λ. This paper describesthe methods of the derivation of the formulas of these kernels. This is the first part of thepaper in which only the toroidal oscillations are considered. They are much simpler thanthose of the spheroidal ones, that we will consider in the second part of the paper. The data of the two types of oscillations are equally important in the solution of theinversion problem, and should be employed simultaneously, and we know that the toroidaloscillations are much simpler than the spheroidal ones, it seems wise to divide the whole programof the inversion problem into steps: first, by employing the toroidal data to correct the twoparameters ρ and μ in the mantle, then by using the spheroidal data to correct the remain-ing parameters, i.e. th     

Modeling and simulation of free radical polymerization of styrene under semibatch reactor conditions

The first part of this approach is concerned with the elaboration of a radical polymerization model of styrenne, based on a kinetic diagram that includes chemical and thermal initiation, propagation, termination by recombination and chain transfer to the monomer. Furthermore, volume contraction during polymerization is considered, as well as the gel and glass effects. The mathematical formalism that describes the model in terms of moments is explored in detail. The model was then used to predict the changes in monomer conversion and molecular weight after intermediate addition of initiator and monomer. The results of this operation are dependent on the conditions of the reaction mass, quantity, and moment of substance addition. Therefore, the simulations were performed at different times with respect to the gel effect; before, during and after this phenomenon, and also with respect to different temperatures and initiators. Increasing the initiator concentration before the gel effect leads to an earlier appearance of the phenomenon and to a decrease in molecular weight. The ratio reveals a polydispersity index smaller for the intermediate addition of initiator. No significant changes take place during or after the gel effect. If along with the initiator, unreacted monomver (used to dissolve the initiator) enters the reactor, a small dip in conversion is observed. The general conclusion of this paper reveals the intermediate addition of initiator as a method to control polymer properties and to prevent the “dead-end” polymerization of styrene.     

Free energies of transfer of alkali halides from water to urea-water mixtures at 25°C from EMF measurements

The standard potentials_s E ^o of M/M⁺ (M=Li, Na, K, Rb, and Cs) electrodes in aqueous urea solutions containing 12, 20, 30 and 37% by weight of urea have been determined at 25°C from emf measurements on the cell M(Hg)/MCl (m), solvent/AgCl–Ag, from the activities of metals in amalgams by use of a similar type of cell in water, and from the values of_s E ^o of the Ag/AgCl electrode determined earlier. The standard free energies of transfer of MCl, G _t ^o (MCl), from water to the mixed solvents, computed by use of these values and those for the Ag–AgCl electrode, rise sharply from Li⁺ to Na⁺ but fall from Na⁺ to K⁺ and rather sharply from K⁺ to Cs⁺ with a maximum at Na⁺ in all the solvent compositions. This has been attributed to the superimposition of soft-soft interactions on the electrostatic interactions between the ions and the negative charge centers of the possible hydrogen-bonded solvent complexes in the mixed solvents. Comparison of G _t ^o (i) values for individual ions, obtained by a simultaneous extrapolation procedure, with those in aqueous mixtures of methanol,t-butanol, and dimethyl sulfoxide leads to the conclusion that the solvation of these ions in all these solvents is chiefly dictated by the acid-base type of ion-solvent interactions.