A numerical study of laminar axisymmetric plumes due to the combined buoyancy of heat and mass diffusion

This paper presents the results of a computational study of laminar axisymmetric plumes generated by the simultaneous diffusion of thermal energy and chemical species. Species concentrations are assumed small. The plume is treated as a boundary layer. Boussinesq approximations are incorporated and the governing conservation equations of mass, momentum, energy and species are suitably non-dimensionalised. These equations are solved using one time-step-forward explicit finite-difference method. Upwind differencing is employed for convective terms. The results thus obtained are explained in terms of the basic physical mechanisms that govern these flows. They show many interesting aspects of the complex interaction of the two buoyant mechanisms.

Phase-shift randomness in one-dimensional disordered conductors in the quasi-metallic domain

An invariant imbedding method yields exact analytical results for the distribution of the phase theta (L) of the reflection amplitude and for low-order resistance moments (pn) for a disordered conductor of length L in the quasi-metallic regime L<momentum. For realistic values of 2k0L, the resistance moments are dominated by the terms obtained previously by arbitrarily assuming uniformly distributed phases. A direct proof of the validity of the random-phase assumption for studying scaling of resistance in one dimension is thus obtained.

Numerical simulation of solidification of liquid aluminum alloy flowing on cooling slope

Preparation of semisolid slurry using a cooling slope is increasingly becoming popular, primarily because of the simplicity in design and ease control of the process. In this process, liquid alloy is poured down an inclined surface which is cooled from underneath. The cooling enables partial solidification and the incline provides the necessary shear for producing semisolid slurry. However, the final microstructure of the ingot depends on several process parameters such as cooling rate, incline angle of the cooling slope, length of the slope and initial melt superheat. In this work, a CFD model using volume of fluid (VOF) method for simulating flow along the cooling slope was presented. Equations for conservation of mass, momentum, energy and species were solved to predict hydrodynamic and thermal behavior, in addition to predicting solid fraction distribution and macrosegregation. Solidification was modeled using an enthalpy approach and a volume averaged technique for the different phases. The mushy region was modeled as a multi-layered porous medium consisting of fixed columnar dendrites and mobile equiaxed/fragmented grains. The alloy chosen for the study was aluminum alloy A356, for which adequate experimental data were available in the literature. The effects of two key process parameters, namely the slope angle and the pouring temperature, on temperature distribution, velocity distribution and macrosegregation were also studied.

Dynamics of sheared inelastic dumbbells

We study the dynamical properties of the homogeneous shear flow of inelastic dumbbells in two dimensions as a first step towards examining the effect of shape on the properties of flowing granular materials. The dumbbells are modelled as smooth fused disks characterized by the ratio of the distance between centres (L) and the disk diameter (D), with an aspect ratio (L/D) varying between 0 and 1 in our simulations. Area fractions studied are in the range 0.1-0.7, while coefficients of normal restitution (e(n)) from 0.99 to 0.7 are considered. The simulations use a modified form of the event-driven methodology for circular disks. The average orientation is characterized by an order parameter S, which varies between 0 (for a perfectly disordered fluid) and 1 (for a fluid with the axes of all dumbbells in the same direction). We investigate power-law fits of S as a function of (L D) and (1 - e(n)(2)) There is a gradual increase in ordering as the area fraction is increased, as the aspect ratio is increased or as the coefficient of restitution is decreased. The order parameter has a maximum value of about 0.5 for the highest area fraction and lowest coefficient of restitution considered here. The mean energy of the velocity fluctuations in the flow direction is higher than that in the gradient direction and the rotational energy, though the difference decreases as the area fraction increases, due to the efficient collisional transfer of energy between the three directions. The distributions of the translational and rotational velocities are Gaussian to a very good approximation. The pressure is found to be remarkably independent of the coefficient of restitution. The pressure and dissipation rate show relatively little variation when scaled by the collision frequency for all the area fractions studied here, indicating that the collision frequency determines the momentum transport and energy dissipation, even at the lowest area fractions studied here. The mean angular velocity of the particles is equal to half the vorticity at low area fractions, but the magnitude systematically decreases to less than half the vorticity as the area fraction is increased, even though the stress tensor is symmetric.

Spiral states in the square-lattice Hubbard model

We present a variety of physical implications of a mean-field theory for spiral spin-density-wave states in the square-lattice Hubbard model for small deviations from half filling. The phase diagram with the paramagnetic metal, two spiral (semimetallic) states, and ferromagnet is calculated. The momentum distribution function and the (quasiparticle) density of states are discussed. There is a significant broadening of the quasiparticle bands when the antiferromagnetic insulator is doped. The evolution of the Fermi surface and the variation of the plasma frequency and a charge-stiffness constant with U/t and δ are calculated. The connection to results based on the Schwinger-boson-slave-fermion formalism is made.

Steady incompressible flow of cohesionless granular materials through a wedge-shaped hopper: Frictional 14kinetic solution to the smooth wall, radial gravity problem

