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.

A finite difference scheme for steady, supersonic, two-dimensional, compressible flow of real gases

A finite difference scheme is presented for the solution of the two-dimensional equations of steady, supersonic, compressible flow of real gases. The scheme incorparates numerical characteristic decomposition, is shock-capturing by design and incorporates space-marching as a result of the assumption that the flow is wholly supersonic in at least one space dimension. Results are shown for problems involving oblique hydraulic jumps and reflection from a wall.

A numerical scheme for two-dimensional, open channel flows with non-rectangular geometries

An algorithm based on flux difference splitting is presented for the solution of two-dimensional, open channel flows. A transformation maps a non-rectangular, physical domain into a rectangular one. The governing equations are then the shallow water equations, including terms of slope and friction, in a generalized coordinate system. A regular mesh on a rectangular computational domain can then be employed. The resulting scheme has good jump capturing properties and the advantage of using boundary/body-fitted meshes. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Riemann solvers with primitive parameter vectors for two-dimensional compressible flows of a real gas

An analysis of averaging procedures is presented for an approximate Riemann solver for the equations governing the compressible flow of a real gas. This study extends earlier work for the Euler equations with ideal gases.

Ultrafast charge separation in a photoreactive rhenium-appended porphyrin assembly monitored by picosecond transient infrared spectroscopy

Rhenium(bipyridine)(tricarbonyl)(picoline) units have been linked covalently to tetraphenylmetalloporphyrins of magnesium and zinc via an amide bond between the bipyridine and one phenyl substituent of the porphyrin. The resulting complexes, abbreviated as [Re(CO)(3)(Pic)Bpy-MgTPP][OTf] and [Re(CO)(3)(Pic)Bpy-ZnTPP][OTf], exhibit no signs of electronic interaction between the Re(CO)(3)(bpy) units and the metalloporphyrin units in their ground states. However, emission spectroscopy reveals solvent-dependent quenching of porphyrin emission on irradiation into the long-wavelength absorption bands localized on the porphyrin. The characteristics of the excited states have been probed by picosecond time-resolved absorption (TRVIS) spectroscopy and time-resolved infrared (TRIR) spectroscopy in nitrile solvents. The presence of the charge-separated state involving electron transfer from MgTPP or ZnTPP to Re(bpy) is signaled in the TRIR spectra by a low-frequency shift in the nu(CO) bands of the Re(CO)(3) moiety similar to that observed by spectroelectrochemical reduction. Long-wavelength excitation of [Re(CO)(3)(Pic)Bpy-MTPP][OTf] results in characteristic TRVIS spectra of the S-1 state of the porphyrin that decay with a time constant of 17 ps (M = Mg) or 24 ps (M = Zn). The IR bands of the CS state appear on a time scale of less than 1 ps (Mg) or ca. 5 ps (Zn) and decay giving way to a vibrationally excited (i.e., hot) ground state via back electron transfer. The IR bands of the precursors recover with a time constant of 35 ps (Mg) or 55 ps (Zn). The short lifetimes of the charge-transfer states carry implications for the mechanism of reaction in the presence of triethylamine.

Switching on a two-dimensional time-harmonic scalar wave in the presence of a diffracting edge

This paper concerns the switching on of two-dimensional time-harmonic scalar waves. We first review the switch-on problem for a point source in free space, then proceed to analyse the analogous problem for the diffraction of a plane wave by a half-line (the ‘Sommerfeld problem’), determining in both cases the conditions under which the field is well-approximated by the solution of the corresponding frequency domain problem. In both cases the rate of convergence to the frequency domain solution is found to be dependent on the strength of the singularity on the leading wavefront. In the case of plane wave diffraction at grazing incidence the frequency domain solution is immediately attained along the shadow boundary after the arrival of the leading wavefront. The case of non-grazing incidence is also considered.

Dynamic stall for a Vertical Axis Wind Turbine in a two-dimensional study

