An Energy Conserving Finite Difference Scheme for Simulation of Collisions

Nonlinear phenomena play an essential role in the sound production process of many musical instruments. A common source of these effects is object collision, the numerical simulation of which is known to give rise to stability
issues. This paper presents a method to construct numerical schemes that conserve the total energy in simulations of one-mass systems involving collisions, with no conditions imposed on any of the physical or numerical parameters.
This facilitates the adaptation of numerical models to experimental data, and allows a more free parameter adjustment in sound synthesis explorations. The energy preservedness of the proposed method is tested and demonstrated though several examples, including a bouncing ball and a non-linear oscillator, and implications regarding the wider applicability are discussed.

Evaluation of a substructured soft Real-time Hybrid Test for performing Seismic analysis of complex structural systems

The hybrid test method is a relatively recently developed dynamic testing technique that uses numerical modelling combined with simultaneous physical testing. The concept of substructuring allows the critical or highly nonlinear part of the structure that is difficult to numerically model with accuracy to be physically tested whilst the remainder of the structure, that has a more predictable response, is numerically modelled. In this paper, a substructured soft-real time hybrid test is evaluated as an accurate means of performing seismic tests of complex structures. The structure analysed is a three-storey, two-by-one bay concentrically braced frame (CBF) steel structure subjected to seismic excitation. A ground storey braced frame substructure whose response is critical to the overall response of the structure is tested, whilst the remainder of the structure is numerically modelled. OpenSees is used for numerical modelling and OpenFresco is used for the communication between the test equipment and numerical model. A novel approach using OpenFresco to define the complex numerical substructure of an X-braced frame within a hybrid test is also presented. The results of the hybrid tests are compared to purely numerical models using OpenSees and a simulated test using a combination of OpenSees and OpenFresco. The comparative results indicate that the test method provides an accurate and cost effective procedure for performing
full scale seismic tests of complex structural systems.

Integrating geophysical data and upscaling techniques to support regional groundwater flow modelling: A practical example of the Lagan Valley, Northern Ireland

When studying heterogeneous aquifer systems, especially at regional scale, a degree of generalization is anticipated. This can be due to sparse sampling regimes, complex depositional environments or lack of accessibility to measure the subsurface. This can lead to an inaccurate conceptualization which can be detrimental when applied to groundwater flow models. It is important that numerical models are based on observed and accurate geological information and do not rely on the distribution of artificial aquifer properties. This can still be problematic as data will be modelled at a different scale to which it was collected. It is proposed here that integrating geophysics and upscaling techniques can assist in a more realistic and deterministic groundwater flow model. In this study, the sedimentary aquifer of the Lagan Valley in Northern Ireland is chosen due to intruding sub-vertical dolerite dykes. These dykes are of a lower permeability than the sandstone aquifer. The use of airborne magnetics allows the delineation of heterogeneities, confirmed by field analysis. Permeability measured at the field scale is then upscaled to different levels using a correlation with the geophysical data, creating equivalent parameters that can be directly imported into numerical groundwater flow models. These parameters include directional equivalent permeabilities and anisotropy. Several stages of upscaling are modelled in finite element. Initial modelling is providing promising results, especially at the intermediate scale, suggesting an accurate distribution of aquifer properties. This deterministic based methodology is being expanded to include stochastic methods of obtaining heterogeneity location based on airborne geophysical data. This is through the Direct Sample method of Multiple-Point Statistics (MPS). This method uses the magnetics as a training image to computationally determine a probabilistic occurrence of heterogeneity. There is also a need to apply the method to alternate geological contexts where the heterogeneity is of a higher permeability than the host rock.

Type Ia supernova bolometric light curves and ejected mass estimates from the nearby supernova factory

We present a sample of normal Type Ia supernovae (SNe Ia) from the Nearby Supernova Factory data set with spectrophotometry at sufficiently late phases to estimate the ejected mass using the bolometric light curve.Wemeasure Ni masses from the peak bolometric luminosity, then compare the luminosity in the Co-decay tail to the expected rate of radioactive energy release from ejecta of a given mass. We infer the ejected mass in a Bayesian context using a semi-analytic model of the ejecta, incorporating constraints from contemporary numerical models as priors on the density structure and distribution of Ni throughout the ejecta. We find a strong correlation between ejected mass and light-curve decline rate, and consequently Ni mass, with ejected masses in our data ranging from 0.9 to 1.4 M. Most fast-declining (SALT2 x <-1) normal SNe Ia have significantly sub-Chandrasekhar ejected masses in our fiducial analysis.

Development of dynamic structural testing to current state-of-the-art hybrid testing

