Modelling of different measures for improving removal in a stormwater pond

The effect of retrofitting an existing pond on removal efficiency and hydraulic performance was modelled using the commercial software Mike21 and compartmental modelling. The Mike21 model had previously been calibrated on the studied pond. Installation of baffles, the addition of culverts under a causeway and removal of an existing island were all studied as possible improvement measures in the pond. The subsequent effect on hydraulic performance and removal of. suspended solids was then evaluated. Copper, cadmium, BOD, nitrogen and phosphorus removal were,also investigated for that specific improvement measure showing the best results. Outcomes of this study reveal that all measures increase the removal efficiency of suspended solids. The hydraulic efficiency is improved for all cases, except for the case where the island is removed. Compartmental modelling was also used to evaluate hydraulic performance and facilitated a better understanding of the way each of the different measures affected the flow pattern and performance. It was concluded that the installation of baffles is the best of the studied measures resulting in a reduction in the annual load on the receiving lake by approximately 8,000 kg of suspended solids (25% reduction of the annual load), 2 kg of copper (10% reduction of the annual load) and 600 kg of BOD (10% reduction of the annual load).

Modelling oxygen diffusion and cell growth in a porous, vascularising scaffold for soft tissue engineering applications

Soft tissue engineering presents significant challenges compared to other tissue engineering disciplines such as bone, cartilage or skin engineering. The very high cell density in most soft tissues, often combined with large implant dimensions, means that the supply of oxygen is a critical factor in the success or failure of a soft tissue scaffold. A model is presented for oxygen diffusion in a 15-60 mm diameter dome-shaped scaffold fed by a blood vessel loop at its base. This model incorporates simple models for vascular growth, cell migration and the effect of cell density on the effective oxygen diffusivity. The model shows that the dynamic, homogeneous cell seeding method often employed in small-scale applications is not applicable in the case of larger scale scaffolds such as these. Instead, we propose the implantation of a small biopsy of tissue close to a blood supply within the scaffold as a technique more likely to be successful. Crown Copyright (c) 2005 Published by Elsevier Ltd. All rights reserved.

Swash overtopping and sediment overwash on a truncated beach

New experimental laboratory data are presented on swash overtopping and sediment overwash on a truncated beach, approximating the conditions at the crest of a beach berm or inter-tidal ridge-runnel. The experiments provide a measure of the uprush sediment transport rate in the swash zone that is unaffected by the difficulties inherent in deploying instrumentation or sediment trapping techniques at laboratory scale. Overtopping flow volumes are compared with an analytical solution for swash flows as well as a simple numerical model, both of which are restricted to individual swash events. The analytical solution underestimates the overtopping volume by an order of magnitude while the model provides good overall agreement with the data and the reason for this difference is discussed. Modelled flow velocities are input to simple sediment transport formulae appropriate to the swash zone in order to predict the overwash sediment transport rates. Calculations performed with traditional expressions for the wave friction factor tend to underestimate the measured transport. Additional sediment transport calculations using standard total load equations are used to derive an optimum constant wave friction factor of f(w)=0.024. This is in good agreement with a broad range of published field and laboratory data. However, the influence of long waves and irregular wave run-up on the overtopping and overwash remains to be assessed. The good agreement between modelled and measured sediment transport rates suggests that the model provides accurate predictions of the uprush sediment transport rates in the swash zone, which has application in predicting the growth and height of beach berms. (c) 2005 Elsevier B.V. All rights reserved.

Pendant drop formation of shear-thinning and yield stress fluids

A volume-of-fluid numerical method is used to predict the dynamics of shear-thinning liquid drop formation in air from a circular orifice. The validity of the numerical calculation is confirmed for a Newtonian liquid by comparison with experimental measurements. For particular values of Weber number and Froude number, predictions show a more rapid pinch-off, and a reduced number of secondary droplets, with increasing shear-thinning. Also a minimum in the limiting drop length occurs for the smallest Weber number as the zero-shear viscosity is varied. At the highest viscosity, the drop length is reduced due to shear-thinning, whereas at lower viscosities there is little effect of shear-thinning. The evolution of predicted drop shape, drop thickness and length, and the configuration at pinch-off are discussed for shear-thinning drops. The evolution of a drop of Bingham yield stress liquid is also considered as a limiting case. In contrast to the shear-thinning cases, it exhibits a plug flow prior to necking, an almost step-change approach to pinch-off of a torpedo shaped drop following the onset of necking, and a much smaller neck length; no secondary drops are formed. The results demonstrate the potential of the numerical model as a design tool in tailoring the fluid rheology for controlling drop formation behaviour. (c) 2006 Elsevier Inc. All rights reserved.

