Parallel Monte Carlo approach for integration of the rendering equation

This paper is addressed to the numerical solving of the rendering equation in realistic image creation. The rendering equation is integral equation describing the light propagation in a scene accordingly to a given illumination model. The used illumination model determines the kernel of the equation under consideration. Nowadays, widely used are the Monte Carlo methods for solving the rendering equation in order to create photorealistic images. In this work we consider the Monte Carlo solving of the rendering equation in the context of the parallel sampling scheme for hemisphere. Our aim is to apply this sampling scheme to stratified Monte Carlo integration method for parallel solving of the rendering equation. The domain for integration of the rendering equation is a hemisphere. We divide the hemispherical domain into a number of equal sub-domains of orthogonal spherical triangles. This domain partitioning allows to solve the rendering equation in parallel. It is known that the Neumann series represent the solution of the integral equation as a infinity sum of integrals. We approximate this sum with a desired truncation error (systematic error) receiving the fixed number of iteration. Then the rendering equation is solved iteratively using Monte Carlo approach. At each iteration we solve multi-dimensional integrals using uniform hemisphere partitioning scheme. An estimate of the rate of convergence is obtained using the stratified Monte Carlo method. This domain partitioning allows easy parallel realization and leads to convergence improvement of the Monte Carlo method. The high performance and Grid computing of the corresponding Monte Carlo scheme are discussed.

Parallel Monte Carlo sampling scheme for sphere and hemisphere

The sampling of certain solid angle is a fundamental operation in realistic image synthesis, where the rendering equation describing the light propagation in closed domains is solved. Monte Carlo methods for solving the rendering equation use sampling of the solid angle subtended by unit hemisphere or unit sphere in order to perform the numerical integration of the rendering equation. In this work we consider the problem for generation of uniformly distributed random samples over hemisphere and sphere. Our aim is to construct and study the parallel sampling scheme for hemisphere and sphere. First we apply the symmetry property for partitioning of hemisphere and sphere. The domain of solid angle subtended by a hemisphere is divided into a number of equal sub-domains. Each sub-domain represents solid angle subtended by orthogonal spherical triangle with fixed vertices and computable parameters. Then we introduce two new algorithms for sampling of orthogonal spherical triangles. Both algorithms are based on a transformation of the unit square. Similarly to the Arvo's algorithm for sampling of arbitrary spherical triangle the suggested algorithms accommodate the stratified sampling. We derive the necessary transformations for the algorithms. The first sampling algorithm generates a sample by mapping of the unit square onto orthogonal spherical triangle. The second algorithm directly compute the unit radius vector of a sampling point inside to the orthogonal spherical triangle. The sampling of total hemisphere and sphere is performed in parallel for all sub-domains simultaneously by using the symmetry property of partitioning. The applicability of the corresponding parallel sampling scheme for Monte Carlo and Quasi-D/lonte Carlo solving of rendering equation is discussed.

Forward modelling to determine the observational signatures of white-light imaging and interplanetary scintillation for the propagation of an interplanetary shock in the ecliptic plane

Recent coordinated observations of interplanetary scintillation (IPS) from the EISCAT, MERLIN, and STELab, and stereoscopic white-light imaging from the two heliospheric imagers (HIs) onboard the twin STEREO spacecraft are significant to continuously track the propagation and evolution of solar eruptions throughout interplanetary space. In order to obtain a better understanding of the observational signatures in these two remote-sensing techniques, the magnetohydrodynamics of the macro-scale interplanetary disturbance and the radio-wave scattering of the micro-scale electron-density fluctuation are coupled and investigated using a newly constructed multi-scale numerical model. This model is then applied to a case of an interplanetary shock propagation within the ecliptic plane. The shock could be nearly invisible to an HI, once entering the Thomson-scattering sphere of the HI. The asymmetry in the optical images between the western and eastern HIs suggests the shock propagation off the Sun–Earth line. Meanwhile, an IPS signal, strongly dependent on the local electron density, is insensitive to the density cavity far downstream of the shock front. When this cavity (or the shock nose) is cut through by an IPS ray-path, a single speed component at the flank (or the nose) of the shock can be recorded; when an IPS ray-path penetrates the sheath between the shock nose and this cavity, two speed components at the sheath and flank can be detected. Moreover, once a shock front touches an IPS ray-path, the derived position and speed at the irregularity source of this IPS signal, together with an assumption of a radial and constant propagation of the shock, can be used to estimate the later appearance of the shock front in the elongation of the HI field of view. The results of synthetic measurements from forward modelling are helpful in inferring the in-situ properties of coronal mass ejection from real observational data via an inverse approach.

