NSInitialCondition.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
// Navier-Stokes includes
#include "NS.h"
#include "NSInitialCondition.h"
 
// FluidProperties includes
#include "IdealGasFluidProperties.h"
 
// MOOSE includes
#include "MooseVariable.h"
 
registerMooseObject("NavierStokesApp", NSInitialCondition);
 
template <>
InputParameters
validParams<NSInitialCondition>()
{
  InputParameters params = validParams<InitialCondition>();
  params.addClassDescription("NSInitialCondition sets intial constant values for all variables.");
  params.addRequiredParam<Real>("initial_pressure",
                                "The initial pressure, assumed constant everywhere");
  params.addRequiredParam<Real>("initial_temperature",
                                "The initial temperature, assumed constant everywhere");
  params.addRequiredParam<RealVectorValue>("initial_velocity",
                                           "The initial velocity, assumed constant everywhere");
  params.addRequiredParam<UserObjectName>("fluid_properties",
                                          "The name of the user object for fluid properties");
 
  return params;
}
 
NSInitialCondition::NSInitialCondition(const InputParameters & parameters)
  : InitialCondition(parameters),
    _initial_pressure(getParam<Real>("initial_pressure")),
    _initial_temperature(getParam<Real>("initial_temperature")),
    _initial_velocity(getParam<RealVectorValue>("initial_velocity")),
    _fp(getUserObject<IdealGasFluidProperties>("fluid_properties"))
{
}
 
Real
NSInitialCondition::value(const Point & /*p*/)
{
  const Real rho_initial = _fp.rho_from_p_T(_initial_pressure, _initial_temperature);
 
  // TODO: The internal energy could be computed by the IdealGasFluidProperties.
  const Real e_initial = _fp.cv() * _initial_temperature;
  const Real E_initial = e_initial + 0.5 * _initial_velocity.norm_sq();
  const Real v_initial = 1. / rho_initial;
 
  if (_var.name() == NS::enthalpy)
    return E_initial + _initial_pressure / rho_initial;
 
  if (_var.name() == NS::internal_energy)
    return e_initial;
 
  if (_var.name() == NS::mach_number)
    return _initial_velocity.norm() / _fp.c_from_v_e(v_initial, e_initial);
 
  if (_var.name() == NS::pressure)
    return _initial_pressure;
 
  if (_var.name() == NS::density)
    return rho_initial;
 
  if (_var.name() == NS::momentum_x)
    return rho_initial * _initial_velocity(0);
 
  if (_var.name() == NS::momentum_y)
    return rho_initial * _initial_velocity(1);
 
  if (_var.name() == NS::momentum_z)
    return rho_initial * _initial_velocity(2);
 
  if (_var.name() == NS::total_energy)
    return rho_initial * E_initial;
 
  if (_var.name() == NS::specific_volume)
    return v_initial;
 
  if (_var.name() == NS::temperature)
    return _initial_temperature;
 
  if (_var.name() == NS::velocity_x)
    return _initial_velocity(0);
 
  if (_var.name() == NS::velocity_y)
    return _initial_velocity(1);
 
  if (_var.name() == NS::velocity_z)
    return _initial_velocity(2);
 
  // If we got here, then the variable name was not one of the ones we know about.
  mooseError("Unrecognized variable: ", _var.name());
  return 0.;
}