A model of an exploding,radiating star in general relativity

In a previous paper Adams, Cary and Cohen (1994) presented a model of a supernova. In that paper the equations of General Relativity describing the evolution of a spherically symmetric, radiating star were solved analytically. The evolution of the star was determined by the application of boundary conditions at the center and at the edge. Due to lmitations in the presupernova model, only the very slow inward motion of an unstable, degenerate core could be considered. The solution was also limited by the need to exclude a runaway term, one that increased exponentially with time. Without the exclusion of the runaway, the luminosity would have increased without bound and the mass would have become negative.This paper presents a completely analytic solution to the equations of General Relativity describing the evolution of a Type II supernova. Professor S.E. Woosley kindly gave us data on the physical variables of a 12M ₀ presupernova star. In our model the core collapses within 1 s, leaving a 1.3M ₀ remnant. Shortly afterward 10.6M ₀ is ejected to infinity, and 0.17M ₀ is radiated away in the form of neutrinos. The distance of the edge from the center increases proportionally to the two-thirds power of the time. The luminosity decreases proportionally to the inverse four-thirds power.Although the runaway solution was modified by the exploding rather than a static envelope, it must still be excluded by adjusting initial conditions. Its character is changed from an exponential to a very large power (55) of time. The removal of a degree of freedom by this exclusion leads to physically non-sensical results such as negative luminosity. The inclusion of a term describing motion of the mantle due to neutrino interactions provides the additional degree of freedom necessary for physically reasonable results.     

An analysis of 1983 observations of facular contrast with an extreme limb photometer

Analysis of facular contrast I/I _qs from Extreme Limb Photometer (ELP) data of the summer of 1983 yield a mean contrast of 0.91 ± 0.19% and 1.57 ± 0.16% for apertures 1 and 2 located at = 0.198 and at = 0.111, respectively. The ratio of the mean contrast in the outer aperture (closer to the limb) to that of the inner one is 1.71 ± 0.40, indicating an increase in the mean facular contrast toward the limb. This result is in agreement with observations made in 1975, 1979, and 1982. The errors are dominated by the random presence of solar active regions. The combined results from all seasons follow an approximately ^–1 curve. Facular excess solar oblateness signals for 1983 are 33.8 ± 6.6 arc ms and 16.5 ± 2.1 arc ms for ELP apertures 1 and 2, respectively, in reasonable agreement with the 1983 excess solar oblateness results of Dicke, Kuhn, and Libbrecht (1985).     

Propagation of MHD body and surface waves in magnetically structured regions of the solar atmosphere

The fact that magnetically structured regions exist in the solar atmosphere has been known for a number of years. It has been suggested that different kinds of magnetohydrodynamic (MHD) waves can be efficiently damped in these regions and that the dissipated wave energy may be responsible for the observed enhancement in radiative losses. From a theoretical point of view, an important task would be to investigate the propagation and dissipation of MHD waves in these highly structured regions of the solar atmosphere. In this paper, we study the behavior of MHD body and surface waves in a medium with either a single or double (slab) magnetic interface by use of a nonlinear, two-dimensional, time-dependent, ideal MHD numerical model constructed on the basis of a Lagrangean grid and semi-implicit scheme. The processes of wave confinement and wave energy leakage are discussed in detail. It is shown that the obtained results depend strongly on the type of perturbations imposed on the interface or slab and on the plasma parameter, . The relevance of the obtained results to the heating problem of the upper parts of the solar atmosphere is also discussed.     

Propagation of solar cosmic rays: Diffusion versus focused diffusion

This paper is designed to bring to the attention the fact that the effect of focusing of solar energetic particles is always essential as compared with scattering, no matter how small the value of the mean free path may be. That is why, an ordinary (focusing-free) diffusion approach can not be applied to the solar cosmic ray transport. In the case of high-energy solar particles, the focused diffusion is demonstrated to lead to a power law decay of energetic particle intensity much like an ordinary diffusion. However, the power law index of the decay is renormalized by the focusing.     

Coronal magnetic fields from microwave polarization observations

The solar active region (AR) 7530 was observed at 6 cm on July 3 and 4, 1993 with the Westerbork Synthesis Radio Telescope, using a multi-channel receiver with very narrow bandwidth. We compare the radio data with Yohkoh SXT observations and with the magnetic field extrapolated from the Marshall vector magnetograms in the force-free and current-free approximations. The comparison with soft X-rays shows that, although a general agreement exists between the shape of the radio intensity map and the X-ray loops, the brightness temperature, T _b, obtained using the parameters derived from the SXT is much lower than that observed. The comparison with the extrapolated photospheric fields shows instead that they account very well for the observed T _b above the main sunspots, if gyroresonance emission is assumed. In the observation of July 4 an inversion and strong suppression of the circular polarization was clearly present above different portions of the AR, which indicates that particular relationships exist between the electron density and the magnetic field in the region where the corresponding lines of sight cross the field quasi-perpendicularly. The extrapolated magnetic field at a much higher level ( 10¹⁰ cm), satisfies the constraints required by the wave propagation theory all over the AR. However, a rather low electron density is derived.     

A wind tunnel study of turbulent flow around single and multiple windbreaks, part I: Velocity fields

