ISAT models

Buildings/Resources/C-Sources/cfdExchangeData.c

@@ -7,8 +7,11 @@
  * \author Wangda Zuo
  *         University of Miami
  *         [email protected]
+ *         Xu Han
+ *         University of Colorado Boulder
+ *         [email protected] 
  *
- * \date   8/3/2013
+ * \date   4/5/2020
  *
  */
 #include "cfdCosimulation.h"
@@ -33,7 +36,53 @@ int cfdExchangeData(double t0, double dt, double *u, size_t nU, size_t nY,
 
   /*check if current modelica time equals to last time*/
   /*if yes, it means cfdExchangeData() was called multiple times at one synchronization point, then directly return*/
-  if(abs(cosim->modelica->lt - t0) < 1E-6){
+  if(cosim->modelica->lt - t0 < 1E-6 && t0 - cosim->modelica->lt < 1E-6){
+ /****************************************************************************
+  | Copy data from CFD
+  ****************************************************************************/
+  /* Get the temperature/heat flux for solid surface*/
+  for(i=0; i<cosim->para->nSur; i++) {
+    y[i] = cosim->ffd->temHea[i];
+  }
+
+  /* Get the averaged room temperature*/
+  y[i] = cosim->ffd->TRoo;
+  i++;
+
+  /* Get the temperature of shading device if there is a shading device*/
+  if(cosim->para->sha==1) {
+    for(j=0; j<cosim->para->nConExtWin; i++, j++) {
+      y[i] = cosim->ffd->TSha[j];
+    }
+  }
+
+  /* Get the temperature fluid at the fluid ports*/
+  for(j=0; j<cosim->para->nPorts; i++, j++) {
+    y[i] = cosim->ffd->TPor[j];
+  }
+
+  /* Get the mass fraction at fluid ports*/
+  for(j=0; j<cosim->para->nPorts; j++)
+    for(k=0; k<cosim->para->nXi; k++, i++) {
+       y[i] = cosim->ffd->XiPor[j][k];
+    }
+
+  /* Get the trace substance at fluid ports*/
+  for(j=0; j<cosim->para->nPorts; j++)
+    for(k=0; k<cosim->para->nC; k++, i++) {
+       y[i] = cosim->ffd->CPor[j][k];
+    }
+
+  /* Get the sensor data*/
+  for(j=0; j<cosim->para->nSen; j++, i++) {
+    y[i] = cosim->ffd->senVal[j];
+  }
+
+  /* Update the data status
+  cosim->ffd->flag = 0;*/
+
+  *t1 = cosim->ffd->t;
+
     return 0;
   }
 
@@ -75,6 +124,13 @@ int cfdExchangeData(double t0, double dt, double *u, size_t nU, size_t nY,
   cosim->modelica->latentHeat = u[i];
   i++;
 
+  if(cosim->para->nSou>0){
+    for(j=0; j<cosim->para->nSou; j++) {
+      cosim->modelica->sourceHeat[j] = u[i+j];
+    }
+  }
+  i = i + cosim->para->nSou;
+
   cosim->modelica->p = u[i];
   i++;
 

Buildings/Resources/C-Sources/cfdSendStopCommand.c

@@ -21,7 +21,7 @@
 void cfdSendStopCommand(void *thread) {
 
   size_t i = 0;
-  size_t imax = 10000;
+  size_t imax = 1000000;
 
   /*send stop command to FFD*/
   cosim->para->flag = 0;
@@ -133,6 +133,16 @@ void cfdSendStopCommand(void *thread) {
   if (cosim->ffd->temHea != NULL){
     free(cosim->ffd->temHea);
   }
+
+  if(cosim->para->nSou>0){
+    if (cosim->para->souName != NULL){
+      free(cosim->para->souName);
+    }
+    if (cosim->modelica->sourceHeat != NULL){
+      free(cosim->modelica->sourceHeat);
+    }
+  }
+
   if (cosim->para != NULL){
     free(cosim->para);
   }

Buildings/Resources/C-Sources/cfdStartCosimulation.c

@@ -44,7 +44,7 @@
 int cfdStartCosimulation(char *cfdFilNam, char **name, double *A, double *til,
                 int *bouCon, int nPorts, char** portName, int haveSensor,
                 char **sensorName, int haveShade, size_t nSur, size_t nSen,
-                size_t nConExtWin, size_t nXi, size_t nC, double rho_start) {
+                size_t nConExtWin, size_t nXi, size_t nC, int haveSource, size_t nSou, char **sourceName, double rho_start) {
   size_t i;
   size_t nBou;
 
@@ -66,6 +66,18 @@ int cfdStartCosimulation(char *cfdFilNam, char **name, double *A, double *til,
   cosim->para->nXi = nXi;
   cosim->para->rho_start = rho_start;
 
+  cosim->para->Sou = haveSource;	
+  cosim->para->nSou = nSou;
+
+  if(cosim->para->nSou>0){
+    cosim->para->souName = (char**) malloc(nSou*sizeof(char *));
+    cosim->modelica->sourceHeat = (REAL *) malloc(nSou*sizeof(REAL));
+    for(i=0; i<nSou; i++) {
+      cosim->para->souName[i] = (char *)malloc(sizeof(char)*(strlen(sourceName[i])+1));
+      strcpy(cosim->para->souName[i], sourceName[i]);
+    }
+  }
+
   nBou = nSur + nPorts;
 
   cosim->para->name = (char**) malloc(nSur*sizeof(char *));

Buildings/Resources/C-Sources/cfdcosim.c

@@ -96,6 +96,10 @@ void *cfdcosim() {
   cosim->modelica->shaConSig = NULL;
   cosim->modelica->shaAbsRad = NULL;
   cosim->ffd->TSha = NULL;
+  cosim->para->Sou = 0;	
+  cosim->para->nSou = 0;
+  cosim->para->souName = NULL;
+	cosim->modelica->sourceHeat = NULL;	
 
   return (void*) cosim;
 } /* End of cfdcosim()*/

Buildings/Resources/C-Sources/isatCosimulation.h

@@ -0,0 +1,42 @@
+/*
+ *
+ * @file   isatCosimulation.h
+ *
+ * @brief  Header file for coupled simulation with ISAT based on cfdCosimulation.h
+ *
+ * @author
+ *         Xu Han
+ *         University of Colorado Boulder
+ *         [email protected]
+ *
+ *         Wangda Zuo
+ *         University of Miami
+ *         [email protected] 
+ *
+ * \date   9/21/2019
+ *
+ */
+
+#if defined(_MSC_VER) || defined(__WIN32__) /* Windows */
+#include <windows.h>
+#else /* Linux*/
+#include <dlfcn.h>  /*For load shared library*/
+#include <unistd.h> /*For Linux function*/
+#define Sleep(x) sleep(x/1000)
+#endif
+
+#ifndef _MODELICA_FFD_COMMON_H
+#define _MODELICA_FFD_COMMON_H
+#include "../src/FastFluidDynamics/modelica_ffd_common.h"
+#include <stdint.h> /* Needed to detect 32 vs. 64 bit using UINTPTR_MAX*/
+#endif
+
+#if UINTPTR_MAX == 0xffffffff
+/* 32-bit */
+#error "*** ISAT is not yet on Windows 32 bit."
+#endif
+
+CosimulationData *cosim;
+
+/*declare the isat_dll function in DLL*/
+void *isat_dll(CosimulationData *cosim);