Numerical simulation of incompressible flow problems using an unstructured staggered mesh method

A number of two dimensional staggered unstructured discretisation schemes for the solution of fluid flow and heat transfer problems have been developed. All schemes store and solve velocity vector components at cell faces with scalar variables solved at cell centres. The velocity is resolved into face-normal and face-parallel components and the various schemes investigated differ in the treatment of the parallel component. Steady-state and time-dependent fluid flow and thermal energy equations are solved with the well known pressure correction scheme, SIMPLE, employed to couple continuity and momentum. The numerical methods developed are tested on well known benchmark cases: the Lid-Driven Cavity, Natural Convection in a Cavity and Melting of Gallium in a rectangular domain. The results obtained are shown to be comparable to benchmark, but with accuracy dependent on scheme selection.

CFD analysis of one-dimensional infiltration in vadose zone

This study presents a CFD analysis constructed around PHYSICA, an open framework for multi-physics computational continuum mechanics modelling, to investigate the water movement in unsaturated porous media. The modelling environment is based on a cell-centred finite-volume discretisation technique. A number of test cases are performed in order to validate the correct implementation of Richard's equation for compressible and incompressible fluids. The pressure head form of the equation is used together with the constitutive relationships between pressure, volumetric water content and hydraulic conductivity described by Haverkamp and Van Genuchten models. The flow problems presented are associated with infiltration into initially dry soils with homogeneous or layered geologic settings. Comparison of results with the problems selected from literature shows a good agreement and validates the approach and the implementation.

Investigating the representation of merging behavior at the floor–stair interface in computer simulations of multi-floor building evacuations

In this article, the representation of the merging process at the floor— stair interface is examined within a comprehensive evacuation model and trends found in experimental data are compared with model predictions. The analysis suggests that the representation of floor—stair merging within the comprehensive model appears to be consistent with trends observed within several published experiments of the merging process. In particular: (a) The floor flow rate onto the stairs decreases as the stair population density increases. (b) For a given stair population density, the floor population's flow rate onto the stairs can be maximized by connecting the floor to the landing adjacent to the incoming stair. (c) In situations where the floor is connected adjacent to the incoming stair, the merging process appears to be biased in favor of the floor population. It is further conjectured that when the floor is connected opposite the incoming stair, the merging process between the stair and floor streams is almost in balance for high stair population densities, with a slight bias in favor of the floor stream at low population densities. A key practical finding of this analysis is that the speed at which a floor can be emptied onto a stair can be enhanced simply by connecting the floor to the landing at a location adjacent to the incoming stair rather than opposite the stair. Configuring the stair in this way, while reducing the floor emptying time, results in a corresponding decrease in the descent flow rate of those already on the stairs. While this is expected to have a negligible impact on the overall time to evacuate the building, the evacuation time for those higher up in the building is extended while those on the lower flows is reduced. It is thus suggested that in high-rise buildings, floors should be connected to the landing on the opposite side to the incoming stair. Information of this type will allow engineers to better design stair—floor interfaces to meet specific design objectives.

Multilevel approaches applied to the capacitated clustering problem

This paper presents two multilevel refinement algorithms for the capacitated clustering problem. Multilevel refinement is a collaborative technique capable of significantly aiding the solution process for optimisation problems. The central methodologies of the technique are filtering solutions from the search space and reducing the level of problem detail to be considered at each level of the solution process. The first multilevel algorithm uses a simple tabu search while the other executes a standard local search procedure. Both algorithms demonstrate that the multilevel technique is capable of aiding the solution process for this combinatorial optimisation problem.

Image denoise by fourth-order PDE based on the changes of laplacian

Fourth-order partial differential equation (PDE) proposed by You and Kaveh (You-Kaveh fourth-order PDE), which replaces the gradient operator in classical second-order nonlinear diffusion methods with a Laplacian operator, is able to avoid blocky effects often caused by second-order nonlinear PDEs. However, the equation brought forward by You and Kaveh tends to leave the processed images with isolated black and white speckles. Although You and Kaveh use median filters to filter these speckles, median filters can blur the processed images to some extent, which weakens the result of You-Kaveh fourth-order PDE. In this paper, the reason why You-Kaveh fourth-order PDE can leave the processed images with isolated black and white speckles is analyzed, and a new fourth-order PDE based on the changes of Laplacian (LC fourth-order PDE) is proposed and tested. The new fourth-order PDE preserves the advantage of You-Kaveh fourth-order PDE and avoids leaving isolated black and white speckles. Moreover, the new fourth-order PDE keeps the boundary from being blurred and preserves the nuance in the processed images, so, the processed images look very natural.

An algorithm for total variation inpainting based on nonlinear multi-grid methods

