A Visualization of Null Geodesics for the Bonnor Massive Dipole

In this work we simulate null geodesics for the Bonnor massive dipole metric by implementing a symbolic-numerical algorithm in Sage and Python. This program is also capable of visualizing in 3D, in principle, the geodesics for any given metric. Geodesics are launched from a common point, collectively forming a cone of light beams, simulating a solid-angle section of a point source in front of a massive object with a magnetic field. Parallel light beams also were considered, and their bending due to the curvature of the space-time was simulated.Comment: 9 pages, 5 figures, Presented at XIX Simposio Internacional de M\'etodos Matem\'aticos Aplicados a las Ciencias (19th International Symposium of Mathematical Methods Applied to Sciences, XIX-SIMMAC

MFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver

MFC is an open-source tool for solving multi-component, multi-phase, and bubbly compressible flows. It is capable of efficiently solving a wide range of flows, including droplet atomization, shock–bubble interaction, and bubble dynamics. We present the 5- and 6-equation thermodynamically-consistent diffuse-interface models we use to handle such flows, which are coupled to high-order interface-capturing methods, HLL-type Riemann solvers, and TVD time-integration schemes that are capable of simulating unsteady flows with strong shocks. The numerical methods are implemented in a flexible, modular framework that is amenable to future development. The methods we employ are validated via comparisons to experimental results for shock–bubble, shock–droplet, and shock–water-cylinder interaction problems and verified to be free of spurious oscillations for material-interface advection and gas–liquid Riemann problems. For smooth solutions, such as the advection of an isentropic vortex, the methods are verified to be high-order accurate. Illustrative examples involving shock–bubble-vessel-wall and acoustic–bubble-net interactions are used to demonstrate the full capabilities of MFC

Improvements to the APBS biomolecular solvation software suite

The Adaptive Poisson-Boltzmann Solver (APBS) software was developed to solve the equations of continuum electrostatics for large biomolecular assemblages that has provided impact in the study of a broad range of chemical, biological, and biomedical applications. APBS addresses three key technology challenges for understanding solvation and electrostatics in biomedical applications: accurate and efficient models for biomolecular solvation and electrostatics, robust and scalable software for applying those theories to biomolecular systems, and mechanisms for sharing and analyzing biomolecular electrostatics data in the scientific community. To address new research applications and advancing computational capabilities, we have continually updated APBS and its suite of accompanying software since its release in 2001. In this manuscript, we discuss the models and capabilities that have recently been implemented within the APBS software package including: a Poisson-Boltzmann analytical and a semi-analytical solver, an optimized boundary element solver, a geometry-based geometric flow solvation model, a graph theory based algorithm for determining p$K_a$ values, and an improved web-based visualization tool for viewing electrostatics

Fast and accurate frequency-dependent radiation transport for hydrodynamics simulations in massive star formation

Context: Radiative feedback plays a crucial role in the formation of massive stars. The implementation of a fast and accurate description of the proceeding thermodynamics in pre-stellar cores and evolving accretion disks is therefore a main effort in current hydrodynamics simulations. Aims: We introduce our newly implemented three-dimensional frequency dependent radiation transport algorithm for hydrodynamics simulations of spatial configurations with a dominant central source. Methods: The module combines the advantage of the speed of an approximate Flux Limited Diffusion (FLD) solver with the high accuracy of a frequency dependent first order ray-tracing routine. Results: We prove the viability of the scheme in a standard radiation benchmark test compared to a full frequency dependent Monte-Carlo based radiative transfer code. The setup includes a central star, a circumstellar flared disk, as well as an envelope. The test is performed for different optical depths. Considering the frequency dependence of the stellar irradiation, the temperature distributions can be described precisely in the optically thin, thick, and irradiated transition regions. Resulting radiative forces onto dust grains are reproduced with high accuracy. The achievable parallel speedup of the method imposes no restriction on further radiative (magneto-) hydrodynamics simulations. Conclusions: The proposed approximate radiation transport method enables frequency dependent radiation hydrodynamics studies of the evolution of pre-stellar cores and circumstellar accretion disks around an evolving massive star in a highly efficient and accurate manner.Comment: 16 pages, 11 figure