Even though computational power used for structural analysis is ever increasing, there is still a fundamental need for testing in structural engineering, either for validation of complex numerical models or to assess material behaviour. In addition to analysis of structures using scale models, many structural engineers are aware to some extent of cyclic and shake-table test methods, but less so of ‘hybrid testing’. The latter is a combination of physical testing (e.g. hydraulic
actuators) and computational modelling (e.g. finite element modelling). Over the past 40 years, hybrid testing of engineering structures has developed from concept through to maturity to become a reliable and accurate dynamic testing technique. The hybrid test method provides users with some additional benefits that standard dynamic testing methods do not, and the method is more cost-effective in comparison to shake-table testing. This article aims to provide the reader with a basic understanding of the hybrid test method, including its contextual development and potential as a dynamic testing technique.

Design of Scaled Bridge Prototype for Energy Harvesting Testing

Energy harvesting from ambient vibration is a promising field, especially for applications in larger infrastructures such as bridges. These structures are more frequently monitored for damage detection because of their extended life, increased traffic load and environmental deterioration. In this regard, the possibility of sourcing the power necessary for the sensors from devices embedded in the structure, thus cutting the cost due to the management of battery replacing over the lifespan of the structure, is particularly attracting. Among others, piezoelectric devices have proven to be especially effective and easy to apply since they can be bonded to existing host structure. For these devices the energy harvesting capacity is achieved directly from the variation in the strain conditions from the surface of the structure. However these systems need to undergo significant research for optimisation of their harvesting capacity and for assessing the feasibility of application to various ranges of bridge span and load. In this regard scaled bridge prototypes can be effectively used not only to assess numerical models and studies in an inexpensive and repeatable way but also to test the electronic devices under realistic field conditions. In this paper the theory of physical similitude is applied to the design of bridge beams with embedded energy harvesting systems and health monitoring sensors. It will show both how bridge beams can be scaled in such a way to apply and test energy harvesting systems and 2) how experimental data from existing bridges can be applied to prototypes in a laboratory environment. The study will be used for assessing the reliability of the system over a train bridge case study undergoing a set load cycles and induced localised damage.

Experimental Investigation of Transient Saltwater Intrusion in Heterogeneous Porous Media

A 2D sandbox style experiment was developed to compare the results of numerical modelling to physical testing for saltwater intrusion in homogeneous and heterogeneous aquifers. The sandbox consisted of a thin central viewing chamber filled with glass beads of varying diameters (780μm, 1090μm and 1325μm) under fully saturated conditions. Dyed saltwater (SW) was introduced at the side boundary and a head difference imposed across the porous media. Images of the SW wedge were recorded at intervals in order to assess the suitability of the numerical models predictions of transient SW intrusion. Numerical modelling of the experimental cases were simulated using SUTRA. Two main parameters were chosen to express the condition of the intruding SW wedge at each recorded time step; the toe penetration length (TL) and the width of the mixing zone (WMZ). The WMZ was larger under transient conditions in the heterogeneous case, while the TL was longer for the homogeneous case. The increased variability in the flow field fo the heterogeneous case resulted in increased dispersion, and thus, increased WMZ.

Overview of laser-driven generation of electron–positron beams

Electron–positron (e–p) plasmas are widely thought to be emitted, in the form of ultra-relativistic winds or collimated jets, by some of the most energetic or powerful objects in the Universe, such as black-holes, pulsars, and quasars. These phenomena represent an unmatched astrophysical laboratory to test physics at its limit and, given their immense distance from Earth (some even farther than several billion light years), they also provide a unique window on the very early stages of our Universe. However, due to such gigantic distances, their properties are only inferred from the indirect interpretation of their radiative signatures and from matching numerical models: their generation mechanism and dynamics still pose complicated enigmas to the scientific community. Small-scale reproductions in the laboratory would represent a fundamental step towards a deeper understanding of this exotic state of matter. Here we present recent experimental results concerning the laser-driven production of ultra-relativistic e–p beams. In particular, we focus on the possibility of generating beams that present charge neutrality and that allow for collective effects in their dynamics, necessary ingredients for the testing pair-plasma physics in the laboratory. A brief discussion of the analytical and numerical modelling of the dynamics of these plasmas is also presented in order to provide a summary of the novel plasma physics that can be accessed with these objects. Finally, general considerations on the scalability of laboratory plasmas up to astrophysical scenarios are given.

Predicting the crushing behaviour of composite material using high-fidelity finite element modelling

