On the computation of two-dimensional compressible flow fields by the modified local stability scheme

The modified local stability scheme is applied to several two-dimensional problems—blunt body flow, regular reflection of a shock and lambda shock. The resolution of the flow features obtained by the modified local stability scheme is found to be better than that achieved by the other first order schemes and almost identical to that achieved by the second order schemes incorporating artificial viscosity. The scheme is easy for coding, consumes moderate amount of computer storage and time. The scheme can be advantageously used in place of second order schemes.

Studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring

Numerical and experimental studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring are performed. The alloy is electromagnetically stirred to produce semisolid slurry in a cylindrical graphite mould placed in the annulus of a linear electromagnetic stirrer. The mould is cooled at the bottom, such that solidification progresses from the bottom to the top of the cylindrical mould. A numerical model is developed for simulating the transport phenomena associated with the solidification process using a set of single-phase governing equations of mass. momentum, energy. and species conservation. The viscosity variation of the slurry, used in the model, is determined experimentally using a rotary viscometer. The set of governing equations is solved using a pressure-based finite volume technique, along with an enthalpy based phase change algorithm. The numerical study involves prediction of temperature, velocity, species and solid fraction distribution in the mould. Corresponding solidification experiments are performed, with time-temperature history recorded at key locations. The microstructures at various temperature measurement locations in the solidified billet are analyzed. The numerical predictions of temperature variations are in good agreement with experiments, and the predicted flow field evolution correlates well with the microstructures observed at various locations.

Multicarrier analysis of magnetotransport data at low and high electric fields

We present some results on multicarrier analysis of magnetotransport data, Both synthetic as well as data from narrow gap Hg0.8Cd0.2Te samples are used to demonstrate applicability of various algorithms vs. nonlinear least square fitting, Quantitative Mobility Spectrum Analysis (QMSA) and Maximum Entropy Mobility Spectrum Analysis (MEMSA). Comments are made from our experience oil these algorithms, and, on the inversion procedure from experimental R/sigma-B to S-mu specifically with least square fitting as an example. Amongst the conclusions drawn are: (i) Experimentally measured resistivity (R-xx, R-xy) should also be used instead of just the inverted conductivity (sigma(xx), sigma(xy)) to fit data to semiclassical expressions for better fits especially at higher B. (ii) High magnetic field is necessary to extract low mobility carrier parameters. (iii) Provided the error in data is not large, better estimates to carrier parameters of remaining carrier species can be obtained at any stage by subtracting highest mobility carrier contribution to sigma from the experimental data and fitting with the remaining carriers. (iv)Even in presence of high electric field, an approximate multicarrier expression can be used to guess the carrier mobilities and their variations before solving the full Boltzmann equation.

Electron energy distributions and transport coefficients in N2 in E times B fields

Using the concept of energy-dependent effective field intensity, electron transport coefficients in nitrogen have been determined in E times B fields (E = electric field intensity, B = magnetic flux density) by the numerical solution of the Boltzmann transport equation for the energy distribution of electrons. It has been observed that as the value of B/p (p = gas pressure) is increased from zero, the perpendicular drift velocity increased linearly at first, reaches a maximum value, and then decreases with increasing B/p. In general, the electron mean energy is found to be a function of Eavet/p( Eavet = averaged effective electric field intensity) only, but the other transport coefficients, such as transverse drift velocity, perpendicular drift velocity, and the Townsend ionization coefficient, are functions of both E/p and B/p.

Mixed mode (I and II) crack tip fields in bulk metallic glasses

Stationary crack tip fields in bulk metallic glasses under mixed mode (I and II) loading are studied through detailed finite element simulations assuming plane strain, small scale yielding conditions. The influence of internal friction or pressure sensitivity on the plastic zones. notch deformation, stress and plastic strain fields is examined for different mode mixities. Under mixed mode loading, the notch deforms into a shape such that one part of its surface sharpens while the other part blunts. Increase in mode If component of loading dramatically enhances the normalized plastic zone size, lowers the stresses but significantly elevates the plastic strain levels near the notch tip. Higher internal friction reduces the peak tangential stress but increases the plastic strain and stretching near the blunted part of the notch. The simulated shear bands are straight and extend over a long distance ahead of the notch tip under mode II dominant loading. The possible variations of fracture toughness with mode mixity corresponding to failure by brittle micro-cracking and ductile shear banding are predicted employing two simple fracture criteria. The salient results from finite element simulations are validated by comparison with those from mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass.