Eastward propagation of a plasma convection enhancement following a southward turning of the interplanetary magnetic field

On October 27th 1984, high-latitude ionospheric convection was observed by the European incoherent scatter (EISCAT) radar. For a nine-hour period, simultaneous observations of the interplanetary magnetic field (IMF) were obtained sunward of the Earth's bow shock. During this period, the IMF abruptly turned southward, having previously been predominantly northward for approximately three hours, and a strong enhancement in convection was observed 11 ± 1 minutes later. Using the very high time resolution of the EISCAT data, it is shown that the convection enhancement propagated eastward, around the afternoon magnetic local time sector, at a speed of the order of 1 kms−1. These results are interpreted in terms of the effects of an onset of steady IMF-geomagnetic field merging and are the first to show how a new pattern of enhanced convection is established in the high latitude ionosphere.

A solar storm observed from the Sun to Venus using the STEREO, Venus Express, and MESSENGER spacecraft

The suite of SECCHI optical imaging instruments on the STEREO-A spacecraft is used to track a solar storm, consisting of several coronal mass ejections (CMEs) and other coronal loops, as it propagates from the Sun into the heliosphere during May 2007. The 3-D propagation path of the largest interplanetary CME (ICME) is determined from the observations made by the SECCHI Heliospheric Imager (HI) on STEREO-A (HI-1/2A). Two parts of the CME are tracked through the SECCHI images, a bright loop and a V-shaped feature located at the rear of the event. We show that these two structures could be the result of line-of-sight integration of the light scattered by electrons located on a single flux rope. In addition to being imaged by HI, the CME is observed simultaneously by the plasma and magnetic field experiments on the Venus Express and MESSENGER spacecraft. The imaged loop and V-shaped structure bound, as expected, the flux rope observed in situ. The SECCHI images reveal that the leading loop-like structure propagated faster than the V-shaped structure, and a decrease in in situ CME speed occurred during the passage of the flux rope.We interpret this as the result of the continuous radial expansion of the flux rope as it progressed outward through the interplanetary medium. An expansion speed in the radial direction of ~30 km s-1 is obtained directly from the SECCHI-HI images and is in agreement with the difference in speed of the two structures observed in situ. This paper shows that the flux rope location can be determined from white light images, which could have important space weather applications.

Changes in light intensity can influence age at sexual maturity in domestic pullets

1. Shaver White and ISA Brown pullets were reared to 140 d in groups of 8 in cages on a 10-h photoperiod of incandescent light and maintained at an illuminance of 3 or 25 lux, or transferred from 3 to 25 lux or from 25 to 3 lux at 63 or 112 d of age. 2. There was no significant difference in sexual maturity, measured as eggs per 100 bird.d at 139 and 140 d, for ISA Brown maintained on 3 or 25 lux, but Shaver White pullets exposed to constant 3 lux matured significantly later than those maintained on 25 lux. 3. In Shaver Whites, sexual maturity was significantly delayed by an increase from 3 to 25 lux at 63 and 112 d, and advanced by a decrease from 25 to 3 lux at 112 d. Sexual maturity of ISA Browns was not significantly affected by a change in illuminance at 63 or 112 d, though responses were in the same direction as for Shaver Whites. 4. In both breeds, total feed consumed to 112 d was higher for birds on 3 lux than 25 lux, but lower between 112 d and 140 d when birds on 25 lux underwent rapid sexual development. In both breeds, body weight at 63 d was higher for birds exposed to 3 lux than 25 lux, but body weight gain thereafter was similar for the two light intensities. 5. In both breeds, plasma luteinising hormone (LH) concentration at 63 and 112 d was lower in birds maintained on 3 lux than 25 lux. At 63 and 112 d, transfers from 25 to 3 lux depressed, whereas transfers from 3 to 25 lux at 63 d, but not at 112 d, increased plasma LH. 6. Advances or delays in sexual maturity induced by changes in illuminance were not correlated with differences in feed intake, body weight gain, or with changes in plasma LH. 7. One possible explanation for the inverse relationship between the direction of change in illuminance at 63 and 112 d in pullets exposed to a 10-h photoperiod and the age at which they became sexually mature is that changes in light intensity and/or spectral composition affect the entrainment of the circadian rhythm of photoinducibility, to effect a phase shift in the photoinducible phase and/or the responsiveness of phototransduction pathways.

The effects of photoperiod and light spectrum on stock plant growth and rooting of cuttings of Cotinus coggygria 'Royal Purple'