The capability to numerically model the crushing behaviour of composite structures will enable the efficient design of structures with high specific energy absorption capacity. This is particularly relevant to the aerospace and automotive industries where cabin structures need to be shown to be crashworthy. In this paper, a three-dimensional damage model is presented, which accurately represents the behaviour of composite laminates under crush loading. Both intralaminar and interlaminar failure mechanisms are taken into account. The crush damage model was implemented in ABAQUS/Explicit as a VUMAT subroutine. Numerical predictions are shown to agree well with experimental results, accurately capturing the intralaminar and interlaminar damage for a range of stacking sequences, triggers and composite materials. The use of measured material parameters required by the numerical models, without the need to ‘calibrate’ this input data, demonstrates this computational tool's predictive capabilities

Overview of Seismic Hybrid Testing

Even though computational power used for structural analysis is ever increasing, there is still a fundamental need for testing in structural engineering, either for validation of complex numerical models or material behaviour. Many structural engineers/researchers are aware of cyclic and shake table test methods, but less so hybrid testing. Over the past 40 years, hybrid testing of engineering structures has developed from concept through to maturity to become a reliable and accurate dynamic testing technique. In particular, the application of hybrid testing as a seismic testing technique in recent years has increased notably. The hybrid test method provides users with some additional benefits that standard dynamic testing methods do not, and the method is much more cost effective in comparison to shake table testing. This paper aims to provide the reader with a basic understanding of the hybrid test method and its potential as a dynamic testing technique.

Investigation of tanpura string vibrations using a two-dimensional time-domain model incorporating coupling and bridge friction

Tanpura string vibrations have been investigated previously using numerical models based on energy conserving schemes derived from a Hamiltonian description in one-dimensional form. Such time-domain models have the property that, for the lossless case, the numerical Hamiltonian (representing total energy of the system) can be proven to be constant from one time step
to the next, irrespective of any of the system parameters; in practice the Hamiltonian can be shown to be conserved within machine precision. Models of this kind can reproduce a jvari effect, which results from the bridge-string interaction. However the one-dimensional formulation has recently been shown to fail to replicate the jvaris strong dependence on the thread placement. As a first step towards simulations which accurately emulate this sensitivity to the thread placement, a twodimensional model is proposed, incorporating coupling of controllable level between the two string polarisations at the string termination opposite from the barrier. In addition, a friction force acting when the string slides across the bridge in horizontal direction is introduced, thus effecting a further damping mechanism. In this preliminary study, the string is terminated at the position of the thread. As in the one-dimensional model, an implicit scheme has to be used to solve the system, employing Newton's method to calculate the updated positions and momentums of each string segment. The two-dimensional model is proven to be energy conserving when the loss parameters are set to zero, irrespective of the coupling constant. Both frequency-dependent and independent losses are then added to the string, so that the model can be compared to analogous instruments. The influence of coupling and the bridge friction are investigated.

Wave Fields around Wave Energy Converter Arrays

Wave energy converters, by their nature, extract large amounts of energy
from incident waves. If the industry is to progress such that wave energy
becomes a significant provider of power in the future, large wave farms will
be required. Presently, consenting for these sites is a long and problematic
process, mainly due to a lack of knowledge of the potential environmental
impacts. Accurate numerical modelling of the effect of wave energy extraction
on the wave field and subsequent evaluation of changes to coastal
processes is therefore required. Modelling the wave field impact is also
necessary to allow optimum wave farm configurations to be determined.
This thesis addresses the need for more accurate representation of wave
energy converters in numerical models so that the effect on the wave field,
and subsequently the coastal processes, may be evaluated. Using a hybrid
of physical and numerical modelling (MIKE21 BW and SW models) the
effect of energy extraction and operation of a WEC array on the local wave
climate has been determined.
The main outcomes of the thesis are: an improved wave basin facility, in
terms of wave climate homogeneity, reducing the standard deviation of wave
amplitude by up to 50%; experimental measurement of the wave field around
WEC arrays, showing that radiated waves account for a significant proportion
of the wave disturbance; a new representation method of WECs for use
with standard numerical modelling tools, validated against experimental
results.
The methodology and procedures developed here allow subsequent evaluation
of changes to coastal processes and sediment transport due to WEC
arrays.

Alguns aspectos da modelação ecológica na costa portuguesa (Aveiro) : efeitos físicos na distribuição de nutrientes e biomassa fitoplanctónica