Stationary crack tip fields in elastic-plastic solids: an overview of recent numerical simulations

In this paper, an overview of some recent numerical simulations of stationary crack tip fields in elastic-plastic solids is presented. First, asymptotic analyses carried out within the framework of 2D plane strain or plane stress conditions in both pressure insensitive and pressure sensitive plastic solids are reviewed. This is followed by discussion of salient results obtained from recent computational studies. These pertain to 3D characteristics of elastic-plastic near-front fields under mixed mode loading, mechanics of fracture and simulation of near-tip shear banding process of amorphous alloys and influence of crack tip constraint on the structure of near-tip fields in ductile single crystals. These results serve to illustrate several important features associated with stress and strain distributions near the crack tip and provide the foundation for understanding the operative failure mechanisms. The paper concludes by highlighting some of the future prospects for this field of study.

Crack Tip Fields in a Single Edge Notched Aluminum Single Crystal Specimen

We report a combined experimental and computational study of a low constraint aluminum single crystal fracture geometry and investigate the near-tip stress and strain fields. To this end, a single edge notched tensile (SENT) specimen is considered. A notch, with a radius of 50 µm, is taken to lie in the (010) plane and its front is aligned along the [101] direction. Experiments are conducted by subjecting the specimen to tensile loading using a special fixture inside a scanning electron microscope chamber. Both SEM micrographs and electron back-scattered diffraction (EBSD) maps are obtained from the near-tip region. The experiments are complemented by performing 3D and 2D plane strain finite element simulations within a continuum crystal plasticity framework assuming an isotropic hardening response characterized by the Pierce–Asaro–Needleman model. The simulations show a distinct slip band forming at about 55 deg with respect to the notch line corresponding to slip on (11-bar 1)[011] system, which corroborates well with experimental data. Furthermore, two kink bands occur at about 45 deg and 90 deg with respect to the notch line within which large rotations in the crystal orientation take place. These predictions are in good agreement with the EBSD observations. Finally, the near-tip angular variations of the 3D stress and plastic strain fields in the low constraint SENT fracture geometry are examined in detail.

Effect of process parameters on transport phenomena during solidification in the presence of electromagnetic stirring

In the present work, a numerical study is performed to predict the effect of process parameters on transport phenomena during solidification of aluminium alloy A356 in the presence of electromagnetic stirring. A set of single-phase governing equations of mass, momentum, energy and species conservation is used to represent the solidification process and the associated fluid flow, heat and mass transfer. In the model, the electromagnetic forces are incorporated using an analytical solution of Maxwell equation in the momentum conservation equations and the slurry rheology during solidification is represented using an experimentally determined variable viscosity function. Finally, the set of governing equations is solved for various process conditions using a pressure based finite volume technique, along with an enthalpy based phase change algorithm. In present work, the effect of stirring intensity and cooling rate are considered. It is found that increasing stirring intensity results in increase of slurry velocity and corresponding increase in the fraction of solid in the slurry. In addition, the increasing stirring intensity results uniform distribution of species and fraction of solid in the slurry. It is also found from the simulation that the distribution of solid fraction and species is dependent on cooling rate conditions. At low cooling rate, the fragmentation of dendrites from the solid/liquid interface is more.

Scattering polarization in the presence of magnetic and electric fields

The polarization of radiation by scattering on an atom embedded in combined external quadrupole electric and uniform magnetic fields is studied theoretically. Limiting cases of scattering under Zeeman effect, and Hanle effect in weak magnetic fields are discussed. The theory is general enough to handle scattering in intermediate magnetic fields (Hanle-Zeeman effect) and for arbitrary orientation of magnetic field. The quadrupolar electric field produces asymmetric line shifts, and causes interesting level-crossing phenomena either in the absence of an ambient magnetic field, or in its presence. It is shown that the quadrupolar electric field produces an additional depolarization in the Q/I profiles and rotation of the plane of polarization in the U/I profile over and above that arising from magnetic field itself. This characteristic may have a diagnostic potential to detect steady-state and time-varying electric fields that surround radiating atoms in solar atmospheric layers. (c) 2007 Elsevier Ltd. All rights reserved.

Frequency characteristics of very fast transient currents in a 245-kV GIS

