fluid_dynamics_inner.h

Go to the documentation of this file.

/* -------------------------------------------------------------------------*
*                               SPHinXsys                                   *
* --------------------------------------------------------------------------*
* SPHinXsys (pronunciation: s'finksis) is an acronym from Smoothed Particle *
* Hydrodynamics for industrial compleX systems. It provides C++ APIs for    *
* physical accurate simulation and aims to model coupled industrial dynamic *
* systems including fluid, solid, multi-body dynamics and beyond with SPH   *
* (smoothed particle hydrodynamics), a meshless computational method using  *
* particle discretization.                                                  *
*                                                                           *
* SPHinXsys is partially funded by German Research Foundation               *
* (Deutsche Forschungsgemeinschaft) DFG HU1527/6-1, HU1527/10-1             *
* and HU1527/12-1.                                                          *
*                                                                           *
* Portions copyright (c) 2017-2020 Technical University of Munich and       *
* the authors' affiliations.                                                *
*                                                                           *
* Licensed under the Apache License, Version 2.0 (the "License"); you may   *
* not use this file except in compliance with the License. You may obtain a *
* copy of the License at http://www.apache.org/licenses/LICENSE-2.0.        *
*                                                                           *
* --------------------------------------------------------------------------*/
#ifndef FLUID_DYNAMICS_INNER_H
#define FLUID_DYNAMICS_INNER_H

#include "all_particle_dynamics.h"
#include "base_kernel.h"
#include "all_body_relations.h"
#include "fluid_body.h"
#include "fluid_particles.h"
#include "weakly_compressible_fluid.h"
#include "riemann_solver.h"

namespace SPH
{
    namespace fluid_dynamics
    {
        typedef DataDelegateSimple<FluidBody, FluidParticles, Fluid> FluidDataSimple;
        typedef DataDelegateInner<FluidBody, FluidParticles, Fluid> FluidDataInner;

        class FluidInitialCondition : public ParticleDynamicsSimple, public FluidDataSimple
        {
        public:
            explicit FluidInitialCondition(FluidBody &fluid_body);
            virtual ~FluidInitialCondition(){};

        protected:
            StdLargeVec<Vecd> &pos_, &vel_;
        };

        class DensitySummationInner : public InteractionDynamicsWithUpdate, public FluidDataInner
        {
        public:
            explicit DensitySummationInner(BaseBodyRelationInner &inner_relation);
            virtual ~DensitySummationInner(){};

        protected:
            Real W0_, rho0_, inv_sigma0_;
            StdLargeVec<Real> &Vol_, &rho_, &mass_, &rho_sum_;

            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
            virtual void Update(size_t index_i, Real dt = 0.0) override;
            virtual Real ReinitializedDensity(Real rho_sum, Real rho_0, Real rho_n) { return rho_sum; };
        };

        class ViscousAccelerationInner : public InteractionDynamics, public FluidDataInner
        {
        public:
            explicit ViscousAccelerationInner(BaseBodyRelationInner &inner_relation);
            virtual ~ViscousAccelerationInner(){};

        protected:
            Real mu_;
            Real smoothing_length_;
            StdLargeVec<Real> &Vol_, &rho_, &p_;
            StdLargeVec<Vecd> &vel_, &acc_prior_;

            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
        };

        class AngularConservativeViscousAccelerationInner : public ViscousAccelerationInner
        {
        public:
            explicit AngularConservativeViscousAccelerationInner(BaseBodyRelationInner &inner_relation)
                : ViscousAccelerationInner(inner_relation){};
            virtual ~AngularConservativeViscousAccelerationInner(){};

        protected:
            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
        };

        class TransportVelocityCorrectionInner : public InteractionDynamics, public FluidDataInner
        {
        public:
            explicit TransportVelocityCorrectionInner(BaseBodyRelationInner &inner_relation);
            virtual ~TransportVelocityCorrectionInner(){};

        protected:
            StdLargeVec<Real> &Vol_, &rho_;
            StdLargeVec<Vecd> &pos_;
            StdLargeVec<int> &surface_indicator_;
            Real p_background_;

            virtual void setupDynamics(Real dt = 0.0) override;
            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
        };

        class AcousticTimeStepSize : public ParticleDynamicsReduce<Real, ReduceMax>, public FluidDataSimple
        {
        public:
            explicit AcousticTimeStepSize(FluidBody &fluid_body);
            virtual ~AcousticTimeStepSize(){};

        protected:
            StdLargeVec<Real> &rho_, &p_;
            StdLargeVec<Vecd> &vel_;
            Real smoothing_length_;
            Real ReduceFunction(size_t index_i, Real dt = 0.0) override;
            Real OutputResult(Real reduced_value) override;
        };

        class AdvectionTimeStepSize : public ParticleDynamicsReduce<Real, ReduceMax>, public FluidDataSimple
        {
        public:
            explicit AdvectionTimeStepSize(FluidBody &fluid_body, Real U_max);
            virtual ~AdvectionTimeStepSize(){};

