Inclusion of turbulence in solar modeling

The general consensus is that in order to reproduce the observed solar p-mode oscillation frequencies, turbulence should be included in solar models. However, until now there has not been any well-tested efficient method to incorporate turbulence into solar modeling. We present here two methods to include turbulence in solar modeling within the framework of the mixing length theory, using the turbulent velocity obtained from numerical simulations of the highly superadiabatic layer of the sun at three stages of its evolution. The first approach is to include the turbulent pressure alone, and the second is to include both the turbulent pressure and the turbulent kinetic energy. The latter is achieved by introducing two variables: the turbulent kinetic energy per unit mass, and the effective ratio of specific heats due to the turbulent perturbation. These are treated as additions to the standard thermodynamic coordinates (e.g. pressure and temperature). We investigate the effects of both treatments of turbulence on the structure variables, the adiabatic sound speed, the structure of the highly superadiabatic layer, and the p-mode frequencies. We find that the second method reproduces the SAL structure obtained in 3D simulations, and produces a p-mode frequency correction an order of magnitude better than the first method.Comment: 10 pages, 12 figure

Dewatering saturated, networked suspensions with a screw press

A model is presented for the dewatering of a saturated two-phase medium in a screw press. The model accounts for the detailed two-phase rheological behaviour of the pressed material and splits the press into two zones, an initial well-mixed constant-pressure region followed by an axial transport region in which the total pressure steadily increases. In this latter region, a slowly varying helical coordinate transformation is introduced to help reduce the dynamics to an annular bi-axial compression of the two-phase medium. Unlike previous modelling, the transition point between the two zones is determined self-consistently, rather than set a priori, and the pressure along the length of the press is deduced from the rheology of the two-phase flow rather than averaging the two-phase dynamics over a cross section of the press. The model is compared to experimental observations of the dewatering of a paper-making fibre suspension and of a clay slurry, and is shown to reproduce operational data

Probing Solar Convection

In the solar convection zone acoustic waves are scattered by turbulent sound speed fluctuations. In this paper the scattering of waves by convective cells is treated using Rytov's technique. Particular care is taken to include diffraction effects which are important especially for high-degree modes that are confined to the surface layers of the Sun. The scattering leads to damping of the waves and causes a phase shift. Damping manifests itself in the width of the spectral peak of p-mode eigenfrequencies. The contribution of scattering to the line widths is estimated and the sensitivity of the results on the assumed spectrum of the turbulence is studied. Finally the theoretical predictions are compared with recently measured line widths of high-degree modes.Comment: 26 pages, 7 figures, accepted by MNRA

Imbibition with swelling: Capillary rise in thin deformable porous media

The imbibition of a liquid into a thin deformable porous substrate driven by capillary suction is considered. The substrate is initially dry and has uniform porosity and thickness. Two-phase flow theory is used to describe how the liquid flows through the pore space behind the wetting front when out-of-plane deformation of the solid matrix is considered. Neglecting gravity and evaporation, standard shallow-layer scalings are used to construct a reduced model of the dynamics. The model predicts convergence to a self-similar behavior in all regions except near the wetting front, where a boundary layer arises whose structure narrows with the advance of the front. Over time, the rise height approaches the similarity scaling of t^1/2, as in the classical Washburn or BCLW law. The results are compared with a series of laboratory experiments using cellulose paper sheets, which provide qualitative agreement.Kruger Products, NSER

Obstructed viscoplastic flow in a Hele–Shaw cell

Experiments are conducted exploring the flow of Carbopol past obstacles in a narrow slot and compared with predictions of a model based on the Herschel-Bulkley constitutive law and the conventional Hele-Shaw approximation. Although Carbopol is often assumed to be a relatively simple yield-stress fluid, the flow pattern around an obstacle markedly lacks the fore-aft symmetry expected theoretically. Such asymmetry has been observed previously for viscoplastic flows past obstacles in unconfined geometries, but the narrowness of the Hele-Shaw cell ensures that the stress state is very different, placing further constraints on the underlying origin. The asymmetry is robust, as demonstrated by varying the shape and number of the obstacles, the surfaces of the cell walls, and the steadiness of the flow rate. The results suggest that rheological hysteresis near the yield point may be the cause of the asymmetry

Cardiac amyloidosis in non-transplant cardiac surgery