Hybrid frictional-kinetic equations are used to predict the velocity, grain temperature, and stress fields in hoppers. A suitable choice of dimensionless variables permits the pseudo-thermal energy balance to be decoupled from the momentum balance. These balances contain a small parameter, which is analogous to a reciprocal Reynolds number. Hence an approximate semi-analytical solution is constructed using perturbation methods. The energy balance is solved using the method of matched asymptotic expansions. The effect of heat conduction is confined to a very thin boundary layer near the exit, where it causes a marginal change in the temperature. Outside this layer, the temperature T increases rapidly as the radial coordinate r decreases. In particular, the conduction-free energy balance yields an asymptotic solution, valid for small values of r, of the form T proportional r-4. There is a corresponding increase in the kinetic stresses, which attain their maximum values at the hopper exit. The momentum balance is solved by a regular perturbation method. The contribution of the kinetic stresses is important only in a small region near the exit, where the frictional stresses tend to zero. Therefore, the discharge rate is only about 2.3% lower than the frictional value, for typical parameter values. As in the frictional case, the discharge rate for deep hoppers is found to be independent of the head of material.

Boundary layer swirling flow through a conical hydrocyclone

The paper deals with the flow and heat-transfer problem of a steady axisymmetric laminar incompressible boundary layer swirling flow of a fluid through a conical hydrocyclone. The implicit finitedifference scheme is used to solve the partial differential equations governing the flow. The effect of swirl is found to be more pronounced on the longitudinal skin friction than on the tangential skin friction and heat transfer. The skin friction and heat transfer increase with swirl or with longitudinal distance. Swirl also gives rise to velocity overshoot in the longitudinal velocity profiles and the magnitude of the velocity overshoot increases as the swirl parameter increases. The results are found to be in good agreement with those of the local nonsimilarity and momentum integral methods but they differ appreciably from those of the local similarity method except for the longitudinal skin friction which is fairly in good agreement with that of the local similarity method.Die Arbeit beschäftigt sich mit der Strömung und dem Wärmeübergang in einem konischen Zyklon unter der Voraussetzung stationärer, achsensymmetrischer, laminarer, inkompressibler Grenzschichtströmung. Ein implizites Differenzenverfahren wird benutzt, um die partiellen Differentialgleichungen zu lösen. Der Einfluß des Dralls ist besonders ausgeprägt auf die longitudinale Komponente der Oberflächenreibung, weniger dagegen bei der tangentialen Komponente und beim Wärmeübergang. Die Oberflächenreibung und der Wärmeübergang nehmen zu mit dem Drall, sowie mit dem longitudinalen Abstand. Der Drall erzeugt ein Überschießen der Geschwindigkeit in der longitudinalen Abstand. Der Drall erzeugt ein Überschießen der Geschwindigkeit in der Längsrichtung. Die Größe des Überschusses nimmt mit wachsendem Drallparameter zu. Die Resultate stimmen gut mit den Ergebnissen der Theorie der lokalen Nichtähnlichkeit und der Impulsintegralmethode überein. Dagegen weichen sie mit Ausnahme der longitudinalen Komponente der Oberflächenreibung beträchtlich von den Resultaten der Theorie der lokalen Ähnlichkeit ab.

Wigner distribution for angle coordinates in quantum mechanics

The method of Wigner distribution functions, and the Weyl correspondence between quantum and classical variables, are extended from the usual kind of canonically conjugate position and momentum operators to the case of an angle and angular momentum operator pair. The sense in which one has a description of quantum mechanics using classical phase‐space language is much clarified by this extension.

Analysis of conjugate laminar mixed convection cooling in a shrouded array of electronic components

The temperature variation in the insulation around an electronic component, mounted on a horizontal circuit board is studied numerically. The flow is assumed to be laminar and fully developed. The effect of mixed convection and two different types of insulation are considered. The mass, momentum and energy conservation equations in the fluid and conduction equation in the insulation are solved using the SIMPLER algorithm. Computations are carried out for liquid Freon and water, for different conductivity ratios, and different Rayleigh numbers. It is demonstrated that the temperature variation within the insulation becomes important when the thermal conductivity of the insulation is less than ten times the thermal conductivity of the cooling medium.

Time-Dependent Dynamics of a Planar Galaxy Model

Time-dependent models of collisionless stellar systems with harmonic potentials allowing for an essentially exact analytic description have recently been described. These include oscillating spheres and spheroids. This paper extends the analysis to time-dependent elliptic discs. Although restricted to two space dimensions, the systems are richer in that their parameters form a 10-dimensional phase space (in contrast to six for the earlier models). Apart from total energy and angular momentum, two additional conserved quantities emerge naturally. These can be chosen as the areas of extremal sections of the ellipsoidal region of phase space occupied by the system (their product gives the conserved volume). The present paper describes the construction of these models. An application to a tidal encounter is given which allows one to go beyond the impulse approximation and demonstrates the effects of rotation of the perturbed system on energy and angular-momentum transfer. The angular-momentum transfer is shown to scale inversely as the cube of the encounter velocity for an initial configuration of the perturbed galaxy with zero quadrupole moment.

