CEA 기후 제어 — 온도, 습도, 센서, 응답, 상호결합

두 가지 브라우저 실험

위 패널은 냉방 응답을 실행합니다. 액추에이터 시간 상수, 조명에 의한 현열, 냉방 용량을 변경할 수 있습니다. 0/300/900/1800초 프리셋은 다른 입력값을 발표된 기준값으로 유지합니다. 온도와 명령/전달된 냉방은 별도로 그려집니다. Python 설치는 필요하지 않으며, 아래의 로컬 재현 명령은 선택 사항입니다.

결합 제어(coupled-control) 패널은 모든 파라미터 변경에 대해 ON/OFF와 PI를 모두 실행합니다. 여섯 개의 개별 막대 그래프가 밴드 위반, 냉방/제습기 명령 중첩, 냉방열, 제거된 수분, 액추에이터 변동을 비교합니다. 각 그래프는 고유한 단위와 척도를 가집니다. A/B는 컨트롤러 결과 쌍 전체를 저장하므로, 컨트롤러 비교와 조건 변경이 서로 구분된 채 유지됩니다. 전체 지표, 수치 시계열, 다운로드, 공유 조건은 두 연구실 모두에서 이용할 수 있습니다.

냉방 응답은 원래의 12시간 모델, 10초 실내 오일러 스텝, 정확한 구간 평균 1차 액추에이터 업데이트를 유지합니다. 조명은 2시간과 8시간에 전환되며, 초기/설정 온도는 24°C, 외부 22°C, 기본 열량 800W, 열용량 12MJ/K, UA 300W/K, 고정 PI 게인은 0.6 및 0.0004입니다. 최종적으로 유지되는 명령을 포함하여 모든 원본 샘플이 그려집니다. 평가 밴드는 23–25°C로 유지됩니다. 오차 및 밴드 시간은 구간 시작값을 사용합니다. 극값에는 양 끝점이 포함됩니다. 이 모델은 수분과 센서 지연을 제외합니다.

결합 제어는 60초 명시적 스텝으로 24시간을 유지하며, 조명은 08–24시, 초기값은 24°C/RH70%, 300m³, 12MJ/K, UA 450W/K, 환기량 0.05m³/s입니다. 외부 온도는 21±6°C, RH60%이며, 조명 중/소등 시 증산량은 각각 2/0.35kg/h로 지정됩니다. ON/OFF 임계값은 23.5/24.5°C와 RH65/75%로 유지되며, PI 목표값은 24°C/RH70%로 원래의 고정 게인과 조건부 적분을 유지합니다. 위반 여부는 갱신된 상태값을 23–25°C 및 RH60–80%와 비교하여 판단합니다. 총 변동량은 절대 명령 변화의 합이며, 전환 횟수가 아닙니다. 원래의 중첩 지표는 실내 재가열이 0으로 설정되어 있더라도 양의 냉방 및 제습기 명령을 계수합니다.

결합 모델은 온도를 업데이트하기 전, 이전 온도에서의 포화 상태를 기준으로 수분을 제한합니다. 냉방 이후 응결을 재평형화하지 않으며 응결열을 반환하지도 않습니다. 냉방은 현열만을 다루며, 건조는 별도의 재가열을 동반하는 지정된 질량 흡수원입니다. 코일 모델, 액추에이터 지연, 센서 동역학, CO₂ 결합은 없습니다. 두 연구실은 서로 교환 가능한 설정이 아니라 별개의 합성 실험입니다. 열적 kWh는 전력이 아니며, 한계값/프리셋은 작물 설정값이나 장비 규격이 아닌 교육적 가정입니다. 원본 커널, 설정, ZIP 파일, 참조 결과는 보존됩니다.

cooling: experiment.py · cooling_kernel.py · config.json

coupled: experiment.py · coupled_kernel.py · config.json

CEA 온도 및 습도 제어 — 냉방, 건조, 환기가 상호작용하는 이유

온도와 RH는 서로 다른 목표값을 가질 수 있지만, 액추에이터가 하나의 값에만 영향을 미치는 경우는 드뭅니다. 냉각 코일은 현열을 제거하며 이슬점 이하에서 수분을 응결시킬 수 있습니다. 일부 제습기는 물을 모으는 동시에 현열을 실내로 되돌립니다. 환기는 외부 조건에 따라 온도, 수증기, CO₂를 함께 이동시킵니다. 따라서 독립된 온도 및 습도 스위치는 서로 충돌할 수 있습니다.

