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.

Automated nonlinear system modelling with multiple neural networks

This article discusses the identification of nonlinear dynamic systems using multi-layer perceptrons (MLPs). It focuses on both structure uncertainty and parameter uncertainty, which have been widely explored in the literature of nonlinear system identification. The main contribution is that an integrated analytic framework is proposed for automated neural network structure selection, parameter identification and hysteresis network switching with guaranteed neural identification performance. First, an automated network structure selection procedure is proposed within a fixed time interval for a given network construction criterion. Then, the network parameter updating algorithm is proposed with guaranteed bounded identification error. To cope with structure uncertainty, a hysteresis strategy is proposed to enable neural identifier switching with guaranteed network performance along the switching process. Both theoretic analysis and a simulation example show the efficacy of the proposed method.

Modelling of the sodium complex of a calixarene tetraester in the 1,3-alternate conformation

This paper describes the use of molecular mechanics to model the geometry of the sodium complex of a calix[4] arene tetraester, in the 1,3-alternate conformation 1. Partial charges were assigned to the calixarene on the basis of semi-empirical (AM1, PM3, MNDO, INDO, CNDO and ZINDO) calculations and the binding of the sodium ion to the calixarene was modelled using molecular mechanics. Agreement between the optimised and X-ray structures of the complex was very good. The effect of placing the cation in different starting positions on the energy-minimised geometry of the complex is described.

Efficient modelling and optimisation of hybrid multilayered plates subject to ballistic impact

Among the key challenges present in the modelling and optimisation of composite structures against impact is the computational expense involved in setting up accurate simulations of the impact event and then performing the iterations required to optimise the designs. It is of more interest to find good designs given the limitations of the resources and time available rather than the best possible design. In this paper, low cost but sufficiently accurate finite element (FE) models were generated in LS Dyna for several experimentally characterised materials by semi-automating the modelling process and using existing material models. These models were then used by an optimisation algorithm to generate new hybrid offspring, leading to minimum weight and/or cost designs from a selection of isotropic metals, polymers and orthotropic fibre-reinforced laminates that countered a specified impact threat. Experimental validation of the optimal designs thus identified was then successfully carried out using a single stage gas gun. With sufficient computational hardware, the techniques developed in this pilot study can further utilise fine meshes, equations of state and sophisticated material models, so that optimal hybrid systems can be identified from a wide range of materials, designs and threats.

Developing an integrated hydrograph separation and lumped modelling approach to quantifying hydrological pathways in Irish river catchments

An appreciation of the quantity of streamflow derived from the main hydrological pathways involved in transporting diffuse contaminants is critical when addressing a wide range of water resource management issues. In order to assess hydrological pathway contributions to streams, it is necessary to provide feasible upper and lower bounds for flows in each pathway. An important first step in this process is to provide reliable estimates of the slower responding groundwater pathways and subsequently the quicker overland and interflow pathways. This paper investigates the effectiveness of a multi-faceted approach applying different hydrograph separation techniques, supplemented by lumped hydrological modelling, for calculating the Baseflow Index (BFI), for the development of an integrated approach to hydrograph separation. A semi-distributed, lumped and deterministic rainfall runoff model known as NAM has been applied to ten catchments (ranging from 5 to 699 km2). While this modelling approach is useful as a validation method, NAM itself is also an important tool for investigation. These separation techniques provide a large variation in BFI, a difference of 0.741 predicted for BFI in a catchment with the less reliable fixed and sliding interval methods and local minima turning point methods included. This variation is reduced to 0.167 with these methods omitted. The Boughton and Eckhardt algorithms, while quite subjective in their use, provide quick and easily implemented approaches for obtaining physically realistic hydrograph separations. It is observed that while the different separation techniques give varying BFI values for each of the catchments, a recharge coefficient approach developed in Ireland, when applied in conjunction with the Master recession Curve Tabulation method, predict estimates in agreement with those obtained using the NAM model, and these estimates are also consistent with the study catchments’ geology. These two separation methods, in conjunction with the NAM model, were selected to form an integrated approach to assessing BFI in catchments.

Modelling the chromosphere and transition region of ε Eri (K2 V)

Measurements of ultraviolet line fluxes from Space Telescope Imaging Spectrograph and Far-Ultraviolet Spectroscopic Explorer spectra of the K2-dwarf e Eri are reported. These are used to develop new emission measure distributions and semi-empirical atmospheric models for the chromosphere and lower transition region of the star. These models are the most detailed constructed to date for a main-sequence star other than the Sun. New ionization balance calculations, which account for the effect of finite density on dielectronic recombination rates, are presented for carbon, nitrogen, oxygen and silicon. The results of these calculations are significantly different from the standard Arnaud & Rothenflug ion balance, particularly for alkali-like ions. The new atmospheric models are used to place constraints on possible first ionization potential (FIP)-related abundance variations in the lower atmosphere and to discuss limitations of single-component models for the interpretation of certain optically thick line fluxes. © 2005 RAS.

