Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userModeling in-ice radio propagation with parabolic equation methods
81
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We investigate the use of parabolic equation (PE) methods for solving radio-wave propagation in polar ice. PE methods provide an approximate solution to Maxwells equations, in contrast to full-field solutions such as finite-difference-time-domain (FDTD) methods, yet provide a more complete model of propagation than simple geometric ray-tracing (RT) methods that are the current state of the art for simulating in-ice radio detection of neutrino-induced cascades. PE are more computationally efficient than FDTD methods, and more flexible than RT methods, allowing for the inclusion of diffractive effects, and modeling of propagation in regions that cannot be modeled with geometric methods. We present a new PE approximation suited to the in-ice case. We conclude that current ray-tracing methods may be too simplistic in their treatment of ice properties, and their continued use could overestimate experimental sensitivity for in-ice neutrino detection experiments. We discuss the implications for current in-ice Askaryan-type detectors and for the upcoming Radar Echo Telescope; two families of experiments for which these results are most relevant. We suggest that PE methods be investigated further for in-ice radio applications.
rate research
Read More
A parabolic equation for the propagation of periodic internal waves over varying bottom topography is derived using the multiple-scale perturbation method. Some computational aspects of the numerical implementation are discussed. The results of numerical experiments on propagation of an incident plane wave over a circular-type shoal are presented in comparison with the analytical result, based on Born approximation.
The Askaryan Radio Array (ARA) experiment at the South Pole is designed to detect high-energy neutrinos which, via in-ice interactions, produce coherent radiation at frequencies up to 1000 MHz. In Dec. 2018, a custom high-amplitude radio-frequency transmitter was lowered into the 1700 m SPICE ice core to provide test sources for ARA receiver stations sensitive to vertical and horizontal polarizations. For these tests, signal geometries correspond to obliquely propagating radio waves from below. The ARA collaboration has recently measured the polarization-dependent time delay variation, and report more significant time delays for trajectories perpendicular to ice flow. Here we use fabric data from the SPICE ice core to construct a bounding model for the ice birefringence and the polarization time delays across ARA. The data-model comparison is consistent with the vertical girdle fabric at the South Pole having the prevailing horizontal crystallographic axis oriented near-perpendicular to ice flow. This study presents the possibility that ice birefringence can be used to constrain the range to a neutrino interaction, and hence aid in neutrino energy reconstruction, for in-ice experiments such as ARA.
In a previous work, we developed the idea to solve Keplers equation with a CORDIC-like algorithm, which does not require any division, but still multiplications in each iteration. Here we overcome this major shortcoming and solve Keplers equation using only bitshifts, additions, and one initial multiplication. We prescale the initial vector with the eccentricity and the scale correction factor. The rotation direction is decided without correction for the changing scale. We find that double CORDIC iterations are self-correcting and compensate possible wrong rotations in subsequent iterations. The algorithm needs 75% more iterations and delivers the eccentric anomaly and its sine and cosine terms times the eccentricity. The algorithm can be adopted for the hyperbolic case, too. The new shift-and-add algorithm brings Keplers equation close to hardware and allows to solve it with cheap and simple hardware components.
We present the recent development of hybridizable and embedded discontinuous Galerkin (DG) methods for wave propagation problems in fluids, solids, and electromagnetism. In each of these areas, we describe the methods, discuss their main features, display numerical results to illustrate their performance, and conclude with bibliography notes. The main ingredients in devising these DG methods are (i) a local Galerkin projection of the underlying partial differential equations at the element level onto spaces of polynomials of degree k to parametrize the numerical solution in terms of the numerical trace; (ii) a judicious choice of the numerical flux to provide stability and consistency; and (iii) a global jump condition that enforces the continuity of the numerical flux to obtain a global system in terms of the numerical trace. These DG methods are termed hybridized DG methods, because they are amenable to hybridization (static condensation) and hence to more efficient implementations. They share many common advantages of DG methods and possess some unique features that make them well-suited to wave propagation problems.
Radiative transfer has a strong impact on the collapse and the fragmentation of prestellar dense cores. We present the radiation-hydrodynamics solver we designed for the RAMSES code. The method is designed for astrophysical purposes, and in particular for protostellar collapse. We present the solver, using the co-moving frame to evaluate the radiative quantities. We use the popular flux limited diffusion approximation, under the grey approximation (one group of photon). The solver is based on the second-order Godunov scheme of RAMSES for its hyperbolic part, and on an implicit scheme for the radiation diffusion and the coupling between radiation and matter. We report in details our methodology to integrate the RHD solver into RAMSES. We test successfully the method against several conventional tests. For validation in 3D, we perform calculations of the collapse of an isolated 1 M_sun prestellar dense core, without rotation. We compare successfully the results with previous studies using different models for radiation and hydrodynamics. We have developed a full radiation hydrodynamics solver in the RAMSES code, that handles adaptive mesh refinement grids. The method is a combination of an explicit scheme and an implicit scheme, accurate to the second-order in space. Our method is well suited for star formation purposes. Results of multidimensional dense core collapse calculations with rotation are presented in a companion paper.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