Diagram 1 · Use the button to switch views

다이어그램: Duskcoil. 특정 장비 배치와 무관한 책임 수준의 개략도입니다.

상태와 명령을 분리하십시오

열수지 연구실은 실내 온도 T를 저장하고, 수분수지 연구실은 수증기 질량 m_v를 저장합니다. RH는 별도의 보존 상태가 되지 않고 계속 유도값으로 남습니다.

미국 에너지부(DOE)의 HVAC 지침은 현열 제어와 잠열 제어를 구분하고 과냉각/재가열 상호작용을 설명합니다. 오클라호마 주립대학교 익스텐션 역시 온실 환기를 온도 관리와 습도 관리 모두의 일부로 다룹니다.

감독 체계는 우발적 충돌을 방지합니다

습도에 의한 냉방 뒤에 저온 난방이 이어지면 냉방과 난방이 동시에 작동할 수 있습니다. 이것이 자동으로 잘못된 것은 아닙니다. 의도적인 재가열은 습도를 낮출 수 있지만, 의도, 허용 지속시간, 에너지 비용을 통해 설계된 모드와 서로를 상쇄하는 컨트롤러를 구분해야 합니다.

최소한의 감독 체계는 우선순위를 명시적으로 만듭니다.

1. 센서 결함, 응결, 장비 보호 제약
2. 엄격한 온도 및 절대 수분 한계
3. 정상 목표 밴드 및 히스테리시스
4. 동시 작동을 허용하거나 금지하는 조건
5. 전력 한계 및 환기로 인한 CO₂ 손실과 같은 부작용

컨트롤러 비교를 위한 수용 조건

On/off, 규칙 기반, PI 비교는 기상 조건, 지정된 증산량, LED 일정, 초기 상태를 공유해야 합니다. 액추에이터 용량, 최소 가동/정지 시간, 지연, 포화, 전력 또는 수집된 물의 회계를 고정하십시오. 온도 및 RH 밴드를 벗어난 시간, 응결, 전환 횟수, 동시 난방/냉방 시간, 에너지를 각각 보고하십시오. 하나의 평균 오차만으로는 충분하지 않습니다.

센서를 소리 없이 클리핑하거나, 누락값을 선언 없이 유지하거나, 한쪽 컨트롤러에만 미래 기상 정보를 제공하는 실행은 거부하십시오. 근본적인 에너지 또는 수분 잔차 검사를 위반하는 실행도 거부하십시오.

이는 제어 설계의 출발점이며, 작물별 설정값이나 장비 권장 사항이 아닙니다. 다음 구현에서는 CO₂를 별도의 모니터링 대상 제약으로 유지하면서 열적 상태와 수분 상태만을 하나의 시뮬레이션으로 결합해야 합니다.

참고 자료

• Grid-interactive Efficient Buildings: HVAC (U.S. DOE)
• The Hobby Greenhouse (Oklahoma State University Extension)
• Assessing the Risk of Disease in Greenhouses (Penn State Extension)

온도 및 습도 센서 — 배치, 응답, 교정

컨트롤러는 어떤 온도를 보는가?

냉방 응답 연구실은 냉방 명령을 전달된 냉방과 분리했습니다. 그 컨트롤러는 여전히 현재 모델 온도를 직접 수신했습니다. 실제 시설로 나아가려면 물리적 환경, 센서 판독값, 제어로 전달되는 값을 별도로 구분해야 합니다.

24°C에서 26°C로의 이상적인 공기 온도 스텝과, 60초 1차 시간 상수로 모델링된 측정 시스템을 가정해 보십시오. 30초 후, 판독값은 24.787°C입니다. 이 판독값을 25°C와 비교하면 입력값에 이미 존재하는 초과분을 놓치게 됩니다. 이는 설명을 위한 해석적 계산이며, 시설 데이터나 제품 사양이 아닙니다.

이 글은 측정 설계 가이드입니다. 이전의 폐루프 연구실을 다시 실행하거나 제품 순위를 매기지 않으며, 예시로 사용된 임계값은 작물 권장 사항이 아닙니다.

각 측정 위치에 목적을 부여하십시오

