On the coupling of Navier–Stokes and linearised Euler equations for aeroacoustic simulation

Aerodynamic generation of sound is governed by the Navier–Stokes equations while acoustic propagation in a non-uniform medium is effectively described by the linearised Euler equations. Different numerical schemes are required for the efficient solution of these two sets of equations, and therefore, coupling techniques become an essential issue. Two types of one-way coupling between the flow solver and the acoustic solver are discussed: (a) for aerodynamic sound generated at solid surfaces, and (b) in the free stream. Test results indicate how the coupling achieves the necessary accuracy so that Computational Fluid Dynamics codes can be used in aeroacoustic simulations.

Effect of adding 0.3wt% Ni into the Sn-0.7wt% Cu solder - Part II: Growth of intermetallic layer with Cu during wetting and aging

The growth behavior of intermetallic layer with or without adding 0.3 wt% Ni into the Sn-0.7Cu solder was studied during the wetting reaction on Cu-substrate and thereafter in solid-state aging condition. The Cu-solder reaction couple was prepared at 255, 275 and 295 °C for 10 s. The samples reacted at 255 °C were then isothermally aged for 2-14 days at 150 °C. The reaction species formed for the Sn-0.7Cu/Cu and Sn-0.7Cu-0.3Ni/Cu soldering systems were Cu6Sn5 and (CuNi)6Sn5, respectively. The thickness of the intermetallic compounds formed at the solder/Cu interfaces and also in the bulk of both solders increased with the increase of reaction temperature. It was found that Ni-containing Sn-0.7Cu solder exhibited lower growth of intermetallic layer during wetting and in the early stage of aging and eventually exceeded the intermetallic layer thickness of Sn-0.7Cu/Cu soldering system after 6 days of aging. As the aging time proceeds, a non-uniform intermetallic layer growth tendency was observed for the case of Sn-0.7Cu-0.3Ni solder. The growth behavior of intermetallic layer during aging for both solders followed the diffusion-controlled mechanism. The intermetallic layer growth rate constants for Sn-0.7Cu and Sn-0.7Cu-0.3Ni solders were calculated as 1.41 × 10-17 and 1.89 × 10-17 m2/s, respectively which indicated that adding 0.3 wt% Ni with Sn-0.7Cu solder contributed to the higher growth of intermetallic layer during aging. © 2006 Elsevier B.V. All rights reserved.

Fracture mechanics analysis of cracks in solder joint intermetallic compounds

The trend towards miniaturization of electronic products leads to the need for very small sized solder joints. Therefore, there is a higher reliability risk that too large a fraction of solder joints will transform into Intermetallic Compounds (IMCs) at the solder interface. In this paper, fracture mechanics study of the IMC layer for SnPb and Pb-free solder joints was carried out using finite element numerical computer modelling method. It is assumed that only one crack is present in the IMC layer. Linear Elastic Fracture Mechanics (LEFM) approach is used for parametric study of the Stress Intensity Factors (SIF, KI and KII), at the predefined crack in the IMC layer of solder butt joint tensile sample. Contrary to intuition, it is revealed that a thicker IMC layer in fact increases the reliability of solder joint for a cracked IMC. Value of KI and KII are found to decrease with the location of the crack further away from the solder interfaces while other parameters are constant. Solder thickness and strain rate were also found to have a significant influence on the SIF values. It has been found that soft solder matrix generates non-uniform plastic deformation across the solder-IMC interface near the crack tip that is responsible to obtain higher KI and KII.

Models for magnetohydrodynamics of aluminium electrolysis cells

The electric current and the associated magnetic field in aluminium electrolysis cells create effects limiting the cell productivity and possibly cause instabilities: surface waving, ‘anode effects’, erosion of pot lining, feed material sedimentation, etc. The instructive analysis is presented via a step by step inclusion of different physical coupling factors affecting the magnetic field, electric current, velocity and wave development in the electrolysis cells. The full time dependent model couples the nonlinear turbulent fluid dynamics and the extended electromagnetic field in the cell, and the whole bus bar circuit with the ferromagnetic effects. Animated examples for the high amperage cells are presented. The theory and numerical model of the electrolysis cell is extended to the cases of variable cell bottom of aluminium layer and the variable thickness of the electrolyte due to the anode non-uniform burn-out process and the presence of the anode channels. The problem of the channel importance is well known Moreau-Evans model) for the stationary interface and the velocity field, and was validated against measurements in commercial cells, particularly with the recently published ‘benchmark’ test for the MHD models of aluminium cells [1]. The presence of electrolyte channels requires also to reconsider the previous magnetohydrodynamic instability theories and the dynamic wave development models. The results indicate the importance of a ‘sloshing’ parametrically excited MHD wave development in the aluminium production cells.