Cardiac amyloidosis is a rare infiltrative cardiomyopathy that portends a poor prognosis. There is a growing recognition of co-existent aortic valve stenosis and transthyretin cardiac amyloidosis, with some studies suggesting that dual pathology may be associated increased risk of complication and mortality during surgical intervention. This review aims to evaluate the available literature on non-transplant cardiac surgical interventions in patients with cardiac amyloidosis, with particular focus on diagnosis, high surgical risk and areas of uncertainty that require further research

Space and Ground Based Pulsation Data of Eta Bootis Explained with Stellar Models Including Turbulence

The space telescope MOST is now providing us with extremely accurate low frequency p-mode oscillation data for the star Eta Boo. We demonstrate in this paper that these data, when combined with ground based measurements of the high frequency p-mode spectrum, can be reproduced with stellar models that include the effects of turbulence in their outer layers. Without turbulence, the l=0 modes of our models deviate from either the ground based or the space data by about 1.5-4.0 micro Hz. This discrepancy can be completely removed by including turbulence in the models and we can exactly match 12 out of 13 MOST frequencies that we identified as l=0 modes in addition to 13 out of 21 ground based frequencies within their observational 2 sigma tolerances. The better agreement between model frequencies and observed ones depends for the most part on the turbulent kinetic energy which was taken from a 3D convection simulation for the Sun.Comment: 13 pages, 7 figures, ApJ in pres

Solar Oscillations and Convection: II. Excitation of Radial Oscillations

Solar p-mode oscillations are excited by the work of stochastic, non-adiabatic, pressure fluctuations on the compressive modes. We evaluate the expression for the radial mode excitation rate derived by Nordlund and Stein (Paper I) using numerical simulations of near surface solar convection. We first apply this expression to the three radial modes of the simulation and obtain good agreement between the predicted excitation rate and the actual mode damping rates as determined from their energies and the widths of their resolved spectral profiles. We then apply this expression for the mode excitation rate to the solar modes and obtain excellent agreement with the low l damping rates determined from GOLF data. Excitation occurs close to the surface, mainly in the intergranular lanes and near the boundaries of granules (where turbulence and radiative cooling are large). The non-adiabatic pressure fluctuations near the surface are produced by small instantaneous local imbalances between the divergence of the radiative and convective fluxes near the solar surface. Below the surface, the non-adiabatic pressure fluctuations are produced primarily by turbulent pressure fluctuations (Reynolds stresses). The frequency dependence of the mode excitation is due to effects of the mode structure and the pressure fluctuation spectrum. Excitation is small at low frequencies due to mode properties -- the mode compression decreases and the mode mass increases at low frequency. Excitation is small at high frequencies due to the pressure fluctuation spectrum -- pressure fluctuations become small at high frequencies because they are due to convection which is a long time scale phenomena compared to the dominant p-mode periods.Comment: Accepted for publication in ApJ (scheduled for Dec 10, 2000 issue). 17 pages, 27 figures, some with reduced resolution -- high resolution versions available at http://www.astro.ku.dk/~aake/astro-ph/0008048

Thermal structure of a gas-permeable lava dome and timescale separation in its response to perturbation

The thermal boundary layer at the surface of a volcanic lava dome is investigated through a continuum model of the thermodynamic advection diffusion processes resulting from magmatic gas flow through the dome matrix. The magmatic gas mass flux, porosity and permeability of the rock are identified as key parameters. New, theoretical, nonlinear steady-state thermal profiles are reported which give a realistic surface temperature of 210 degC for a region of lava dome surface through which a gas flux of 3.5 x 10-3 kg s-1 m-2 passes. This contrasts favourably with earlier purely diffusive thermal models, which cool too quickly. Results are presented for time-dependent perturbations of the steady states as a response to: changes in surface pressure, a sudden rockfall from the lava dome surface, and a change in the magmatic gas mass flux at depth. Together with a generalized analysis using the method of multiple scales, this identifies two characteristic time scales associated with the thermal evolution of a dome carapace: a short time scale of several minutes, over which the magmatic gas mass flux, density, and pressure change to a new quasi-steady-state, and a longer time scale of several days, over which the thermal profile changes to a new equilibrium distribution. Over the longer time scale the dynamic properties of the dome continue to evolve, but only in slavish response to the ongoing temperature evolution. In the light of this time scale separation, the use of surface temperature measurements to infer changes in the magmatic gas flux for use in volcanic hazard prediction is discussed