Cardea modellezés

Aszteroida becsapódás szimuláció

--------------------

import math

# ============================================================

# Cardea aszteroida - egyszerű becsapódási szimuláció

# Oktatási célú modell

# ============================================================

# ---------- Aszteroida paraméterei ----------

ATMERO = 163              # méter

SURUSEG = 2700            # kg/m³ (kőzet)

SEBESSEG = 20000          # m/s (20 km/s)

# ---------- Becsapódási paraméterek ----------

BECSAPODASI_SZOG = 45     # fok

# ---------- New York becsült népsűrűsége ----------

NEPSURUSEG = 11000        # fő/km²

# ------------------------------------------------

sugar = ATMERO / 2

terfogat = (4/3) * math.pi * sugar**3

tomeg = terfogat * SURUSEG

energia = 0.5 * tomeg * SEBESSEG**2

MEGATONNA = energia / 4.184e15

# ------------------------------------------------

# Egyszerű empirikus becslések

# ------------------------------------------------

if MEGATONNA < 10:

    pusztulas = 2

elif MEGATONNA < 50:

    pusztulas = 3

elif MEGATONNA < 150:

    pusztulas = 4

elif MEGATONNA < 300:

    pusztulas = 5

elif MEGATONNA < 600:

    pusztulas = 7

else:

    pusztulas = 10

# lökéshullám

lokeshullam = pusztulas * 2

# hősugárzás

hosugar = pusztulas * 1.5

# kráter (nagyon durva közelítés)

krater = ATMERO * 10 / 1000

# ------------------------------------------------

# Áldozatok becslése

# ------------------------------------------------

terulet = math.pi * pusztulas**2

lakossag = terulet * NEPSURUSEG

# a teljes pusztulási zónában

halalozas = lakossag * 0.95

# sérültek

serultek = lakossag * 1.8

# ------------------------------------------------

print("=" * 60)

print("CARDEA ASZTEROIDA BECSAPÓDÁSI SZIMULÁCIÓ")

print("=" * 60)

print(f"Aszteroida átmérője:      {ATMERO} m")

print(f"Tömege:                   {tomeg:,.0f} kg")

print(f"Sebessége:                {SEBESSEG/1000:.1f} km/s")

print(f"Becsapódási energia:      {MEGATONNA:.1f} Mt TNT")

print()

print(f"Kráter átmérő:            {krater:.2f} km")

print(f"Teljes pusztulási sugár:  {pusztulas:.1f} km")

print(f"Lökéshullám sugara:       {lokeshullam:.1f} km")

print(f"Hősugárzás sugara:        {hosugar:.1f} km")

print()

print("----- New York hipotetikus becslés -----")

print(f"Teljes pusztulási terület: {terulet:.1f} km²")

print(f"Becsült halálos áldozatok: {halalozas:,.0f} fő")

print(f"Becsült sérültek:          {serultek:,.0f} fő")

print("=" * 60)

print("FIGYELEM!")

print("Ez kizárólag oktatási célú, egyszerűsített modell.")

print("Nem valódi veszély-előrejelzés.")

print("=" * 60)

----------------

Eredmény;

============================================================

CARDEA ASZTEROIDA BECSAPÓDÁSI SZIMULÁCIÓ

============================================================

Aszteroida átmérője:      163 m

Tömege:                   6,122,449,332 kg

Sebessége:                20.0 km/s

Becsapódási energia:      292.7 Mt TNT

Kráter átmérő:            1.63 km

Teljes pusztulási sugár:  5.0 km

Lökéshullám sugara:       10.0 km

Hősugárzás sugara:        7.5 km

----- New York hipotetikus becslés -----

Teljes pusztulási terület: 78.5 km²

Becsült halálos áldozatok: 820,741 fő

Becsült sérültek:          1,555,088 fő

============================================================

Közvetett áldozatok; 8 millió fő

============================================

Futtás; https://onecompiler.com/python#draft-8v69