Modelling of saline intrusion in marine outfalls

Since their introduction in the 1950s, marine outfalls with diffusers have been prone to saline intrusion, a process in which seawater ingresses into the outfall. This can greatly reduce the dilution and subsequent dispersion of wastewater discharged, sometimes resulting in serious deterioration of coastal water quality. Although long aware of the difficulties posed by saline intrusion, engineers still lack satisfactory methods for its prediction and robust design methods for its alleviation. However, with recent developments in numerical methods and computer power, it has been suggested that commercially available computational fluid dynamics (CFD) software may be a useful aid in combating this phenomenon by improving understanding through synthesising likely behaviour. This document reviews current knowledge on saline intrusion and its implications and then outlines a model-scale investigation of the process undertaken at Queen's University Belfast, using both physical and CFD methods. Results are presented for a simple outfall configuration, incorporating several outlets. The features observed agree with general observations from full-scale marine outfalls, and quantify the intricate internal flow mechanisms associated with saline intrusion. The two-dimensional numerical model was found to represent saline intrusion, but in a qualitative manner, not yet adequate for design purposes. Specific areas requiring further development were identified. The ultimate aim is to provide a reliable, practical and cost effective means by which engineers can minimise saline intrusion through optimised outfall design.

Numerical simulation of separation factor of NaCl in surface force pore flow model

This work deals with modelling and experimental verification of desalination theory (surface force pore flow) . The work has direct application in desalination of sea water.

Numerical analysis of bodyworn UHF antenna systems

Bodyworn antennas are found in a wide range of medical, military and personal communication applications, yet reliable communication from the surface of the human body still presents a range of engineering challenges. At UHF and microwave frequencies, bodyworn antennas can suffer from reduced efficiency due to electromagnetic absorption in tissue, radiation pattern fragmentation and variations in feed-point impedance. The significance and nature of these effects are system specific and depend on the operating frequency, propagation environment and physical constraints on the antenna itself. This paper describes how numerical electromagnetic modelling techniques such as FDTD (finite-difference time-domain) can be used in the design of bodyworn antennas. Examples are presented for 418 MHz, 916 .5 MHz and 2 . 45 GHz, in the context of both biomedical signalling and wireless personal-area networking applications such as the Bluetooth(TM)* wireless technology.

Diagnostic based modelling of radio-frequency driven atmospheric pressure plasmas

Diagnostic based modelling (DBM) actively combines complementary advantages of numerical plasma simulations and relatively simple optical emission spectroscopy (OES). DBM is employed to determine absolute atomic oxygen ground state densities in a helium–oxygen radio-frequency driven atmospheric pressure plasma jet. A comparatively simple one-dimensional simulation yields detailed information on electron properties governing the population dynamics of excited states. Important characteristics of the electron dynamics are found to be largely insensitive to details of the chemical composition and to be in very good agreement with space and phase-resolved OES. Benchmarking the time and space resolved simulation allows us to subsequently derive effective excitation rates as the basis for DBM with simple space and time integrated OES. The population dynamics of the upper O 3p 3P (? = 844 nm) atomic oxygen state is governed by direct electron impact excitation, dissociative excitation, radiation losses and collisional induced quenching. Absolute values for atomic oxygen densities are obtained through tracer comparison with the upper Ar 2p1 (? = 750.4 nm) state. The presented results for the atomic oxygen density show excellent quantitative agreement with independent two-photon laser-induced fluorescence measurements.

The challenge of revealing and tailoring the dynamics of radio-frequency plasmas

The complex dynamics of ionization and excitation mechanisms in capacitively coupled radio-frequency plasmas is discussed for single- and dual-frequency operations in low-pressure and atmospheric pressure plasmas. Electrons are energized through the dynamics of electric fields in the vicinity of the plasma boundary sheaths. Distinctly different power dissipation mechanisms can either co-exist or initiate mode transitions exhibiting characteristic spatio-temporal ionization structures. Phase resolved optical emission spectroscopy, in combination with adequate modelling of the population dynamics of excited states, and numerical simulations reveal dissipation associated with sheath expansion, sheath collapse, transient electron avalanches and wave–particle interactions. In dual-frequency systems the relative phase between the two frequency components provides additional strategies to tailor the plasma dynamics.