재배 구역, 공급 공기, 반환 공기 판독값을 단순히 평균하면 변화가 발생한 위치에 대한 정보가 사라집니다. 집계 규칙을 선택하기 전에 관심 대상 수량과 위치를 정의하십시오. 다음의 조사 계획은 보편적인 센서 개수 요건이 아니라 편집상의 제안입니다.

Sensirion의 설계 가이드는 대표성 있는 배치, 인클로저 공기 교환, 열적 영향을 시스템 수준의 고려 사항으로 규정합니다. 센서 단독의 정확도가 설치된 상태의 정확도를 보장하지는 않습니다.

조사를 위해서는 조명, 팬 작동, 작물 생장 단계를 기록한 뒤 해당 운영 조건에서 위치들을 비교하십시오. 휴대용 센서가 위치를 이동하는 경우, 고정된 기준점을 유지하십시오. 그렇지 않으면 조사 중 실내 전체의 변화가 공간적 차이처럼 보일 수 있습니다. 두 계측기를 위치 간에 교환하는 것도 도움이 될 수 있지만, 교환 중 조건이 변할 수 있습니다. 교환과 공동 배치를 함께 사용하여 원인을 좁히십시오.

응답 시간과 기록 간격을 분리하십시오

1초마다 값을 읽는다고 해서 1초의 물리적 응답이 성립하는 것은 아닙니다. 온도와 RH도 서로 다르게 반응할 수 있습니다. 정의와 설치된 구성을 확인하십시오. 설계 가이드는 스텝 변화의 약 63%에 도달하는 시간으로 응답을 논의합니다.

이 예시에서만 목표 온도를 T, 측정 온도를 T_m, 측정 시간 상수를 \tau_s라고 합시다.

24°C에서의 초기 평형, 26°C의 일정한 입력, \tau_s=60초일 때, 해는 T_m(t)=26-2\exp(-t/60)입니다.

판독값은 약 41.6초 후에 예시로 사용된 25°C 임계값에 도달합니다. 이산 기록은 감지를 다음 샘플까지 지연시킬 수 있습니다. 이 모델은 즉시 응답을 시작하며, 아무 응답도 없는 순수한 데드타임과는 다릅니다.

다음 표준 라이브러리 Python 코드를 sensor_example.py로 저장하고 python3 sensor_example.py를 실행하여 표와 임계값 교차 시간을 재현하십시오. 외부 파일을 읽거나 쓰지 않습니다.

from math import exp, log

initial, final, tau = 24.0, 26.0, 60.0
for seconds in (30, 60, 180):
    measured = final + (initial - final) * exp(-seconds / tau)
    print(f"{seconds:3d} s: {measured:.3f} C, gap {final-measured:.3f} C")
threshold = 25.0
crossing = -tau * log((final - threshold) / (final - initial))
print(f"threshold crossing: {crossing:.1f} s")

실제 실내는 스텝 형태로 변하지 않을 수 있으며, 설치된 측정 시스템도 단일 시간 상수를 갖지 않을 수 있습니다. 로그로부터 응답을 조사할 때는 기준 계측기의 동역학과 타임스탬프 정렬을 고려하십시오. 냉방 명령과 변화된 판독값 사이의 간격에는 액추에이터, 공기 이동, 측정 효과가 결합되어 있으므로, 그 간격만으로는 각 기여도를 식별할 수 없습니다.

필터링 전후의 값을 보존하십시오

제어에서 이동평균이나 다른 필터를 적용하는 경우, 획득된 계측기 출력값과 처리된 제어 입력값을 모두 기록하십시오. 획득된 출력값에는 이미 내부 처리가 포함되어 있을 수 있습니다. 필터 설정 역시 보존하십시오.

이전 연구실을 확장하려면, 평가를 위해 모델의 실제 온도 T를 유지하면서 PI 오차로 T_m-T_{sp}를 사용하십시오. 밴드를 벗어난 시간을 측정 온도와 실제 온도에 대해 각각 계산하십시오. 여기서의 계산에는 피드백 루프가 없으므로, 결과적인 진동이나 추가적인 오버슈트에 대해 어떠한 수치적 주장도 하지 않습니다.

시설에서는 실제 온도를 직접 얻을 수 없습니다. 독립적인 기준 측정값에도 불확실성이 있으므로, 비교 결과를 보고할 때는 기준 계측기, 위치, 교정 정보를 문서화하십시오.

