Effects of ions on whistler-mode ray tracing.
Smith subprotonospheric whistlers and digital computation of ion effects on whistler-mode ray tracing
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Smith subprotonospheric whistlers and digital computation of ion effects on whistler-mode ray tracing
Ray-tracing (RT) of Lunar Digital Elevation Models (DEM)'s is performed to virtually derive the degree of radiation incident to terrain as a function of time, orbital and ephemeris constraints [I- 4]. This process is an integral modeling process in lunar polar research and exploration due to the present paucity of terrain information at the poles and mission planning activities for the anticipated spring 2009 launch of the Lunar Reconnaissance Orbiter (LRO). As part of the Lunar Exploration Neutron Detector (LEND) and Lunar Crater Observation and Sensing Satellite (LCROSS) preparations RI methods are used to estimate the critical conditions presented by the combined effects of high latitude, terrain and the moons low obliquity [5-7]. These factors yield low incident solar illumination and subsequently extreme thermal, and radiation conditions. The presented research uses RT methods both for radiation transport modeling in space and regolith related research as well as to derive permanently shadowed regions (PSR)'s in high latitude topographic minima, e.g craters. These regions are of scientific and human exploration interest due to the near constant low temperatures in PSRs, inferred to be < 100 K. Hydrogen is thought to have accumulated in PSR's through the combined effects of periodic cometary bombardment and/or solar wind processes, and the extreme cold which minimizes hydrogen sublimation [8-9]. RT methods are also of use in surface position optimization for future illumination dependent on surface resources e.g. power and communications equipment.
A ray-trace simulation code for the Hard X-ray Telescope (HXT) on board the Hitomi (ASTRO-H) satellite is beingdeveloped. The half power diameter and effective area simulated based on the code are consistent with groundmeasurements within 10%. The HXT observed the pulsar wind nebula G21.5-0.9 for 105 ksec. We confirmed that theencircled energy function and the half power diameter obtained from the data are consistent with the ground measurements.
The Root-Sum-Squared, or “RSS” wavefront error model is a simple, scalar tool, commonly used for space telescope error budgeting. At the same time, much more detailed models, combining ray-trace and Fourier optics with optical alignments and wavefront controls, can provide accurate, high -resolution simulations for detailed system and subsystem design. This paper makes a connection between the two modeling approaches by deriving RSS model coefficients from ray-trace models, including the effects of wavefront controls, for computing system performance from component error statistics. It is shown that, properly constructed, the simple RSS error budget is a covariance analysis, and can be as accurate as high-resolution wavefront models for statistical wavefront error prediction. A notional segmented-aperture space telescope is used to illustrate this error modeling process.
The Voyager 1 and 2 Planetary Radio Astronomy observations of Io-dependent decametric (DAM) radiation originating from the Southern Hemisphere of Jupiter were compared with the results of three-dimensional model ray tracing calculations of the DAM radiation. The ray trajectories for sources located at constant sub-Io longitudes of 260 and 300 deg were computed for both the Northern and the Southern Jovian Hemisphere sources. The model results of wave propagation agree with the Voyager observations obtained with Io located at 260 and 300 deg in Jovian system III longitude. The agreement between the Voyager observations and the model ray tracings allows identification of the origin of several of the emission components.
Analysis of energy emitted from simple or complex cavity designs can lead to intricate solutions due to nonuniform radiosity and irradiation within a cavity. A numerical ray tracing technique was applied to simulate radiation propagating within and from various cavity designs. To obtain the energy balance relationships between isothermal and nonisothermal cavity surfaces and space, the computer code NEVADA was utilized for its statistical technique applied to numerical ray tracing. The analysis method was validated by comparing results with known theoretical and limiting solutions, and the electrical resistance network method. In general, for nonisothermal cavities the performance (apparent emissivity) is a function of cylinder length-to-diameter ratio, surface emissivity, and cylinder surface temperatures. The extent of nonisothermal conditions in a cylindrical cavity significantly affects the overall cavity performance. Results are presented over a wide range of parametric variables for use as a possible design reference.
An MHD simulation was performed to obtain a self-consistent model of magnetic field and plasma density near the X point reconnection region. The MHD model was used to perform extensive ray tracing calculations in order to clarify the propagation characteristics of the plasma waves near the X point reconnection region. The dynamic wave spectra possibly observed by the Geotail spacecraft during a typical cross-tail trajectory are reconstructed. By comparing the extensive ray tracing calculations with the plasma wave data from Geotail, it is possible to perform a kind of 'remote sensing' to identify the location and structure of potential X point reconnection regions.
This paper provides an overview of the computational techniques used to characterize the viability of different lunar architectures and their ability to provide communication services to the lunar surface. This analysis was done with modern ray tracing techniques that allow for the computations to be done on Graphics Processing Unit (GPU) clusters for a high level of parallelism and severe reduction in computation time. The ray tracing computations were done with the GPU platform Compute Unified Device Architecture (CUDA) provided by NVIDIA which utilizes general-purpose computing on graphics processing units (GPGPU). This new method provides the advantage of being able to characterize a much larger portion of the lunar surface due to its computational efficiency as well as providing a more accurate representation of elevation angle limits instead of the typical and often inaccurate elevation angle mask. The Lunar surface can now be characterized with metrics such as contact time, outage time, and received data rate. With these metrics, different proposed Lunar architectures can be rapidly evaluated. This reduction in computation time not only leads to more accurate results but allows these results to be obtained in a time frame that allows for the complete characterization of the trade space. It is expected that these different architecture comparisons will lead to a conclusive determination of the optimal Lunar architecture and will allow for future Lunar missions to operate as close to real time as possible. In addition, this computation method can be used to recreate visibility figures generated by previous methods but with an increased level of accuracy.
A computationally efficient 3DHZETRN code with enhanced neutron and light ion (Z is less than or equal to 2) propagation was recently developed for complex, inhomogeneous shield geometry described by combinatorial objects. Comparisons were made between 3DHZETRN results and Monte Carlo (MC) simulations at locations within the combinatorial geometry, and it was shown that 3DHZETRN agrees with the MC codes to the extent they agree with each other. In the present report, the 3DHZETRN code is extended to enable analysis in ray-trace geometry. This latest extension enables the code to be used within current engineering design practices utilizing fully detailed vehicle and habitat geometries. Through convergence testing, it is shown that fidelity in an actual shield geometry can be maintained in the discrete ray-trace description by systematically increasing the number of discrete rays used. It is also shown that this fidelity is carried into transport procedures and resulting exposure quantities without sacrificing computational efficiency.
Global Positioning System (GPS) Radio Occultations (RO) bending angles and refractivity data characterize mostly the vertical structure of the Earth's atmosphere. We answer the question whether proper simulation of GPS RO data for data assimilation can be obtained with one-dimensional vertical operators, or if accounting also for horizontal atmospheric structures via ray-tracing makes a positive difference when compared with real data. We present a detailed implementation of a geometrical optics multi-plane two-dimensional (2D) ray-tracing as an observation operator to simulate GPS RO bending angles and refractivities within the Finite Volume Data Assimilation System (FVDAS). Comparisons of the outputs of that 2D observation operator with those of simpler ID observation operators are used to generate estimates of errors induced by neglecting tangent point drift (TPD) and horizontal gradients (HG). These error estimates are then confronted with errors estimates derived using 6335 real CHAMP and SAC-C occultations. The agreement for TPD-induced (HG-induced) errors is remarkably positive at altitudes 10-30 km (below 10 km). Comparisons in bending angles O - B STD of the outputs of the multi-plane 2D ray-tracer with those of a vertical Abel transform show reductions of about 8% of the usual O - B bending angle STD due to TPD in the stratosphere (3% due to HG, in the troposphere only). In terms of refractivity, the O - B STD reductions are about 1520% for TPD and 3-5% for HG in the same regions. These reductions are obtained using either 6-hour forecasts or analyses as backgrounds, and using Geometrical Optics (GO) or Canonical Transform (CT) data.
Light Detection And Ranging (LiDAR) provides unique data on the 3-D structure of atmosphere constituents and the Earth's surface. Simulating LiDAR returns for different laser technologies and Earth scenes is fundamental for evaluating and interpreting signal and noise in LiDAR data. Different types of models are capable of simulating LiDAR waveforms of Earth surfaces. Semi-empirical and geometric models can be imprecise because they rely on simplified simulations of Earth surfaces and light interaction mechanisms. On the other hand, Monte Carlo ray tracing (MCRT) models are potentially accurate but require long computational time. Here, we present a new LiDAR waveform simulation tool that is based on the introduction of a quasi-Monte Carlo ray tracing approach in the Discrete Anisotropic Radiative Transfer (DART) model. Two new approaches, the so-called "box method" and "Ray Carlo method", are implemented to provide robust and accurate simulations of LiDAR waveforms for any landscape, atmosphere and LiDAR sensor configuration (view direction, footprint size, pulse characteristics, etc.). The box method accelerates the selection of the scattering direction of a photon in the presence of scatterers with non-invertible phase function. The Ray Carlo method brings traditional ray-tracking into MCRT simulation, which makes computational time independent of LiDAR field of view (FOV) and reception solid angle. Both methods are fast enough for simulating multi-pulse acquisition. Sensitivity studies with various landscapes and atmosphere constituents are presented, and the simulated LiDAR signals compare favorably with their associated reflectance images and Laser Vegetation Imaging Sensor (LVIS) waveforms. The LiDAR module is fully integrated into DART, enabling more detailed simulations of LiDAR sensitivity to specific scene elements (e.g., atmospheric aerosols, leaf area, branches, or topography) and sensor configuration for airborne or satellite LiDAR sensors.
Ray tracing digital computer program for magnetospheric duct propagation studies
The properties of a pencil of light as defined approximately in the geometric optics ray tracing method are investigated. The vector Kirchhoff integral is utilized to accurately compute the electromagnetic near field in and around the pencil of light with various beam base sizes, shapes, propagation directions and medium refractive indices. If a pencil of light has geometric mean cross section size of the order p times the wavelength, it can propagate independently to a distance p2 times the wavelength, where most of the beam energy diffuses out of the beam region. This is consistent with a statement that van de Hulst made in a classical text on light scattering. The electromagnetic near fields in the pencil of light are not uniform, have complicated patterns within short distances from the beam base, and the fields tend to converge to Fraunhofer diffraction fields far away from the base.
A 'shooting and bouncing ray' (SBR) formulation is presented for treating the electromagnetic scattering from electrically large, inhomogeneous objects. A dense grid of rays representing the incident plane wave is shot toward the inhomogeneous object. At the scatterer boundary, reflected rays and refracted rays are generated due to discontinuity of the medium parameters. The trajectory, amplitude, phase and polarization of the rays inside the inhomogeneous object are traced based on geometrical optics. Whenever the rays cross the scatterer surface, additional reflected/refracted rays are generated and are tracked. This process is repeated until the intensities of the refracted/reflected rays become negligible. The contributions of the exiting rays to the total scattered field are calculated by using the equivalence principle in conjunction with a ray-tube integration scheme. The ray formulation is applied to calculate the backscattering from cylinders and spheres and good agreement with the exact series solutions is observed in the high frequency range. In addition, the backscattering mechanisms in penetrable objects are interpreted in terms of simple ray pictures.
Assessing and optimizing polarization performance in the context of ray-based optical design can be challenging. We describe an approach to this problem that decouples polarization effects from optical system geometry for reflective systems. Each surface's polarization properties are parameterized in terms of their impact on retardance and diattenuation in the small angle-of-incidence limit, separating polarization assessment from the task of coating design. A low-resolution ray trace of the system is adequate to determine ray geometry at each interface, which can then be interpolated to rapidly evaluate net Jones Matrix pupil functions. Coating behavior can be easily varied using the ellipsometric parameters to investigate impacts and compensation. Desired values of these parameters can then be specified as constraints in coating design. Investigation with candidate telescope optical designs for LUVOIR show baseline root-mean-square wavefront errors in the nm range for the on-diagonal Jones matrix component, and throughputs of tens of parts per million. Promising possibilities for compensation using a purpose-designed coating on the secondary mirror are discussed, which reduce the on-diagonal wavefront error by a factor 20, with accompanying but more modest reductions in coupling into off-diagonal terms.
Line-of-sight deviations and ray trace analyses of Apollo side window for optical experiments
High frequency radio wave propagation research using Alouette 2 topside sounder data and digital ray tracing techniques
The ray tracing computer program described by Menietti and Lin (1985) was used to investigate the propagation of Z mode emissions from sources in the auroral zone. The calculation results indicate that the propagation is primarily perpendicular to the magnetic field line (horizontal propagation) for almost all sources considered, except when the wave source is assumed to be unrestricted by the cyclotron maser mechanism. While the propagation is initially downward, refraction rapidly increases the wave normal angle, and only propagation at near-constant altitude occurs. The frequency bandwidth estimated from the ray paths was found to be less than that observed by the DE 1 satellite. Large horizontal magnetic field gradients associated with field-aligned currents produce downward propagation but are still not sufficient to produce the observed bandwidths. The results suggest that wave growth may be due to mechanisms in addition to cyclotron maser resonance.