Convergence and tracking analysis of a variable normalised LMF (XE-NLMF) algorithm

The least-mean-fourth (LMF) algorithm is known for its fast convergence and lower steady state error, especially in sub-Gaussian noise environments. Recent work on normalised versions of the LMF algorithm has further enhanced its stability and performance in both Gaussian and sub-Gaussian noise environments. For example, the recently developed normalised LMF (XE-NLMF) algorithm is normalised by the mixed signal and error powers, and weighted by a fixed mixed-power parameter. Unfortunately, this algorithm depends on the selection of this mixing parameter. In this work, a time-varying mixed-power parameter technique is introduced to overcome this dependency. A convergence analysis, transient analysis, and steady-state behaviour of the proposed algorithm are derived and verified through simulations. An enhancement in performance is obtained through the use of this technique in two different scenarios. Moreover, the tracking analysis of the proposed algorithm is carried out in the presence of two sources of nonstationarities: (1) carrier frequency offset between transmitter and receiver and (2) random variations in the environment. Close agreement between analysis and simulation results is obtained. The results show that, unlike in the stationary case, the steady-state excess mean-square error is not a monotonically increasing function of the step size. (c) 2007 Elsevier B.V. All rights reserved.

Inelastic quantum transport in nanostructures: The self-consistent Born approximation and correlated electron-ion dynamics

A dynamical method for inelastic transport simulations in nanostructures is compared to a steady-state method based on nonequilibrium Green's functions. A simplified form of the dynamical method produces, in the steady state in the weak-coupling limit, effective self-energies analogous to those in the Born approximation due to electron-phonon coupling. The two methods are then compared numerically on a resonant system consisting of a linear trimer weakly embedded between metal electrodes. This system exhibits an enhanced heating at high biases and long phonon equilibration times. Despite the differences in their formulation, the static and dynamical methods capture local current-induced heating and inelastic corrections to the current with good agreement over a wide range of conditions, except in the limit of very high vibrational excitations where differences begin to emerge.

Dynamics of Gompertzian tumour growth under environmental fluctuations

We investigate the effect of correlated additive and multiplicative Gaussian white noise oil the Gompertzian growth of tumours. Our results are obtained by Solving numerically the time-dependent Fokker-Planck equation (FPE) associated with the stochastic dynamics. In Our numerical approach we have adopted B-spline functions as a truncated basis to expand the approximated eigenfunctions. The eigenfunctions and eigenvalues obtained using this method are used to derive approximate solutions of the dynamics under Study. We perform simulations to analyze various aspects, of the probability distribution. of the tumour cell populations in the transient- and steady-state regimes. More precisely, we are concerned mainly with the behaviour of the relaxation time (tau) to the steady-state distribution as a function of (i) of the correlation strength (lambda) between the additive noise and the multiplicative noise and (ii) as a function of the multiplicative noise intensity (D) and additive noise intensity (alpha). It is observed that both the correlation strength and the intensities of additive and multiplicative noise, affect the relaxation time.

Testing Tidal Turbines - Part 1: Steady Towing Tests vs Tidal Moored Tests

Tidal turbines have been tested extensively at many scales in steady state flow. Testing medium- or full-scale devices in turbulent flow has been less thoroughly examined. The differences between turbine performances in these two different states are needed for testing method verification and numerical model validation. The work in this paper documents the performance of a 1/10 scale turbine in steady state pushing tests and tidal moored tests. The overall performance of the device appears to decrease with turbulent flow, though there is increased data scatter and therefore, reduced uncertainty. At maximum power performance, as velocity increases the mechanical power and electrical power reduction from steady to unsteady flow increases. The drive train conversion efficiency also decreases. This infers that the performance for this turbine design is affected by the presence of turbulent flow.

Nonconservative dynamics in long atomic wires

The effect of nonconservative current-induced forces on the ions in a defect-free metallic nanowire is investigated using both steady-state calculations and dynamical simulations. Nonconservative forces were found to have a major influence on the ion dynamics in these systems, but their role in increasing the kinetic energy of the ions decreases with increasing system length. The results illustrate the importance of nonconservative effects in short nanowires and the scaling of these effects with system size. The dependence on bias and ion mass can be understood with the help of a simple pen and paper model. This material highlights the benefit of simple preliminary steady-state calculations in anticipating aspects of brute-force dynamical simulations, and provides rule of thumb criteria for the design of stable quantum wires.

ELM-induced transient tungsten melting in the JET divertor

