Stochastic fluctuations of the solar dynamo

Attempts in the past to model the irregularities of the solar cycle (such as the Maunder minimum) were based on studies of the nonlinear feedback of magnetic fields on the dynamo source terms. Since the alpha-coefficient is obtained by averaging over the turbulence, it is expected to have stochastic fluctuations, and we show that these fluctuations can explain the irregularities of the solar cycle in a more satisfactory way. We solve the dynamo equations in a slab with a single mode, taking the alpha-coefficient to be constant in space but fluctuating stochastically in time with some given amplitude and given correlation time. The same level of percentile fluctuations (about 10 %) produces no effect on an alpha-omega dynamo, but makes an alpha-2 dynamo completely chaotic. The level of irregularities in an alpha-2-omega dynamo qualitatively agrees with the solar behavior, reinforcing the conclusion of Choudhuri (1990a) that the solar dynamo is of the alpha-2-omega-type. The irregularities are found to increase on increasing either the amplitude or the correlation time of the stochastic fluctuations. The alpha-quenching mechanism tends to make the system stable against the irregularities and hence it is inferred that the alpha-quenching should not be too strong so that the irregularities are not completely suppressed. We also present a simple-minded analysis to understand why the stochastic fluctuations in the alpha-omega, alpha-2-omega and alpha-2 regimes have such different outcomes.

Estimation of surface temperature from MONTBLEX data

It is observed that the daily mean temperature of the soil is linear with depth and the variation of the temperature is sinusoidal with a period of a day. Based on these observations the one-dimensional heat conduction equation for the soil can be solved which gives the amplitude and phase variation of the temperature wave with depth. Given the temperature data at three levels below the surface, the amplitude and phase variation and hence the surface temperature variation over the day are estimated. The daily mean temperature of the surface is estimated from linear extrapolation of the daily means at the three levels below the surface. Estimated values of soil thermal diffusivity show a subtantial change after sudden and heavy rains.

Peristaltic transport of a biofluid in a pipe of elliptic cross section

Peristaltic transport of two fluids occupying the peripheral layer and the core in an elliptic tube is, investigated in elliptic cylindrical co-ordinate system, under long wavelength and low Reynolds number approximations. The effect of peripheral-layer viscosity on the flow rate and the frictional force for a slightly elliptic tube is discussed. The limiting results for the one-fluid model are obtained for different eccentricities of the undisturbed tube cross sections with the same area. As a result of non-uniformity of the peristaltic wave, two different amplitude ratios are defined and the time-averaged flux and mechanical efficiency are studied for different eccentricities. It is observed that the time-averaged flux is not affected significantly by the pressure drop when the eccentricity is large. For the peristaltic waves with same area variation, the pumping seems to improve with the eccentricity.

Simulation of ENSO-Related Surface Winds in the Tropical Pacific by an Atmospheric General Circulation Model Forced by Observed Sea Surface Temperatures

