Ejemplo n.º 1
0
double
LatticeTransportMaterial :: giveCharacteristicValue(MatResponseMode mode,
                                                    GaussPoint *gp,
                                                    TimeStep *tStep)
{
    LatticeTransportMaterialStatus *status = static_cast< LatticeTransportMaterialStatus * >( this->giveStatus(gp) );
    double suction = status->giveTempField().at(1);

    if ( mode == Capacity ) {
        return computeCapacity(suction, gp);
    } else if ( mode == Conductivity ) {
        return computeConductivity(suction, gp, tStep);
    } else {
        OOFEM_ERROR("unknown mode");
    }

    return 0; // to make compiler happy
}
Ejemplo n.º 2
0
Real
DiffTensorKernel::computeQpJacobian()
{
  RealTensorValue k = computeConductivity(_t, _qp);
  return k * _grad_test[_i][_qp] * _grad_phi[_j][_qp];
}
Ejemplo n.º 3
0
Real
DiffTensorKernel::computeQpResidual()
{
  RealTensorValue k = computeConductivity(_t, _qp);
  return k * _grad_test[_i][_qp] * _grad_u[_qp] - 4.0;
}