The original goals of the JET ITER-like wall included the study of the impact of an all W divertor on plasma operation (Coenen et al 2013 Nucl. Fusion 53 073043) and fuel retention (Brezinsek et al 2013 Nucl. Fusion 53 083023). ITER has recently decided to install a full-tungsten (W) divertor from the start of operations. One of the key inputs required in support of this decision was the study of the possibility of W melting and melt splashing during transients. Damage of this type can lead to modifications of surface topology which could lead to higher disruption frequency or compromise subsequent plasma operation. Although every effort will be made to avoid leading edges, ITER plasma stored energies are sufficient that transients can drive shallow melting on the top surfaces of components. JET is able to produce ELMs large enough to allow access to transient melting in a regime of relevance to ITER.

Transient W melt experiments were performed in JET using a dedicated divertor module and a sequence of I-P = 3.0 MA/B-T = 2.9 T H-mode pulses with an input power of P-IN = 23 MW, a stored energy of similar to 6 MJ and regular type I ELMs at Delta W-ELM = 0.3 MJ and f(ELM) similar to 30 Hz. By moving the outer strike point onto a dedicated leading edge in the W divertor the base temperature was raised within similar to 1 s to a level allowing transient, ELM-driven melting during the subsequent 0.5 s. Such ELMs (delta W similar to 300 kJ per ELM) are comparable to mitigated ELMs expected in ITER (Pitts et al 2011 J. Nucl. Mater. 415 (Suppl.) S957-64).

Although significant material losses in terms of ejections into the plasma were not observed, there is indirect evidence that some small droplets (similar to 80 mu m) were released. Almost 1 mm (similar to 6 mm(3)) of W was moved by similar to 150 ELMs within 7 subsequent discharges. The impact on the main plasma parameters was minor and no disruptions occurred. The W-melt gradually moved along the leading edge towards the high-field side, driven by j x B forces. The evaporation rate determined from spectroscopy is 100 times less than expected from steady state melting and is thus consistent only with transient melting during the individual ELMs. Analysis of IR data and spectroscopy together with modelling using the MEMOS code Bazylev et al 2009 J. Nucl. Mater. 390-391 810-13 point to transient melting as the main process. 3D MEMOS simulations on the consequences of multiple ELMs on damage of tungsten castellated armour have been performed.

These experiments provide the first experimental evidence for the absence of significant melt splashing at transient events resembling mitigated ELMs on ITER and establish a key experimental benchmark for the MEMOS code.

Unsteady tip leakage characteristics and heat transfer on turbine blade tip and casing

An unsteady numerical investigation was performed to examine time dependent behaviors of the tip leakage flow structures and heat transfer on the rotor blade tip and casing in a single stage gas turbine engine. A transonic, high-pressure
turbine stage was modeled and simulated using a stage pressure ratio of 3.2. The rotor’s tip clearance was 1.2 mm in height (3% of the rotor span) and its speed was set at 9500 rpm. Periodic flow is observed for each vane passing period. Tip leakage flow as well as heat transfer data showed highly time dependent behaviors. A stator trailing edge shock appears as the turbine stage is operating at transonic conditions. The shock alters the flow condition in the rotor section, namely, the tip leakage flow structures and heat transfer rate distributions. The instantaneous Nusselt number distributions are compared to the time averaged and steady-state results. The same patterns in tip leakage flow
structures and heat transfer rate distributions were observed in both unsteady and steady simulations. However, the unsteady simulation captured the locally time-dependent high heat transfer phenomena caused by the unsteady interaction with the upstream vane trailing-edge shock and the passing wake.

Field test of multiple 1/10 scale tidal turbines in steady flows

Queen's University Belfast and Wave Barrier Ltd have developed a tidal testing platform to test hydrokinetic turbines at medium scale. Multiple turbines can be pushed through still water conditions, in steady-state pushing tests. Experiments were conducted to evaluate the interactions between two identical, mono-strut, horizontal axis tidal turbines (HATTs) of 1.5 m diameter (D) rotor. Their relative performance when located individually, in-plane and in-line are investigated. The data shows a high consistency in the power curves at different flow speeds, which indicates high repeatability in this Reynolds range. For an individual turbine, there is no performance difference when the rotor is mounted either upstream or downstream of the supporting structure. When placed in-plane, the turbines have no adverse effect on one another. When spaced in-line with 2D separation, there is a 63% reduction in the performance of the downstream turbine. At 6D downstream this performance reduction is still 59%, indicating some wake recovery between 2D and 6D, though the influence from the upstream rotor persists to at least 6D downstream of the first device. In contrast the performance of the downstream turbine when placed at 1.5D offset of the upstream device at 6D downstream is approximately recovered to the individual turbine performance.