This paper describes wind-tunnel experiments on the flow around single and multiple porous windbreaks (height H), sheltering a model plant canopy (height H/3). The mean wind is normal to the windbreaks, which span the width of the wind tunnel. The incident turbulent flow simulates the adiabatic atmospheric surface layer. Five configurations are examined: single breaks of three solidities (low, medium, high; solidity = 1 - porosity), and medium-solidity multiple breaks of streamwise spacing 12H and 6H. The experimental emphases are on the interactions of the windbreak flow with the underlying plant canopy; the effects of solidity; the differences in shelter between single and multiple windbreaks; and the scaling properties of the flow. Principal results are: (1) the "quiet zones" behind each windbreak are smaller in multiple than single arrays, because of the higher turbulence level in the very rough-wall internal boundary layer which develops over the multiple arrays. Nevertheless, the overall shelter effectiveness is higher for multiple arrays than single windbreaks because of the "nonlocal shelter" induced by the array as a whole. (2) The flow approaching the windbreak decelerates above the canopy but accelerates within the canopy, particularly when the windbreak solidity is high. (3) A strong mixing layer forms just downwind of the top of each windbreak, showing some of the turbulence and scaling properties of the classical mixing layer formed between uniform, coflowing streams. (4) No dramatic increase in turbulence levels in the canopy is evident at the point where the deepening mixing layer contacts the canopy (around x/H = 3) but the characteristic inflection in the canopy wind profile is eliminated at this point.     

Conditional statistics of vertical heat fluxes in local advection conditions

Turbulent fluctuations have been investigated in the internal boundary layer (IBL) which forms after a dry-to-wet surface transition. The IBL is defined as that part of the atmospheric surface layer where the influence of the downstream surface is noticeable. The results of the application of three different quadrant analysis techniques are presented. The three techniques, in increasing order of the amount of information supplied, provide:

1. the diurnal variation of quadrant contribution (C _i), number fraction (T _i) and conditional average (% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% GaeyykJeUabm4DayaafaGabm4CayaafaGaeyOkJe-aaSbaaSqaaiaa% dMgaaeqaaaaa!4215!\[\langle w's'\rangle _i \], with s = T or q) of vertical sensible and latent heat fluxes,
2. the quadrant contribution and number of samples of different sizes depending on the relative magnitude of each sample, and
3. the distribution of the nondimensional probability density function.

The results show that in the IBL the vertical flux of sensible heat is maintained by (i) a small fraction of large samples with warm air carried upwards, and (ii) a larger fraction of small samples with cool air carried downwards. Both processes are almost equal in importance. In the morning and near the top of the IBL negative temperature fluctuations are limited by the near-uniform temperature conditions upstream and above the IBL. This limitation reduces, at that location, the conditional average of the sinking motions of cool air. Closer to the wet surface the negative temperature fluctuations are less susceptible to the above mentioned limitation. As a consequence contributions from all four quadrants are almost equal leading to a very small vertical heat flux. In the presence of a temperature inversion over both the upstream and the downstream terrain, shear-generated turbulence appears to be the cause of the relative abundance of sinking motions of warm air and rising motions of cool air, leading to a reversal of the sensible heat flux. The latent heat flux is positive (i.e. directed away from the surface) at all times and is maintained in almost equal amount by (i) a small number of large magnitude samples with moist air carried upwards, and (ii) small magnitude samples with sinking motions of dry air. These sinking motions of dry air are far more numerous, especially in the morning, but their conditional average is very small. The abundance of sinking motions of dry air is attributed to the fact that over the downstream terrain evaporation is greatly enhanced, leading to a skewed w′q′ signal. This skewness is clearly visible in the w′q′-probability density distribution of the morning runs. In the evening the asymmetry between these two different contributions disappears. This is because evaporation is greatly reduced and large positive humidity fluctuations no longer occur.     

Simulation of large-scale gas structures formed in the interaction between galaxies. I. Method and preliminary results

Computer programs developed to study large-scale, transient gas structures in galaxies are described and test results are given. Gas-dynamic quantities are determined on the basis of a three-dimensional algorithm using so-called "smoothed particle hydrodynamics" (SPH). Preliminary calculations were made to simulate the formation of a gas ring around a spheroidal galaxy when it absorbs a low-mass, gas-rich companion, as well due to gas accretion during the flyby of a spiral galaxy of comparable mass. The evolution of tidal gas tails of disk galaxies is investigated.Translated from Astrofizika, Vol. 39, No. 2, pp. 265–284, April–June, 1996.     

On solving Kepler's equation for nearly parabolic orbits

We deal here with the efficient starting points for Kepler's equation in the special case of nearly parabolic orbits. Our approach provides with very simple formulas that allow calculating these points on a scientific vest-pocket calculator. Moreover, srtarting with these points in the Newton's method we can calculate a root of Kepler's equation with an accuracy greater than 0.001 in 0–2 iterations. This accuracy holds for the true anomaly || 135° and |e – 1| 0.01. We explain the reason for this effect also.Dedicated to the memory of Professor G.N. Duboshin (1903–1986).     

A filament eruption and accompanying coronal field changes on November 5, 1992

An H filament eruption on November 5, 1992 was fully observed in H with the Hida Flare Monitoring Telescope, while Yohkoh's Soft X-ray Telescope observed the pre- and post-eruption evolution of the coronal magnetic fields. From the H data, including the red and blue wings, we have reconstructed the rise of the filament, including trajectory, velocity, and acceleration. In combination with the Yohkoh data this reconstruction suggests that the filament had several interactions with other coronal magnetic fields during the eruption. The Yohkoh data also shows pre-eruption changes in the coronal fields and several post-eruption bright coronal structures. The pre-eruption changes are interpreted as a partial opening of the corona, indicating that it is not necessary to have a complete opening of the corona in order for a filament to erupt and we discuss the several possible contributions from emerging flux. The post-event bright coronal structures are compared with theory and with a cleaner filament eruption event on July 31, 1992. These comparisons suggest that, although there are many similarities, it is hard to completely reconcile the observations with the existing theory.