The last few years have proved that Vertical Axis Wind Turbines (VAWTs) are more suitable for urban areas than Horizontal Axis Wind Turbines (HAWTs). To date, very little has been published in this area to assess good performance and lifetime of VAWTs either in open or urban areas. At low tip speed ratios (TSRs<5), VAWTs are subjected to a phenomenon called 'dynamic stall'. This can really affect the fatigue life of a VAWT if it is not well understood. The purpose of this paper is to investigate how CFD is able to simulate the dynamic stall for 2-D flow around VAWT blades. During the numerical simulations different turbulence models were used and compared with the data available on the subject. In this numerical analysis the Shear Stress Transport (SST) turbulence model seems to predict the dynamic stall better than the other turbulence models available. The limitations of the study are that the simulations are based on a 2-D case with constant wind and rotational speeds instead of considering a 3-D case with variable wind speeds. This approach was necessary for having a numerical analysis at low computational cost and time. Consequently, in the future it is strongly suggested to develop a more sophisticated model that is a more realistic simulation of a dynamic stall in a three-dimensional VAWT.

Prediction of coagulation properties, titratable acidity and pH of bovine milk using mid-infrared spectroscopy

This study investigated the potential application of mid-infrared spectroscopy (MIR 4,000–900 cm−1) for the determination of milk coagulation properties (MCP), titratable acidity (TA), and pH in Brown Swiss milk samples (n = 1,064). Because MCP directly influence the efficiency of the cheese-making process, there is strong industrial interest in developing a rapid method for their assessment. Currently, the determination of MCP involves time-consuming laboratory-based measurements, and it is not feasible to carry out these measurements on the large numbers of milk samples associated with milk recording programs. Mid-infrared spectroscopy is an objective and nondestructive technique providing rapid real-time analysis of food compositional and quality parameters. Analysis of milk rennet coagulation time (RCT, min), curd firmness (a30, mm), TA (SH°/50 mL; SH° = Soxhlet-Henkel degree), and pH was carried out, and MIR data were recorded over the spectral range of 4,000 to 900 cm−1. Models were developed by partial least squares regression using untreated and pretreated spectra. The MCP, TA, and pH prediction models were improved by using the combined spectral ranges of 1,600 to 900 cm−1, 3,040 to 1,700 cm−1, and 4,000 to 3,470 cm−1. The root mean square errors of cross-validation for the developed models were 2.36 min (RCT, range 24.9 min), 6.86 mm (a30, range 58 mm), 0.25 SH°/50 mL (TA, range 3.58 SH°/50 mL), and 0.07 (pH, range 1.15). The most successfully predicted attributes were TA, RCT, and pH. The model for the prediction of TA provided approximate prediction (R2 = 0.66), whereas the predictive models developed for RCT and pH could discriminate between high and low values (R2 = 0.59 to 0.62). It was concluded that, although the models require further development to improve their accuracy before their application in industry, MIR spectroscopy has potential application for the assessment of RCT, TA, and pH during routine milk analysis in the dairy industry. The implementation of such models could be a means of improving MCP through phenotypic-based selection programs and to amend milk payment systems to incorporate MCP into their payment criteria.

Visible-near infrared spectroscopy sensor for predicting curd and whey composition during cheese processing

The potential of visible-near infrared spectra, obtained using a light backscatter sensor, in conjunction with chemometrics, to predict curd moisture and whey fat content in a cheese vat was examined. A three-factor (renneting temperature, calcium chloride, cutting time), central composite design was carried out in triplicate. Spectra (300–1,100 nm) of the product in the cheese vat were captured during syneresis using a prototype light backscatter sensor. Stirring followed upon cutting the gel, and samples of curd and whey were removed at 10 min intervals and analyzed for curd moisture and whey fat content. Spectral data were used to develop models for predicting curd moisture and whey fat contents using partial least squares regression. Subjecting the spectral data set to Jack-knifing improved the accuracy of the models. The whey fat models (R = 0.91, 0.95) and curd moisture model (R = 0.86, 0.89) provided good and approximate predictions, respectively. Visible-near infrared spectroscopy was found to have potential for the prediction of important syneresis indices in stirred cheese vats.