Modelling of Organic Rankine Cycle for Waste Heat Recovery Process in Supercritical Condition

Organic Rankine Cycle (ORC) is the most commonly used method for recovering energy from small sources of heat. The investigation of the ORC in supercritical condition is a new research area as it has a potential to generate high power and thermal efficiency in a waste heat recovery system. This paper presents a steady state ORC model in supercritical condition and its simulations with a real engine’s exhaust data. The key component of ORC, evaporator, is modelled using finite volume method, modelling of all other components of the waste heat recovery system such as pump, expander and condenser are also presented. The aim of this paper is to investigate the effects of mass flow rate and evaporator outlet temperature on the efficiency of the waste heat recovery process. Additionally, the necessity of maintaining an optimum evaporator outlet temperature is also investigated. Simulation results show that modification of mass flow rate is the key to changing the operating temperature at the evaporator outlet.

Magnetospheric Structure and Atmospheric Joule Heating of Habitable Planets Orbiting M-dwarf Stars

We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvénic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvénic sectors, while no bow shock forms in the sub-Alfvénic sectors. The planets reside most of the time in the sub-Alfvénic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the intense stellar wind. For the steady-state solution, the heating is about 0.1%-3% of the total incoming stellar irradiation, and it is enhanced by 50% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport. Here we assume constant ionospheric Pedersen conductance similar to that of the Earth. The conductance could be greater due to the intense EUV radiation leading to smaller heating rates. We plan to quantify the ionospheric conductance in future study.

Length matters:Keeping atomic wires in check

Dynamical effects of non-conservative forces in long, defect free atomic wires are investigated. Current flow through these wires is simulated and we find that during the initial transient, the kinetic energies of the ions are contained in a small number of phonon modes, closely clustered in frequency. These phonon modes correspond to the waterwheel modes determined from preliminary static calculations. The static calculations allow one to predict the appearance of non-conservative effects in advance of the more expensive real-time simulations. The ion kinetic energy redistributes across the band as non-conservative forces reach a steady state with electronic factional forces. The typical ion kinetic energy is found to decrease with system length, increase with atomic mass, and its dependence on bias, mass and length is supported with a pen and paper model. This paper highlights the importance of non-conservative forces in current carrying devices and provides criteria for the design of stable atomic wires.

Nonequilibrium processes from generalized Langevin equations: Realistic nanoscale systems connected to two thermal baths

We extend the generalized Langevin equation (GLE) method [L. Stella, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nosé-Hoover thermostats). We concentrate on the steady-state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e., ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient.

Assessment of turbulence model predictions for a centrifugal compressor simulation

Steady-state computational fluid dynamics (CFD) simulations are an essential tool in the design process of centrifugal compressors. Whilst global parameters, such as pressure ratio and efficiency, can be predicted with reasonable accuracy, the accurate prediction of detailed compressor flow fields is a much more significant challenge. Much of the inaccuracy is associated with the incorrect selection of turbulence model. The need for a quick turnaround in simulations during the design optimisation process, also demands that the turbulence model selected be robust and numerically stable with short simulation times.
In order to assess the accuracy of a number of turbulence model predictions, the current study used an exemplar open CFD test case, the centrifugal compressor ‘Radiver’, to compare the results of three eddy viscosity models and two Reynolds stress type models. The turbulence models investigated in this study were (i) Spalart-Allmaras (SA) model, (ii) the Shear Stress Transport (SST) model, (iii) a modification to the SST model denoted the SST-curvature correction (SST-CC), (iv) Reynolds stress model of Speziale, Sarkar and Gatski (RSM-SSG), and (v) the turbulence frequency formulated Reynolds stress model (RSM-ω). Each was found to be in good agreement with the experiments (below 2% discrepancy), with respect to total-to-total parameters at three different operating conditions. However, for the off-design conditions, local flow field differences were observed between the models, with the SA model showing particularly poor prediction of local flow structures. The SST-CC showed better prediction of curved rotating flows in the impeller. The RSM-ω was better for the wake and separated flow in the diffuser. The SST model showed reasonably stable, robust and time efficient capability to predict global and local flow features.