The authors present the simulation of the tropical Pacific surface wind variability by a low-resolution (R15 horizontal resolution and 18 vertical levels) version of the Center for Ocean-Land-Atmosphere Interactions, Maryland, general circulation model (GCM) when forced by observed global sea surface temperature. The authors have examined the monthly mean surface winds acid precipitation simulated by the model that was integrated from January 1979 to March 1992. Analyses of the climatological annual cycle and interannual variability over the Pacific are presented. The annual means of the simulated zonal and meridional winds agree well with observations. The only appreciable difference is in the region of strong trade winds where the simulated zonal winds are about 15%-20% weaker than observed, The amplitude of the annual harmonics are weaker than observed over the intertropical convergence zone and the South Pacific convergence zone regions. The amplitudes of the interannual variation of the simulated zonal and meridional winds are close to those of the observed variation. The first few dominant empirical orthogonal functions (EOF) of the simulated, as well as the observed, monthly mean winds are found to contain a targe amount of high-frequency intraseasonal variations, While the statistical properties of the high-frequency modes, such as their amplitude and geographical locations, agree with observations, their detailed time evolution does not. When the data are subjected to a 5-month running-mean filter, the first two dominant EOFs of the simulated winds representing the low-frequency EI Nino-Southern Oscillation fluctuations compare quite well with observations. However, the location of the center of the westerly anomalies associated with the warm episodes is simulated about 15 degrees west of the observed locations. The model simulates well the progress of the westerly anomalies toward the eastern Pacific during the evolution of a warm event. The simulated equatorial wind anomalies are comparable in magnitude to the observed anomalies. An intercomparison of the simulation of the interannual variability by a few other GCMs with comparable resolution is also presented. The success in simulation of the large-scale low-frequency part of the tropical surface winds by the atmospheric GCM seems to be related to the model's ability to simulate the large-scale low-frequency part of the precipitation. Good correspondence between the simulated precipitation and the highly reflective cloud anomalies is seen in the first two EOFs of the 5-month running means. Moreover, the strong correlation found between the simulated precipitation and the simulated winds in the first two principal components indicates the primary role of model precipitation in driving the surface winds. The surface winds simulated by a linear model forced by the GCM-simulated precipitation show good resemblance to the GCM-simulated winds in the equatorial region. This result supports the recent findings that the large-scale part of the tropical surface winds is primarily linear.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

The effect of grain size on the high-strain, high-strain-rate behavior of copper

Copper with four widely differing grain sizes was subjected to high-strain-rate plastic deformation in a special experimental arrangement in which high shear strains of approximately 2 to 7 were generated. The adiabatic plastic deformation produced temperature rises in excess of 300 K, creating conditions favorable for dynamic recrystallization, with an attendant change in the mechanical response. Preshocking of the specimens to an amplitude of 50 GPa generated a high dislocation density; twinning was highly dependent on grain size, being profuse for the 117- and 315-mu m grain-size specimens and virtually absent for the 9.5-mu m grain-size specimens. This has a profound effect on the subsequent mechanical response of the specimens, with the smaller grain-size material undergoing considerably more hardening than the larger grain-size material. A rationale is proposed which leads to a prediction of the shock threshold stress for twinning as a function of grain size. The strain required for localization of plastic deformation was dependent on the combined grain size/shock-induced microstructure, with the large grain-size specimens localizing more readily. The experimental results obtained are rationalized in terms of dynamic recrystallization, and a constitutive equation is applied to the experimental results; it correctly predicts the earlier onset of localization for the large grain-size specimens. It is suggested that the grain-size dependence of shock response can significantly affect the performance of shaped charges.

Dynamic Behavior of Poly(2-vinylpyridine) in Solution Studied by 13C Nuclear Magnetic Relaxation

Abstract: The dynamics of poly(2-vinylpyridine) in chloroform solution has been examined by C-13 spin-lattice relaxation time and NOE measurements as a function of temperature. The experiments were performed at 50.3 and 100.6 MHz. The backbone carbon relaxation data have been analyzed in terms of six motional models. Among these models, the models which consider conformational transitions and bond librations for the backbone were found to be more successful. Pyridyl ring motion has been modeled as a restricted rotation with the rotational amplitude varying with temperature. The activation energy parameters obtained from the relaxation data of the pyridyl ring carbon have been compared with the energy barrier for ring rotation estimated from conformational energy calculations using the AM1 semiempirical quantum chemical method. The results of the conformational energy calculations support the description of pyridyl ring motion as a restricted rotation.

Lopsidedness in WHISP galaxies II. Morphological lopsidedness