Dynamic rupture in a 3-D particle-based simulation of a rough planar fault

An appreciation of the physical mechanisms which cause observed seismicity complexity is fundamental to the understanding of the temporal behaviour of faults and single slip events. Numerical simulation of fault slip can provide insights into fault processes by allowing exploration of parameter spaces which influence microscopic and macroscopic physics of processes which may lead towards an answer to those questions. Particle-based models such as the Lattice Solid Model have been used previously for the simulation of stick-slip dynamics of faults, although mainly in two dimensions. Recent increases in the power of computers and the ability to use the power of parallel computer systems have made it possible to extend particle-based fault simulations to three dimensions. In this paper a particle-based numerical model of a rough planar fault embedded between two elastic blocks in three dimensions is presented. A very simple friction law without any rate dependency and no spatial heterogeneity in the intrinsic coefficient of friction is used in the model. To simulate earthquake dynamics the model is sheared in a direction parallel to the fault plane with a constant velocity at the driving edges. Spontaneous slip occurs on the fault when the shear stress is large enough to overcome the frictional forces on the fault. Slip events with a wide range of event sizes are observed. Investigation of the temporal evolution and spatial distribution of slip during each event shows a high degree of variability between the events. In some of the larger events highly complex slip patterns are observed.

Measurement and modelling of controlled beach groundwater levels under wave action

New laboratory scale experimental data are presented on the forcing of beach groundwater levels by wave run-up. The experimental setup simulates a coastal barrier dividing the ocean from a relatively constant back beach water level, conditions approximating a closed off lagoon system or beach aquifer. The data are critically compared to an advanced numerical model for simulating wave and beach groundwater interaction in the coastal zone, and provide the first experimental verification of such a model. Overall model-data comparisons are good, but some systematic discrepancies are apparent, and reasons for these are discussed.

Indagine modellistica sulla tenuta idraulica delle difese arginali deputate alla protezione idraulica del territorio

La prima parte di questo lavoro di tesi tratta dell’interazione tra un bacino di laminazione e il sottostante acquifero: è in fase di progetto, infatti, la costruzione di una cassa di espansione sul torrente Baganza, a monte della città di Parma. L’obiettivo di tale intervento è di ridurre il rischio di esondazione immagazzinando temporaneamente, in un serbatoio artificiale, la parte più pericolosa del volume di piena che verrebbe rilasciata successivamente con portate che possono essere agevolmente contenute nel tratto cittadino del torrente. L’acquifero è stato preliminarmente indagato e monitorato permettendone la caratterizzazione litostratigrafica. La stratigrafia si può riassumere in una sequenza di strati ghiaioso-sabbiosi con successione di lenti d’argilla più o meno spesse e continue, distinguendo due acquiferi differenti (uno freatico ed uno confinato). Nel presente studio si fa riferimento al solo acquifero superficiale che è stato modellato numericamente, alle differenze finite, per mezzo del software MODFLOW_2005. L'obiettivo del presente lavoro è di rappresentare il sistema acquifero nelle condizioni attuali (in assenza di alcuna opera) e di progetto. La calibrazione è stata condotta in condizioni stazionarie utilizzando i livelli piezometrici raccolti nei punti d’osservazione durante la primavera del 2013. I valori di conducibilità idraulica sono stati stimati per mezzo di un approccio geostatistico Bayesiano. Il codice utilizzato per la stima è il bgaPEST, un software gratuito per la soluzione di problemi inversi fortemente parametrizzati, sviluppato sulla base dei protocolli del software PEST. La metodologia inversa stima il campo di conducibilità idraulica combinando osservazioni sullo stato del sistema (livelli piezometrici nel caso in esame) e informazioni a-priori sulla struttura dei parametri incogniti. La procedura inversa richiede il calcolo della sensitività di ciascuna osservazione a ciascuno dei parametri stimati; questa è stata valutata in maniera efficiente facendo ricorso ad una formulazione agli stati aggiunti del codice in avanti MODFLOW_2005_Adjoint. I risultati della metodologia sono coerenti con la natura alluvionale dell'acquifero indagato e con le informazioni raccolte nei punti di osservazione. Il modello calibrato può quindi essere utilizzato come supporto alla progettazione e gestione dell’opera di laminazione. La seconda parte di questa tesi tratta l'analisi delle sollecitazioni indotte dai percorsi di flusso preferenziali causati da fenomeni di piping all’interno dei rilevati arginali. Tali percorsi preferenziali possono essere dovuti alla presenza di gallerie scavate da animali selvatici. Questo studio è stato ispirato dal crollo del rilevato arginale del Fiume Secchia (Modena), che si è verificato in gennaio 2014 a seguito di un evento alluvionale, durante il quale il livello dell'acqua non ha mai raggiunto la sommità arginale. La commissione scientifica, la cui relazione finale fornisce i dati utilizzati per questo studio, ha attribuito, con molta probabilità, il crollo del rilevato alla presenza di tane di animali. Con lo scopo di analizzare il comportamento del rilevato in condizioni integre e in condizioni modificate dall'esistenza di un tunnel che attraversa il manufatto arginale, è stato realizzato un modello numerico 3D dell’argine mediante i noti software Femwater e Feflow. I modelli descrivono le infiltrazioni all'interno del rilevato considerando il terreno in entrambe le porzioni sature ed insature, adottando la tecnica agli elementi finiti. La tana è stata rappresentata da elementi con elevata permeabilità e porosità, i cui valori sono stati modificati al fine di valutare le diverse influenze sui flussi e sui contenuti idrici. Per valutare se le situazioni analizzate presentino o meno il verificarsi del fenomeno di erosione, sono stati calcolati i valori del fattore di sicurezza. Questo è stato valutato in differenti modi, tra cui quello recentemente proposto da Richards e Reddy (2014), che si riferisce al criterio di energia cinetica critica. In ultima analisi è stato utilizzato il modello di Bonelli (2007) per calcolare il tempo di erosione ed il tempo rimanente al collasso del rilevato.