Determining surface tension using DC magnetic levitation

The values of material physical properties are vital for the successful use of numerical simulations for electromagnetic processing of materials. The surface tension of materials can be determined from the experimental measurement of the surface oscillation frequency of liquid droplets. In order for this technique to be used, a positioning field is required that results in a modification to the oscillation frequency. A number of previous analytical models have been developed that mainly focus on electrically conducting droplets positioned using an A.C. electromagnetic field, but due to the turbulent flow resulting from the high electromagnetic fields required to balance gravity, reliable measurements have largely been limited to microgravity. In this work axisymmetric analytical and numerical models are developed, which allow the surface tension of a diamagnetic droplet positioned in a high DC magnetic field to be determined from the surface oscillations. In the case of D.C. levitation there is no internal electric currents with resulting Joule heating, Marangoni flow and other effects that introduce additional physics that complicates the measurement process. The analytical solution uses the linearised Navier-Stokes equations in the inviscid case. The body force from a DC field is potential, in contrast to the AC case, and it can be derived from Maxwell equations giving a solution for the magnetic field in the form of a series expansion of Legendre polynomials. The first few terms in this expansion represent a constant and gradient magnetic field valid close to the origin, which can be used to position the droplet. Initially the mathematical model is verified in microgravity conditions using a numerical model developed to solve the transient electromagnetics, fluid flow and thermodynamic equations. In the numerical model (as in experiment) the magnetic field is obtained using electrical current carrying coils, which provides the confinement force for a liquid droplet. The model incorporates free surface deformation to accurately model the oscillations that result from the interaction between the droplet and the non-uniform external magnetic field. A comparison is made between the analytical perturbation theory and the numerical pseudo spectral approximation solutions for small amplitude oscillations.

Thermoelectric MHD effects on equiaxed crystal morphology

The effects of a constant uniform magnetic field on a growing equiaxed crystal are investigated using a 3-dimensional enthalpy based numerical model. Two cases are considered: The first case looks at unconstrained growth, where the current density is generated through the thermo-electric effect and the current circulates between the tips and roots of the dendrite, the second represents an imposed potential difference across the domain. A jump in the electrical conductivity between the liquid and solid causes the current density to be non uniform. In both cases the resulting Lorentz force drives fluid flow in the liquid phase, this in turn causes advection of the thermal and solute field altering the free energy close to the interface and changing the morphology of the dendrite. In the first case the flow field is complex comprising of many circulations, the morphological changes are modelled using a 2D model with a quasi 3D approximation. The second case is comparable to classic problems involving a constant velocity boundary.

Tephra-linked peat humification records from Irish ombrotrophic bogs question nature of solar forcing at 850 cal. yr BC.

This paper investigates evidence for palaeoclimatic changes during the period ca. 1500-500 cal. yr BC through peat humification studies on seven Irish ombrotrophic bogs. The sites are well-correlated by the identification of three mid-first millennium BC tephras, which enable the humification records at specific points in time to be directly compared. Phases of temporarily increased wetness are suggested at ca. 1300-1250 cal. yr BC, ca. 1150-1050 cal. yr BC, ca. 940 cal. yr BC and ca. 740 cal. yr BC. The last of these is confirmed to be synchronous at five sites, suggesting external forcing on a regional scale. The timing of this wet-shift is constrained by two closely dated tephras and is demonstrated to be distinct from the widely reported changes to cooler/wetter conditions associated with a solar minimum at 850-760 cal. yr BC, at which time the Irish sites appear instead to experience drier conditions. The results suggest the possibility of either non-uniform responses to solar forcing in northwest Europe at this time, or the existence of unrelated climate events in the early first millennium BC. The findings caution against the correlation of loosely dated palaeoclimate data if the effects of forcing mechanisms are to be understood.

Plane wave and Coulomb asymptotics

A simple plane wave solution of the Schrodinger-Helmholtz equation is a quantum eigenfunction obeying both energy and linear momentum correspondence principles. Inclusion of the outgoing wave with scattering amplitude f asymptotic development of the plane wave, we show that there is a problem with angular momentum when we consider forward scattering at the point of closest approach and at large impact parameter given semiclassically by (l + 1/2)/k where l is the azimuthal quantum number and may be large (J. Leech et al., Phys. Rev. Lett. 88. 257901 (2002)). The problem is resolved via non- uniform, non-standard analysis involving the Heaviside step function, unifying classical, semiclassical and quantum mechanics, and the treatment is extended to the case of pure Coulomb scattering.