Radio wave heating of astrophysical plasmas

The coherent plasma process such as parametric decay instability (PDI) has been applied to a homogeneous and unmagnetized plasma. These instabilities cause anomalous absorption of strong electromagnetic radiation under specific conditions of energy and momentum conservation and thus cause anomalous heating of the plasma. The maximum plasma temperatures reached are functions of luminosity of the radio radiation and plasma parameters. We believe that these processes may be taking place in many astrophysical objects. Here, the conditions in the sources 3C 273, 3C 48 and Crab Nebula are shown to be conducive to the excitation of PDI. These processes also contribute towards the absorption of 21cm radiation

Effect of Turbulence on Emerging Magnetic Flux Tubes in the Convection Zone

Working under the hypothesis that magnetic flux in the sun is generated at the bottom of the convection zone, Choudhuri and Gilman (1987; Astrophys. J. 316, 788) found that a magnetic flux tube symmetric around the rotation axis, when released at the bottom of the convection zone, gets deflected by the Coriolis force and tends to move parallel to the rotation axis as it rises in the convection zone. As a result, all the flux emerges at rather high latitudes and the flux observed at the typical sunspot latitudes remains unexplained. Choudhuri (1989; Solar Physics, in press) finds that non-axisymmetric perturbations too cannot subdue the Coriolis force. In this paper, we no longer treat the convection zone to be passive as in the previous papers, but we consider the role of turbulence in the convection zone in inhibiting the Coriolis force. The interaction of the flux tubes with the turbulence is treated in a phenomenological way as follows: (1) Large scale turbulence on the scale of giant cells can physically drag the tubes outwards, thus pulling the flux towards lower latitudes by dominating over the Coriolis force. (2) Small scale turbulence of the size of the tubes can exchange angular momentum with the tube, thus suppressing the growth of the Coriolis force and making the tubes emerge at lower latitudes. Numerical simulations show that the giant cells can drag the tubes and make them emerge at lower latitufes only if the velocities within the giant cells are unrealistically large of if the radii of the flux tubes are as small as 10 km. However, small scale turbulence can successfully suppress the growth of the Coriolis force if the tubes have radii smaller than about 300 km which may not be unreasonable. Such flux tubes can then emerge at low latitudes where sunspots are seen.

Can equipartition fields produce the tilts of bipolar magnetic regions?

The effect of turbulence on the nonaxisymmetric flux rings of equipartition field strength in bipolar magnetic regions is studied on the basis of the small-scale momentum exchange mechanism and the giant cell drag combined with the Kelvin-Helmholtz drag mechanism. It is shown that the giant cell drag and small-scale momentum exchange mechanism can make equipartition flux loops emerge at low latitudes, in addition to making them exhibit the observed tilts. However, the sizes of the flux tubes have to be restricted to a couple of hundred kilometers. An ad hoc constraint on the footpoints of the flux loops is introduced by not letting them move in the phi direction, and it is found that equipartition fields of any size can be made to emerge at sunspot latitudes with the observed tilts by suitably adjusting the footpoint separations.

Kaon B parameter with Wilson fermions

We calculate the kaon B parameter in quenched lattice QCD at beta=6.0 using Wilson fermions at kappa=0.154 and 0.155. We use two kinds of nonlocal (''smeared'') sources for quark propagators to calculate the matrix elements between states of definite momentum. The use of smeared sources yields results with much smaller errors than obtained in previous calculations with Wilson fermions. By combining results for p=(0,0,0) and p=(0,0,1), we show that one can carry out the noperturbative subtraction necessary to remove the dominant lattice artifacts induced by the chiral-symmetry-breaking term in the Wilson action. Our final results are in good agreement with those obtained using staggered fermions. We also present results for B parameters of the DELTAI = 3/2 part of the electromagnetic penguin operators, and preliminary results for B(K) in the presence of two flavors of dynamical quarks.

Effect of uncertainty on helicopter performance predictions

The effect of uncertainties on performance predictions of a helicopter is studied in this article. The aeroelastic parameters such as the air density, blade profile drag coefficient, main rotor angular velocity, main rotor radius, and blade chord are considered as uncertain variables. The propagation of these uncertainties in the performance parameters such as thrust coefficient, figure of merit, induced velocity, and power required are studied using Monte Carlo simulation and the first-order reliability method. The Rankine-Froude momentum theory is used for performance prediction in hover, axial climb, and forward flight. The propagation of uncertainty causes large deviations from the baseline deterministic predictions, which undoubtedly affect both the achievable performance and the safety of the helicopter. The numerical results in this article provide useful bounds on helicopter power requirements.