Coupler: single-block hyperbolic-elliptic coupler (Poisson -> aux -> advance loop). More...

#include <pops/core/foundation/types.hpp>
#include <pops/coupling/base/aux_fill.hpp>
#include <pops/coupling/base/coupling_policy.hpp>
#include <pops/coupling/base/elliptic_rhs.hpp>
#include <pops/numerics/elliptic/interface/elliptic_problem.hpp>
#include <pops/numerics/elliptic/interface/elliptic_solver.hpp>
#include <pops/numerics/elliptic/mg/geometric_mg.hpp>
#include <pops/numerics/time/integrators/time_integrator.hpp>
#include <pops/numerics/time/integrators/time_steppers.hpp>
#include <pops/mesh/layout/box_array.hpp>
#include <pops/mesh/layout/distribution_mapping.hpp>
#include <pops/mesh/storage/fab2d.hpp>
#include <pops/mesh/execution/for_each.hpp>
#include <pops/mesh/geometry/geometry.hpp>
#include <pops/mesh/storage/mf_arith.hpp>
#include <pops/mesh/storage/multifab.hpp>
#include <pops/mesh/boundary/physical_bc.hpp>
#include <pops/numerics/fv/reconstruction.hpp>
#include <pops/numerics/fv/spatial_discretisation.hpp>
#include <pops/numerics/spatial_operator.hpp>
#include <pops/parallel/comm.hpp>
#include <functional>
#include <type_traits>
#include <utility>

Include dependency graph for coupler.hpp:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	pops::Coupler< Model, Elliptic >
	Single-block hyperbolic-elliptic coupler. More...

Namespaces
namespace	pops

namespace	pops::detail

Functions
template<class Model >
void	pops::detail::coupler_eval_rhs (const MultiFab &state, MultiFab &rhs, const Model &model)
	Assemble the single-model elliptic RHS: rhs = model.elliptic_rhs(U) on valid cells (delegated to SingleModelEllipticRhs).

void	pops::detail::coupler_grad_phi (const MultiFab &phi, MultiFab &aux, Real cx, Real cy)
	Set aux = (phi, d phi/dx, d phi/dy) by centered differences (factors cx, cy = 1/(2 dx), 1/(2 dy)).

Detailed Description

Coupler: single-block hyperbolic-elliptic coupler (Poisson -> aux -> advance loop).

At each integrator stage (stage-by-stage coupling): (1) RHS f = elliptic_rhs(model, U); (2) solve lap(phi) = f with the elliptic backend (warm start); (3) aux = (phi, grad phi) by centered differences; (4) assemble the hyperbolic residual with this aux. For drift transport aux enters through the FLUX (E x B); for a self-gravitating fluid through the SOURCE. Three orthogonal axes, all template parameters: Limiter (reconstruction), Policy (PerStage vs OncePerStep), NumericalFlux (Rusanov by default). Compatible with a SINGLE model; multi-species goes through SystemCoupler. The detail:: helpers are at namespace scope (an POPS_HD extended lambda cannot live in a private method, an nvcc restriction).

include/pops/coupling/single/coupler.hpp File Reference¶

Classes

Namespaces

Functions

Detailed Description