LaserParticleContainer.H

Go to the documentation of this file.

/* Copyright 2019 Andrew Myers, Axel Huebl, David Grote
 * Luca Fedeli, Maxence Thevenet, Remi Lehe
 * Weiqun Zhang
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef WARPX_LaserParticleContainer_H_
#define WARPX_LaserParticleContainer_H_
 
#include "Evolve/WarpXDtType.H"
#include "Laser/LaserProfiles.H"
#include "WarpXParticleContainer.H"
 
#include <AMReX_Extension.H>
#include <AMReX_Random.H>
#include <AMReX_REAL.H>
#include <AMReX_RealBox.H>
#include <AMReX_Vector.H>
 
#include <AMReX_BaseFwd.H>
#include <AMReX_AmrCoreFwd.H>
 
#include <limits>
#include <memory>
#include <string>
 
class LaserParticleContainer
    : public WarpXParticleContainer
{
public:
    LaserParticleContainer (amrex::AmrCore* amr_core, int ispecies, const std::string& name);
    virtual ~LaserParticleContainer () {}
 
    virtual void InitData () final;
 
    virtual void DefaultInitializeRuntimeAttributes (
                        amrex::ParticleTile<amrex::Particle<NStructReal, NStructInt>,
                                            NArrayReal, NArrayInt, amrex::PinnedArenaAllocator>& /*pinned_tile*/,
                        const int /*n_external_attr_real*/,
                        const int /*n_external_attr_int*/,
                        const amrex::RandomEngine& /*engine*/) override final {}
 
    virtual void ReadHeader (std::istream& is) final;
 
    virtual void WriteHeader (std::ostream& os) const final;
 
    virtual void Evolve (int lev,
                         const amrex::MultiFab&, const amrex::MultiFab&, const amrex::MultiFab&,
                         const amrex::MultiFab&, const amrex::MultiFab&, const amrex::MultiFab&,
                         amrex::MultiFab& jx, amrex::MultiFab& jy, amrex::MultiFab& jz,
                         amrex::MultiFab*, amrex::MultiFab*, amrex::MultiFab*,
                         amrex::MultiFab* rho, amrex::MultiFab* crho,
                         const amrex::MultiFab*, const amrex::MultiFab*, const amrex::MultiFab*,
                         const amrex::MultiFab*, const amrex::MultiFab*, const amrex::MultiFab*,
                         amrex::Real t, amrex::Real dt, DtType a_dt_type=DtType::Full,
                         bool skip_deposition=false) final;
 
    virtual void PushP (int lev, amrex::Real dt,
                        const amrex::MultiFab& ,
                        const amrex::MultiFab& ,
                        const amrex::MultiFab& ,
                        const amrex::MultiFab& ,
                        const amrex::MultiFab& ,
                        const amrex::MultiFab& ) final;
 
    virtual void PostRestart () final;
 
    void calculate_laser_plane_coordinates (const WarpXParIter& pti, const int np,
                                            amrex::Real * AMREX_RESTRICT const pplane_Xp,
                                            amrex::Real * AMREX_RESTRICT const pplane_Yp);
 
    void update_laser_particle (WarpXParIter& pti, const int np, amrex::ParticleReal * AMREX_RESTRICT const puxp,
                                amrex::ParticleReal * AMREX_RESTRICT const puyp,
                                amrex::ParticleReal * AMREX_RESTRICT const puzp,
                                amrex::ParticleReal const * AMREX_RESTRICT const pwp,
                                amrex::Real const * AMREX_RESTRICT const amplitude,
                                const amrex::Real dt);
 
protected:
 
    std::string m_laser_name;
 
private:
    // runtime paramters
    amrex::Vector<amrex::Real> m_position; 
    amrex::Vector<amrex::Real> m_nvec; 
    amrex::Vector<amrex::Real> m_p_X;// ! Polarization
 
    amrex::Real               m_e_max       = std::numeric_limits<amrex::Real>::quiet_NaN();
    amrex::Real               m_wavelength  = std::numeric_limits<amrex::Real>::quiet_NaN();
 
    amrex::Real               m_Z0_lab = 0; // Position of the antenna in the lab frame
 
    long m_min_particles_per_mode = 4;
 
    // computed using runtime parameters
    amrex::Vector<amrex::Real> m_p_Y;
    amrex::Vector<amrex::Real> m_u_X;
    amrex::Vector<amrex::Real> m_u_Y;
    amrex::Real m_weight   = std::numeric_limits<amrex::Real>::quiet_NaN();
    amrex::Real m_mobility = std::numeric_limits<amrex::Real>::quiet_NaN();
 
 
    // laser particle domain
    amrex::RealBox m_laser_injection_box;
    // Theoretical position of the antenna. Used if do_continuous_injection=1.
    // Track the position of the antenna until it enters the simulation domain.
    amrex::Vector<amrex::Real> m_updated_position;
 
    void ComputeSpacing (int lev, amrex::Real& Sx, amrex::Real& Sy) const;
    void ComputeWeightMobility (amrex::Real Sx, amrex::Real Sy);
    void InitData (int lev);
    // Inject the laser antenna during the simulation, if it started
    // outside of the simulation domain and enters it.
    void ContinuousInjection(const amrex::RealBox& injection_box) override;
    // Update position of the antenna
    void UpdateContinuousInjectionPosition(amrex::Real dt) override;
 
    // Unique (smart) pointer to the laser profile
    std::unique_ptr<WarpXLaserProfiles::ILaserProfile> m_up_laser_profile;
 
    // Flag to disable the laser (e.g., if e_max is 0)
    bool m_enabled = true;
};
 
#endif