Numerical simulations of fluid-structure interactions in single-reed mouthpieces

Most single-reed woodwind instrument models rely on a quasistationary approximation to describe the relationship between the volume flow and. the pressure difference across the reed channel. Semiempirical models based on the quasistationary approximation are very useful in explaining the fundamental characteristics of this family of instruments such as self-sustained oscillations and threshold of blowing pressure. However, they fail at explaining more complex phenomena associated with the fluid-structure interaction during dynamic flow regimes, such as the transient and steady-state behavior of the system as a function. of the mouthpiece geometry. Previous studies have discussed the accuracy of the quasistationary approximation but the amount of literature on the subject is sparse, mainly due to the difficulties involved in the measurement of dynamic flows in channels with an oscillating reed. In this paper, a numerical technique based on the lattice Boltzmann method and a finite difference scheme is proposed in order to investigate the characteristics of fully coupled fluid-structure interaction in single-reed mouthpieces with different channel configurations. Results obtained for a stationary simulation with a static reed agree very well with those predicted by the literature based on the quasistationary approximation. However, simulations carried out for a dynamic regime with dn oscillating reed show that the phenomenon associated with flow detachment and reattachment diverges considerably frorn the theoretical assumptions. Furthermore, in the case of long reed channels, the results obtained for the vena contracta factor are in significant disagreement with those predicted by theory. For short channels, the assumption of constant vena contracta was found to be valid for only 40% of the duty cycle. (c) 2007 Acoustical Society of America.

NUMERICAL SIMULATIONS OF BURST PROCESSES IN FIBER-BUNDLES

We discuss a very effective numerical method for simulating fibre-bundle models with equal load-sharing and with local load-sharing. Particular attention is paid to the case of the local load-sharing model, in which the critical load x(c) is defined as the average load per fibre that causes the final complete failure. It is shown that x(c) --> 0 when the size of the system N --> infinity. We also show analytically that the power law of the burst size distribution, D(Delta) alpha Delta(-xi), is approximately correct. The exponent xi in the local load-sharing case is not universal, since it depends on the strength distribution as well on as the size of the system.

A Review of Wave Makers for 3D Numerical Simulations

Numerical methods have enabled the simulation of complex problems in off-shore and marine engineering. A significant challenge in these simulations is the creation of a realistic wave field. A good numerical tank requires wave creation and absorption of waves at various locations. Several numerical wavemakers with these capabilities have been presented in the past. This paper reviews four different wave-maker methods and discusses limitations, computational efficiency, requirements on the mesh and preprocessing and complexity of implementation.

Optimising power take-off of an oscillating wave surge converter using high fidelity numerical simulations

Oscillating wave surge converters are a promising technology to harvest ocean wave energy in the near shore region. Although research has been going on for many years, the characteristics of the wave action on the structure and especially the phase relation between the driving force and wave quantities like velocity or surface elevation have not been investigated in detail. The main reason for this is the lack of suitable methods. Experimental investigations using tank tests do not give direct access to overall hydrodynamic loads, only damping torque of a power take off system can be measured directly. Non-linear computational fluid dynamics methods have only recently been applied in the research of this type of devices. This paper presents a new metric named wave torque, which is the total hydrodynamic torque minus the still water pitch stiffness at any given angle of rotation. Changes in characteristics of that metric over a wave cycle and for different power take off settings are investigated using computational fluid dynamics methods. Firstly, it is shown that linearised methods cannot predict optimum damping in typical operating states of OWSCs. We then present phase relationships between main kinetic parameters for different damping levels. Although the flap seems to operate close to resonance, as predicted by linear theory, no obvious condition defining optimum damping is found.

Propagating Wave Phenomena Detected in Observations and Simulations of the Lower Solar Atmosphere

