Cabaret finite-difference schemes for the one-dimensional Euler equations

In the present paper we consider second order compact upwind schemes with a space split time derivative (CABARET) applied to one-dimensional compressible gas flows. As opposed to the conventional approach associated with incorporating adjacent space cells we use information from adjacent time layer to improve the solution accuracy. Taking the first order Roe scheme as the basis we develop a few higher (i.e. second within regions of smooth solutions) order accurate difference schemes. One of them (CABARET3) is formulated in a two-time-layer form, which makes it most simple and robust. Supersonic and subsonic shock-tube tests are used to compare the new schemes with several well-known second-order TVD schemes. In particular, it is shown that CABARET3 is notably more accurate than the standard second-order Roe scheme with MUSCL flux splitting.

The post-processed Galerkin method applied to non-linear shell vibrations

We present the results of a computational study of the post-processed Galerkin methods put forward by Garcia-Archilla et al. applied to the non-linear von Karman equations governing the dynamic response of a thin cylindrical panel periodically forced by a transverse point load. We spatially discretize the shell using finite differences to produce a large system of ordinary differential equations (ODEs). By analogy with spectral non-linear Galerkin methods we split this large system into a 'slowly' contracting subsystem and a 'quickly' contracting subsystem. We then compare the accuracy and efficiency of (i) ignoring the dynamics of the 'quick' system (analogous to a traditional spectral Galerkin truncation and sometimes referred to as 'subspace dynamics' in the finite element community when applied to numerical eigenvectors), (ii) slaving the dynamics of the quick system to the slow system during numerical integration (analogous to a non-linear Galerkin method), and (iii) ignoring the influence of the dynamics of the quick system on the evolution of the slow system until we require some output, when we 'lift' the variables from the slow system to the quick using the same slaving rule as in (ii). This corresponds to the post-processing of Garcia-Archilla et al. We find that method (iii) produces essentially the same accuracy as method (ii) but requires only the computational power of method (i) and is thus more efficient than either. In contrast with spectral methods, this type of finite-difference technique can be applied to irregularly shaped domains. We feel that post-processing of this form is a valuable method that can be implemented in computational schemes for a wide variety of partial differential equations (PDEs) of practical importance.

Can frequency diversity provide performance gains for WSNs at 2.4GHz for the fire hydrant to above ground channel

Wireless Sensor Networks (WSNs) which utilise IEEE 802.15.4 technology offer the potential for low cost deployment and maintenance compared with conventional wired sensor networks, enabling effective and efficient condition monitoring of aged civil engineering infrastructure. We will address wireless propagation for a below to above ground scenario where one of the wireless nodes is located in a below ground fire hydrant chamber to permit monitoring of the local water distribution network. Frequency Diversity (FD) is one method that can be used to combat the damaging effects of multipath fading and so improve the reliability of radio links. However, no quantitative investigation concerning the potential performance gains from the use of FD at 2.4GHz is available for the outlined scenario. In this paper, we try to answer this question by performing accurate propagation measurements using modified and calibrated off-the-shelf 802.15.4 based sensor nodes. These measurement results are also compared with those obtained from simulations that employ our Modified 2D Finite-Difference Time-Domain (FDTD) approach. ©2009 IEEE.

Frequency diversity measurements at 2.4 GHz for wireless sensor networks deployed in tunnels

Wireless Sensor Networks (WSNs) which utilise IEEE 802.15.4 technology operate primarily in the 2.4 GHz globally compatible ISM band. However, the wireless propagation channel in this crowded band is notoriously variable and unpredictable, and it has a significant impact on the coverage range and quality of the radio links between the wireless nodes. Therefore, the use of Frequency Diversity (FD) has potential to ameliorate this situation. In this paper, the possible benefits of using FD in a tunnel environment have been quantified by performing accurate propagation measurements using modified and calibrated off-the-shelf 802.15.4 based sensor motes in the disused Aldwych underground railway tunnel. The objective of this investigation is to characterise the performance of FD in this confined environment. Cross correlation coefficients are calculated from samples of the received power on a number of frequency channels gathered during the field measurements. The low measured values of the cross correlation coefficients indicate that applying FD at 2.4 GHz will improve link performance in a WSN deployed in a tunnel. This finding closely matches results obtained by running a computational simulation of the tunnel radio propagation using a 2D Finite-Difference Time-Domain (FDTD) method. ©2009 IEEE.

The effect of an insoluble surfactant on the skin friction of a bubble

In this paper, we develop a novel moving mesh method suitable for solving axisymmetric free-boundary problems, including the Marangoni effect induced by surfactant or temperature variation. This method employs a body-fitted grid system where the gas-liquid interface is one line of the grid system. We model the surfactant equation of state with a non-linear Langmuir law, and, for simplicity, we limit ourselves to the situation of an insoluble surfactant. We solve complicated dynamic boundary conditions accurately on the gas-liquid interface in the framework of finite-volume methods. Our method is used to study the effect of a surfactant on the skin friction of a bubble in a uniaxial flow. For the limiting case where the surface diffusivity is zero, the effect of a tangential stress generated by the surface tension gradient, allows us to explain a new phenomenon in high concentration regimes: larger surface tension, but also larger deformation. Furthermore, this condition leads to the formation of boundary layers and flow separation at high Reynolds numbers. The influence of these complex flow patterns is examined. © 2005 Elsevier SAS. All rights reserved.

Coupled deformation-flow analysis for methane hydrate production by depressurized wells

A formulation for coupled flow-deformation analysis of methane-hydrate extraction problems is presented. By assuming that the hydrate does not flow, a two phase flow formulation is considered, based on Darcy's law and capillary pressure relation. The formulation is implemented in the finite difference code FLAC. The code was used to investigate the stability of a methane extraction well by depressurizing the well. © 2005 Taylor & Francis Group, London.

Transient thermal analysis of power devices based on Fourier-series thermal model

A new thermal model based on Fourier series expansion method has been presented for dynamic thermal analysis on power devices. The thermal model based on the Fourier series method has been programmed in MATLAB SIMULINK and integrated with a physics-based electrical model previously reported. The model was verified for accuracy using a two-dimensional Fourier model and a two-dimensional finite difference model for comparison. To validate this thermal model, experiments using a 600V 50A IGBT module switching an inductive load, has been completed under high frequency operation. The result of the thermal measurement shows an excellent match with the simulated temperature variations and temperature time-response within the power module. ©2008 IEEE.

Building flexible user interfaces for solving PDEs

FEniCS is a collection of software tools for the automated solution of differential equations by finite element methods. In this note, we describe how FEniCS can be used to solve a simple nonlinear model problem with varying levels of automation. At one extreme, FEniCS provides tools for the fully automated and adaptive solution of nonlinear partial differential equations. At the other extreme, FEniCS provides a range of tools that allow the computational scientist to experiment with novel solution algorithms. © 2010 American Institute of Physics.

A FEniCS-based programming framework for modeling turbulent flow by the Reynolds-averaged Navier-Stokes equations

Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model. © 2011 Elsevier Ltd.