eulerian_weakly_compressible_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.        *
*                                                                           *
* --------------------------------------------------------------------------*/
#pragma once

#include "fluid_dynamics_inner.h"

#include "all_particle_dynamics.h"
#include "base_kernel.h"
#include "external_force.h"
#include "riemann_solver.h"

namespace SPH
{
    namespace eulerian_weakly_compressible_fluid_dynamics
    {
        typedef DataDelegateSimple<EulerianFluidBody, WeaklyCompressibleFluidParticles, Fluid> EulerianWeaklyCompressibleFluidDataSimple;
        typedef DataDelegateInner<EulerianFluidBody, WeaklyCompressibleFluidParticles, Fluid> EulerianWeaklyCompressibleFluidDataInner;

        class EulerianFlowTimeStepInitialization
            : public ParticleDynamicsSimple,
            public EulerianWeaklyCompressibleFluidDataSimple
        {
        private:
            UniquePtrKeeper<Gravity> gravity_ptr_keeper_;

        public:
            explicit EulerianFlowTimeStepInitialization(SPHBody &sph_body);
            EulerianFlowTimeStepInitialization(SPHBody &sph_body, Gravity &gravity);
            virtual ~EulerianFlowTimeStepInitialization() {};

        protected:
            StdLargeVec<Real> &rho_, &mass_;
            StdLargeVec<Vecd> &pos_, &vel_, &dmom_dt_prior_;
            Gravity *gravity_;
            virtual void setupDynamics(Real dt = 0.0) override;
            virtual void Update(size_t index_i, Real dt = 0.0) override;
        };

        class FreeSurfaceIndicationInner
            : public InteractionDynamicsWithUpdate, public EulerianWeaklyCompressibleFluidDataInner
        {
        public:
            explicit FreeSurfaceIndicationInner(BaseBodyRelationInner &inner_relation, Real thereshold = 0.75);
            virtual ~FreeSurfaceIndicationInner() {};

        protected:
            Real thereshold_by_dimensions_;
            StdLargeVec<Real>& Vol_;
            StdLargeVec<int>& surface_indicator_;
            StdLargeVec<Real> pos_div_;
            Real smoothing_length_;

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

        class ViscousAccelerationInner
            : public InteractionDynamics, public EulerianWeaklyCompressibleFluidDataInner
        {
        public:
            explicit ViscousAccelerationInner(BaseBodyRelationInner &inner_relation);
            virtual ~ViscousAccelerationInner() {};
        protected:
            Real mu_;
            Real smoothing_length_;
            StdLargeVec<Real> &Vol_, &rho_, &p_;
            StdLargeVec<Vecd> &vel_, &dmom_dt_prior_;

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

        class AcousticTimeStepSize :
            public ParticleDynamicsReduce<Real, ReduceMax>, public EulerianWeaklyCompressibleFluidDataSimple
        {
        public:
            explicit AcousticTimeStepSize(EulerianFluidBody &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 VorticityInner
            : public InteractionDynamics, public EulerianWeaklyCompressibleFluidDataInner
        {
        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 EulerianWeaklyCompressibleFluidDataInner
        {
        public:
            explicit BaseRelaxation(BaseBodyRelationInner &inner_relation);
            virtual ~BaseRelaxation() {};
        protected:
            StdLargeVec<Real>& Vol_, &mass_, &rho_, &p_, &drho_dt_;
            StdLargeVec<Vecd>&vel_, &mom_, &dmom_dt_, &dmom_dt_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;
        };

        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 PressureRelaxationAcousticRiemannInner = BasePressureRelaxationInner<AcousticRiemannSolver>;
        using PressureRelaxationHLLCRiemannInner = BasePressureRelaxationInner<HLLCRiemannSolverInWeaklyCompressibleFluid>;
        using PressureRelaxationHLLCWithLimiterRiemannInner = BasePressureRelaxationInner<HLLCRiemannSolverWithLimiterInWeaklyCompressibleFluid>;

        class BaseDensityAndEnergyRelaxation : public BaseRelaxation
        {
        public:
            explicit BaseDensityAndEnergyRelaxation(BaseBodyRelationInner &inner_relation);
            virtual ~BaseDensityAndEnergyRelaxation() {};
        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 BaseDensityAndEnergyRelaxationInner : public BaseDensityAndEnergyRelaxation
        {
        public:
            explicit BaseDensityAndEnergyRelaxationInner(BaseBodyRelationInner &inner_relation);
            virtual ~BaseDensityAndEnergyRelaxationInner() {};
            RiemannSolverType riemann_solver_;
        protected:
            virtual void Interaction(size_t index_i, Real dt = 0.0) override;
        };
        using DensityAndEnergyRelaxationInner = BaseDensityAndEnergyRelaxationInner<NoRiemannSolver>;
        using DensityAndEnergyRelaxationAcousticRiemannInner = BaseDensityAndEnergyRelaxationInner<AcousticRiemannSolver>;
        using DensityAndEnergyRelaxationHLLCRiemannInner = BaseDensityAndEnergyRelaxationInner<HLLCRiemannSolverInWeaklyCompressibleFluid>;
        using DensityAndEnergyRelaxationHLLCWithLimiterRiemannInner = BaseDensityAndEnergyRelaxationInner<HLLCRiemannSolverWithLimiterInWeaklyCompressibleFluid>;

        class NonReflectiveBoundaryVariableCorrection :public InteractionDynamics, public EulerianWeaklyCompressibleFluidDataInner
        {
        public:
            NonReflectiveBoundaryVariableCorrection(BaseBodyRelationInner &inner_relation) :
                InteractionDynamics(*inner_relation.sph_body_), EulerianWeaklyCompressibleFluidDataInner(inner_relation),
                rho_(particles_->rho_), p_(particles_->p_), vel_(particles_->vel_),
                mom_(particles_->mom_), pos_(particles_->pos_), mass_(particles_->mass_), Vol_(particles_->Vol_),
                surface_indicator_(particles_->surface_indicator_)
            {
                particles_->registerVariable(n_, "NormalDirection");
            };
            virtual ~NonReflectiveBoundaryVariableCorrection() {};
        protected:
            Real p_farfield_, rho_farfield_, gamma_, sound_speed_;
            Vecd vel_farfield_;
            StdLargeVec<Real>& rho_, &p_, &mass_, &Vol_;
            StdLargeVec<Vecd>& vel_, &mom_, &pos_;
            StdLargeVec<Vecd> n_;
            StdLargeVec<int>& surface_indicator_;

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