Questions tagged [astronomy]

Questions on the astronomical applications of Mathematica.

Filter by
Sorted by
Tagged with
2
votes
1answer
63 views

Charting sun location throughout the year from a specified geographic location

I'd like to chart the location of the sun (azimuth and elevation) from a given location (e.g., Chicago) throughout the year for a stated year (e.g., 2010). Here is the basic extraction of the current ...
1
vote
1answer
44 views

Units Question for PlanetData[..]

I asked a similar question yesterday, which was answered splendidly, but I missed adding a detail. Attached is a command, its output, and the issue I have now. This outputs g/cm^3, but I'd rather have ...
2
votes
2answers
75 views

Converting Units [duplicate]

I'm using MMA 12.2. Trying to put together some astronomy pages, but the output of the command below is in miles. I would rather have it in SI units instead. On my system evaluating ...
2
votes
1answer
55 views

Getting error Solve::units [closed]

I'm trying to solve the following formula for i, given values for all other variables. There are two equations, the first one is a formula to show the transit ...
22
votes
1answer
342 views

Travelling around Jupiter

Since we have soon historical conjunction of Jupiter and Saturn I have prepared code with visualization planet Jupiter and Saturn with moons. This code first been prepared for version 11.3 and with ...
0
votes
1answer
47 views

Quick question about “shifting x axis”

Imagine I have two lists: x={x1,x2,x3,x4,...} y={y1,y2,y3,y4,...} And I want a listplot ListPlot[Transpose[{x,y}]] Would do ...
5
votes
0answers
113 views

Fast N-Body Simulation with CUDA [closed]

There is a Mathematica code for 3 body simulation ...
3
votes
1answer
320 views

2 body orbits around a barycenter

I'm attempting to plot a barycentric orbit of the sun and Jupiter. Initial Conditions ...
4
votes
1answer
386 views

gravitational two body problem for the orbit of the sun and jupiter around their barycenter

This is a simulation of the sun and jupiter orbiting the respective barycenter. Initial conditions ...
2
votes
1answer
318 views

N body simulation

I'm trying to solve a 3 body problem consisting of the earth moon and sun. I'm wanting to show the earth orbiting the sun and the moon orbiting the earth as a result of the gravitational forces. ...
7
votes
1answer
114 views

Detect sub figures and split (automatically?)

I would like to take a large figure showing a sequence of individual plots and automatically split the into series of images showing the individual sub-plots (to be exported as animated GIF later). ...
2
votes
1answer
106 views

NDSolve Maxsteps Reached Error

I'm trying to solve a system of Equations. The code is commented below, I am mainly conserved with the last few lines, specifically the NDSolve line. No matter what I do to the Maxsteps option, the ...
3
votes
2answers
183 views

Solar Power Calculations

Good Afternoon, I have calculated the Incident beam radiation for a given locations local sunlight hours. ...
4
votes
1answer
59 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: <...
7
votes
1answer
472 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,...
3
votes
1answer
73 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, ...
7
votes
1answer
167 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 ...
0
votes
2answers
158 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; ...
1
vote
0answers
200 views
0
votes
1answer
489 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 ...
5
votes
1answer
66 views

Select Specific StarData

The help menu for StarData has many examples such as StarData[EntityClass[ "Star", {EntityProperty["Star", "DistanceFromEarth"] -> TakeSmallest[10]}]] ...
-1
votes
2answers
106 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 ...
0
votes
1answer
87 views
0
votes
1answer
408 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....
6
votes
1answer
1k 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 (...
1
vote
0answers
313 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?
4
votes
1answer
160 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 ...
1
vote
0answers
558 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: ...
0
votes
2answers
157 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 ...
0
votes
1answer
155 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:...
15
votes
1answer
1k 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 ...
6
votes
1answer
142 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, ...
17
votes
2answers
922 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’...
2
votes
1answer
79 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. ...
1
vote
2answers
102 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 ...
1
vote
0answers
105 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 ...
4
votes
1answer
133 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 ...
1
vote
1answer
191 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 ...
1
vote
1answer
88 views

TrueAnomaly of planet orbit as a function?

It is possible to obtain the true anomaly of a planet this way: ...
1
vote
1answer
2k 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 ...
2
votes
1answer
391 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 ...
4
votes
1answer
444 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 ...
2
votes
0answers
406 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 ...
4
votes
1answer
553 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 ...
6
votes
1answer
188 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 ...
11
votes
1answer
449 views

Calculate moonrise

Mathematica has Sunrise[], but no Moonrise[]. But it has MoonPosition[], so I thought it ...
12
votes
3answers
605 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)$, ...
10
votes
1answer
451 views

Plotting how sunrise times change over the year

I created a somewhat lengthy Python code to show how the time of sunrise varies over a year, resulting in the following figure (daylight savings were disregarded): The dataset can be found here. The ...
1
vote
0answers
211 views

Orbit recognition error in planet motion [closed]

...
4
votes
1answer
123 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: ...