A ray tracing digital computer program for the study of magnetospheric duct propagation
Ray tracing digital computer program for magnetospheric duct propagation studies
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.
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
Ocean surface gravity waves are an important component of air-sea interaction, influencing energy, momentum, and gas exchanges across the ocean-atmosphere interface. In specific applications such as refraction by ocean currents or bathymetry, ray tracing provides a computationally efficient way to gain insight into wave propagation. In this paper, we introduce Mantaray, an open-source software package implemented in Rust, with a Python interface, that solves the ray equations for ocean surface gravity waves. Mantaray is designed for performance, robustness, and ease of use. The package is modular to facilitate further development and can currently be applied to both idealized and realistic wave propagation problems (Fig. 1).
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.
We introduce Mahakala, a Python-based, modular, radiative ray-tracing code for curved spacetimes. We employ Google's JAX framework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr spacetimes. JAX also enables Mahakala to run in parallel on both CPUs and GPUs. Mahakala natively uses the Cartesian Kerr–Schild coordinate system, which avoids numerical issues caused by the pole in spherical coordinate systems. We demonstrate Mahakala's capabilities by simulating 1.3 mm wavelength images (the wavelength of Event Horizon Telescope observations) of general relativistic magnetohydrodynamic simulations of low-accretion rate supermassive black holes. The modular nature of Mahakala allows us to quantitatively explore how different regions of the flow influence different image features. We show that most of the emission seen in 1.3 mm images originates close to the black hole and peaks near the photon orbit. We also quantify the relative contribution of the disk, forward jet, and counterjet to 1.3 mm images.
Nonvertical propagation effects on high altitude topside ionosphere sounder data reduction, using ray tracing
Two-dimensional ray tracing of Z-mode radiation from sources in the auroral zone has been performed. The Z-mode waves are assumed to be excited by the cyclotron maser mechanism at frequencies near the electron gyrofrequency and at wave normal angles near 90 deg. The emission is found to propagate primarily perpendicular to the magnetic field line (horizontal propagation) for all sources considered, and thus not to substantially lower altitudes. The frequency bandwidth estimated from the ray paths is much less than that observed by DE-1. It is suggested that the wave source needs to be in the region of field-aligned current in order for the cyclotron maser mechanism to explain the broad frequency bandwidth of Z-mode emissions.
Using a model of Jupiter's magnetized plasma environment, radiation raypaths were calculated with a three-dimension ray tracing program. It is assumed that energetic particles produce the emission in the planet's auroral zone at frequencies just above the electron gyrofrequencies. This radiation is generated in narrow sheets defined by the angle of a ray with respect to the magnetic field line. By specifying the source position: latitude, longitude, and radial distance from the planet, signatures in the spectrum of frequency versus time seen by Voyager 1 and 2 were duplicated. The frequency range and the curvature of the decametric arcs in these dynamic spectra are the result of the geometry of the radiation sheets (imposed by the plasma and by the B-field) and illumination of Voyager 1 and 2 as the rotating magnetosphere mimics a pulsar.
Laser melting is a core component of the ongoing industrial revolution, dubbed Industry 4.0, as lasers facilitate fast and precise melting and fusion in advanced manufacturing. There is a strong need to optimize the laser process using simulations. However, this has proven challenging as high fidelity simulations are needed for predictive modeling and this is currently prohibitively expensive even when run on hundreds of processors on high performance computers. The challenge is capturing complex physics of laser material interaction, fluid dynamics, thermal physics and material phase transformations at various length and time scales. To close this technological gap, we modified a squeeze U-net to accelerate the laser ray tracing component of such high fidelity models by ~4x–40x while preserving the core physics principle of conservation of energy with 97% accuracy. This approach enables the accurate modeling of global laser energy absorption as a function of local surface temperatures and complex surface topologies, which govern the reflection directions and energy losses of laser rays upon interacting with the material surface.
We use an analytical fit to an emission lobe profile together with three-dimensional ray tracing to model the broad-banded smooth Uranian kilometric radiation (UKR). We assume the radiation is gyroemission from sources along magnetic field lines. Using an iterative technique that modifies the lobe function and source region, the results are compared to observations at a frequency of 481 kHz. The best-fit calculations are compared to previously published models and to recent ultraviolet (UV) observations.
We present SHADOW4, a new version of the popular ray tracing code. The SHADOW kernel has been completely rewritten in Python applying modern concepts of software engineering. A new user interface is available in the OASYS ecosystem. The new tool has been designed and implemented preparing the future needs both in computing (cloud computing, AI integration) and in the transit to fourth generation sources and beyond.
Radio wave propagation predictions are described in which modern comprehensive theoretical ionospheric models are coupled with ray-tracing programs. In the computer code described, a network of electron density and collision frequency parameters along a band about the great circle path is calculated by specifying the transmitter and receiver geographic coordinates, time, the day number, and the 2800-MHz solar flux. The ray paths are calculated on specifying the frequency, mode, range of elevation angles, and range of azimuth angles from the great circle direction. The current program uses a combination of the Penn State MKI E and F region models and the Mitra-Rowe D and E region model. Application of the technique to the prediction of satellite to ground propagation and calculation of oblique incidence propagation paths and absorption are described. The implications of the study to the development of the next generation of ionospheric models are discussed.
Raypaths for decametric wavelength radiation in Jupiter's magnetosphere were calculated. The model-dependent raypaths with the Voyager observations were compared. Characteristics of the source regions and the influence of propagation effects were deduced. A three dimensional ray tracing program was employed to calculate the raypaths. Families of rays were launched at particular angles with respect to the magnetic field lines to generate conical sheets of radiation for various frequencies and various source locations. As the planet's magnetic field rotates, these warped sheets of radiation sweep past the observer, producing signatures in frequency versus time plots. These signatures match some of those found in the Voyager data. The greatest propagation effects occur in and around the source regions in the Io auroral oval.
Errors may be introduced in satellite laser ranging data by atmospheric refractivity. Ray tracing data have indicated that horizontal refractivity gradients may introduce nearly 3-cm rms error when satellites are near 10-degree elevation. A correction formula to compensate for the horizontal gradients has been developed. Its accuracy is evaluated by comparing it to refractivity profiles. It is found that if both spherical and gradient correction formulas are employed in conjunction with meteorological measurements, a range resolution of one cm or less is feasible for satellite elevation angles above 10 degrees.