Questions tagged [astronomy]

Questions on the astronomical applications of Mathematica.

Filter by
Sorted by
Tagged with
19
votes
3answers
2k views

Only final result from NDSolve

Finally, I started to play with differential equations in Mathematica. And I have faced the problem, which seems to me so basic that I'm afraid this question is going to be closed soon. However, I'...
49
votes
4answers
4k views

Mollweide maps in Mathematica

Context In my field of research, many people use the following package: healpix (for Hierarchical Equal Area isoLatitude Pixelization) which has been ported to a few different languages (F90, C,C++, ...
26
votes
4answers
7k views

Composition: how to make a day and night world map?

Given the following world images: ...
15
votes
3answers
3k views

How do I compute sidereal time and JD?

There was once a package that provided a wide range of functions for computing astronomical information. This package now appears to be out of date; but as near as I can tell, many of the useful ...
19
votes
2answers
1k views

Using map projections with astronomical data

I noticed that all important "Geoprojections" are available in projections for a spherical reference models: GeoProjectionData function; 1 - I am trying using the sinusoidal projection for ...
14
votes
1answer
684 views

Speed up NDSolve compared to Python (calls to LSODA)

I migrated a numerical model code from Python to Mathematica and am surprised how much faster the Python version runs. Profiling of the Python version tells me that it is about 100 times faster (120 ...
10
votes
1answer
614 views

How can I get this Graphics3D to rotate around a particular point?

I made this Graphics3D of the Sun with the plotted orbits of a few astronomical objects. ...
4
votes
1answer
415 views

Orbit followed by a particle around Schwarzschild Black Hole

The following is a equation which describes various possible orbits of a particle around the Schwarzschild black hole spacetime in general relativity. I want to solve it from ...
36
votes
2answers
3k views

The Orbit and Perigee of the Flamsteed comet

Historical context This year we have the 330-th anniversary of the Battle of Vienna - one of the great formative events of European history, it took place on September 12, 1683. Kara Mustafa, Grand ...
7
votes
3answers
1k views

RK4 Gravity Simulator

I have the following RK4 solver which splits the two 2nd order ODEs, used to calculate x and y positions under the influence of a gravitating body where $$x''(t)=\frac{G m x(t)}{(x(t)^2+y(t)^2)^{3/2}}$...
16
votes
1answer
403 views

Is there a way to get AbsoluteTime to correctly report absolute time differences?

AbsoluteTime supposedly reports the "total number of seconds" between two dates, but it does not. For example, because of the leap second applied at 2012-06-30T23:...
14
votes
1answer
664 views

GeoProjection for astronomical data - wrong ticks

Here is an excellent answer to what I wanted to do, and it works like a charm. Unfortunately, when I want to change it a bit, it fails. For setting the stage, here is the datafile, and this is my ...
11
votes
3answers
488 views

How to speed up plotting with PlotPoints?

I'm plotting an ellipse 1/(1 - e Cos[x]) with e=0.9 using PolarPlot. Taking $x\in[0,2\pi)$, ...
9
votes
1answer
333 views

How do I determine astronomical transit times?

Mathematica, through the "legacy" Scientific Astronomer package, used to have the ability to easily determine accurate transit times of astronomical objects simply, with ...
26
votes
1answer
950 views

Interstellar Image Effects

I'd like to know how to create this type of gravitational lensing effect in an image (here are some star fields from the hubble) "I wrote down the equations and I tested them in Mathematica by ...
29
votes
2answers
9k views

Creating a simulation of our Solar System

As the title says, my objective is to create a simulation of the motion of the planets in our Solar System using Mathematica. All the theoretical background regarding the equations of motion of ...
21
votes
1answer
1k views

Right Ascension and Declination in astronomical functions

In the process of answering this question, I was forced to confront the various astronomical coordinate systems used by Mathematica. Background In astronomy, positions of celestial objects (stars, ...
13
votes
1answer
1k views

Geocentric, Earth-fixed, coordinates of Sun, Moon (and others objects eventually)

I know I can get "HelioCoordinates" of the Earth and the Moon with ...
8
votes
2answers
645 views

Extracting and displaying data from GW150914 hdf5 files

Has anyone looked at the LIGO data using Mathematica or know how to do this? I'm not familiar with the hdf5 format but I've been able to import the data, consisting primarily of 131,072 floating ...
13
votes
2answers
444 views

How is AstronomicalData interpreting its date arguments?

It is easy to confirm that the AstronomicalData is not interpreting its arguments as UTC (or "GMT", which is, in any case, ambiguous) dates by verifying that the ...
5
votes
1answer
218 views

Inaccurate Heliocentric coordinates of PlanetaryMoonData

I have noticed that the heliocentric coordinates obtained with PlanetaryMoonData are very inaccurate. Fore example, consider the Galilean moons of Jupiter: ...
7
votes
1answer
186 views

Mass distribution in NDSolve

I'm trying to model the solar system with NDSolve, basically a numerical n-body solution. I model the initial positions of bodies from a data set thusly: ...
5
votes
1answer
590 views

Using ParametricNDSolve and NMinimize to find accurate interplanetary trajectory arrival periapse radius

I have an interplanetary trajectory simulation that calculates the required velocity (using a Lambert solver) and escape angle for a spacecraft to travel from Earth to Mars. This model does not use a ...
2
votes
1answer
240 views

Using DeleteCases with multiple AstronomicalData properties

I am attempting to use AstronomicalData and DeleteCases together to remove invalid cases of two properties while listing the ...
8
votes
2answers
745 views

Problems getting clean sunrise, sunset, civil, nautical and astronomical twilight data

If in a notebook you enter: ...
8
votes
3answers
503 views

Why is AstronomicalData so slow?

I'm having trouble getting some pretty simple functions using AstronomicalData to perform at anywhere near the speed I need them to. For example, I have many ...
10
votes
2answers
530 views

Impact Crater Simulation

I have been trying to model crater formation on a given planetary surface - $500\,\textrm{km}^2$. The locations of impacts are random, however, if an impact is within $30\,\textrm{km}$ of another, the ...
7
votes
2answers
970 views

Combining three R, G, B matrices into a single matrix of 3-tuples?

Working with Mathematica 8 on OSX here. I have three $1024 \times 1024$ matrices representing the counts from consecutive exposures onto a CCD with $R$, $G$, and $B$ filters on it, and am looking to ...
6
votes
1answer
365 views

Can NDSolve address spherical gravitational collapse?

As far as I know spherical gravitational collapse -of central importance to theoretical physics- is thought to be out of the scope of Wolfram Mathematica. However there are cases, as in this article,...
1
vote
0answers
940 views

Solar System Orbital Parameters [duplicate]

I recently made an n-body simulation and thought it would be interesting to try and model it on the interior of our solar system (Sun to Mars), but I cannot find initial conditions for such an ...
7
votes
0answers
161 views

MoonPosition and TimeZone support

Bug introduced in 10.0 and fixed in 10.2 MoonPosition is new in 10.0 With this question, I noticed some some weird results of ...
3
votes
1answer
193 views

Computing planet conjunctions with 2D circular orbits still hard?

All methods I've seen for computing planetary/stellar conjuctions (when two planets or a planet and a star appear close together in the sky because their angular separation, as measured from Earth, ...
2
votes
2answers
166 views

Finding Minimium Periapse Distance Between Spacecraft and Planet

I have the following simple model that simulates a spacecraft orbiting Mars. What I'm hoping to do is be able to find the periapse radius of the spacecraft as I did below, except instead of using ...
0
votes
1answer
89 views

Plotting the orbit of the Earth around the Sun. Error: step size is effectively zero; singularity or stiff system suspected

I've been asked to plot the orbit of the Earth around the sun over a period of three years. I've copied my exact code below. When I try to run it I get a series of errors, the first being: NDSolve:...
0
votes
1answer
146 views

Plotting elliptical orbits using Verlet method

I'm looking to solve a 2 body system consisting of the Earth and a body of mass in orbit. I'm using the Verlet method to numerically integrate Newton's equation of motion in order to plot an ...