Python 프로그래밍/Chemistry

3. 열역학, scipy로 적분 하기

palefaceman 2021. 3. 24. 18:19

1. 적분하기

예시 1

quad(적분 함수(), start, end, 적분함수 인자)

quad 함수는 값과, 오차값을 배열(array)로 반환한다.

1.6667은 적분 값, 1.85e-14는 오차값

예시 2

예시 3 가우스 함수(종 함수) 구하기

a 값은 종 함수의 width로 크기가 클수록 넓어지고, 크기가 작으면 좁아진다.

2. PV-그래프(PV-diagram)

PV-그래프를 그린다.

아래에 면적은 일(work)이다. PV = W [J] 

2.1 등온 과정(isothermal process)

2.1.1. 적분해서 면적(일)을 구한다.

부피 10~40 사이에 등온(isothermal) 팽창한 기체의 200,300,400K에서 시스템 밖으로 한 일을 구한다.

2.1.2. 열역학 1법칙으로 열(heat)를 구한다.

W= Q에서 열(heat)를 구한다.

2.2. 등온(isothermal)과 단열(adiabatic)을 비교해본다.

2.2.1.  PV그래프를 그린다.

2.2.2. 적분해서 일(work)을 구한다.

2.2.3. 단열(adiabtic) 과정에 T2와 P2를 구한다.