cs_cdofb_ac.h

Go to the documentation of this file.

#ifndef CS_CDOFB_AC_H
#define CS_CDOFB_AC_H
 
/*============================================================================
 * Build an algebraic CDO face-based system for the Navier-Stokes equations
 * and solved it with an artificial compressibility algorithm
 *============================================================================*/
 
/*
  This file is part of code_saturne, a general-purpose CFD tool.
 
  Copyright (C) 1998-2026 EDF S.A.
 
  This program is free software; you can redistribute it and/or modify it under
  the terms of the GNU General Public License as published by the Free Software
  Foundation; either version 2 of the License, or (at your option) any later
  version.
 
  This program is distributed in the hope that it will be useful, but WITHOUT
  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
  details.
 
  You should have received a copy of the GNU General Public License along with
  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
  Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
 
/*----------------------------------------------------------------------------*/
 
#include "base/cs_defs.h"
 
/*----------------------------------------------------------------------------
 * Standard C library headers
 *----------------------------------------------------------------------------*/
 
/*----------------------------------------------------------------------------
 *  Local headers
 *----------------------------------------------------------------------------*/
 
#include "base/cs_base.h"
#include "base/cs_time_step.h"
#include "cdo/cs_cdo_connect.h"
#include "cdo/cs_cdo_quantities.h"
#include "cdo/cs_cdofb_navsto.h"
#include "cdo/cs_equation.h"
#include "cdo/cs_navsto_context.h"
#include "cdo/cs_navsto_coupling.h"
#include "cdo/cs_navsto_param.h"
#include "cdo/cs_source_term.h"
#include "mesh/cs_mesh.h"
 
/*============================================================================
 * Macro definitions
 *============================================================================*/
 
/*============================================================================
 * Type definitions
 *============================================================================*/
 
struct cs_cdofb_ac_t : public cs::cdo_navsto_ac_ctx_t {
  cs_cdofb_navsto_source_t *add_gravity_term;
 
};
 
/*============================================================================
 * Inline static function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
inline static cs_real_t *
cs_cdofb_ac_get_face_velocity(bool   previous)
{
  return cs_equation_get_face_values(cs_equation_by_name("momentum"),
                                     previous);
}
 
/*============================================================================
 * Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_cdofb_ac_init_common(const cs_cdo_quantities_t     *quant,
                        const cs_cdo_connect_t        *connect,
                        const cs_time_step_t          *time_step);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs::cdo_navsto_ctx_t *
cs_cdofb_ac_init_scheme_context(const cs_navsto_param_t *nsp,
                                cs_adv_field_t          *adv_field,
                                cs_real_t               *mflux,
                                cs_real_t               *mflux_pre,
                                cs_boundary_type_t      *fb_type,
                                void                    *nsc_input);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void *
cs_cdofb_ac_free_scheme_context(cs_cdofb_ac_t *sc);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_cdofb_ac_compute_implicit(const cs_mesh_t         *mesh,
                             const cs_navsto_param_t *nsp,
                             cs_cdofb_ac_t           *sc);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_cdofb_ac_compute_implicit_nl(const cs_mesh_t         *mesh,
                                const cs_navsto_param_t *nsp,
                                cs_cdofb_ac_t           *sc);
 
#endif /* CS_CDOFB_AC_H */