Experiments to assess the transport of fibrous insulation debris

Loss of coolant accidents (LOCA) in the primary cooling circuit of a nuclear reactor may result in damage to insulation materials that are located near to the leak. The insulation materials released may compromise the operation of the emergency core cooling system (ECCS). Insulation material in the form of mineral wool fibre agglomerates (MWFA) maybe transported to the containment sump strainers mounted at the inlet of the emergency cooling pumps, where the insulation fibres may block or penetrate the strainers. In addition to the impact of MWFA on the pressure drop across the strainers, corrosion products formed over time may also accumulate in the fibre cakes on the strainers, which can lead to a significant increase in the strainer pressure drop and result in cavitation in the ECCS. Thus, knowledge of transport characteristics of the damaged insulation materials in various scenarios is required to help plan for the long-term operability of nuclear reactors, which undergo LOCA. An experimental and theoretical study performed by the Helmholtz-Zentrum Dresden-Rossendorf and the Hochschule Zittau/Görlitz1 is investigating the phenomena that maybe observed in the containment vessel during a LOCA. The study entails the generation of fibre agglomerates, the determination of their transport properties in single and multi-effect experiments and the long-term effect that corrosion of the containment internals by the coolant has on the strainer pressure drop. The focus of this presentation is on the experiments performed that characterize the horizontal transport of MWFA, whereas the corresponding CFD simulations are described in an accompanying contribution (see abstract of Cartland Glover et al.). The experiments were performed a racetrack type channel that provided a near uniform horizontal flow. The channel is 0.1 wide by 1.2 m high with a straight length of 5 m and two bends of 0.5 m. The measurement techniques include particle imaging (both wide-angle and macro lens), concurrent particle image velocimetry, ultravelocimetry, laser detection sensors to sense the presence of absence of MWFA and pertinent measurements of the MWFA concentration and quiescent settling characteristics. The transport of the MWFA was observed at velocities of 0.1 and 0.25 m s-1 to verify numerical model behaviour in and just beyond expected velocities in the containment sump of a nuclear reactor.

Nondestructive index profiling of long period gratings in photonic crystal fibres

A numerical model of a long period grating in photonic crystal fibre fabricated by an electric arc is proposed that allows for the spectral characterisation of the grating. In the combination with the suggested model of the photonic crystal and the experimentally recorded grating growth it is used to find the index change induced by the electric arc.

Distributed temperature sensing using a chirped fibre Bragg grating

A fully distributed temperature sensor consisting of a chirped fibre Bragg grating has been demonstrated. By fitting a numerical model of the grating response including temperature change, position and width of localized heating applied to the grating, we achieve measurements of these parameters to within 2.2 K, 149 µm and 306 µm of applied values, respectively. Assuming that deviation from linearity is accounted for in making measurement, much higher precision is achievable and the standard deviations for these measurements are 0.6 K, 28.5 µm and 56.0 µm, respectively.

Hydrodynamic modeling of mineral wool fiber suspensions in a two-dimensional channel flow

