An efficient numerical scheme for the two-dimensional shallow water equations using arithmetic averaging

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gas dynamics is defined, and a scheme, based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem, and incorporates the technique of operator splitting. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. An extension to the two-dimensional equations with source terms is included. The scheme is applied to the one-dimensional problems of a breaking dam and reflection of a bore, and in each case the approximate solution is compared to the exact solution of ideal fluid flow. The scheme is also applied to a problem of stationary bore generation in a channel of variable cross-section. Finally, the scheme is applied to two other dam-break problems, this time in two dimensions with one having cylindrical symmetry. Each approximate solution compares well with those given by other authors.

Solutions of a two-dimensional dam-break problem

Solutions of a two-dimensional dam break problem are presented for two tailwater/reservoir height ratios. The numerical scheme used is an extension of one previously given by the author [J. Hyd. Res. 26(3), 293–306 (1988)], and is based on numerical characteristic decomposition. Thus approximate solutions are obtained via linearised problems, and the method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids non-physical, spurious oscillations.

Prediction of supercritical flow in open channels

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow water equations in open channels. A linearised problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Flux difference splitting for open channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.

Error estimates for a fully discrete spectral scheme for a class of nonlinear, nonlocal dispersive wave equations

We analyze a fully discrete spectral method for the numerical solution of the initial- and periodic boundary-value problem for two nonlinear, nonlocal, dispersive wave equations, the Benjamin–Ono and the Intermediate Long Wave equations. The equations are discretized in space by the standard Fourier–Galerkin spectral method and in time by the explicit leap-frog scheme. For the resulting fully discrete, conditionally stable scheme we prove an L2-error bound of spectral accuracy in space and of second-order accuracy in time.

Time-resolved gas-phase kinetic study of the germylene addition reaction, GeH2 + C2D2

Time-resolved studies of germylene, GeH2, generated by laser. ash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate coefficients for its bimolecular reaction with C2D2. The reaction was studied in the gas phase, mainly at a total pressure of 1.3 kPa (in SF6 bath gas) at five temperatures in the range 298-558 K. Pressure variation measurements over the range 0.13-13 kPa ( SF6) at 298, 397 and 558 K revealed a small pressure dependence but only at 558 K. After correction for this, the second-order rate coefficients gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.96 +/- 0.05) + ( 6.16 +/- 0.37 kJ mol(-1))/RT ln 10 Comparison with the reaction of GeH2 + C2H2 (studied earlier) showed a similar behaviour with almost identical rate coefficients. The lack of a significant isotope effect is consistent with a rate-determining addition process and is explained by irreversible decomposition of the reaction intermediate to give Ge(P-3) + C2H4. This result contrasts with that for GeH2 + C2H4/C2D4 and those for the analogous silylene reactions. The underlying reasons for this are discussed.

Gas-phase kinetics of chlorosilylene reactions II. ClSiH + C2H4: absolute rate measurements and quantum chemical and RRKM calculations for the prototype pi addition reaction

Time-resolved studies of chlorosilylene, ClSiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, are carried out to obtain rate constants for its bimolecular reaction with ethene, C2H4, in the gas-phase. The reaction is studied over the pressure range 0.13-13.3 kPa (with added SF6) at five temperatures in the range 296-562 K. The second order rate constants, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1))=(-10.55 +/- 0.10) + (3.86 +/- 0.70) kJ mol(-1)/RT ln10. The Arrhenius parameters correspond to a loose transition state and the rate constant at room temperature is 43% of that for SiH2 + C2H4, showing that the deactivating effect of Cl-for-H substitution in the silylene is not large. Quantum chemical calculations of the potential energy surface for this reaction at the G3MP2//B3LYP level show that, as well as 1-chlorosilirane, ethylchlorosilylene is a viable product. The calculations reveal how the added effect of the Cl atom on the divalent state stabilisation of ClSiH influences the course of this reaction. RRKM calculations of the reaction pressure dependence suggest that ethylchlorosilylene should be the main product. The results are compared and contrasted with those of SiH2 and SiCl2 with C2H4.

Time-resolved gas-phase kinetic, quantum chemical and RRKM studies of reactions of silylene with cyclic ethers

