Line data Source code
1 : //* This file is part of the MOOSE framework 2 : //* https://mooseframework.inl.gov 3 : //* 4 : //* All rights reserved, see COPYRIGHT for full restrictions 5 : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT 6 : //* 7 : //* Licensed under LGPL 2.1, please see LICENSE for details 8 : //* https://www.gnu.org/licenses/lgpl-2.1.html 9 : 10 : #ifdef MOOSE_MFEM_ENABLED 11 : 12 : #include "MFEMConvectiveHeatFluxBC.h" 13 : #include "MFEMProblem.h" 14 : 15 : registerMooseObject("MooseApp", MFEMConvectiveHeatFluxBC); 16 : 17 : InputParameters 18 8660 : MFEMConvectiveHeatFluxBC::validParams() 19 : { 20 8660 : InputParameters params = MFEMIntegratedBC::validParams(); 21 : // FIXME: Should these really be specified via properties? T_infinity in particular? Use functions 22 : // instead? 23 17320 : params.addClassDescription( 24 : "Convective heat transfer boundary condition with temperature and heat " 25 : "transfer coefficent given by material properties to add to MFEM problems."); 26 34640 : params.addParam<MFEMScalarCoefficientName>( 27 : "T_infinity", "0.", "Name of a coefficient specifying the far-field temperature"); 28 25980 : params.addParam<MFEMScalarCoefficientName>( 29 : "heat_transfer_coefficient", 30 : "1.", 31 : "Name of the coefficient specifying the heat transfer coefficient"); 32 8660 : return params; 33 0 : } 34 : 35 15 : MFEMConvectiveHeatFluxBC::MFEMConvectiveHeatFluxBC(const InputParameters & parameters) 36 : : MFEMIntegratedBC(parameters), 37 15 : _heat_transfer_coef(getScalarCoefficient("heat_transfer_coefficient")), 38 30 : _T_inf_coef(getScalarCoefficient("T_infinity")), 39 15 : _external_heat_flux_coef( 40 30 : getMFEMProblem().getCoefficients().declareScalar<mfem::ProductCoefficient>( 41 30 : "__ConvectiveHeatFluxBC_" + parameters.get<std::string>(MooseBase::unique_name_param), 42 : _heat_transfer_coef, 43 15 : _T_inf_coef)) 44 : { 45 15 : } 46 : 47 : // Create a new MFEM integrator to apply to the RHS of the weak form. Ownership managed by the 48 : // caller. 49 : mfem::LinearFormIntegrator * 50 57 : MFEMConvectiveHeatFluxBC::createLFIntegrator() 51 : { 52 57 : return new mfem::BoundaryLFIntegrator(_external_heat_flux_coef); 53 : } 54 : 55 : // Create a new MFEM integrator to apply to LHS of the weak form. Ownership managed by the caller. 56 : mfem::BilinearFormIntegrator * 57 113 : MFEMConvectiveHeatFluxBC::createBFIntegrator() 58 : { 59 113 : return new mfem::BoundaryMassIntegrator(_heat_transfer_coef); 60 : } 61 : 62 : #endif