Plasma dynamics in an inductively coupled magnetic neutral loop discharge

Magnetic neutral loop discharges (NLDs) can be operated at significantly lower pressures than conventional radio-frequency (rf) inductively coupled plasmas (ICPs). These low pressure conditions are favourable for technological applications, in particular anisotropic etching. An ICP–NLD has been designed providing excellent diagnostics access for detailed investigations of fundamental mechanisms. Spatially resolved Langmuir probe measurements have been performed in the plasma production region (NL region) as well as in the remote application region downstream from the NL region. Depending on the NL gradient two different operation modes have been observed exhibiting different opportunities for control of plasma uniformity. The efficient operation at comparatively low pressures results in ionization degrees exceeding 1%. In this regime neutral dynamics has to be considered and can influence neutral gas and process uniformity. Neutral gas depletion through elevated gas temperatures and high ionization rates have been quantified. At pressures above 0.1 Pa, gas heating is the dominant depletion mechanism. At lower pressures neutral gas is predominantly depleted through high ionization rates and rapid transport of ions by ambipolar diffusion along the magnetic field lines. Non-uniform profiles of the ionization rate can, therefore, result in localized neutral gas depletion and non-uniform processing. We have also investigated the electron dynamics within the radio-frequency cycle using phase resolved optical emission spectroscopy and Thomson scattering. In these measurements electron drift phenomena along the NL torus have been identified.

Neutral gas depletion mechanisms in dense low-temperature argon plasmas