온도와 RH를 한 쌍으로 취급하십시오

상대습도는 해당 온도에서의 포화 증기압에 대한 증기압의 비율을 나타냅니다. VPD 계산에는 동일한 공기를 나타내는 온도와 RH를 사용하고 그 타임스탬프를 정렬하십시오. 선반 온도와 공급 공기 RH를 결합한다고 해서 선반의 VPD를 측정하는 것은 아닙니다.

가열식 습도 프로브는 가열되는 감지 영역과 주변 공기를 구분해야 합니다. Vaisala의 가열 프로브 노트는 주변 RH를 얻기 위해 별도의 주변 온도 측정값을 사용하는 방법을 설명합니다. 모든 가열식 프로브 출력을 재배 구역의 온도 및 RH로 취급하지 말고, 선택한 계측기의 출력 사양을 따르십시오.

VPD 연구실에 값을 제공할 때는, 공기 온도 기반 계산과 잎 온도를 사용한 잎-공기 계산을 구분하십시오. 잎 온도 측정값이 없다면 공기 기반 결과에 측정된 잎 상태라는 라벨을 붙이지 마십시오.

교정, 조정, 배치 조사를 분리하여 유지하십시오

Vaisala의 교정 노트는 기준값과의 비교를 계측기 조정과 구분합니다. 한 지점에서의 일치가 전체 운용 범위를 검증하지는 않습니다. 운용 범위와 기준 불확실성을 염두에 두고 비교 지점과 허용 한계를 선택하십시오.

계측기 ID, 조정 전 판독값, 기준 판독값, 안정화 조건, 날짜를 기록하십시오. 조정이 이루어지면 설정 변경 사항과 이후의 비교 결과를 추가하십시오. 원래의 불일치를 폐기하면 이전 재배 기록을 평가하기 어려워집니다.

두 계측기가 일치한다고 해서 공통된 편향이 배제되는 것은 아닙니다. 교정 역시 설치된 위치가 재배 구역을 대표한다는 것을 입증하지 않습니다. 계측기 비교와 위치 조사를 별개의 작업으로 완료하십시오.

제어로 전달되는 기록을 정의하십시오

다음 필드는 제안된 기록 설계이며, 측정된 데이터셋이 아닙니다.

통신 장애 중 마지막 판독값을 유지하면 겉보기에 안심되는 평탄한 표시가 만들어질 수 있습니다. 값과 함께 샘플 경과 시간과 유효성을 전달하십시오. 실제 장비와 작물에 대해 오래된 데이터의 한계와 그에 따른 제어 조치를 정의한 뒤, 커미셔닝 중 그 동작을 검증하십시오.

다음 실험 전에, 공동 배치를 통한 계측기 차이, 위치 조사를 통한 공간적 차이, 획득에서 제어 입력까지의 지연을 조사하십시오. 이러한 기록들이 함께 있으면 컨트롤러를 변경할지 측정 프로세스를 개선할지 결정하는 데 도움이 됩니다.

냉방 응답 연구실 — PI 명령을 전달된 냉방과 분리

동일한 냉방 용량이 동일한 응답을 의미하지는 않습니다

이전의 결합 제어 비교는 명령의 변화가 즉시 냉방 및 건조 용량을 변화시킨다고 가정했습니다. 이 실험은 현열 냉방만을 분리하고 전달된 용량이 명령을 점진적으로 따라가도록 합니다. 이는 측정된 재배 시설이 아니라 합성 실내 모델입니다.

동일한 12시간 동안, 액추에이터 시간 상수를 0에서 900초로 늘리면 적분된 절대 온도 오차는 2.287에서 4.351 K·h로, 최저 온도는 23.419에서 23.044°C로 변합니다. 두 경우 모두 선택된 23–25°C 밴드 내에 머무릅니다. 따라서 이진법적인 합격/불합격 밴드 검사만으로는 응답의 상당한 차이가 가려집니다.

Diagram 2 · Use the button to switch views

그림: Duskcoil, 배포된 Python 시뮬레이션으로부터 생성. 빨간색 1800초 궤적은 스트레스 사례입니다. 밴드는 작물 권장 사항이 아니라 실험적 가정입니다.

의도적으로 축소된 실험