Time-resolved kinetic studies of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reactions with oxirane, oxetane, and tetrahydrofuran (THF). The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at four or five temperatures in the range 294-605 K. All three reactions showed pressure dependences characteristic of third-body-assisted association reactions with, surprisingly, SiH2 + oxirane showing the least and SiH2 + THF showing the most pressure dependence. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equations where the error limits are single standard deviations: log(k(oxirane)(infinity)/cm(3) molecule(-1) s(-1)) = (-11.03 +/- 0.07) + (5.70 +/- 0.51) kJ mol(-1)/RT In 10 log(k(oxetane)(infinity)/cm(3) molecule(-1) s(-1)) = (-11.17 +/- 0.11) + (9.04 +/- 0.78) kJ mol(-1)/RT In 10 log(k(THF)(infinity)/cm(3) molecule(-1) s(-1)) = (-10.59 +/- 0.10) + (5.76 +/- 0.65) kJ mol(-1)/RT In 10 Binding-energy values of 77, 97, and 92 kJ mol(-1) have been obtained for the donor-acceptor complexes of SiH2 with oxirane, oxetane, and THF, respectively, by means of quantum chemical (ab initio) calculations carried Out at the G3 level. The use of these values to model the pressure dependences of these reactions, via RRKM theory, provided a good fit only in the case of SiH2 + THF. The lack of fit in the other two cases is attributed to further reaction pathways for the association complexes of SiH2 with oxirane and oxetane. The finding of ethene as a product of the SiH2 + oxirane reaction supports a pathway leading to H2Si=O + C2H4 predicted by the theoretical calculations of Apeloig and Sklenak.

Time-resolved gas-phase kinetic study of the germylene addition reaction, CH3C CCH3

Time-resolved kinetic studies of the reaction of germylene, GeH2, generated by laser. ash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reaction with 2-butyne, CH3C CCH3. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 300-556 K. The second order rate constants obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.46 +/- 10.06) + (5.16 +/- 10.47) kJ mol(-1)/ RT ln 10 Calculations of the energy surface of the GeC4H8 reaction system were carried out employing the additivity principle, by combining previous quantum chemical calculations of related reaction systems. These support formation of 1,2-dimethylvinylgermylene (rather than 2,3-dimethylgermirene) as the end product. RRKM calculations of the pressure dependence of the reaction are in reasonable agreement with this finding. The reactions of GeH2 with C2H2 and with CH3CRCCH3 are compared and contrasted.

The gas-phase reaction between silylene and 2-butyne: kinetics, isotope studies, pressure dependence studies and quantum chemical calculations

Time-resolved kinetic studies of the reactions of silylene, SiH2, and dideutero-silylene, SiD2, generated by laser. ash photolysis of phenylsilane and phenylsilane-d(3), respectively, have been carried out to obtain rate coefficients for their bimolecular reactions with 2-butyne, CH3C CCH3. The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 294-612 K. The second-order rate coefficients, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equations where the error limits are single standard deviations: log(k(H)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.67 +/- 0.04) + (1.71 +/- 0.33) kJ mol(-1)/RTln10 log(k(D)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.65 +/- 0.01) + (1.92 +/- 0.13) kJ mol(-1)/RTln10 Additionally, pressure-dependent rate coefficients for the reaction of SiH2 with 2-butyne in the presence of He (1-100 Torr) were obtained at 301, 429 and 613 K. Quantum chemical (ab initio) calculations of the SiC4H8 reaction system at the G3 level support the formation of 2,3-dimethylsilirene [cyclo-SiH2C(CH3)=C(CH3)-] as the sole end product. However, reversible formation of 2,3-dimethylvinylsilylene [CH3CH=C(CH3)SiH] is also an important process. The calculations also indicate the probable involvement of several other intermediates, and possible products. RRKM calculations are in reasonable agreement with the pressure dependences at an enthalpy value for 2,3-dimethylsilirene fairly close to that suggested by the ab initio calculations. The experimental isotope effects deviate significantly from those predicted by RRKM theory. The differences can be explained by an isotopic scrambling mechanism, involving H - D exchange between the hydrogens of the methyl groups and the D-atoms in the ring in 2,3-dimethylsilirene-1,1-d(2). A detailed mechanism involving several intermediate species, which is consistent with the G3 energy surface, is proposed to account for this.

Motion estimation through approximated densities

Many techniques are currently used for motion estimation. In the block-based approaches the most common procedure applied is the block-matching based on various algorithms. To refine the motion estimates resulting from the full search or any coarse search algorithm, one can find few applications of Kalman filtering, mainly in the intraframe scheme. The Kalman filtering technique applicability for block-based motion estimation is rather limited due to discontinuities in the dynamic behaviour of the motion vectors. Therefore, we propose an application of the concept of the filtering by approximated densities (FAD). The FAD, originally introduced to alleviate limitations due to conventional Kalman modelling, is applied to interframe block-motion estimation. This application uses a simple form of FAD involving statistical characteristics of multi-modal distributions up to second order.