Neutral gas depletion mechanisms are investigated in a dense low-temperature argon plasma-an inductively coupled magnetic neutral loop (NL) discharge. Gas temperatures are deduced from the Doppler profile of the 772.38 nm line absorbed by argon metastable atoms. Electron density and temperature measurements reveal that at pressures below 0.1 Pa, relatively high degrees of ionization (exceeding 1%) result in electron pressures, p(e) = kT(e)n(e), exceeding the neutral gas pressure. In this regime, neutral dynamics has to be taken into account and depletion through comparatively high ionization rates becomes important. This additional depletion mechanism can be spatially separated due to non-uniform electron temperature and density profiles (non-uniform ionization rate), while the gas temperature is rather uniform within the discharge region. Spatial profiles of the depletion of metastable argon atoms in the NL region are observed by laser induced fluorescence spectroscopy. In this region, the depletion of ground state argon atoms is expected to be even more pronounced since in the investigated high electron density regime the ratio of metastable and ground state argon atom densities is governed by the electron temperature, which peaks in the NL region. This neutral gas depletion is attributed to a high ionization rate in the NL zone and fast ion loss through ambipolar diffusion along the magnetic field lines. This is totally different from what is observed at pressures above 10 Pa where the degree of ionization is relatively low (

A study of the biological effects of modulated 6 MV radiation fields

The delivery of spatially modulated radiation fields has been shown to impact on in vitro cell survival responses. To study the effect of modulated fields on cell survival, dose response curves were determined for human DU-145 prostate, T98G glioma tumour cells and normal primary AGO-1552 fibroblast cells exposed to modulated and non-modulated field configurations delivered using a 6 MV Linac with multi-leaf collimator. When exposed to uniform fields delivered as a non-modulated or modulated configuration, no significant differences in survival were observed with the exception of DU-145 cells at a dose of 8 Gy (p = 0.024). Survival responses were determined for exposure to non-uniform-modulated beams in DU-145 and T98G and showed no deviation from the survival response observed following uniform non-modulated exposures. The results of these experiments indicate no major deviation in response to modulated fields compared to uniform exposures.

Modelling the mechanical response of the Reed-mouthpiece-lip system of a clarinet. Part I. A one-dimensional distributed model

The motion of a clarinet reed that is clamped to a mouthpiece and supported by a lip is simulated in the time-domain using finite difference methods. The reed is modelled as a bar with non-uniform cross section, and is described using a one-dimensional, fourth-order partial differential equation. The interactions with the mouthpiece Jay and the player's lip are taken into account by incorporating conditional contact forces in the bar equation. The model is completed by clamped-free boundary conditions for the reed. An implicit finite difference method is used for discretising the system, and values for the physical parameters are chosen both from laboratory measurements and by accurate tuning of the numerical simulations. The accuracy of the numerical system is assessed through analysis of frequency warping effects and of resonance estimation. Finally, the mechanical properties of the system are studied by analysing its response to external driving forces. In particular, the effects of reed curling are investigated.

Single-phase fluid flow distribution and heat transfer in microstructured reactors

Single-phase microreactors and micro-heat-exchangers have been widely used in industrial and scientific applications over the last decade. In several cases, operation of microreactors has shown that their expected efficiency cannot be reached either due to non-uniform distribution of reactants between different channels or due to flow maldistribution between individual microreactors working in parallel. The latter problem can result in substantial temperature deviations between different microreactors resulting in thermal run away which could arise from an exothermicreaction. Thus advances in the understanding of heat transfer and fluid flow distribution continue to be crucial in achieving improved performance, efficiency and safety in microstructured reactors used for different applications. This paper presents a review of the experimental and numerical results on fluid flow distribution, heat transfer and combination thereof, available in the open literature. Heat transfer in microchannels can be suitably described by standard theory and correlations, but scaling effects (entrance effects, conjugate heat transfer, viscous heating, and temperature-dependent properties) have often to be accounted for in microsystems. Experiments with single channels are in good agreement with predictions from the published correlations. The accuracy of multichannel experiments is lower due to flow maldistribution. Special attention is devoted to theoretical and experimental studies on the effect of a flow maldistribution on the thermal and conversion response of catalytic microreactors. There view concludes with a set of design recommendations aimed at improving the reactor performance. (C) 2010 Elsevier Ltd. All rights reserved.

Relativistic laser pulse compression in plasmas with a linear axial density gradient

The self-compression of a relativistic Gaussian laser pulse propagating in a non-uniform plasma is investigated. A linear density inhomogeneity (density ramp) is assumed in the axial direction. The nonlinear Schrodinger equation is first solved within a one-dimensional geometry by using the paraxial formalism to demonstrate the occurrence of longitudinal pulse compression and the associated increase in intensity. Both longitudinal and transverse self-compression in plasma is examined for a finite extent Gaussian laser pulse. A pair of appropriate trial functions, for the beam width parameter (in space) and the pulse width parameter (in time) are defined and the corresponding equations of space and time evolution are derived. A numerical investigation shows that inhomogeneity in the plasma can further boost the compression mechanism and localize the pulse intensity, in comparison with a homogeneous plasma. A 100 fs pulse is compressed in an inhomogeneous plasma medium by more than ten times. Our findings indicate the possibility for the generation of particularly intense and short pulses, with relevance to the future development of tabletop high-power ultrashort laser pulse based particle acceleration devices and associated high harmonic generation. An extension of the model is proposed to investigate relativistic laser pulse compression in magnetized plasmas.

Optimization of heat transfer characteristics, flow distribution, and reaction processing for a microstructured reactor/heat-exchanger for optimal performance in platinum catalyzed ammonia oxidation

The present work is focused on the demonstration of the advantages of miniaturized reactor systems which are essential for processes where potential for considerable heat transfer intensification exists as well as for kinetic studies of highly exothermic reactions at near-isothermal conditions. The heat transfer characteristics of four different cross-flow designs of a microstructured reactor/heat-exchanger (MRHE) were studied by CFD simulation using ammonia oxidation on a platinum catalyst as a model reaction. An appropriate distribution of the nitrogen flow used as a coolant can decrease drastically the axial temperature gradient in the reaction channels. In case of a microreactor made of a highly conductive material, the temperature non-uniformity in the reactor is strongly dependent on the distance between the reaction and cooling channels. Appropriate design of a single periodic reactor/heat-exchanger unit, combined with a non-uniform inlet coolant distribution, reduces the temperature gradients in the complete reactor to less than 4degreesC, even at conditions corresponding to an adiabatic temperature rise of about 1400degreesC, which are generally not accessible in conventional reactors because of the danger of runaway reactions. To obtain the required coolant flow distribution, an optimization study was performed to acquire the particular geometry of the inlet and outlet chambers in the microreactor/heat-exchanger. The predicted temperature profiles are in good agreement with experimental data from temperature sensors located along the reactant and coolant flows. The results demonstrate the clear potential of microstructured devices as reliable instruments for kinetic research as well as for proper heat management in the case of highly exothermic reactions. (C) 2002 Elsevier Science B.V. All rights reserved.