이는 이전의 결합 모델을 한 개의 파라미터만 바꿔 재실행한 것이 아니라, 액추에이터 응답을 분리하기 위한 별도의 실험입니다. 습도, 증산, 제습기 재가열, CO₂를 제외합니다. 모든 경우에서 외부 공기 22°C, 배경 열량 800W, 그리고 2시간에서 8시간 사이 조명에 의한 추가적인 현열 8,000W를 사용합니다. 이 8,000W는 조명의 전력 수요 추정치가 아닙니다.

DOE의 현열/잠열 분리 프로젝트는 온도와 수분 부하를 별도로 처리하는 기술을 설명합니다. 이 실내 모델은 그 장비를 재현하지 않습니다. 수분 제어와 전력에는 추가적인 모델과 평가가 필요합니다.

명령과 전달된 출력값에 별도의 변수를 부여하십시오

설정값보다 높은 온도가 냉방을 증가시키도록 오차를 정의합니다. 0–1로 고정된 PI 출력이 명령 u입니다. 최대 용량 대비 비율로 전달된 냉방은 a입니다.

이는 1차 지연입니다. 일정한 명령 하에서, 하나의 시간 상수는 요청된 변화의 약 63.2%를 담당합니다. 900초의 시간 상수가 900초 동안 아무 일도 일어나지 않는다는 뜻은 아닙니다. 그것은 데드타임에 해당할 것입니다. 전달된 출력값은 0에서 시작하며, 시간 상수가 0인 경우에만 명령을 즉시 따릅니다.

실내는 현열 수지를 사용합니다.

코드는 각 10초 구간마다 명령을 고정하고, 액추에이터의 종점값과 구간 평균값을 해석적으로 계산합니다. 실내 온도는 구간 평균 냉방을 사용한 전진 오일러법을 사용합니다. 이는 종점 냉방률이 구간 전체에 걸쳐 전달된 것처럼 계산되는 것을 방지합니다.

PI 컨트롤러는 조건부 적분을 사용합니다. 제약이 없는 명령이 상한을 초과하고 오차가 포화 상태를 더 심화시킬 경우 적분을 정지하지만, 포화 상태를 해제하는 적분은 허용합니다. MathWorks는 이 방법을 실제 액추에이터 출력을 추적하는 것과 구분합니다. 여기서는 조건부 적분만 구현되어 있으며, 그 자체로 액추에이터 동역학을 보상하지는 않습니다.

동일한 구간에서 지표를 비교하십시오

적분된 절대 오차는 |T-T_{sp}|를 누적하므로, 따뜻한 방향과 차가운 방향의 이탈이 서로 상쇄되지 않습니다. 밴드 지속시간과 오차는 각 구간 시작값의 온도를 사용합니다. 즉 10초 단위의 근사값입니다. 최저값과 최고값에는 초기 및 최종 상태를 포함한 모든 격자 끝점이 포함됩니다.

조명이 꺼지는 시점에 현열 입력은 즉시 감소하지만 전달된 냉방은 점진적으로 감쇠합니다. 명령이 감소한 이후 하단 패널에 남아 있는 출력값을 상단 패널의 저온 이탈과 대조해 보십시오. 냉방열의 증가를 전력이나 비용의 증가로 직접 해석할 수는 없습니다. 냉방 COP, 펌프, 팬 모델이 없기 때문입니다.

실행한 뒤 조건 하나만 변경하십시오

재현 가능한 ZIP 파일을 빈 디렉토리에 압축 해제하십시오. 계산과 테스트에는 Python 표준 라이브러리만 필요하며, 검증은 Python 3.12로 수행했습니다.

python3 reproduce.py
python3 -m unittest -v
python3 experiment.py config.json --tau 900 --csv trace-local.csv
python3 experiment.py config.json --tau 1800

reproduce.py는 배포된 파일 해시를 검증하고 저장된 0초, 300초, 900초 사례를 다시 계산합니다. CSV 내보내기는 기존 파일을 덮어쓰지 않습니다. 마지막 CSV 행은 이전에 유지된 명령값을 담고 있으며, 종점에서의 새로운 제어 업데이트를 나타내지 않습니다.

개별 다운로드: 설정, 계산 소스, 참조 결과, 1800초 결과, 900초 CSV, 매니페스트.

