Numerical simulation of tides and tidal currents in Liaodong Bay with POM

A three-dimensional ocean circulation model, called Princeton Ocean Model (POM), is employed to simulate tides and tidal currents in Liaodong Bay. The nested grid technique is adopted to improve the computation precision. Computed harmonic constants of M-1, M-2 tides at five tidal gauge stations and surface elevations at two oil platforms are compared with those observed, and show good agreements with them. Based on the calculated results, the co-amplitude and co-phase tag lines of nil and M-2 tidal constituents, the residual current field of M-2 constituent, tidal form, tidal Current ellipse and the moving style of tidal current are given. It is found that diurnal tidal constituents have no amphidromic point whereas semi-diurnal constituents have one in the region of interest. Meanwhile, some meaningful results are concluded and presented, which are conducive to a thorough knowledge of the characteristics of tides and tidal currents in the Liaodong Bay.

Simulation of barotropic and baroclinic tides in the South China Sea

The four leading tidal constituents M-2, S-2, K-1 and O-1 in the South China Sea are simulated by using POM. The model is forced with tide-generating potential and four leading tidal constituents at the open boundary. In order to simulate more exactly, TOPEX/Poseidon altimeter data are assimilated into the model and the open boundary is optimized. The computed co-tidal charts for M-2 and K-1 constituents are generally consistent with previous results in this region. The numerical simulation shows that energetic internal tides are generated over the bottom topography such as the Dongsha Islands, the Xisha Islands, the Zhongsha Islands, the Nansba Islands and the Luzon Strait.

Numerical simulation of the mass loss process in pyrolizing char materials

We present here a decoupling technique to tackle the entanglement of the nonlinear boundary condition and the movement of the char/virgin front for a thermal pyrolysis model for charring materials. Standard numerical techniques to solve moving front problems — often referred to as Stefan problems — encounter difficulties when dealing with nonlinear boundaries. While special integral methods have been developed to solve this problem, they suffer from several limitations which the technique described here overcomes. The newly developed technique is compared with the exact analytical solutions for some simple ideal situations which demonstrate that the numerical method is capable of producing accurate numerical solutions. The pyrolysis model is also used to simulate the mass loss process from a white pine sample exposed to a constant radiative flux in a nitrogen atmosphere. Comparison with experimental results demonstrates that the predictions of mass loss rates and temperature profile within the solid material are in good agreement with the experiment.

Flight simulation

Dr Fuchen Jia, Dr Mayer Patel and Professor Edwin Galea explain how advanced fire models were used to unravel the secrets of Swissair Flight 111, which crashed off the coast of Canada in 1998.

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.

Numerical simulation of microwave thawing using a coupled solver approach

Thawing of a frozen food product in a domestic microwave oven is numerically simulated using a coupled solver approach. The approach consists of a dedicated electromagnetic FDTD solver and a closely coupled UFVM multi-physics package. Two overlapping numerical meshes are defined; the food material and container were meshed for heat transfer and phase change solution, whilst the microwave oven cavity and waveguide were meshed for the microwave irradiation. The two solution domains were linked using a cross-mapping routine. This approach allowed the rotation of the food load to be captured. Power densities obtained on the structured FDTD mesh were interpolated onto the UFVM mesh for each timestep/turntable position. The UFVM solver utilised the power density data to advance the temperature and phase distribution solution. The temperature-dependant dielectric and thermo-physical properties of the food load were updated prior to revising the electromagnetic solution. Changes in thermal/electric properties associated with the phase transition were fully accounted for as well as heat losses from product to cavity. Two scenarios were investigated: a centric and eccentric placement on the turntable. Developing temperature fields predicted by the numerical solution are validated against experimentally obtained data. Presented results indicate the feasibility of fully coupled simulations of the microwave heating of a frozen product. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Low SNR Capacity for MIMO Rician and Rayleigh-Product Fading Channels with Single Co-channel Interferer and Noise