The distribution of stars and gas in many galaxies is asymmetric. This so-called lopsidedness is expected to significantly affect the dynamics and evolution of the disc, including the star formation activity. Here, we measure the degree of lopsidedness for the gas distribution in a selected sample of 70 galaxies from the Westerbork Hi Survey of Spiral and Irregular Galaxies. This complements our earlier work (Paper I) where the kinematic lopsidedness was derived for the same galaxies. The morphological lopsidedness is measured by performing a harmonic decomposition of the surface density maps. The amplitude of lopsidedness A(1), the fractional value of the first Fourier component, is typically quite high (about 0.1) within the optical disc and has a constant phase. Thus, lopsidedness is a common feature in galaxies and indicates a global mode. We measure A(1) out to typically one to four optical radii, sometimes even further. This is, on average, four times larger than the distance to which lopsidedness was measured in the past using near-IR as a tracer of the old stellar component, and therefore provides a new, more stringent constraint on the mechanism for the origin of lopsidedness. Interestingly, the value of A(1) saturates beyond the optical radius. Furthermore, the plot of A(1) versus radius shows fluctuations that we argue are due to local spiral features. We also try to explain the physical origin of this observed disc lopsidedness. No clear trend is found when the degree of lopsidedness is compared to a measure of the isolation or interaction probability of the sample galaxies. However, this does not rule out a tidal origin if the lopsidedness is long-lived. In addition, we find that the early-type galaxies tend to be more morphologically lopsided than the late-type galaxies. Both results together indicate that lopsidedness has a tidal origin.

Solvation dynamics in monohydroxy alcohols: Agreement between theory and different experiments

Recently three different experimental studies on ultrafast solvation dynamics in monohydroxy straight-chain alcohols (C-1-C-4) have been carried out, with an aim to quantify the time constant (and the amplitude) of the ultrafast component. The results reported are, however, rather different from different experiments. In order to understand the reason for these differences, we have carried out a detailed theoretical study to investigate the time dependent progress of solvation of both an ionic and a dipolar solute probe in these alcohols. For methanol, the agreement between the theoretical predictions and the experimental results [Bingemann and Ernsting J. Chem. Phys. 1995, 102, 2691 and Horng et al. J: Phys, Chern, 1995, 99, 17311] is excellent. For ethanol, propanol, and butanol, we find no ultrafast component of the time constant of 70 fs or so. For these three liquids, the theoretical results are in almost complete agreement with the experimental results of Horng et al. For ethanol and propanol, the theoretical prediction for ionic solvation is not significantly different from that of dipolar solvation. Thus, the theory suggests that the experiments of Bingemann and Ernsting and those of Horng et al. studied essentially the polar solvation dynamics. The theoretical studies also suggest that the experimental investigations of Joo et al. which report a much faster and larger ultrafast component in the same series of solvents (J. Chem. Phys. 1996, 104, 6089) might have been more sensitive to the nonpolar part of solvation dynamics than the polar part. In addition, a discussion on the validity of the present theoretical approach is presented. In this theory the ultrafast component arises from almost frictionless inertial motion of the individual solvent molecules in the force field of its neighbors.

Spin and charge density waves in the extended Hubbard model: A slave-boson approach

We use the extended Hubbard model to investigate the properties of the charge- and spin-density-wave phases in the presence of a nearest-neighbors repulsion term in the framework of the slave-boson technique. We show that, contrary to Hartree-Fock results, an instablity may occur for sufficiently high values of the Hubbard repulsion, both in the spin- and charge-density-wave phase, which makes the system discontinuously jump to a phase with a smaller or zero wave amplitude. The limits of applicability of our approach are discussed and our results are compared with previous numerical analysis. The phase diagram of the model at half-filling is determined.

Isomerization dynamics in viscous liquids: Microscopic investigation of the coupling and decoupling of the rate to and from solvent viscosity and dependence on the intermolecular potential

