Ejemplo n.º 1
0
//Methods
void CFDRunClass::Run()
{
SetInitialConditions();
bool velcheck;
UpdateFluidFlag();
SetPreliminaryGuessValues();
float t=0;
int m=0;
do
{
 	RunNavierStokes();
	velcheck=UpdateVelocity();
	RunFreeSurfaceCalc();
	t+=dt;
	m+=1;
	vector<float> timevec;
	timevec.push_back(0.02);timevec.push_back(0.05);timevec.push_back(0.1);timevec.push_back(0.2);timevec.push_back(0.3);timevec.push_back(0.4);timevec.push_back(0.5);timevec.push_back(0.6);timevec.push_back(0.7);timevec.push_back(0.8);timevec.push_back(0.9);timevec.push_back(1.0);
	//timevec.push_back(0.01);timevec.push_back(0.02);timevec.push_back(0.05);timevec.push_back(0.1);timevec.push_back(0.2);timevec.push_back(0.3);timevec.push_back(0.4);
	for(int mit=0;mit<=11;mit++)
	{
	if(t==timevec[mit] || ((t-dt)<timevec[mit] && t>timevec[mit]))
	{
	PrintData(m,t);
	PrintFreeSurfaceData(m);
	}
	}
}while((t<tmax)&&(velcheck==true));
PrintData(m,t);
}
Ejemplo n.º 2
0
Archivo: lb.hpp Proyecto: UCL/hemelb
    void LBM<LatticeType>::Initialise(vis::Control* iControl,
                                      iolets::BoundaryValues* iInletValues,
                                      iolets::BoundaryValues* iOutletValues,
                                      const util::UnitConverter* iUnits)
    {
      mInletValues = iInletValues;
      mOutletValues = iOutletValues;
      mUnits = iUnits;

      InitCollisions();

      SetInitialConditions();

      mVisControl = iControl;
    }
Ejemplo n.º 3
0
int main()
{
    InitRCC();
    InitVTimers();
    InitControllerPeripheral();
    SetInitialConditions();
    InitLCD();
    
    while(1)
    {
        ControllerDisplayDataTask();
    }
    
    return 0;
}
Ejemplo n.º 4
0
void InitialConditions(struct DATA *params, struct GRIDINFO *grid, struct FIELDCONTAINER *field){
	
	grid->SetTime(0,grid);

	for(int i = grid->imin; i < grid->imax; i++){
		grid->GetPos(i,grid,0);	
		for(int j = grid->jmin; j < grid->jmax; j++){
			grid->GetPos(j,grid,1);
			for(int k = grid->kmin; k < grid->kmax; k++){
				grid->GetPos(k,grid,2);
								
				SetInitialConditions(params,grid,field);
				
			} // END k-loop
		} // END j-loop
	} // END i-loop

} // END InitialConditions
Ejemplo n.º 5
0
void GremlinsWorker::Setup()
{
  Register_Stream(stGremlins);

  SetInitialConditions();

  Sweeps = 0;

  if (!Online_GetGremlinsLoop()) 
  {
    LastSweep = 1;

    if (Online_GetGremlinsRoundtrip())
    {
      LastSweep ++;
    }
  }
  else
  {
    LastSweep = -1;
  }
}
Ejemplo n.º 6
0
void APE::v_DoInitialise()
{
    SetInitialConditions();
}