We present high-cadence observations and simulations of the solar photosphere, obtained using the Rapid Oscillations in the Solar Atmosphere imaging system and the MuRAM magnetohydrodynamic (MHD) code, respectively. Each data set demonstrates a wealth of magnetoacoustic oscillatory behavior, visible as periodic intensity fluctuations with periods in the range 110–600 s. Almost no propagating waves with periods less than 140 s and 110 s are detected in the observational and simulated data sets, respectively. High concentrations of power are found in highly magnetized regions, such as magnetic bright points and intergranular lanes. Radiative diagnostics of the photospheric simulations replicate our observational results, confirming that the current breed of MHD simulations are able to accurately represent the lower solar atmosphere. All observed oscillations are generated as a result of naturally occurring magnetoconvective processes, with no specific input driver present. Using contribution functions extracted from our numerical simulations, we estimate minimum G-band and 4170 Å continuum formation heights of 100 km and 25 km, respectively. Detected magnetoacoustic oscillations exhibit a dominant phase delay of −8◦ between the G-band and 4170 Å continuum observations, suggesting the presence of upwardly propagating waves.More than 73% of MBPs (73% from observations and 96% from simulations) display upwardly propagating wave phenomena, suggesting the abundant nature of oscillatory behavior detected higher in the solar atmosphere may be traced back to magnetoconvective processes occurring in the upper layers of the Sun’s convection zone.

Numerical and Experimental Investigation of Aircraft Panel Deformations During Riveting Process

In collaboration with Airbus-UK, the dimensional growth of aircraft panels while being riveted with stiffeners is investigated. Small panels are used in this investigation. The stiffeners have been fastened to the panels with rivets and it has been observed that during this operation the panels expand in the longitudinal and transverse directions. It has been observed that the growth is variable and the challenge is to control the riveting process to minimize this variability. In this investigation, the assembly of the small panels and longitudinal stiffeners has been simulated using static stress and nonlinear explicit finite element models. The models have been validated against a limited set of experimental measurements; it was found that more accurate predictions of the riveting process are achieved using explicit finite element models. Yet, the static stress finite element model is more time efficient, and more practical to simulate hundreds of rivets and the stochastic nature of the process. Furthermore, through a series of numerical simulations and probabilistic analyses, the manufacturing process control parameters that influence panel growth have been identified. Alternative fastening approaches were examined and it was found that dimensional growth can be controlled by changing the design of the dies used for forming the rivets.

Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. It has been proposed that some of these outflows are driven off gas infalling towards the central supermassive black hole. We perform radiative transfer calculations to compute the gas ionization state and the emergent X-ray spectra for both two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of relatively dense (10 -10 g cm ) outflowing gas. This outflow causes strong absorption for observer orientation angles of ?35°. Particularly in 3D, the strength of this absorption varies significantly for different lines of sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features in the X-ray spectra since the ionization state is predicted to be very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a Compton-scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow are predicted to give rise to several emission features including a multicomponent Fe Ka emission complex for all observer orientations. For the hydrodynamical simulations considered, we predict both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. Higher resolution hydrodynamic simulations are needed to determine whether the outflows fragment on scales unresolved in our current study, which may yield the denser lower ionization material that could reconcile the models and the observations. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Vlasov-kinetic computer simulations of electrostatic waves in dusty plasmas: an overview of recent results

A series of numerical simulations based on a recurrence-free Vlasov kinetic model using kinetic phase point trajectories are presented. Electron-ion plasmas and three-component (electron-ion-dust) dusty or complex plasmas are considered, via independent simulations. Considering all plasma components modeled through a kinetic approach, the linear and nonlinear behavior of ion-acoustic excitations is investigated. Maxwellian and kappa-type (superthermal) distribution functions are assumed, as initial conditions, in separate simulations for the sake of comparison. The focus is on the parametric dependence of ion-acoustic waves on the electron-to-ion temperature ratio and on the dust concentration. © 2014 EDP Sciences, SIF, Springer-Verlag Berlin Heidelberg.

High resolution numerical simulation of triple point collision and origin of unburned gas pockets in turbulent detonations

A key issue in pulse detonation engine development is better understanding of the detonation structure and its propagation mechanism. Thus, in the present work the turbulent structure of an irregular detonation is studied through very high resolution numerical simulations of 600 points per half reaction length. The aim is to explore the nature of the transverse waves during the collision and reflection processes of the triple point with the channel walls. Consequently the formation and consumption mechanism of unreacted gas pockets is studied. Results show that the triple point and the transverse wave collide simultaneously with the wall. The strong transverse wave switches from a primary triple point before collision to a new one after reflection. Due to simultaneous interaction of the triple point and the transverse wave with the wall in the second half of the detonation cell, a larger high-pressurised region appears on the wall. During the reflection the reaction zone detaches from the shock front and produces a pocket of unburned gas. Three mechanisms found to be of significance in the re-initiation mechanism of detonation at the end of the detonation cell; i: energy resealed via consumption of unburned pockets by turbulent mixing ii: compression waves arise due to collision of the triple point on the wall which helps the shock to jump abruptly to an overdriven detonation iii: drastic growth of the Richtmyer–Meshkov instability causing a part of the front to accelerate with respect to the neighbouring portions.