A consequence of a loss of coolant accident is the damage of adjacent insulation materials (IM). IM may then be transported to the containment sump strainers where water is drawn into the ECCS (emergency core cooling system). Blockage of the strainers by IM lead to an increased pressure drop acting on the operating ECCS pumps. IM can also penetrate the strainers, enter the reactor coolant system and then accumulate in the reactor pressure vessel. An experimental and theoretical study that concentrates on mineral wool fiber transport in the containment sump and the ECCS is being performed. The study entails fiber generation and the assessment of fiber transport in single and multi-effect experiments. The experiments include measurement of the terminal settling velocity, the strainer pressure drop, fiber sedimentation and resuspension in a channel flow and jet flow in a rectangular tank. An integrated test facility is also operated to assess the compounded effects. Each experimental facility is used to provide data for the validation of equivalent computational fluid dynamic models. The channel flow facility allows the determination of the steady state distribution of the fibers at different flow velocities. The fibers are modeled in the Eulerian-Eulerian reference frame as spherical wetted agglomerates. The fiber agglomerate size, density, the relative viscosity of the fluid-fiber mixture and the turbulent dispersion of the fibers all affect the steady state accumulation of fibers at the channel base. In the current simulations, two fiber phases are separately considered. The particle size is kept constant while the density is modified, which affects both the terminal velocity and volume fraction. The relative viscosity is only significant at higher concentrations. The numerical model finds that the fibers accumulate at the channel base even at high velocities; therefore, modifications to the drag and turbulent dispersion forces can be made to reduce fiber accumulation.

Development of multiphase and multiscale mathematical models for thermal spray process

High velocity oxyfuel (HVOF) thermal spraying is one of the most significant developments in the thermal spray industry since the development of the original plasma spray technique. The first investigation deals with the combustion and discrete particle models within the general purpose commercial CFD code FLUENT to solve the combustion of kerosene and couple the motion of fuel droplets with the gas flow dynamics in a Lagrangian fashion. The effects of liquid fuel droplets on the thermodynamics of the combusting gas flow are examined thoroughly showing that combustion process of kerosene is independent on the initial fuel droplet sizes. The second analysis copes with the full water cooling numerical model, which can assist on thermal performance optimisation or to determine the best method for heat removal without the cost of building physical prototypes. The numerical results indicate that the water flow rate and direction has noticeable influence on the cooling efficiency but no noticeable effect on the gas flow dynamics within the thermal spraying gun. The third investigation deals with the development and implementation of discrete phase particle models. The results indicate that most powder particles are not melted upon hitting the substrate to be coated. The oxidation model confirms that HVOF guns can produce metallic coating with low oxidation within the typical standing-off distance about 30cm. Physical properties such as porosity, microstructure, surface roughness and adhesion strength of coatings produced by droplet deposition in a thermal spray process are determined to a large extent by the dynamics of deformation and solidification of the particles impinging on the substrate. Therefore, is one of the objectives of this study to present a complete numerical model of droplet impact and solidification. The modelling results show that solidification of droplets is significantly affected by the thermal contact resistance/substrate surface roughness.

Distributed temperature sensing using a chirped fibre Bragg grating

A fully distributed temperature sensor consisting of a chirped fibre Bragg grating has been demonstrated. By fitting a numerical model of the grating response including temperature change, position and width of localized heating applied to the grating, we achieve measurements of these parameters to within 2.2 K, 149 μm and 306 μm of applied values, respectively. Assuming that deviation from linearity is accounted for in making measurement, much higher precision is achievable and the standard deviations for these measurements are 0.6 K, 28.5 μm and 56.0 μm, respectively. © 2004 IOP Publishing Ltd.

Nondestructive index profiling of long period gratings in photonic crystal fibres

A numerical model of a long period grating in photonic crystal fibre fabricated by an electric arc is proposed that allows for the spectral characterisation of the grating. In the combination with the suggested model of the photonic crystal and the experimentally recorded grating growth it is used to find the index change induced by the electric arc. © Springer Science+Business Media B.V. 2007.

Modelling vortex formation in an unbaffled stirred tank reactors

Agitating liquids in unbaffled stirred tank leads to the formation of a vortex in the region of the impeller shaft when operating in the turbulent flow regime. A numerical model is presented here that captures such a vortex. The volume of fluid model, a multiphase flow model was employed in conjunction with a multiple reference frame model and the shear stress turbulence model. The dimensions of the tank considered here, were 0.585 m for the liquid depth and tank diameter with a 0.2925 m diameter impeller at a height of 0.2925 m. The impeller considered was an eight-bladed paddle type agitator that was rotating with an angular velocity of 7.54 rad s (72 rpm) giving a Reynolds number of 10 and Froude number of 0.043. Preliminary results of a second investigation into the effect of liquid phase properties on the vortex formed are also presented. © 2006 Elsevier B.V. All rights reserved.