This paper studies the ergodic capacity of multiple-input multiple-output (MIMO) systems with a single co-channel interferer in the low signal-to-noise-ratio (SNR) regime. Two MIMO models namely Rician and Rayleigh-product channels are investigated. Exact analytical expressions for the minimum energy per information bit, Eb/N0min, and wideband slope, S0, are derived for both channels. Our results show that the minimum energy per information bit is the same for both channels while their wideband slopes differ significantly. Further, the impact of the numbers of transmit and receive antennas, the Rician K factor, the channel mean matrix and the interference-to-noise-ratio (INR) on the capacity, is addressed. Results indicate that interference degrades the capacity by increasing the required minimum energy per information bit and reducing the wideband slope. Simulation results validate our analytical results.

Conventional 3D staging PET/CT in CT simulation for lung cancer: impact of rigid and deformable target volume alignments for radiotherapy treatment planning

Objective: Positron emission tomography (PET)/CT scans can improve target definition in radiotherapy for non-small cell lung cancer (NSCLC). As staging PET/CT scans are increasingly available, we evaluated different methods for co-registration of staging PET/CT data to radiotherapy simulation (RTP) scans.

Methods: 10 patients underwent staging PET/CT followed by RTP PET/CT. On both scans, gross tumour volumes (GTVs) were delineated using CT (GTVCT) and PET display settings. Four PET-based contours (manual delineation, two threshold methods and a source-to-background ratio method) were delineated. The CT component of the staging scan was co-registered using both rigid and deformable techniques to the CT component of RTP PET/CT. Subsequently rigid registration and deformation warps were used to transfer PET and CT contours from the staging scan to the RTP scan. Dice’s similarity coefficient (DSC) was used to assess the registration accuracy of staging-based GTVs following both registration methods with the GTVs delineated on the RTP PET/CT scan.

Results: When the GTVCT delineated on the staging scan after both rigid registration and deformation was compared with the GTVCT on the RTP scan, a significant improvement in overlap (registration) using deformation was observed (mean DSC 0.66 for rigid registration and 0.82 for deformable registration, p50.008). A similar comparison for PET contours revealed no significant improvement in overlap with the use of deformable registration.

Conclusions: No consistent improvements in similarity measures were observed when deformable registration was used for transferring PET-based contours from a staging PET/CT. This suggests that currently the use of rigid registration remains the most appropriate method for RTP in NSCLC.

High pressure CO2 absorption studies on imidazolium-based ionic liquids:Experimental and simulation approaches

A combined experimental-computational study on the CO absorption on 1-butyl-3-methylimidazolium hexafluophosphate, 1-ethyl-3-methylimidazolium bis[trifluoromethylsulfonyl]imide, and 1-butyl-3-methylimidazolium bis[trifluoromethylsulfonyl]imide ionic liquids is reported. The reported results allowed to infer a detailed nanoscopic vision of the absorption phenomena as a function of pressure and temperature. Absorption isotherms were measured at 318 and 338K for pressures up to 20MPa for ultrapure samples using a state-of-the-art magnetic suspension densimeter, for which measurement procedures are developed. A remarkable swelling effect upon CO absorption was observed for pressures higher than 10MPa, which was corrected using a method based on experimental volumetric data. The experimental data reported in this work are in good agreement with available literature isotherms. Soave-Redlich-Kwong and Peng-Robinson equations of state coupled with bi-parametric van der Waals mixing rule were used for successful correlations of experimental high pressure absorption data. Molecular dynamics results allowed to infer structural, energetic and dynamic properties of the studied CO+ionic liquids mixed fluids, showing the relevant role of the strength of anion-cation interactions on fluid volumetric properties and CO absorption. © 2012 Elsevier B.V.

A co-operative formation-based collision avoidance approach for a group of autonomous vehicles (Special Issue on Adaptive Control and Signal Processing in Marine System)

