cs_navsto_coupling.h

Go to the documentation of this file.

#ifndef CS_NAVSTO_COUPLING_H
#define CS_NAVSTO_COUPLING_H
 
/*============================================================================
 * Functions and structures to handle the settings of different
 * velocity/pressure coupling algorithms
 *  - Artificial Compressibility algorithm
 *  - Monolithic algorithm
 *  - Projection 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.
*/
 
/*----------------------------------------------------------------------------
 *  Local headers
 *----------------------------------------------------------------------------*/
 
#include "base/cs_field.h"
#include "cdo/cs_equation.h"
#include "cdo/cs_navsto_param.h"
 
/*============================================================================
 * Macro definitions
 *============================================================================*/
 
/*============================================================================
 * Type definitions
 *============================================================================*/
 
/* Predefined context structures depending on the settings */
/* ======================================================= */
 
typedef struct {
 
  cs_equation_t   *momentum; 
  cs_property_t   *zeta;
 
} cs_navsto_ac_t;
 
 
typedef struct {
 
  cs_equation_t   *momentum;  
} cs_navsto_monolithic_t;
 
 
typedef struct {
 
  cs_equation_t  *prediction; 
  cs_equation_t  *correction; 
  cs_real_t      *div_st;
 
  cs_real_t      *bdy_pressure_incr;
 
  cs_field_t     *predicted_velocity;
 
  cs_field_t     *phi;
 
  cs_field_t    *pressure_incr_gradient;
 
} cs_navsto_projection_t;
 
/*============================================================================
 * Public function prototypes
 *============================================================================*/
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_equation_param_t *
cs_navsto_coupling_get_momentum_eqp(const cs_navsto_param_t    *nsp,
                                    void                       *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void *
cs_navsto_ac_create_context(cs_param_bc_type_t    bc,
                            cs_navsto_param_t    *nsp);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void *
cs_navsto_ac_free_context(void     *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_navsto_ac_init_setup(const cs_navsto_param_t    *nsp,
                        cs_adv_field_t             *adv_field,
                        void                       *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_navsto_ac_last_setup(const cs_navsto_param_t     *nsp,
                        void                        *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_equation_t *
cs_navsto_ac_get_momentum_eq(void       *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void *
cs_navsto_monolithic_create_context(cs_param_bc_type_t    bc,
                                    cs_navsto_param_t    *nsp);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void *
cs_navsto_monolithic_free_context(void             *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_navsto_monolithic_init_setup(const cs_navsto_param_t    *nsp,
                                cs_adv_field_t             *adv_field,
                                void                       *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_navsto_monolithic_last_setup(const cs_navsto_param_t     *nsp,
                                void                        *context);
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_equation_t *
cs_navsto_monolithic_get_momentum_eq(void       *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void *
cs_navsto_projection_create_context(cs_param_bc_type_t    bc,
                                    cs_navsto_param_t    *nsp);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void *
cs_navsto_projection_free_context(void           *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_navsto_projection_init_setup(const cs_navsto_param_t    *nsp,
                                cs_adv_field_t             *adv_field,
                                bool                        has_previous,
                                void                       *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_navsto_projection_last_setup(const cs_cdo_quantities_t  *quant,
                                const cs_navsto_param_t    *nsp,
                                void                       *context);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_equation_t *
cs_navsto_projection_get_momentum_eq(void       *context);
 
#endif /* CS_NAVSTO_COUPLING_H */