O presente trabalho propõe estudar alguns aspectos da modelação ecológica na Costa Portuguesa, nomeadamente o efeito físico na distribuição de nutrientes e biomassa fitoplanctónica. O principal propósito foi implementar e validar um modelo acoplado tridimensional físico e ecológico, para a costa portuguesa, e aplicá-lo numa área limitada, a norte, por Vila do Conde e a sul pela Figueira da Foz centrada na região de Aveiro, para estudar a distribuição, vertical e horizontal, de temperatura, nutrientes e biomassa fitoplanctónica. A região em estudo está situada na costa oeste da Península Ibérica e faz parte da Região de Afloramento do Atlântico Norte. É caracterizada por condições meteorológicas onde a predominância e prevalência, durante uma grande parte do ano, de ventos de norte/noroeste, constitui um dos principais elementos forçadores do transporte para o largo das águas costeiras e consequente subida das águas mais frias e profundas, ricas em nutrientes. A estas condições juntam-se as boas condições de luminosidade e temperatura essenciais ao desenvolvimento fitoplanctónico, que servirá de alimento às espécies marinhas. Este facto, torna esta região, uma zona de elevada riqueza biológica favorável ao desenvolvimento de várias espécies marinhas, transformando este local num ecossistema de forte produtividade. O modelo foi calibrado e validado para a área em estudo e simulou com sucesso, a resposta do sistema à situação de ventos favoráveis ao afloramento costeiro para a região em estudo. Quando comparada com dados observados, os resultados mostram que o modelo é capaz de prever satisfatoriamente as distribuições superficiais e na coluna de água: da temperatura, dos nutrientes, do oxigénio e da clorofila-a. Os resultados evidenciam o crucial papel desempenhado pelos processos físicos no aumento de fitoplâncton que ocorre ao longo de uma estreita área da costa norte portuguesa, mostrando a estreita ligação entre a distribuição costeira de fitoplâncton e a distribuição costeira de temperatura, à superfície. A produtividade do fitoplâncton que ocorre nas águas costeiras não só é atribuída à disponibilidade de nutrientes mas também à intensidade luminosa. A luz e a intensidade dos ventos de norte/noroeste são os factores chave no controlo dos blooms de fitoplâncton observados nesta região de afloramento, sobretudo no Verão, nomeadamente na camada de mistura pouco profunda e junto à nutriclina. Deste modo os modelos numéricos em associação com dados in situ e imagens de satélite poderão ser considerados uma excelente ferramenta para a análise e previsão de cenários, presentes e futuros, de acções praticadas sobre o meio ambiente.

Sensores ópticos para monitorização dinâmica de estruturas

O presente trabalho centra-se no desenvolvimento e aplicação de sensores de aceleração ópticos, baseados em redes de Bragg gravadas em fibras ópticas, para monitorização da integridade estrutural de estruturas de engenharia civil. Foram implementados dois acelerómetros uniaxiais e um acelerómetro biaxial. Recorreu-se a uma ferramenta de simulação baseada no método dos elementos finitos que permitiu optimizar, sem custos de produção, as características dos sensores, nomeadamente a sua frequência de ressonância. A caracterização dos sensores foi realizada em ambiente laboratorial e a sua resposta comparada com os resultados de simulação, de modo a validar os modelos numéricos. A aplicabilidade e demonstração de conceito foram realizadas na monitorização de estruturas com testes de campo. Foi monitorizado um teste destrutivo de uma parede de adobe, construída à escala real no Departamento de Engenharia Civil da Universidade de Aveiro, onde foram utilizados sensores estáticos e dinâmicos baseados em redes de Bragg gravadas em fibra óptica. Foram realizadas medidas dinâmicas na ponte pedonal do Campus Universitário da Universidade de Aveiro, onde se compararam os resultados obtidos com os sensores ópticos com resultados de sensores electrónicos comerciais. O acelerómetro biaxial foi testado na monitorização de estruturas esbeltas, nomeadamente na monitorização de duas torres de telecomunicações móveis, onde se demonstrou a possibilidade da utilização de sensores ópticos na caracterização dinâmica deste tipo de estruturas. Outro tipo de estruturas de engenharia civil onde foi demonstrada a aplicabilidade dos sensores ópticos desenvolvidos na monitorização estrutural foram os reservatórios de água elevados. Foi realizada a monitorização dinâmica de um exemplo deste tipo de estruturas, localizado no Campus Universitário da Universidade de Aveiro. A monitorização foi realizada recorrendo ao sensor biaxial desenvolvido e a um sismógrafo, ficando o sensor óptico instalado na estrutura de modo a permitir futuras leituras e assim a monitorização periódica da estrutura. Foi ainda desenvolvido um sensor de humidade relativa do ar, com um material sol-gel, que permitiu registar o nível de humidade relativa no interior de blocos de betão, durante um ano. Este sensor pode ser incluído numa rede de sensores multiplexados, na caracterização e monitorização da integridade estrutural de certas estruturas de engenharia civil.