A detailed investigation of viscosity dependence of the isomerization rate is carried out for continuous potentials by using a fully microscopic, self-consistent mode-coupling theory calculation of both the friction on the reactant and the viscosity of the medium. In this calculation we avoid approximating the short time response by the Enskog limit, which overestimates the friction at high frequencies. The isomerization rate is obtained by using the Grote-Hynes formula. The viscosity dependence of the rate has been investigated for a large number of thermodynamic state points. Since the activated barrier crossing dynamics probes the high-frequency frictional response of the liquid, the barrier crossing rate is found to be sensitive to the nature of the reactant-solvent interaction potential. When the solute-solvent interaction is modeled by a 6-12 Lennard-Jones potential, we find that over a large variation of viscosity (eta), the rate (k) can indeed be fitted very well to a fractional viscosity dependence: (k similar to eta(-alpha)), with the exponent alpha in the range 1 greater than or equal to alpha >0. The calculated values of the exponent appear to be in very good agreement with many experimental results. In particular, the theory, for the first time, explains the experimentally observed high value of alpha even at the barrier frequency, omega(b). similar or equal to 9 X 10(12) s(-1) for the isomerization reaction of 2-(2'-propenyl)anthracene in liquid eta-alkanes. The present study can also explain the reason for the very low value of vb observed in another study for the isomerization reaction of trans-stilbene in liquid n-alkanes. For omega(b) greater than or equal to 2.0 X 10(13) s(-1), we obtain alpha similar or equal to 0, which implies that the barrier crossing rate becomes identical to the transition-state theory predictions. A careful analysis of isomerization reaction dynamics involving large amplitude motion suggests that the barrier crossing dynamics itself may become irrelevant in highly viscous liquids and the rate might again be coupled directly to the viscosity. This crossover is predicted to be strongly temperature dependent and could be studied by changing the solvent viscosity by the application of pressure. (C) 1999 American Institute of Physics. [S0021-9606(9950514-X].

Nonequilibrium stationary states of a particle in a gravitational field driven by a vibrating surface

Stationary velocity distribution functions are determined for a particle in a gravitational field driven by a vibrating surface in the limit of small dissipation. It is found that the form of the distribution function is sensitive to the mechanism of energy dissipation, inelastic collisions or viscous drag, and also to the form of the amplitude function of the vibrating surface. The velocity distributions obtained analytically are found to be in excellent agreement with the results of computer simulations in the limit of low dissipation. [S0031-9007(99)08898-5].

Rossby waves in May and the Indian summer monsoon rainfall

Large amplitude stationary Rossby wave trains with wavelength in the range 50 degrees to 60 degrees longitude have been identified in the upper troposphere during May, through the analysis of 200 hPa wind anomalies. The spatial phase of these waves has been shown to differ by about 20 degrees of longitude between the dry and wet Indian monsoon years. It has been shown empirically that the Rossby waves are induced by the heat sources in the ITCZ. These heat sources appear in the Bay of Bengal and adjoining regions in May just prior to the onset of the Indian summer monsoon. The inter-annual spatial phase shift of the Rossby waves has been shown to be related to the shift in the deep convection in the zonal direction.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

A checking method for probabilistic seismic-hazard assessment: case studies on three cities

The conventional Cornell's source-based approach of probabilistic seismic-hazard assessment (PSHA) has been employed all around the world, whilst many studies often rely on the use of computer packages such as FRISK (McGuire FRISK-a computer program for seismic risk analysis. Open-File Report 78-1007, United States Geological Survey, Department of Interior, Washington 1978) and SEISRISK III (Bender and Perkins SEISRISK III-a computer program for seismic hazard estimation, Bulletin 1772. United States Geological Survey, Department of Interior, Washington 1987). A ``black-box'' syndrome may be resulted if the user of the software does not have another simple and robust PSHA method that can be used to make comparisons. An alternative method for PSHA, namely direct amplitude-based (DAB) approach, has been developed as a heuristic and efficient method enabling users to undertake their own sanity checks on outputs from computer packages. This paper experiments the application of the DAB approach for three cities in China, Iran, and India, respectively, and compares with documented results computed by the source-based approach. Several insights regarding the procedure of conducting PSHA have also been obtained, which could be useful for future seismic-hazard studies.