Understanding macroalgal dispersal in a complex hydrodynamic environment: a combined population genetic and physical modelling approach

Gene flow in macroalgal populations can be strongly influenced by spore or gamete dispersal. This, in turn, is influenced by a convolution of the effects of current flow and specific plant reproductive strategies. Although several studies have demonstrated genetic variability in macroalgal populations over a wide range of spatial scales, the associated current data have generally been poorly resolved spatially and temporally. In this study, we used a combination of population genetic analyses and high-resolution hydrodynamic modelling to investigate potential connectivity between populations of the kelp Laminaria digitata in the Strangford Narrows, a narrow channel characterized by strong currents linking the large semi-enclosed sea lough, Strangford Lough, to the Irish Sea. Levels of genetic structuring based on six microsatellite markers were very low, indicating high levels of gene flow and a pattern of isolation-by-distance, where populations are more likely to exchange migrants with geographically proximal populations, but with occasional long-distance dispersal. This was confirmed by the particle tracking model, which showed that, while the majority of spores settle near the release site, there is potential for dispersal over several kilometres. This combined population genetic and modelling approach suggests that the complex hydrodynamic environment at the entrance to Strangford Lough can facilitate dispersal on a scale exceeding that proposed for L. digitata in particular, and the majority of macroalgae in general. The study demonstrates the potential of integrated physical–biological approaches for the prediction of ecological changes resulting from factors such as anthropogenically induced coastal zone changes.

Experimental and discrete element modelling of cohesive iron ore fines

In this study, the behaviour of iron ore fines with varying levels of adhesion was investigated using a confined compression test and a uniaxial test. The uniaxial test was conducted using the semi-automated uniaxial EPT tester in which the cohesive strength of a bulk solid is evaluated from an unconfined compression test following a period of consolidation to a pre-defined vertical stress. The iron ore fines were also tested by measuring both the vertical and circumferential strains on the cylindrical container walls under vertical loading in a separate confined compression tester - the K0 tester, to determine the lateral pressure ratio. Discrete Element Method simulations of both experiments were carried out and the predictions were compared with the experimental observations. A recently developed DEM contact model for cohesive solids, an Elasto-Plastic Adhesive model, was used. This particle contact model uses hysteretic non-linear loading and unloading paths and an adhesion parameter which is a function of the maximum contact overlap. The model parameters for the simulations are phenomenologically based to reproduce the key bulk characteristics exhibited by the solid. The simulation results show a good agreement in capturing the stress history dependent behaviour depicted by the flow function of the cohesive iron ore fines while also providing a reasonably good match for the lateral pressure ratio observed during the confined compression K0 tests. This demonstrates the potential for the DEM model to be used in the simulation of bulk handling applications.

Interactive management control via analytic hierarchy process (AHP). An empirical study in a public university hospital

Management control in public university hospitals is a challenging task because of continuous changes due to external pressures (e.g. economic pressures, stakeholder focuses and scientific progress) and internal complexities (top management turnover, shared leadership, technological evolution, and researcher oriented mission). Interactive budgeting contributed to improving vertical and horizontal communication between hospital and stakeholders and between different organizational levels. This paper describes an application of Analytic Hierarchy Process (AHP) to enhance interactive budgeting in one of the biggest public university hospital in Italy. AHP improved budget allocation facilitating elicitation and formalization of units' needs. Furthermore, AHP facilitated vertical communication among manager and stakeholders, as it allowed multilevel hierarchical representation of hospital needs, and horizontal communication among staff of the same hospital, as it allowed units' need prioritization and standardization, with a scientific multi-criteria approach, without using complex mathematics. Finally, AHP allowed traceability of a complex decision making processes (as budget allocation), this aspect being of paramount importance in public sectors, where managers are called to respond to many different stakeholders about their choices.

A Review of Numerical Modelling of Wave Energy Converter Arrays

