File:Meteor 60m air burst pressure wind sound 1.png

From Wikimedia Commons, the free media repository
Jump to navigation Jump to search

Original file ‎(1,267 × 719 pixels, file size: 125 KB, MIME type: image/png)

Captions

Captions

Pressure, wind and sound of diameter 60 m stony meteor

Summary

[edit]
Description
English: Pressure, wind and sound of diameter 60 m stony meteor
Date
Source Own work
Author ChatGPT 4o and 3 AI

Source of data is "Impact earth". Code by ChatGPT

Impact Earth! Robert Marcus,J Jay Melosh,Gareth Collins https://www.purdue.edu/impactearth/

Properties of asteroid

60m 2750 kgm/m3, angle 45, velocity 20 ms sedimentary rock

Python3 source code

import numpy as np import matplotlib import matplotlib.pyplot as plt from scipy.interpolate import interp1d

  1. Alkuperäinen data

data = np.array([ 

   [1, 17.2, 35.1, 330, np.nan],
   [5, 22.7, 21.1, 232, 103],
   [7, 27.1, 17.5, 202, 102],
   [10, 34.7, 12.1, 152, 99],
   [12, 40.1, 8.46, 114, 96],
   [15, 48.5, 5.43, 78.2, 92],
   [20, 1.05, 3.24, 49.2, 87],
   [25, 1.29, 2.22, 34.6, 84],
   [30, 1.54, 1.66, 26.3, 81],
   [35, 1.79, 1.31, 20.9, 79],
   [40, 2.04, 1.08, 17.3, 78],
   [45, 2.29, 0.909, 14.7, 76],
   [50, 2.54, 0.785, 12.7, 75],
   [80, 4.05, 0.427, 7.00, 70],
   [100, 5.06, 0.326, 5.36, 67],
   [150, 7.58, 0.205, 3.38, 63],
   [200, 10.1, 0.15, 2.47, 60],
   [300, 15.2, 0.0973, 1.61, 57],
   [400, 20.2, 0.0721, 1.19, 54],
   [500, 25.3, 0.0573, 0.95, 52],
   [600, 30.3, 0.0476, 0.788, 51],
   [1000, 50.5, 0.0283, 0.470, 46]

])

  1. Etäisyysarvot (km), joille haluamme interpoloinnin

distances_to_interpolate = np.array([10,20,30,40,50,60,80,100])

  1. Dataa eroteltu eri muuttujiin

r_km = data[:, 0] peak_blast_psi = data[:, 2] max_wind_ms = data[:, 3] loud_db = data[:, 4]

  1. Luodaan logaritminen interpolointifunktio

log_interp_peak_blast_psi = interp1d(np.log(r_km), np.log(peak_blast_psi), fill_value="extrapolate") log_interp_max_wind_ms = interp1d(np.log(r_km), np.log(max_wind_ms), fill_value="extrapolate") log_interp_loud_db = interp1d(np.log(r_km), np.log(loud_db), fill_value="extrapolate")

  1. Lasketaan interpoloidut arvot ja muunnetaan takaisin eksponentiaaliselle asteikolle

interp_peak_blast_psi = np.exp(log_interp_peak_blast_psi(np.log(distances_to_interpolate))) interp_max_wind_ms = np.exp(log_interp_max_wind_ms(np.log(distances_to_interpolate))) interp_loud_db = np.exp(log_interp_loud_db(np.log(distances_to_interpolate)))

  1. Pyöristetään tulokset kahden desimaalin tarkkuuteen

interp_peak_blast_psi = np.round(interp_peak_blast_psi, 2) interp_max_wind_ms = np.round(interp_max_wind_ms, 2) interp_loud_db = np.round(interp_loud_db, 2)

  1. Piirretään kuvaajat

fig, ax1 = plt.subplots(figsize=(10, 6))

matplotlib.rcParams.update({'font.size': 16})

ax1.xaxis.set_tick_params(labelsize=16) ax1.yaxis.set_tick_params(labelsize=16)

  1. Ensimmäinen pystyakseli: tuulen nopeus ja äänen voimakkuus

ax1.set_xlabel("Distance from ground zero (km)", fontsize=18) ax1.set_ylabel("Wind (m/s) ja Sound (dB)", color='black', fontsize=18)

  1. Piirretään tuulen nopeuden käyrä