This paper employs a unique decentralised cooperative control method to realise a formation-based collision avoidance strategy for a group of autonomous vehicles. In this approach, the vehicles' role in the formation and their alert and danger areas are first defined, and the formation-based intra-group and external collision avoidance methods are then proposed to translate the collision avoidance problem into the formation stability problem. The extension–decomposition–aggregation formation control method is next employed to stabilise the original and modified formations, whilst manoeuvring, and subsequently solve their collision avoidance problem indirectly. Simulation study verifies the feasibility and effectiveness of the intra-group and external collision avoidance strategy. It is demonstrated that both formation control and collision avoidance problems can be simultaneously solved if the stability of the expanded formation including external obstacles can be satisfied.

Co-operative tensegrity-based formation control algorithm for a multi-aircraft system

This paper presents a tensegrity-based co-operative control algorithm for an aircraft formation. The 6 degrees-of-freedom model of the well-known Aerosonde unmanned aerial vehicle (UAV), is integrated with the model of the tensegrity structure and a decentralised control scheme is proposed. The strategy is shown to be scalable for 2n number of UAVs and is able to maintain a firm geometry whilst allowing flexible shape transformations. Simulation results demonstrate the effectiveness and stability of the proposed tensegrity-based formation control algorithm in 3D.

Multi-agent electricity market simulation with dynamic strategies & virtual power producers

Distributed energy resources will provide a significant amount of the electricity generation and will be a normal profitable business. In the new decentralized grid, customers will be among the many decentralized players and may even help to co-produce the required energy services such as demand-side management and load shedding. So, they will gain the opportunity to be more active market players. The aggregation of DG plants gives place to a new concept: the Virtual Power Producer (VPP). VPPs can reinforce the importance of these generation technologies making them valuable in electricity markets. In this paper we propose the improvement of MASCEM, a multi-agent simulation tool to study negotiations in electricity spot markets based on different market mechanisms and behavior strategies, in order to take account of decentralized players such as VPP.

Microsatellites can be misleading: an empirical and simulation study.

It has been long recognized that highly polymorphic genetic markers can lead to underestimation of divergence between populations when migration is low. Microsatellite loci, which are characterized by extremely high mutation rates, are particularly likely to be affected. Here, we report genetic differentiation estimates in a contact zone between two chromosome races of the common shrew (Sorex araneus), based on 10 autosomal microsatellites, a newly developed Y-chromosome microsatellite, and mitochondrial DNA. These results are compared to previous data on proteins and karyotypes. Estimates of genetic differentiation based on F- and R-statistics are much lower for autosomal microsatellites than for all other genetic markers. We show by simulations that this discrepancy stems mainly from the high mutation rate of microsatellite markers for F-statistics and from deviations from a single-step mutation model for R-statistics. The sex-linked genetic markers show that all gene exchange between races is mediated by females. The absence of male-mediated gene flow most likely results from male hybrid sterility.

Genetic programming for the RoboCup Rescue Simulation System

The Robocup Rescue Simulation System (RCRSS) is a dynamic system of multi-agent interaction, simulating a large-scale urban disaster scenario. Teams of rescue agents are charged with the tasks of minimizing civilian casualties and infrastructure damage while competing against limitations on time, communication, and awareness. This thesis provides the first known attempt of applying Genetic Programming (GP) to the development of behaviours necessary to perform well in the RCRSS. Specifically, this thesis studies the suitability of GP to evolve the operational behaviours required of each type of rescue agent in the RCRSS. The system developed is evaluated in terms of the consistency with which expected solutions are the target of convergence as well as by comparison to previous competition results. The results indicate that GP is capable of converging to some forms of expected behaviour, but that additional evolution in strategizing behaviours must be performed in order to become competitive. An enhancement to the standard GP algorithm is proposed which is shown to simplify the initial search space allowing evolution to occur much quicker. In addition, two forms of population are employed and compared in terms of their apparent effects on the evolution of control structures for intelligent rescue agents. The first is a single population in which each individual is comprised of three distinct trees for the respective control of three types of agents, the second is a set of three co-evolving subpopulations one for each type of agent. Multiple populations of cooperating individuals appear to achieve higher proficiencies in training, but testing on unseen instances raises the issue of overfitting.