Extending the season of production and improving the scheduling of ornamental crops are key commercial objectives for nurserymen. In some woody species, the period in which cuttings can be rooted successfully is transient, thus limiting the opportunities for scheduled production. Optimum rooting often occurs in early- to mid-summer coinciding with periods of active shoot growth. The relationship between this shoot activity and root initiation was investigated in Cotinus coggygria 'Royal Purple'. Shoot growth on stock plants was manipulated by altering the photoperiod or light quality. Results indicated there were seasonal effects on rooting, but the importance of shoot activity varied with harvest time. Cuttings harvested in August had high rooting percentages, irrespective of photoperiod, and despite shoot growth terminating in response to the short-day treatment. In contrast, by September, rooting percentage was highest in cuttings from plants under long-days, which had maintained greatest shoot growth activity. Cotinus shoots grown in vitro under 16 h days showed reduced shoot growth and increased rooting competence compared with shoots grown under 8 h days. Growing stock plants under polythene films, which altered the amount and quality of the incident light, influenced the rooting of cuttings harvested in August, but no consistent relationship with shoot activity was apparent. From a practical viewpoint, maintaining shoot activity late in the season may prolong the period for propagation by cuttings; but, from a scientific viewpoint, processes associated with an active shoot apex do not provide a complete explanation of seasonal variation in rooting.

Whistler mode wave growth and propagation in the prenoon magnetosphere

Pitch-angle scattering of electrons can limit the stably trapped particle flux in the magnetosphere and precipitate energetic electrons into the ionosphere. Whistler-mode waves generated by a temperature anisotropy can mediate this pitch-angle scattering over a wide range of radial distances and latitudes, but in order to correctly predict the phase-space diffusion, it is important to characterise the whistler-mode wave distributions that result from the instability. We use previously-published observations of number density, pitch-angle anisotropy and phase space density to model the plasma in the quiet pre-noon magnetosphere (defined as periods when AE<100nT). We investigate the global propagation and growth of whistler-mode waves by studying millions of growing ray paths and demonstrate that the wave distribution at any one location is a superposition of many waves at different points along their trajectories and with different histories. We show that for observed electron plasma properties, very few raypaths undergo magnetospheric reflection, most rays grow and decay within 30 degrees of the magnetic equator. The frequency range of the wave distribution at large L can be adequately described by the solutions of the local dispersion relation, but the range of wavenormal angle is different. The wave distribution is asymmetric with respect to the wavenormal angle. The numerical results suggest that it is important to determine the variation of magnetospheric parameters as a function of latitude, as well as local time and L-shell.

Using coordinated observations in polarized white light and Faraday rotation to probe the spatial position and magnetic field of an interplanetary sheath

Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons and therefore depend on both viewing geometry and electron density. The Faraday rotation (FR) of radio waves from extragalactic pulsars and quasars, which arises due to the presence of such solar wind features, depends on the line-of-sight magnetic field component B ∥ and the electron density. To understand coordinated WL and FR observations of CMEs, we perform forward magnetohydrodynamic modeling of an Earth-directed shock and synthesize the signatures that would be remotely sensed at a number of widely distributed vantage points in the inner heliosphere. Removal of the background solar wind contribution reveals the shock-associated enhancements in WL and FR. While the efficiency of Thomson scattering depends on scattering angle, WL radiance I decreases with heliocentric distance r roughly according to the expression Ir –3. The sheath region downstream of the Earth-directed shock is well viewed from the L4 and L5 Lagrangian points, demonstrating the benefits of these points in terms of space weather forecasting. The spatial position of the main scattering site r sheath and the mass of plasma at that position M sheath can be inferred from the polarization of the shock-associated enhancement in WL radiance. From the FR measurements, the local B ∥sheath at r sheath can then be estimated. Simultaneous observations in polarized WL and FR can not only be used to detect CMEs, but also to diagnose their plasma and magnetic field properties.

Modelling the high-latitude ionosphere for time-varying plasma convection

The paper discusses how variations in the pattern of convective plasma flows should beincluded in self-consistent time-dependent models of the coupled ionosphere-thermosphere system. The author shows how these variations depend upon the mechanism by which the solar wind flow excites the convection. The modelling of these effects is not just of relevance to the polar ionosphere. This is because the influence of convection is not confined to high latitudes: the resultant heating and composition changes in the thermosphere are communicated to lower latitudes by the winds which are also greatly modified by the plasma convection. These thermospheric changes alter the global distribution of plasma by modulatingthe rates of the chemical reactions which areresponsible for the loss of plasma. Hence the modelling of these high-latitude processes is of relevanceto the design and operation of HF communication, radar and navigation systems worldwide.