Modelling the mechanical response of the Reed-mouthpiece-lip system of a clarinet. Part I. A one-dimensional distributed model

The motion of a clarinet reed that is clamped to a mouthpiece and supported by a lip is simulated in the time-domain using finite difference methods. The reed is modelled as a bar with non-uniform cross section, and is described using a one-dimensional, fourth-order partial differential equation. The interactions with the mouthpiece Jay and the player's lip are taken into account by incorporating conditional contact forces in the bar equation. The model is completed by clamped-free boundary conditions for the reed. An implicit finite difference method is used for discretising the system, and values for the physical parameters are chosen both from laboratory measurements and by accurate tuning of the numerical simulations. The accuracy of the numerical system is assessed through analysis of frequency warping effects and of resonance estimation. Finally, the mechanical properties of the system are studied by analysing its response to external driving forces. In particular, the effects of reed curling are investigated.

Modelling the mechanical response of the reed-mouthpiece-lip system of a clarinet. Part II: A lumped model approximation

A non-linear lumped model of the reed-mouthpiece-lip system of a clarinet is formulated, in which the lumped parameters are derived from numerical experiments with a finite-difference simulation based on a distributed reed model. The effective stiffness per unit area is formulated as a function of the pressure signal driving the reed, in order to simulate the effects of the reed bending against the lay, and mass and damping terms are added as a first approximation to the dynamic behaviour of the reed. A discrete-time formulation is presented, and its response is compared to that of the distributed model. In addition, the lumped model is applied in the simulation of clarinet tones, enabling the analysis of the effects of using a pressure-dependent stiffness per unit area on sustained oscillations. The analysed effects and features are in qualitative agreement with players' experiences and experimental results obtained in prior studies.

Formulation of Locally Reacting Surfaces in FDTD/K-DWM Modelling of Acoustic Spaces

In this paper, we present new methods for constructing and analysing formulations of locally reacting surfaces that can be used in finite difference time domain (FDTD) simulations of acoustic spaces. Novel FDTD formulations of frequency-independent and simple frequency-dependent impedance boundaries are proposed for 2D and 3D acoustic systems, including a full treatment of corners and boundary edges. The proposed boundary formulations are designed for virtual acoustics applications using the standard leapfrog scheme based on a rectilinear grid, and apply to FDTD as well as Kirchhoff variable digital waveguide mesh (K-DWM) methods. In addition, new analytic evaluation methods that accurately predict the reflectance of numerical boundary formulations are proposed. numerical experiments and numerical boundary analysis (NBA) are analysed in time and frequency domains in terms of the pressure wave reflectance for different angles of incidence and various impedances. The results show that the proposed boundary formulations structurally adhere well to the theoretical reflectance. In particular, both reflectance magnitude and phase are closely approximated even at high angles of incidence and low impedances. Furthermore, excellent agreement was found between the numerical boundary analysis and the experimental results, validating both as tools for researching FDTD boundary formulations.

A numerical investigation of the flow structures and losses for turbulent flow in 90 degrees elbow bends

Computational fluid dynamic modelling was carried out on a series of pipe bends having R/r values of 1.3, 5, and 20, with the purpose of determining the accuracy of numerical models in predicting pressure loss data from which to inform one-dimensional loss models. Four separate turbulence models were studied: the standard k-epsilon model, realizable k-epsilon model, k-omega model, and a Reynolds stress model (RSM). The results are presented for each bend in the form of upstream and downstream pressure profiles, pressure distributions along the inner and outer walls, detailed pressure and velocity fields as well as overall loss values. In each case, measured data were presented to evaluate the predictive ability of each model. The RSM was found to perform the best, producing accurate pressure loss data for bends with R/r values of 5 and 20. For the tightest bend with an R/r value of 1.3, however, predictions were significantly worse due to the presence of flow separation, stronger pressure gradients, and high streamline curvature.

Modelling and Control of Shape Memory Alloy Actuators by using Preisach Model, Genetic Algorithm and Fuzzy Logic