The conducted as well as the induced voltages on control cables and control circuits due to transient electromagnetic (EM) fields generated during switching operations in a gas-insulated substation (GIS) depend on the waveshape of the very fast transient overvoltages and the associated very-fast transient currents (VFTCs). The aim of this paper is to build a basis for characterizing the VFTC generated in gas-insulated switchgear and the,associated equipment during switching operations for the study of transient coupling phenomena. The peak magnitudes of VFTC and their dominant frequency content at various locations have been computed in a 245-kV GIS for different switching operations as well as substation configurations. Finally, the influence of the substation layout on the frequency spectrum, dominant frequencies, and the highest possible frequency component of the VFTC at various distances from the switch have been reported.

Instability of Coupled Longitudinal and Transverse Electromagnetic Modes of Wave‐Propagation in Bounded Streaming Plasmas

The instability of coupled longitudinal and transverse electromagnetic modes associated with long wavelengths is studied in bounded streaming plasmas. The main conclusions are as follows: (i) For long waves for which O (k 2)=0, in the absence of relative streaming motion of electrons and ions and aωp/c<0.66, the whole spectrum of harmonic waves is excited due to finite temperature and boundary effects consisting of two subseries. One of these subseries can be identified with Tonks-Dattner resonance oscillations for the electrons, and arises primarily due to the electrons with frequencies greater than the electrostatic plasma frequency corresponding to the electron density in the midplane in the undisturbed state. The other series arises primarily due to ion motion. When aωp/c>0.66, in addition to the above spectrum of harmonic waves, the system admits an infinite number of growing and decaying waves. The instability associated with these modes is found to arise due to the interaction of the waves inside the plasma with the external electromagnetic field. (ii) For modes with comparatively shorter wavelengths for which O (k3)=0, the coupling due to finite temperature sets in, and it is found that the two series of harmonic waves obtained in (i) deriving energy from the transverse modes also become unstable. Thus, for these wavelengths the system admits three sets of growing and decaying modes, first two for all values of aωp/c and the third for (aωp/c) > 0.66. (iii) The presence of streaming velocities introduces various other coupling mechanisms, and we find that even for the wavelengths for which O (k2)=0, we get three sets of growing and decaying waves. The numerical values for the growth rates show that the streaming velocities enhance the growth rates of instability significantly.

Lightning-Induced Current and Voltage on a Rocket in the Presence of Its Trailing Exhaust Plume

This paper presents time-domain characteristics of induced current and voltage on a rocket in the presence of its exhaust plume when an electromagnetic (EM) wave generated by a nearby lightning discharge is incident on it. For the EM-field interaction with the rocket, the finite-difference time-domain technique has been used. The distributed electrical parameters, such as capacitance and inductance of the rocket and its exhaust plume, are computed using the method of moments technique. For the electrical characterization of the exhaust plume, the computational fluid dynamics technique has been used. The computed peak value of the electrical conductivity of the exhaust plume is 0.12 S/m near the exit plane and it reduces to 0.02 S/m at the downstream end. The relative permittivity varies from 0.91 to 0.99. The exhaust plume behaves as a good conductor for EM fields with frequencies less than 2.285 GHz. It has been observed that the peak value of the induced current on the rocket gets enhanced significantly in the presence of the conducting exhaust plume for the rocket and exhaust plume dimensions and parameters studied. The magnitude of the time-varying induced current at the tail is much more than that of any other section of the rocket.

Lightning-Induced Current and Voltage on a Rocket in the Presence of Its Trailing Exhaust Plume

This paper presents time-domain characteristics of induced current and voltage on a rocket in the presence of its exhaust plume when an electromagnetic (EM) wave generated by a nearby lightning discharge is incident on it. For the EM-field interaction with the rocket, the finite-difference time-domain technique has been used. The distributed electrical parameters, such as capacitance and inductance of the rocket and its exhaust plume, are computed using the method of moments technique. For the electrical characterization of the exhaust plume, the computational fluid dynamics technique has been used. The computed peak value of the electrical conductivity of the exhaust plume is 0.12 S/m near the exit plane and it reduces to 0.02 S/m at the downstream end. The relative permittivity varies from 0.91 to 0.99. The exhaust plume behaves as a good conductor for EM fields with frequencies less than 2.285 GHz. It has been observed that the peak value of the induced current on the rocket gets enhanced significantly in the presence of the conducting exhaust plume for the rocket and exhaust plume dimensions and parameters studied. The magnitude of the time-varying induced current at the tail is much more than that of any other section of the rocket.