Two-fluid and population balance models for subcooled boiling flow

Population balance equations combined with a three-dimensional two-fluid model are employed to predict subcooled boiling flow at low pressure in a vertical annular channel. The MUSIG (MUltiple-SIze-Group) model implemented in the computer code CFX4.4 is further developed to accommodate the wall nucleation at the heated wall and condensation in the subcooled boiling regime. Comparison of model predictions against local measurements is made for the void fraction, bubble Sauter mean diameter and gas and liquid velocities covering a range of different mass and heat fluxes and inlet subcooling temperatures. Additional comparison using empirical relationships for the active nucelation site density and local bubble diameter is also investigated. Good agreement is achieved with the local radial void fraction, bubble Sauter diameter and liquid velocity profiles against measurements. However, significant weakness of the model is evidenced in the prediction of the vapour velocity. Work is in progress to circumvent the deficiency through the consideration of additional momentum equations or developing an algebraic slip model to account for bubble separation.     

Effect of crack tip sharpness on the strength of vulcanized rubbers

The effect of crack tip sharpness on crack propagation in vulcanized rubbers has been studied. For very sharp cracks, tearing is found to occur on a small scale at very low energies not far above the threshold required for the onset of mechanical crack growth. The “small-scale” tearing energies show relatively little variation for rubbers that differ widely in tear strength as normally measured. Thus the latter property appears to be strongly influenced by variations in the ability of rubbers to promote tip blunting. The small-scale tear behavior is of relevance to other fracture phenomena, including cutting by sharp objects and tensile failure. Natural variations in tip sharpness occur during cyclic or time-dependent mechanical crack growth and influence the form of the crack growth characteristics.     

The intersection marker method for 3D interface tracking of deformable surfaces in finite volumes

Currently, the majority of computational fluid dynamics (CFD) codes use the finite volume method to spatially discretise the computational domain, sometimes as an array of cubic control volumes. The Finite volume method works well with single‐phase flow simulations, but two‐phase flow simulations are more challenging because of the need to track the surface interface traversing and deforming within the 3D grid. Surface area and volume fraction details of each interface cell must be accurately accounted for, in order to calculate for the momentum exchange and rates of heat and mass transfer across the interface. To attain a higher accuracy in two‐phase flow CFD calculations, the intersection marker (ISM) method is developed. The ISM method is a hybrid Lagrangian–Eulerian front‐tracking algorithm that can model an arbitrary 3D surface within an array of cubic control volumes. The ISM method has a cell‐by‐cell remeshing capability that is volume conservative and is suitable for the tracking of complex interface deformation in transient two‐phase CFD simulations. Copyright © 2015 John Wiley & Sons, Ltd.     

Automated determination of size and morphology information from soot transmission electron microscope (TEM)-generated images

The thermophoretic sampling of particulates from hot media, coupled with transmission electron microscope (TEM) imaging, is a combined approach that is widely used to derive morphological information. The identification and the measurement of the particulates, however, can be complex when the TEM images are of low contrast, noisy, and have non-uniform background signal level. The image processing method can also be challenging and time consuming, when the samples collected have large variability in shape and size, or have some degree of overlapping. In this work, a three-stage image processing sequence is presented to facilitate time-efficient automated identification and measurement of particulates from the TEM grids. The proposed processing sequence is first applied to soot samples that were thermophoretically sampled from a laminar non-premixed ethylene-air flame. The parameter values that are required to be set to facilitate the automated process are identified, and sensitivity of the results to these parameters is assessed. The same analysis process is also applied to soot samples that were acquired from an externally irradiated laminar non-premixed ethylene-air flame, which have different geometrical characteristics, to assess the morphological dependence of the proposed image processing sequence. Using the optimized parameter values, statistical assessments of the automated results reveal that the largest discrepancies that are associated with the estimated values of primary particle diameter, fractal dimension, and prefactor values of the aggregates for the tested cases, are approximately 3, 1, and 10 %, respectively, when compared with the manual measurements.     

The synthesis and self‐assembly of ABA amphiphilic block copolymers containing styrene and oligo(ethylene glycol) methyl ether methacrylate in dilute aqueous solutions: Elevated cloud point temperatures for thermoresponsive micelles

