달이 멈춘다면 ?!

달이 멈춘다면?!

만약 달이 멈춘다면 어떻게 될까??

달과 지구가 서로 끌어다니는 만유인력으로 인해 상대적으로 질량이 큰 지구쪽으로 끌려간다.

결국 달은 지구와 충돌하게 된다!!

 

이번 시간은 무엇보다 물리코딩의 장점이 드러나는 시간이었다.

물리코딩을 통해 우리가 평소 관측하기 힘든 물리 현상을  쉽게 살펴볼 수 있다!!

비록 작용하는 힘의 요소는 너무 많기 때문에, 실제와 100% 같을 순 없겠지만, 

그럼에도 물리 코딩을 통해 많은 물리 현상을 직접적으로 눈으로 확인할 수 있어서 너무 유익하다~~!

 

GlowScript 3.0 VPython
#Creating Objects
Earth = sphere(pos = vector(0,0,0), radius = 6400000, texture = textures.earth)
Moon = sphere(pos = vector(385000e3,0,0), radius = 1737000, make_trail =True)
sf = 6 #scailing factor
Earth.radius = sf*Earth.radius
Moon.radius = sf*Moon.radius

#Physical Properties
G = 6.67e-11
Earth.mass = 5.972e24
Moon.mass = 7.347e22
Earth.v = vec(0,0,0)
Moon.v = vec(0,0,0)

#time
t = 0
dt = 60*5

scene.waitfor('click')

#Simulation Loop
while True:
    rate(1000)
    #Forces
    r = Moon.pos-Earth.pos
    Moon.f = -G*Earth.mass*Moon.mass/mag(r)**2*norm(r)
    Earth.f= -Moon.f
    #Time Integration
    Moon.v = Moon.v + Moon.f/Moon.mass*dt
    Earth.v = Earth.v + Earth.f/Earth.mass*dt
    Moon.pos = Moon.pos + Moon.v*dt
    Earth.pos = Earth.pos + Earth.v*dt
    t = t + dt
    #Collision Check
    if Earth.radius + Moon.radius >= mag(r):
        print("Collision!")
        print( t/60/60/24, "days")
        break

 

https://www.glowscript.org/#/user/emilyjiminroh/folder/physicscoding/program/%EB%8B%AC%EC%9D%B4%EB%A9%88%EC%B6%98%EB%8B%A4%EB%A9%B4? 

GlowScript IDE