# Questions tagged [astronomy]

Questions on the astronomical applications of Mathematica.

118 questions
Filter by
Sorted by
Tagged with
83 views

### Solar Power Calculations

Good Afternoon, I have calculated the Incident beam radiation for a given locations local sunlight hours. ...
376 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,...
144 views

### How to plot the hour of the sunsets

i have this code: ...
32 views

### StarData, Altitude/Azimuth Anomaly with TimeZones

I am trying to compute lunar parallax using the position of the star Regulus on February 19th of this year. I get different answers depending on how I specify the date. Here are the particulars: <...
114 views

### Is MoonPosition working as expected?

I have noticed some surprising behavior when executing the MoonPosition function. When I provide a range of dates, I get a TemporalData answer: ...
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++, ...
56 views

### Finding the next total SolarEclipse over a given latitude and longitude?

What is a reasonable way to code the title question in Mathematica? Trying the following gave Polygons over GeoPositions over Lists of pairs, but just digging into these seemed inelegant. Besides, ...
122 views

### StarData: proper motion of a star in HelioCoordinates

I am not an astronomer, but am trying to calculate the trajectory of an interstellar space probe, for example to Proxima Centauri. The HelioCoordinates seem appropriate, as typically used within the ...
402 views

### Interpreting Mathematica code on black holes

I am trying to understand the code written down on page 7 of this document (code is in Mathematica) I understand pretty much all of the code on the previous page needed to setup the page 7 code (...
111 views

### Plotting a parametric solution from NDSolve

Im trying to animate multiple points around multiple parametric plots which depicts orbital motion of body around a planet. Initial conditons are; ...
3k views

### Circular Orbit Velocity in Vector Components

I'm trying to find out how to calculate the initial velocity needed to keep a spacecraft in a circular orbit around the Earth given its initial position and an initial angle theta in a Cartesian ...
59 views

...
154 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 ...
697 views

### Can I create a radar chart of sunrise/sunset times for a year

Can I create a radar chart (or just a simple curve) of sunrise/sunset times, for a given location thougout a year? It would be something like this, but with a smooth curve: (source: hassanhameed.com)...
60 views

### Integrating the Planck Function over certain filter ranges [closed]

ClearAll["Global`*"] ...
58 views

### Select Specific StarData

The help menu for StarData has many examples such as StarData[EntityClass[ "Star", {EntityProperty["Star", "DistanceFromEarth"] -> TakeSmallest[10]}]] ...
45 views

### a ListAnimate problem

Thanks to @AlexTrounev, I learned how to use ListAnimate, but I still get stuck using it. Before ListAnimate everything is correct. I just cannot successfully use ListAnimate to make the plot. I have ...
122 views

### Animate particle's trajectory

Does anyone know how to animate the trajectory of a particle, but make sure before the particle reaching a certain spot, there's nothing. I mean I want to make an animation that look like https://en....
129 views

### Solving the values of Christoffel symbols from a given metric tensor

Is there any easy mathematica package by which I can solve the Christoffel symbols from a given metric tensor?
98 views

### Creating the celestial sphere for astronomical demonstrations

Quick question: has anybody produced (and/or willing to produce and share) a nice celestial sphere with some stars and/or constellations? I am traveling , so I cannot do it myself right now, but I ...
54 views

### Calculating Christoffel Symbols [duplicate]

I am totally new in the astrophysics field and I need to calculate Christoffel symbols for a model. The line element for the FLRW metric and the geodesic equation is given below: ...
109 views

### Using WhenEvent to calculate the period of an orbit

I'm having trouble trying to have Mathematica tell me when the period of the orbit of my system is. Namely, I'm trying to use Mathematica to find the orbital period of the Earth around the Sun. Here's ...
131 views

### Does SolarEclipse only cover a limited date range?

If I use the example from the documentation for SolarEclipse and simply replace the years there with older ones, say from around Shakespeare's time, ...
95 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:...
77 views

### InverseFunction's domain is not the range of the function that was inverted

I wrote a function to calculate how long of a period has past when a planet traverses an angle theta from perihelion in an orbit of a given eccentricity. ...
200 views

...
704 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 ...
77 views

### How do I pull a multiple pieces of data from a while loop and put into a table

I am wanting to pull multiple pieces of data from a while loop and put it in a list. The example I have added below is a table with inputted data that I found. Im pretty sure the syntax is unrelated ...
562 views

### Processing images from NASA's Juno Jupiter Orbiter

Question Can you help me improve what I’ve done to create an enhanced-color picture from an image provided by NASA’s JunoCam? [My work is shown in the penultimate image of Jupiter in this post.] I’...
110 views

### Typeset the standard “Moon” symbol

As it says on the tin, I want a way to typeset the astronomical symbol for the moon, which is a crescent. I also want this in such a way that I could integrate it into a button as a label, like the ...
82 views

### SatelliteData for a hypothetical satellite

SatelliteData appears to have some sort of orbit propagator behind it. It works for existing satellites, but is there any way to use it to predict the orbit of a ...
599 views

### Series solution of the Lane-Emden equation

I have the following problem: I'd like to show how the Lane-Emden equation would look like if someone would solve it as a series. Here's the code: ...
116 views

### list of dates for vernal equinox

I'm trying to get a list of dates and times of the vernal equinox for the last 30 years (or some other time interval). I'm flailing about and really have gotten nowhere. When I've tried typing ...
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 ...
1k 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 ...
1k views

### How to convert astronomical RA and Dec from sexagesimal to degree?

Astronomers use right ascension (RA) and declination (Dec) to refer to the sky coordinates of celestial objects. RA and Dec can take be represented degrees or in sexagesimal (base 60). How can I ...
78 views

### TrueAnomaly of planet orbit as a function?

It is possible to obtain the true anomaly of a planet this way: ...
208 views

### Calculating with scientific notation

I want to calculate the gravitational force between the Earth and its Moon using Newton's gravitational formula: (G M1 M2)/d^2 and its appropriate physical units. I'd include the code I wrote for ...
340 views

### Lagrangian form of Earth-Sun Kepler blows up NDSolve?

I've looked at this 1-D SE question. That question seems to have difficulties around diverging potential functions, and I don't have that. I've looked at this 2-D 3-body demonstration, but that starts ...
256 views

### Plotting orbits of trajectories of bodies sharing the earth's orbit in coordinates centered at the barycenter of the solar system [closed]

I am attempting to try and plot the trajectories of a trio of satellites that are orbiting with the earth-moon system around the solar system's barycenter. Originally, I had the satellites in orbit ...
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 ...
5k views

### What is the closest star to Earth other than the sun?

How to find the Sun's nearest neighbor star using Mathematica? I tried but it didn't work
150 views

### TimeZone Anomaly

I ran into the following situation when I was trying to determine the time of sunrise (using the built-in Sunrise function) for a variable location on a variable ...
618 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. ...
419 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 ...
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 ...
343 views

### 3D orbits of moons around their respective planets

I have seen a lot of solar system models (both animated and not), but I am trying to replicate the same thing with moons. Either just one planet and its lunar system or the our whole solar system with ...
670 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 ...
I'm plotting an ellipse 1/(1 - e Cos[x]) with e=0.9 using PolarPlot. Taking $x\in[0,2\pi)$, ...