A series of ABA amphiphilic triblock copolymers possessing polystyrene (PS) central hydrophobic blocks, one group with “short” PS blocks (DP = 54–86) and one with “long” PS blocks (DP = 183–204) were synthesized by atom transfer radical polymerization. The outer hydrophilic blocks were various lengths of poly(oligoethylene glycol methyl ether) methacrylate, a comb‐like polymer. The critical aggregation concentrations were recorded for certain block copolymer samples and were found to be in the range circa 10⁻⁹ mol L⁻¹ for short PS blocks and circa 10⁻¹² mol L⁻¹ for long PS blocks. Dilute aqueous solutions were analyzed by transmission electron microscopy (TEM) and demonstrated that the short PS block copolymers formed spherical micelles and the long PS block copolymers formed predominantly spherical micelles with smaller proportions of cylindrical and Y‐branched cylindrical micelles. Dynamic light scattering analysis results agreed with the TEM observations demonstrating variations in micelle size with PS and POEGMA chain length: the hydrodynamic diameters (D_H) of the shorter PS block copolymer micelles increased with increasing POEGMA block lengths while maintaining similar PS micellar core diameters (D_C); in contrast the values of D_H and D_C for the longer PS block copolymer micelles decreased. Surface‐pressure isotherms were recorded for two of the samples and these indicated close packing of a short PS block copolymer at the air–water interface. The aggregate solutions were demonstrated to be stable over a 38‐day period with no change in aggregate size or noticeable precipitation. The cloud point temperatures of certain block copolymer aggregate solutions were measured and found to be in the range 76–93 °C; significantly these were ∼11 °C higher in temperature than those of POEGMA homopolymer samples with similar chain lengths. © 2008 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 46: 7739–7756, 2008     

Characterization of Acacia mangium polyflavonoid tannins by MALDI-TOF mass spectrometry and CP-MAS C NMR

The MALDI-TOF mass spectrometry (MS) and solid state CP-MAS ¹³C Nuclear Magnetic Resonance (NMR) spectroscopic technique were introduced to characterize Acacia mangium tannin (condensed tannins). The MALDI-TOF MS illustrated a series of peaks corresponding to oligomers of condensed tannins of up to 11 flavonoid units (3200 Da). A. mangium condensed tannins were found to consist predominantly of prorobinetinidin combined with profisetinidin and prodelphinidin. Both the MALDI-TOF mass spectra and the solid state CP-MAS ¹³C NMR indicated that the A. mangium tannins obtained from Kudat, had an almost completely linear structure; In addition, Lembah Beringin, consist of “angular” polymer structure; and Tawau, has included “twice-angular” polymer structures present in oligomers type of up to 7 flavonoid units. The high degree of polymerization of linear, angular type, twice-angular structures and longer oligomer (3200 Da) chains have not been observed in previous studies of condensed tannins. The spectra also indicated that A. mangium tannins are more heavily branched and have higher degree of polymerization (>7.0) compared to commercial mimosa (A. mearnsii) tannin (4.9). Because tannins are phenolic, it was expected that they can be used to replace phenol-formaldehyde (PF) adhesives.     

Numerical modelling of bubbly flows with and without heat and mass transfer

In this paper, modelling gas–liquid bubbly flows is achieved by the introduction of a population balance equation combined with the three-dimensional two-fluid model. For gas–liquid bubbly flows without heat and mass transfer, an average bubble number density transport equation has been incorporated in the commercial code CFX5.7 to better describe the temporal and spatial evolution of the geometrical structure of the gas bubbles. The coalescence and breakage effects of the gas bubbles are modelled according to the coalescence by the random collisions driven by turbulence and wake entrainment while for bubble breakage by the impact of turbulent eddies. Local radial distributions of the void fraction, interfacial area concentration, bubble Sauter mean diameter, and gas and liquid velocities, are compared against experimental data in a vertical pipe flow. Satisfactory agreements for the local distributions are achieved between the predictions and measurements. For gas–liquid bubbly flows with heat and mass transfer, boiling flows at subcooled conditions are considered. Based on the formulation of the MUSIG (multiple-size-group) boiling model and a model considering the forces acting on departing bubbles at the heated surface implemented in the computer code CFX4.4, comparison of model predictions against local measurements is made for the void fraction, bubble Sauter mean diameter, interfacial area concentration, and gas and liquid velocities covering a range of different mass and heat fluxes and inlet subcooling temperatures. Good agreement is achieved with the local radial void fraction, bubble Sauter mean diameter, interfacial area concentration and liquid velocity profiles against measurements. However, significant weakness of the model is evidenced in the prediction of the vapour velocity. Work is in progress through the consideration of additional momentum equations or developing an algebraic slip model to account for the effects of bubble separation.     

Comparative Large Eddy Simulation study of a large-scale buoyant fire

A fully-coupled Large Eddy Simulation model which incorporates all essential combustion, radiation and soot chemistry considerations have been developed to simulate the temporal vortical structure of a large-scale buoyant fire. Numerical results are validated and compared against a full-scale fire measurements and predictions from other LES models. Quantitative comparisons against experimental data suggested that the present model successfully captured the vortical structures and the puffing behaviour of a buoyant fire.