Image inpainting refers to restoring a damaged image with missing information. The total variation (TV) inpainting model is one such method that simultaneously fills in the regions with available information from their surroundings and eliminates noises. The method works well with small narrow inpainting domains. However there remains an urgent need to develop fast iterative solvers, as the underlying problem sizes are large. In addition one needs to tackle the imbalance of results between inpainting and denoising. When the inpainting regions are thick and large, the procedure of inpainting works quite slowly and usually requires a significant number of iterations and leads inevitably to oversmoothing in the outside of the inpainting domain. To overcome these difficulties, we propose a solution for TV inpainting method based on the nonlinear multi-grid algorithm.

Application of particle swarm optimisation to evaluation of polymer cure kinetics models

A particle swarm optimisation approach is used to determine the accuracy and experimental relevance of six disparate cure kinetics models. The cure processes of two commercially available thermosetting polymer materials utilised in microelectronics manufacturing applications have been studied using a differential scanning calorimetry system. Numerical models have been fitted to the experimental data using a particle swarm optimisation algorithm which enables the ultimate accuracy of each of the models to be determined. The particle swarm optimisation approach to model fitting proves to be relatively rapid and effective in determining the optimal coefficient set for the cure kinetics models. Results indicate that the singlestep autocatalytic model is able to represent the curing process more accurately than more complex model, with ultimate accuracy likely to be limited by inaccuracies in the processing of the experimental data.

HPLC-ESI-MS method development for the elucidation of nicardipine degradation products

Purpose: Nicardipine is a member of a family of calcium channel blockers named dihydropiridines that are known to be photolabile and may cause phototoxicity. It is therefore vital to develop analytical method which can study the photodegradation of nicardipine. Method: Forced acid degradation of nicardipine was conducted by heating 12 ml of 1 mg/ml nicardipine with 3 ml of 2.5 M HCl for two hours. A gradient HPLC medthod was developed using Agilent Technologies 1200 series quaternary system. Separation was achieved with a Hichrome (250 x 4.6 mm) 5 μm C18 reversed phase column and mobile phase composition of 70% A(100%v/v water) and 30% B(99%v/v acetonitrile + 1%v/v formic acid) at time zero, composition of A and B was then charged to 60%v/v A;40%v/v B at 10minutes, 50%v/v A; 50%v/v B at 30minutes and 70%v/v A; 30%v/v B at 35minutes. 20μl of 0.8mg/ml of nicardipine degradation was injected at room temperature (25oC). The gradient method was transferred onto a HPLC-ESI-MS system (HP 1050 series - AQUAMAX mass detector) and analysis conducted with an acid degradation concentration of 0.25mg/ml and 20μl injection volume. ESI spectra were acquired in positive ionisation mode with MRM 0-600 m/z. Results: Eleven nicardipine degradation products were detected in the HPLC analysis and the resolution (RS) between the respective degradants where 1.0, 1.2, 6.0, 0.4, 1.7, 3.7, 1.8, 1.0, and 1.7 respectively. Nine degradation products were identified in the ESI spectra with the respective m/z ratio; 171.0, 166.1, 441.2, 423.2, 455.2, 455.2, 331.1, 273.1, and 290.1. The possible molecular formulae for each degradants were ambiguously determined. Conclusion: A sensitive and specific method was developed for the analysis of nicardipine degradants. Method enables detection and quantification of nicardipine degradation products that can be used for the study of the kinetics of nicardipine degradation processes.

Composite front maps for improved visibility of dynamic sea-surface features on cloudy SeaWiFS and AVHRR data

Novel techniques have been developed for increasing the value of cloud-affected sequences of Advanced Very High Resolution Radiometer (AVHRR) sea-surface temperature (SST) data and Sea-viewing Wide Field-of-view Sensor (SeaWiFS) ocean colour data for visualising dynamic physical and biological oceanic processes such as fronts, eddies and blooms. The proposed composite front map approach is to combine the location, strength and persistence of all fronts observed over several days into a single map, which allows intuitive interpretation of mesoscale structures. This method achieves a synoptic view without blurring dynamic features, an inherent problem with conventional time-averaging compositing methods. Objective validation confirms a significant improvement in feature visibility on composite maps compared to individual front maps. A further novel aspect is the automated detection of ocean colour fronts, correctly locating 96% of chlorophyll fronts in a test data set. A sizeable data set of 13,000 AVHRR and 1200 SeaWiFS scenes automatically processed using this technique is applied to the study of dynamic processes off the Iberian Peninsula such as mesoscale eddy generation, and many additional applications are identified. Front map animations provide a unique insight into the evolution of upwelling and eddies.