Large-scale commercial exploitation of wave energy is certain to require the deployment of wave energy converters (WECs) in arrays, creating ‘WEC farms’. An understanding of the hydrodynamic interactions in such arrays is essential for determining optimum layouts of WECs, as well as calculating the area of ocean that the farms will require. It is equally important to consider the potential impact of wave farms on the local and distal wave climates and coastal processes; a poor understanding of the resulting environmental impact may hamper progress, as it would make planning consents more difficult to obtain. It is therefore clear that an understanding the interactions between WECs within a farm is vital for the continued development of the wave energy industry.To support WEC farm design, a range of different numerical models have been developed, with both wave phase-resolving and wave phase-averaging models now available. Phase-resolving methods are primarily based on potential flow models and include semi-analytical techniques, boundary element methods and methods involving the mild-slope equations. Phase-averaging methods are all based around spectral wave models, with supra-grid and sub-grid wave farm models available as alternative implementations.The aims, underlying principles, strengths, weaknesses and obtained results of the main numerical methods currently used for modelling wave energy converter arrays are described in this paper, using a common framework. This allows a qualitative comparative analysis of the different methods to be performed at the end of the paper. This includes consideration of the conditions under which the models may be applied, the output of the models and the relationship between array size and computational effort. Guidance for developers is also presented on the most suitable numerical method to use for given aspects of WEC farm design. For instance, certain models are more suitable for studying near-field effects, whilst others are preferable for investigating far-field effects of the WEC farms. Furthermore, the analysis presented in this paper identifies areas in which the numerical modelling of WEC arrays is relatively weak and thus highlights those in which future developments are required.

Characterization of the unfolding of a weak focus and modulus of analytic classification

La thèse présente une description géométrique d’un germe de famille générique déployant un champ de vecteurs réel analytique avec un foyer faible à l’origine et son complexifié : le feuilletage holomorphe singulier associé. On montre que deux germes de telles familles sont orbitalement analytiquement équivalents si et seulement si les germes de familles de difféomorphismes déployant la complexification de leurs fonctions de retour de Poincaré sont conjuguées par une conjugaison analytique réelle. Le “caractère réel” de la famille correspond à sa Z2-équivariance dans R^4, et cela s’exprime comme l’invariance du plan réel sous le flot du système laquelle, à son tour, entraîne que l’expansion asymptotique de la fonction de Poincaré est réelle quand le paramètre est réel. Le pullback du plan réel après éclatement par la projection monoidal standard intersecte le feuilletage en une bande de Möbius réelle. La technique d’éclatement des singularités permet aussi de donner une réponse à la question de la “réalisation” d’un germe de famille déployant un germe de difféomorphisme avec un point fixe de multiplicateur égal à −1 et de codimension un comme application de semi-monodromie d’une famille générique déployant un foyer faible d’ordre un. Afin d’étudier l’espace des orbites de l’application de Poincaré, nous utilisons le point de vue de Glutsyuk, puisque la dynamique est linéarisable auprès des points singuliers : pour les valeurs réels du paramètre, notre démarche, classique, utilise une méthode géométrique, soit un changement de coordonée (coordonée “déroulante”) dans lequel la dynamique devient beaucoup plus simple. Mais le prix à payer est que la géométrie locale du plan complexe ambiante devient une surface de Riemann, sur laquelle deux notions de translation sont définies. Après avoir pris le quotient par le relèvement de la dynamique nous obtenons l’espace des orbites, ce qui s’avère être l’union de trois tores complexes plus les points singuliers (l’espace résultant est non-Hausdorff). Les translations, le caractère réel de l’application de Poincaré et le fait que cette application est un carré relient les différentes composantes du “module de Glutsyuk”. Cette propriété implique donc le fait qu’une seule composante de l’invariant Glutsyuk est indépendante.

Creation and Annihilation of Fluxons in ac‐driven Semi-annular Josephson Junction

A new geometry (semiannular) for Josephson junction has been proposed and theoretical studies have shown that the new geometry is useful for electronic applications [1, 2]. In this work we study the voltage‐current response of the junction with a periodic modulation. The fluxon experiences an oscillating potential in the presence of the ac‐bias which increases the depinning current value. We show that in a system with periodic boundary conditions, average progressive motion of fluxon commences after the amplitude of the ac drive exceeds a certain threshold value. The analytic studies are justified by simulating the equation using finite‐difference method. We observe creation and annihilation of fluxons in semiannular Josephson junction with an ac‐bias in the presence of an external magnetic field.

A Study of Single Stage Semi‐Dry Anaerobic Digestion of Organic Fraction of Municipal Solid Waste