        protected:
            Real smoothing_length_;
            StdLargeVec<Vecd> &vel_;
            Real ReduceFunction(size_t index_i, Real dt = 0.0) override;
            Real OutputResult(Real reduced_value) override;
        };

        class AdvectionTimeStepSizeForImplicitViscosity : public AdvectionTimeStepSize
        {
        public:
            explicit AdvectionTimeStepSizeForImplicitViscosity(FluidBody &fluid_body, Real U_max);
            virtual ~AdvectionTimeStepSizeForImplicitViscosity(){};
        };

        class VorticityInner : public InteractionDynamics, public FluidDataInner
        {
        public:
            explicit VorticityInner(BaseBodyRelationInner &inner_relation);
            virtual ~VorticityInner(){};

        protected:
            StdLargeVec<Real> &Vol_;
            StdLargeVec<Vecd> &vel_;
            StdLargeVec<AngularVecd> vorticity_;
            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
        };

        class BaseRelaxation : public ParticleDynamics1Level, public FluidDataInner
        {
        public:
            explicit BaseRelaxation(BaseBodyRelationInner &inner_relation);
            virtual ~BaseRelaxation(){};

        protected:
            StdLargeVec<Real> &Vol_, &mass_, &rho_, &p_, &drho_dt_;
            StdLargeVec<Vecd> &pos_, &vel_, &acc_, &acc_prior_;
        };

        class BasePressureRelaxation : public BaseRelaxation
        {
        public:
            explicit BasePressureRelaxation(BaseBodyRelationInner &inner_relation);
            virtual ~BasePressureRelaxation(){};

        protected:
            virtual void Initialization(size_t index_i, Real dt = 0.0) override;
            virtual void Update(size_t index_i, Real dt = 0.0) override;
            virtual Vecd computeNonConservativeAcceleration(size_t index_i);
        };

        template <class RiemannSolverType>
        class BasePressureRelaxationInner : public BasePressureRelaxation
        {
        public:
            explicit BasePressureRelaxationInner(BaseBodyRelationInner &inner_relation);
            virtual ~BasePressureRelaxationInner(){};
            RiemannSolverType riemann_solver_;

        protected:
            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
        };
        using PressureRelaxationInner = BasePressureRelaxationInner<NoRiemannSolver>;
        using PressureRelaxationRiemannInner = BasePressureRelaxationInner<AcousticRiemannSolver>;
        using PressureRelaxationDissipativeRiemannInner = BasePressureRelaxationInner<DissipativeRiemannSolver>;

        class BaseDensityRelaxation : public BaseRelaxation
        {
        public:
            explicit BaseDensityRelaxation(BaseBodyRelationInner &inner_relation);
            virtual ~BaseDensityRelaxation(){};

        protected:
            virtual void Initialization(size_t index_i, Real dt = 0.0) override;
            virtual void Update(size_t index_i, Real dt = 0.0) override;
        };

        template <class RiemannSolverType>
        class BaseDensityRelaxationInner : public BaseDensityRelaxation
        {
        public:
            explicit BaseDensityRelaxationInner(BaseBodyRelationInner &inner_relation);
            virtual ~BaseDensityRelaxationInner(){};
            RiemannSolverType riemann_solver_;

        protected:
            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
        };
        using DensityRelaxationInner = BaseDensityRelaxationInner<NoRiemannSolver>;
        using DensityRelaxationRiemannInner = BaseDensityRelaxationInner<AcousticRiemannSolver>;
        using DensityRelaxationDissipativeRiemannInner = BaseDensityRelaxationInner<DissipativeRiemannSolver>;

        class Oldroyd_B_FluidInitialCondition : public ParticleDynamicsSimple, public FluidDataSimple
        {
        public:
            explicit Oldroyd_B_FluidInitialCondition(FluidBody &fluid_body)
                : ParticleDynamicsSimple(fluid_body), FluidDataSimple(fluid_body){};
            virtual ~Oldroyd_B_FluidInitialCondition(){};
        };

        class PressureRelaxationInnerOldroyd_B : public PressureRelaxationDissipativeRiemannInner
        {
        public:
            explicit PressureRelaxationInnerOldroyd_B(BaseBodyRelationInner &inner_relation);
            virtual ~PressureRelaxationInnerOldroyd_B(){};

        protected:
            StdLargeVec<Matd> &tau_, &dtau_dt_;
            virtual void Initialization(size_t index_i, Real dt = 0.0) override;
            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
        };

        class DensityRelaxationInnerOldroyd_B : public DensityRelaxationDissipativeRiemannInner
        {
        public:
            explicit DensityRelaxationInnerOldroyd_B(BaseBodyRelationInner &inner_relation);
            virtual ~DensityRelaxationInnerOldroyd_B(){};

        protected:
            StdLargeVec<Matd> &tau_, &dtau_dt_;
            Real mu_p_, lambda_;

            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
            virtual void Update(size_t index_i, Real dt = 0.0) override;
        };
    }
}
#endif //FLUID_DYNAMICS_INNER_H