ax1.plot(distances_to_interpolate, interp_max_wind_ms, 'x--', label="Wind (m/s)", color='green')

  1. Piirretään äänen voimakkuuden käyrä

ax1.plot(distances_to_interpolate, interp_loud_db, 'x:', label="Sound (dB)", color='red')

  1. Lisätään merkinnät tuulen nopeuden ja äänen voimakkuuden kohdalle

for i, dist in enumerate(distances_to_interpolate): 

   ax1.text(dist, interp_max_wind_ms[i], f"{interp_max_wind_ms[i]}", ha='right', color='green')
   ax1.text(dist, interp_loud_db[i], f"{interp_loud_db[i]}", ha='right', color='red')
  1. Toinen pystyakseli: paine (psi)

ax2 = ax1.twinx() ax2.xaxis.set_tick_params(labelsize=16) ax2.yaxis.set_tick_params(labelsize=16)

ax2.set_ylabel("Pressure (PSI)", color='blue', fontsize=18)

  1. Piirretään paineen käyrä

ax2.plot(distances_to_interpolate, interp_peak_blast_psi, 'x-', label="Pressure (psi)", color='blue')

  1. Lisätään merkinnät paineen kohdalle

for i, dist in enumerate(distances_to_interpolate): 

   ax2.text(dist, interp_peak_blast_psi[i], f"{interp_peak_blast_psi[i]}", ha='right', color='blue')
  1. Asetetaan molemmille akselille samat x-arvot ja lisätään ruudukko

ax1.set_xticks(distances_to_interpolate) ax1.grid(True) # Normaali ruudukko vasemmalle akselille

  1. Lisätään katkoviivallinen ruudukko oikealle akselille (paineelle)

ax2.grid(True, axis='y', linestyle='--') # Ruudukko vain y-akselille, katkoviivalla

  1. Lisätään otsikko ja legenda

fig.suptitle("Air burst of meteor, that has kinetic E=1.49 MT burst E=1.36 MT", fontsize=18) plt.title("60 m stone, 2750 kg m-3, 20 ms-1, entry angle 45 degrees, fragmentation H=5.58 km", fontsize=14) ax1.legend(loc='upper left', fontsize=16) ax2.legend(loc='upper right', fontsize=16)

  1. Näytetään kuvaaja

plt.show()

Licensing

[edit]
I, the copyright holder of this work, hereby publish it under the following license:
w:en:Creative Commons
attribution
This file is licensed under the Creative Commons Attribution 4.0 International license.
You are free:
  • to share – to copy, distribute and transmit the work
  • to remix – to adapt the work
Under the following conditions:
  • attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.


Public domain
This file is in the public domain because it is the work of a computer algorithm or artificial intelligence and does not contain sufficient human authorship to support a copyright claim.
The United Kingdom and Hong Kong provide a limited term of copyright protection for computer-generated works of 50 years from creation. [1] [2]
AI derivative works Legal disclaimer
Most image-generating AI models were trained using works that are protected by copyright. In some cases, such assets and models can produce images that contain major copyrightable elements of those copyrighted training images, making these outputs derivative works. Accordingly, there is a risk that AI-generated art uploaded on Commons may violate the rights of the authors of the original works. See Commons:AI-generated media for additional details.

azərbaycanca ∙ Deutsch ∙ English ∙ español ∙ français ∙ galego ∙ हिन्दी ∙ 日本語 ∙ português do Brasil ∙ русский ∙ slovenščina ∙ Türkçe ∙ Tiếng Việt ∙ 中文 ∙ 中文（简体） ∙ 中文（繁體） ∙ +/−

File history

Click on a date/time to view the file as it appeared at that time.

Date/TimeThumbnailDimensionsUserComment
current13:21, 26 September 2024Thumbnail for version as of 13:21, 26 September 20241,267 × 719 (125 KB)Merikanto (talk | contribs)Uploaded own work with UploadWizard

You cannot overwrite this file.

There are no pages that use this file.

Metadata

This file contains additional information such as Exif metadata which may have been added by the digital camera, scanner, or software program used to create or digitize it. If the file has been modified from its original state, some details such as the timestamp may not fully reflect those of the original file. The timestamp is only as accurate as the clock in the camera, and it may be completely wrong.

Retrieved from "https://commons.wikimedia.org/w/index.php?title=File:Meteor_60m_air_burst_pressure_wind_sound_1.png&oldid=929207318"