Solid waste generation is a natural consequence of human activity and is increasing along with population growth, urbanization and industrialization. Improper disposal of the huge amount of solid waste seriously affects the environment and contributes to climate change by the release of greenhouse gases. Practicing anaerobic digestion (AD) for the organic fraction of municipal solid waste (OFMSW) can reduce emissions to environment and thereby alleviate the environmental problems together with production of biogas, an energy source, and digestate, a soil amendment. The amenability of substrate for biogasification varies from substrate to substrate and different environmental and operating conditions such as pH, temperature, type and quality of substrate, mixing, retention time etc. Therefore, the purpose of this research work is to develop feasible semi-dry anaerobic digestion process for the treatment of OFMSW from Kerala, India for potential energy recovery and sustainable waste management. This study was carried out in three phases in order to reach the research purpose. In the first phase, batch study of anaerobic digestion of OFMSW was carried out for 100 days at 32°C (mesophilic digestion) for varying substrate concentrations. The aim of this study was to obtain the optimal conditions for biogas production using response surface methodology (RSM). The parameters studied were initial pH, substrate concentration and total organic carbon (TOC). The experimental results showed that the linear model terms of initial pH and substrate concentration and the quadratic model terms of the substrate concentration and TOC had significant individual effect (p < 0.05) on biogas yield. However, there was no interactive effect between these variables (p > 0.05). The optimum conditions for maximizing the biogas yield were a substrate concentration of 99 g/l, an initial pH of 6.5 and TOC of 20.32 g/l. AD of OFMSW with optimized substrate concentration of 99 g/l [Total Solid (TS)-10.5%] is a semi-dry digestion system .Under the optimized condition, the maximum biogas yield was 53.4 L/kg VS (volatile solid).. In the second phase, semi-dry anaerobic digestion of organic solid wastes was conducted for 45 days in a lab-scale batch experiment for substrate concentration of 100 g/l (TS-11.2%) for investigating the start-up performances under thermophilic condition (50°C). The performance of the reactor was evaluated by measuring the daily biogas production and calculating the degradation of total solids and the total volatile solids. The biogas yield at the end of the digestion was 52.9 L/kg VS for the substrate concentration of 100 g/l. About 66.7% of volatile solid degradation was obtained during the digestion. A first order model based on the availability of substrate as the limiting factor was used to perform the kinetic studies of batch anaerobic digestion system. The value of reaction rate constant, k, obtained was 0.0249 day-1. A laboratory bench scale reactor with a capacity of 36.8 litres was designed and fabricated to carry out the continuous anaerobic digestion of OFMSW in the third phase. The purpose of this study was to evaluate the performance of the digester at total solid concentration of 12% (semi-dry) under mesophlic condition (32°C). The digester was operated with different organic loading rates (OLRs) and constant retention time. The performance of the reactor was evaluated using parameters such as pH, volatile fatty acid (VFA), alkalinity, chemical oxygen demand (COD), TOC and ammonia-N as well as biogas yield. During the reactor’s start-up period, the process is stable and there is no inhibition occurred and the average biogas production was 14.7 L/day. The reactor was fed in continuous mode with different OLRs (3.1,4.2 and 5.65 kg VS/m3/d) at constant retention time of 30 days. The highest volatile solid degradation of 65.9%, with specific biogas production of 368 L/kg VS fed was achieved with OLR of 3.1 kg VS/m3/d. Modelling and simulation of anaerobic digestion of OFMSW in continuous operation is done using adapted Anaerobic Digestion Model No 1 (ADM1).The proposed model, which has 34 dynamic state variables, considers both biochemical and physicochemical processes and contains several inhibition factors including three gas components. The number of processes considered is 28. The model is implemented in Matlab® version 7.11.0.584(R2010b). The model based on adapted ADM1 was tested to simulate the behaviour of a bioreactor for the mesophilic anaerobic digestion of OFMSW at OLR of 3.1 kg VS/m3/d. ADM1 showed acceptable simulating results.

Assessing nitrogen dynamics in European ecosystems, integrating measurement and modelling: conclusions

This contribution closes this special issue of Hydrology and Earth System Sciences concerning the assessment of nitrogen dynamics in catchments across Europe within a semi-distributed Integrated Nitrogen model for multiple source assessment in Catchments (INCA). New developments in the understanding of the factors and processes determining the concentrations and loads of nitrogen are outlined. The ability of the INCA model to simulate the hydrological and nitrogen dynamics of different European ecosystems is assessed and the results of the first scenario analyses investigating the impacts of deposition, climatic and land-use change on the nitrogen dynamics are summarised. Consideration is given as to how well the model has performed as a generic too] for describing the nitrogen dynamics of European ecosystems across Arctic, Maritime. Continental and Mediterranean climates, its role in new research initiatives and future research requirements.