Equalisation with adaptive time lag

In an adaptive equaliser, the time lag is an important parameter that significantly influences the performance. Only with the optimum time lag that corresponds to the best minimum-mean-square-error (MMSE) performance, can there be best use of the available resources. Many designs, however, choose the time lag either based on preassumption of the channel or simply based on average experience. The relation between the MMSE performance and the time lag is investigated using a new interpretation of the MMSE equaliser, and then a novel adaptive time lag algorithm is proposed based on gradient search. The proposed algorithm can converge to the optimum time lag in the mean and is verified by the numerical simulations provided.

Laser-driven proton scaling laws and new paths towards energy increase

The past few years have seen remarkable progress in the development of laser-based particle accelerators. The ability to produce ultrabright beams of multi-megaelectronvolt protons routinely has many potential uses from engineering to medicine, but for this potential to be realized substantial improvements in the performances of these devices must be made. Here we show that in the laser-driven accelerator that has been demonstrated experimentally to produce the highest energy protons, scaling laws derived from fluid models and supported by numerical simulations can be used to accurately describe the acceleration of proton beams for a large range of laser and target parameters. This enables us to evaluate the laser parameters needed to produce high-energy and high-quality proton beams of interest for radiography of dense objects or proton therapy of deep-seated tumours.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Plasma boundary sheath in the afterglow of a pulsed inductively coupled RF plasma

The sheath dynamics in the afterglow of a pulsed inductively coupled plasma, operated in hydrogen, is investigated. It is found that the sheath potential does not fully collapse in the early post-discharge. Time resolved measurements of the positive ion flux in a hydrogen plasma, using a mass resolved ion energy analyser, reveal that a constant 2 eV mean ion energy persists for several hundred micro-seconds in the afterglow. The presence of a finite sheath potential is explained by super-elastic collisions between vibrationally excited hydrogen molecules and electrons in the afterglow, leading to an electron temperature of about 0.5 eV. Plasma density decay times measured using both the mass resolved energy analyser and a Langmuir probe are in good agreement. Vibrational temperatures measured using optical emission spectroscopy support the theory of electron heating through super-elastic collisions with vibrationally excited hydrogen molecules. Measurements are also supported by numerical simulations and modelling results.

Idealised operation of zero-voltage-switching series-L/parallel-tuned Class-E power amplifier

An analysis of a modified series-L/parallel-tuned Class-E power amplifier is presented, which includes the effects that a shunt capacitance placed across the switching device will have on Class-E behaviour. In the original series L/parallel-tuned topology in which the output transistor capacitance is not inherently included in the circuit, zero-current switching (ZCS) and zero-current derivative switching (ZCDS) conditions should be applied to obtain optimum Class-E operation. On the other hand, when the output transistor capacitance is incorporated in the circuit, i.e. in the modified series-L/parallel-tuned topology, the ZCS and ZCDS would not give optimum operation and therefore zero-voltage-switching (ZVS) and zero-voltage-derivative switching (ZVDS) conditions should be applied instead. In the modified series-L/parallel-tuned Class-E configuration, the output-device inductance and the output-device output capacitance, both of which can significantly affect the amplifier's performance at microwave frequencies, furnish part, if not all, of the series inductance L and the shunt capacitance COUT, respectively. Further, when compared with the classic shunt-C/series-tuned topology, the proposed Class-E configuration offers some advantages in terms of 44% higher maximum operating frequency (fMAX) and 4% higher power-output capability (PMAX). As in the classic topology, the fMAX of the proposed amplifier circuit is reached when the output-device output capacitance furnishes all of the capacitance COUT, for a given combination of frequency, output power and DC supply voltage. It is also shown that numerical simulations agree well with theoretical predictions.