• Home
  • About Us
  • Contact Us
  • Privacy Policy
  • Special Offers
Business Intelligence Info
  • Business Intelligence
    • BI News and Info
    • Big Data
    • Mobile and Cloud
    • Self-Service BI
  • CRM
    • CRM News and Info
    • InfusionSoft
    • Microsoft Dynamics CRM
    • NetSuite
    • OnContact
    • Salesforce
    • Workbooks
  • Data Mining
    • Pentaho
    • Sisense
    • Tableau
    • TIBCO Spotfire
  • Data Warehousing
    • DWH News and Info
    • IBM DB2
    • Microsoft SQL Server
    • Oracle
    • Teradata
  • Predictive Analytics
    • FICO
    • KNIME
    • Mathematica
    • Matlab
    • Minitab
    • RapidMiner
    • Revolution
    • SAP
    • SAS/SPSS
  • Humor

Make cool animations of an accretion disk/planetary disk

February 27, 2020   BI News and Info

I am trying to make an animation of an accretion disk or a planetary disk.
For example, a cool disk showing below (created by NAHKS TR’EHNL):
w9rv9 Make cool animations of an accretion disk/planetary disk
or in this link

Here is what I have been tried so far:

step1:define some basic parameters.

{rmin, rmax, scale} = {0.1, 1.0, 0.2};

where rmin is the inner boundary of the disk;
rmax is the outer boundary of the disk;
scale is the scale height of the disk

Step 2: Define a function that gives the shape of the intersection/crosssection of the disk:

intersec[r_] := scale*(r - rmin)^0.5/(rmax - rmin)^0.5

You can see that the function =0 at the inner boundary and =scale (in this case 0.2) at the outer boundary of the disk. This will outline the height of the disk. Below is how this function looks like when plotted by the following code:

Plot[intersec[r], {r, 0, rmax}, PlotRange -> All, Axes -> False, Frame -> True, FrameLabel -> {"radius", "disk height"}]

bu3hW Make cool animations of an accretion disk/planetary disk

Step 3 Making a parametric plot of the disk by using the color function:

figdisk = ParametricPlot3D[{r Cos[\[Theta]], r Sin[\[Theta]], 
intersec[r]}, {r, rmin, rmax}, {\[Theta], 0, 2 \[Pi]}, PlotPoints -> 100, Mesh -> False, Axes -> False, Boxed -> False, BoundaryStyle -> Opacity[0], Lighting -> {{"Ambient", White}}, Background -> Black, ColorFunctionScaling -> False, ColorFunction -> (ColorData[{"SolarColors", "Reverse"}][Abs@#3/scale + 0.1 RandomReal[{-1, 1}]] &)];

The trick here is set the color of the disk as a funtion of the radius, or in terms of the intersec funtion, as a function of the disk height. So it shows temperature gradient effect. Then random values are assigned to the colorfunction, so to mimic the random fluctuations.

For black hole, a simple sphere:

figBH = Graphics3D[{Black, Sphere[{0, 0, 0}, rmin]}];

for jet I use a cone:

figjet = ParametricPlot3D[{0.05 z Cos[t], 0.05 z Sin[t], 
z}, {z, -rmin, 5 rmin}, {t, 0, 2 \[Pi]}, Mesh -> False, PlotStyle -> Glow[Blend[{White, Lighter@Blue}]]];

Then

Show[figdisk, figBH, figjet]

kr4F7 Make cool animations of an accretion disk/planetary disk

Now the problem is that:

  1. The way I choose to plot the disk only gives a surface, which is not symmetric above and below the equatorial plane. you can see this by the following figure with viewpoint adjusted and white background;
    Um1OJ Make cool animations of an accretion disk/planetary disk
  2. The appearance of the example given at the top of this post seems to be natural. It is more like a real flow with fluctuations and even some random magnetic field lines on it.
  3. Same for the jet, the example’s jet is more realistic with fluid-like materials ejected outward while mine jet….. I have seen this post that simulated realistic jet engine flame, but it does not look like the fragmented thing in the example.

    Can anyone help me to improve this piece of art? Thank you so much in advance!

Let’s block ads! (Why?)

Recent Questions – Mathematica Stack Exchange

accretion, animations, Cool, Disk, disk/planetary
  • Recent Posts

    • Now make soup!
    • Attach2Dynamics Or SharePoint Security Sync – Choose your smart app for effective document management in Dynamics 365 CRM/Power Apps.
    • 5 jobs that you should apply for this week (before it’s too late)
    • SQL Server authentication methods, logins, and database users
    • DAE solver fails for system of coupled partial differential equations
  • Categories

  • Archives

    • March 2021
    • February 2021
    • January 2021
    • December 2020
    • November 2020
    • October 2020
    • September 2020
    • August 2020
    • July 2020
    • June 2020
    • May 2020
    • April 2020
    • March 2020
    • February 2020
    • January 2020
    • December 2019
    • November 2019
    • October 2019
    • September 2019
    • August 2019
    • July 2019
    • June 2019
    • May 2019
    • April 2019
    • March 2019
    • February 2019
    • January 2019
    • December 2018
    • November 2018
    • October 2018
    • September 2018
    • August 2018
    • July 2018
    • June 2018
    • May 2018
    • April 2018
    • March 2018
    • February 2018
    • January 2018
    • December 2017
    • November 2017
    • October 2017
    • September 2017
    • August 2017
    • July 2017
    • June 2017
    • May 2017
    • April 2017
    • March 2017
    • February 2017
    • January 2017
    • December 2016
    • November 2016
    • October 2016
    • September 2016
    • August 2016
    • July 2016
    • June 2016
    • May 2016
    • April 2016
    • March 2016
    • February 2016
    • January 2016
    • December 2015
    • November 2015
    • October 2015
    • September 2015
    • August 2015
    • July 2015
    • June 2015
    • May 2015
    • April 2015
    • March 2015
    • February 2015
    • January 2015
    • December 2014
    • November 2014
© 2021 Business Intelligence Info
Power BI Training | G Com Solutions Limited