먼저 게인을 변경하지 않고 1800초를 시도해 보십시오. 밴드를 벗어난 시간은 3.650시간입니다. 이는 선택된 수용 기준의 실패를 의미할 뿐, 발산이나 입증된 장비 결함을 의미하지 않습니다. 그런 다음 900초로 돌아가 dt_s를 5로 변경하여 이산화 효과를 확인하십시오. 모든 조건을 한꺼번에 변경하면 개선의 원인이 가려집니다.

수치적 정확성과 시설 적합성은 서로 다른 검사입니다

10개의 테스트는 해석적 1차 응답, 구간 세분화, 액추에이터 경계값, 적분 정지/해제, 일정 열량에서의 알려진 해, 열수지 잔차, 참조 출력값, 잘못된 입력값, 절반 스텝 수렴, 스트레스 사례의 밴드 위반을 다룹니다. 세 가지 표준 사례 전반에서, 스텝을 10초에서 5초로 절반으로 줄이면 극값, 절대 오차, 냉방열은 각각 0.03°C, 0.03 K·h, 0.03kWh 미만으로 변합니다. 수치적 수렴이 합성 파라미터가 실제 시설을 대표한다는 것을 입증하지는 않습니다.

실제 설치의 경우, 전달된 냉방을 추정할 수 있는 명령 이력과 측정값을 수집한 뒤 식별용 데이터와 검증용 데이터를 분리하십시오. 센서 지연, 최소 가동 시간, 단계별 작동, 데드타임, 습도 결합은 추가적인 확장 사항입니다. 작물 설정값, 질병 위험, 에너지 절감 비율은 이 실험의 산출물이 아닙니다.

계속 읽기

온도 및 습도 센서는 판독값을 제어로 전달하기 전의 배치, 측정 응답, 교정을 다룹니다.

결합 제어 연구실 — 동일한 교란 조건에서 on/off와 PI 비교

결과부터 시작

동일한 합성 24시간 동안, PI는 총 액추에이터 변동량을 34.0에서 2.853으로, 냉방 제거량을 182.233에서 178.854kWh로 줄입니다. 그러나 온도 위반 시간은 7.95에서 10.02시간으로 증가하고, 동시 냉방 및 제습기 재가열 시간은 9.37에서 14.65시간으로 증가합니다. 이 실행이 PI의 우월성을 입증하지는 않습니다.

Diagram 3 · Use the button to switch views

다이어그램: Duskcoil. 컨트롤러 비교를 위한 공유 시스템 경계.

공유 조건

열적 수지와 수분 수지는 실내 온도 T와 수증기 질량 m_v라는 단 두 개의 상태로 결합됩니다. 기상 조건, LED, 지정된 증산량, 초기 상태, 60초 스텝, 14kW 냉방, 2.5kg/h 건조는 동일합니다. 합성 제습기는 실내로 1.8kW의 현열을 되돌립니다. CO₂는 별도로 유지됩니다.

히스테리시스 스위치는 24±0.5°C와 65/75% RH를 기준으로 작동합니다. PI 출력은 0–1로 포화되며 조건부 안티와인드업을 사용합니다. 게인은 조정된 장비 설정값이 아니라 합성값입니다.

DOE 지침이 설명하듯, 현열 및 잠열 제어에는 과냉각과 재가열이 수반될 수 있습니다. 위반이나 동시 작동이 증가하는 경우, 낮은 에너지 사용량만으로는 설계를 결정할 수 없습니다.

실행 및 검증

ZIP 파일을 압축 해제하십시오.

python3 experiment.py config.json > result.json
python3 -m unittest -v

config.json, expected.json, experiment.py를 참조하십시오. 여덟 개의 테스트가 두 컨트롤러, 결정론성, 에너지/수분 잔차, 잘못된 입력값을 다룹니다.

게인을 조정하기 전에, 수용 밴드, 지표 가중치, 동시 작동 한계, 전력 환산, 액추에이터 지연을 고정하십시오. 실제 배치에서는 식별용 데이터와 검증용 데이터를 분리해야 합니다. 이 글은 작물 설정값, 실제 장비의 절감 효과, PI 우월성에 대한 보장을 제공하지 않습니다.

냉방 응답 후속편은 명령을 전달된 냉방과 분리하고 액추에이터 시간 상수만을 변화시킵니다.

참고 자료

• Grid-interactive Efficient Buildings: HVAC (U.S. DOE)
• The Hobby Greenhouse (Oklahoma State University Extension)