Velocity-based moving mesh methods for nonlinear partial differential equations

This article describes a number of velocity-based moving mesh numerical methods formultidimensional nonlinear time-dependent partial differential equations (PDEs). It consists of a short historical review followed by a detailed description of a recently developed multidimensional moving mesh finite element method based on conservation. Finite element algorithms are derived for both mass-conserving and non mass-conserving problems, and results shown for a number of multidimensional nonlinear test problems, including the second order porous medium equation and the fourth order thin film equation as well as a two-phase problem. Further applications and extensions are referenced.

Survival of Lactobacillus plantarum in model solutions and fruit juices

The aim of the work was to study the survival of Lactobacillus plantarum NCIMB 8826 in model solutions and develop a mathematical model describing its dependence on pH, citric acid and ascorbic acid. A Central Composite Design (CCD) was developed studying each of the three factors at five levels within the following ranges, i.e., pH (3.0-4.2), citric acid (6-40 g/L), and ascorbic acid (100-1000 mg/L). In total, 17 experimental runs were carried out. The initial cell concentration in the model solutions was approximately 1 × 10(8)CFU/mL; the solutions were stored at 4°C for 6 weeks. Analysis of variance (ANOVA) of the stepwise regression demonstrated that a second order polynomial model fits well the data. The results demonstrated that high pH and citric acid concentration enhanced cell survival; one the other hand, ascorbic acid did not have an effect. Cell survival during storage was also investigated in various types of juices, including orange, grapefruit, blackcurrant, pineapple, pomegranate, cranberry and lemon juice. The model predicted well the cell survival in orange, blackcurrant and pineapple, however it failed to predict cell survival in grapefruit and pomegranate, indicating the influence of additional factors, besides pH and citric acid, on cell survival. Very good cell survival (less than 0.4 log decrease) was observed after 6 weeks of storage in orange, blackcurrant and pineapple juice, all of which had a pH of about 3.8. Cell survival in cranberry and pomegranate decreased very quickly, whereas in the case of lemon juice, the cell concentration decreased approximately 1.1 logs after 6 weeks of storage, albeit the fact that lemon juice had the lowest pH (pH~2.5) among all the juices tested. Taking into account the results from the compositional analysis of the juices and the model, it was deduced that in certain juices, other compounds seemed to protect the cells during storage; these were likely to be proteins and dietary fibre In contrast, in certain juices, such as pomegranate, cell survival was much lower than expected; this could be due to the presence of antimicrobial compounds, such as phenolic compounds.

A comprehensive investigation of the synthesis of prebiotic galactooligosaccharides by whole cells of Bifidobacterium bifidum NCIMB 41171

The synthesis of galactooligosaccharides (GOS) by whole cells of Bifidobacterium bifidum NCIMB 41171 was investigated by developing a set of mathematical models. These were second order polynomial equations, which described responses related to the production of GOS constituents, the selectivity of lactose conversion into GOS, and the relative composition of the produced GOS mixture, as a function of the amount of biocatalyst, temperature, initial lactose concentration, and time. The synthesis reactions were followed for up to 36 h. Samples were withdrawn every 4 h, tested for β-galactosidase activity, and analysed for their carbohydrate content. GOS synthesis was well explained by the models, which were all significant (P < 0.001). The GOS yield increased as temperature increased from 40 °C to 60 °C, as transgalactosylation became more pronounced compared to hydrolysis. The relative composition of GOS produced changed significantly with the initial lactose concentration (P < 0.001); higher ratios of tri-, tetra-, and penta-galactooligosaccharides to transgalactosylated disaccharides were obtained as lactose concentration increased. Time was a critical factor, as a balanced state between GOS synthesis and hydrolysis was roughly attained in most cases between 12 and 20 h, and was followed by more pronounced GOS hydrolysis than synthesis.

Conceptualising and measuring the equity of online brands

Although brand equity is an important source of competitive advantage online, previous conceptualisations and measures overlook the unique characteristics of the internet that render consumers co-creators of brand value. In view of this, a threephased research programme was undertaken to identify the facets of online retail/service (ORS) brand equity and then develop and validate a scale for its measurement. ORS brand equity was found to be a second order construct with five correlated yet distinct dimensions: emotional connection, online experience, responsive service nature, trust, and fulfilment. A series of tests showed that the ensuing 12-item scale has strong psychometric properties. The implications of this research for marketing researchers and practitioners are discussed.