C++ (Cpp) upwind示例

示例#1

文件： UpwindFactory.cpp 项目： fogelson/CFDResearch

void UpwindFactory::produce(Grid * g){
	if(contains(g)){
		return;
	}
	Divergence div(g);
	CellToFaceOperator F;
	F.setGrid(g);
	vector<Face*>::iterator it;
	for(it = g->faces.begin(); it != g->faces.end(); it++){
		if((*it) != 0 && (*it)->isBoundary()){
			FaceIndex faceIndex = (*it)->getIndex();
			Coord pos = (*it)->getCentroid();
			TinyVector<double,2> aSpeed = (*c)(pos);
			double a = sum(aSpeed*((*it)->getNormal()));
			if(a > 0){
				CellIndex upwind((*it)->getFrom()->getI(),(*it)->getFrom()->getJ());
				F.coefficients[faceIndex][upwind] = -a;
			}
			if(a < 0){
				F.constantTerm(faceIndex) = -a;
			}
		}
		else if((*it) != 0 && (*it)->isInterior()){
			FaceIndex faceIndex = (*it)->getIndex();
			Coord pos = (*it)->getCentroid();
			TinyVector<double,2> aSpeed = (*c)(pos);
			double a = sum(aSpeed*((*it)->getNormal()));
			if(a > 0){
				CellIndex upwind((*it)->getFrom()->getI(),(*it)->getFrom()->getJ());
				F.coefficients[faceIndex][upwind] = -a;
			}
			else if(a < 0){
				//cout << "Negative speed through face" << endl;
				CellIndex upwind((*it)->getTo()->getI(),(*it)->getTo()->getJ());
				F.coefficients[faceIndex][upwind] = -a;
			}
		}
	}
	operators[g] = new CellToCellOperator(div(F));
}

示例#2

文件： Q2PPorepressureFlux.C 项目： aeslaughter/moose

void
Q2PPorepressureFlux::computeOffDiagJacobian(unsigned int jvar)
{
  upwind(false, true, jvar);
}

示例#3

文件： RichardsFullyUpwindFlux.C 项目： aeslaughter/moose

void
RichardsFullyUpwindFlux::computeOffDiagJacobian(unsigned int jvar)
{
  upwind(false, true, jvar);
  return;
}

示例#4

文件： RichardsFullyUpwindFlux.C 项目： aeslaughter/moose

void
RichardsFullyUpwindFlux::computeResidual()
{
  upwind(true, false, 0);
  return;
}

示例#5

文件： Q2PSaturationFlux.C 项目： Biyss/moose

void
Q2PSaturationFlux::computeJacobian()
{
  upwind(false, true, _var.number());
}

示例#6

文件： Q2PSaturationFlux.C 项目： Biyss/moose

void
Q2PSaturationFlux::computeResidual()
{
  upwind(true, false, 0);
}

示例#7

文件： Q2PSaturationFlux.C 项目： Biyss/moose

void
Q2PSaturationFlux::computeOffDiagJacobian(unsigned int jvar)
{
  upwind(false, true, jvar);
}

示例#8

文件： Q2PSaturationFlux.C 项目： idaholab/moose

void
Q2PSaturationFlux::computeOffDiagJacobian(MooseVariableFEBase & jvar)
{
  upwind(false, true, jvar.number());
}

示例#9

文件： PorousFlowDarcyBase.C 项目： mangerij/moose

void
PorousFlowDarcyBase::computeResidual()
{
  upwind(JacRes::CALCULATE_RESIDUAL, 0.0);
}

示例#10

文件： PorousFlowDarcyBase.C 项目： mangerij/moose

void
PorousFlowDarcyBase::computeOffDiagJacobian(unsigned int jvar)
{
  upwind(JacRes::CALCULATE_JACOBIAN, jvar);
}

示例#11

文件： PorousFlowAdvectiveFlux.C 项目： Biyss/moose

void
PorousFlowAdvectiveFlux::computeOffDiagJacobian(unsigned int jvar)
{
   upwind(CALCULATE_JACOBIAN, jvar);
}

示例#12

文件： PorousFlowAdvectiveFlux.C 项目： Biyss/moose

void
PorousFlowAdvectiveFlux::computeJacobian()
{
   upwind(CALCULATE_JACOBIAN, _var.number());
}

示例#13

文件： PorousFlowAdvectiveFlux.C 项目： Biyss/moose

void
PorousFlowAdvectiveFlux::computeResidual()
{
  upwind(CALCULATE_RESIDUAL, 0.0);
}

示例#14

文件： RichardsFullyUpwindFlux.C 项目： AhmedAly83/moose

void
RichardsFullyUpwindFlux::computeJacobian()
{
  upwind(false, true, _var.number());
  return;
}

示例#15

文件： Q2PPorepressureFlux.C 项目： aeslaughter/moose

void
Q2PPorepressureFlux::computeResidual()
{
  upwind(true, false, 0);
}

示例#16

文件： Q2PPorepressureFlux.C 项目： aeslaughter/moose

void
Q2PPorepressureFlux::computeJacobian()
{
  upwind(false, true, _var.number());
}

示例#17

文件： PorousFlowDarcyBase.C 项目： mangerij/moose

void
PorousFlowDarcyBase::computeJacobian()
{
  upwind(JacRes::CALCULATE_JACOBIAN, _var.number());
}