Shapememoryalloy (SMA) actuators, which have the ability to return to a predetermined shape when heated, have many potential applications in aeronautics, surgical tools, robotics and so on. Nonlinearity hysteresis effects existing in SMA actuators present a problem in the motion control of these smart actuators. This paper investigates the control problem of SMA actuators in both simulation and experiment. In the simulation, the numerical Preisachmodel with geometrical interpretation is used for hysteresis modeling of SMA actuators. This model is then incorporated in a closed loop PID control strategy. The optimal values of PID parameters are determined by using geneticalgorithm to minimize the mean squared error between desired output displacement and simulated output. However, the control performance is not good compared with the simulation results when these parameters are applied to the real SMA control since the system is disturbed by unknown factors and changes in the surrounding environment of the system. A further automated readjustment of the PID parameters using fuzzylogic is proposed for compensating the limitation. To demonstrate the effectiveness of the proposed controller, real time control experiment results are presented.

Modelling road traffic accidents using macroscopic second-order models of traffic flow

In this paper, we introduce a macroscopic model for road traffic accidents along highway sections. We discuss the motivation and the derivation of such a model, and we present its mathematical properties. The results are presented by means of examples where a section of a crowded one-way highway contains in the middle a cluster of drivers whose dynamics are prone to road traffic accidents. We discuss the coupling conditions and present some existence results of weak solutions to the associated Riemann Problems. Furthermore, we illustrate some features of the proposed model through some numerical simulations. © The authors 2012.

Modelling solute and microorganism transport through a peri-alpine gravel aquifer by considering the aquifer as a dual porosity fractured medium.

Approach:
In-situ passive gradient comparative artificial tracer testing, undertaken using solutes (Uranine and Iodide), Bacteria (E.coli and P.putida) and bacteriophage (H40/1), permitted comparison of the mobility of different sized microorganisms relative to solutes in the sand and gravel aquifer underlying Dornach, Germany.
Tracer breakthrough curves reveal that even though uranine initially arrived at observation wells at the same time as microbiological tracers, maximum relative concentrations were sometimes less than those of microbiological tracers, while solute breakthrough curves proved more disperse.
Monitoring uranine breakthrough with depth suggested tracers arrived in observation wells in discrete 0.5m-1m thick intervals, over the aquifer’s 12m saturated thickness. Nearby exposures of aquifer material suggested that the aquifer consisted of sandy gravels enveloping sequences of open framework (OW) gravel up to 1m thick. Detailed examination of OW units revealed that they contained lenses of silty sand up to 1m long x 30cm thick., while granulometric data suggested that the gravel was two to three orders of magnitude more permeable than the enveloping sandy gravel.
Solute and microorganism tracer responses could not be simulated using conventional advective-dispersive equation solutions employing the same velocity and dispersion terms. By contrast solute tracer responses, modelled using a dual porosity approach for fractured media (DP-1D) corresponded well to observed field data. Simulating microorganism responses using the same transport terms, but no dual porosity term, generated good model fits and explained the higher relative concentration of the bacteria, compared to the non-reactive solute, even with first order removal to account for lower RR. Geologically, model results indicate that the silty units within open framework gravels are accessible to solute tracers, but not to microorganisms.
Importance:
Results highlight the benefits of geological observations developing appropriate conceptual models of solute and micro organism transport and in developing suitable numerical approaches to quantifying microorganism mobility at scales appropriate for the development of groundwater supply (wellhead) protection zones.

Analytical approximations to numerical solutions of theoretical emission measure distributions

Emission line fluxes from cool stars are widely used to establish an apparent emission measure distribution, EmdApp(Te), between temperatures characteristic of the low transition region and the low corona. The true emission measure distribution, EmdTrue(Te), is determined by the energy balance and geometry adopted and, with a numerical model, can be used to predict EmdApp(Te), to guide further modelling. The scaling laws that exist between coronal parameters arise from the dimensions of the terms in the energy balance equation. Here, analytical approximations to numerical solutions for EmdTrue(Te) are presented, which show how the constants in the coronal scaling laws are determined. The apparent emission measure distributions show a minimum value at some T0 and a maximum at the mean coronal temperature Tc (although in some stars, emission from active regions can contribute). It is shown that, for the energy balance and geometry adopted, the analytical values of the emission measure and electron pressure at T0 and Tc depend on only three parameters: the stellar surface gravity and the values of T0 and Tc. The results are tested against full numerical solutions for e Eri (K2 V) and are applied to Procyon (a CMi, F5 IV/V). The analytical approximations can be used to restrict the required range of full numerical solutions, to check the assumed geometry and to show where the adopted energy balance may not be appropriate. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.