Skip to content

Commit

Permalink
fixed Ambipolar Ion Surface Flux
Browse files Browse the repository at this point in the history
  • Loading branch information
tietz-r committed Apr 19, 2022
1 parent 82461d9 commit e9a0f00
Show file tree
Hide file tree
Showing 10 changed files with 520 additions and 5 deletions.
17 changes: 17 additions & 0 deletions regressioncheck/NIG_DSMC/Ambipolar_Diffusion_SF/DSMC.ini
Original file line number Diff line number Diff line change
@@ -0,0 +1,17 @@
! =============================================================================== !
! Species1, N+
! =============================================================================== !
Part-Species1-SpeciesName=NIon1
Part-Species1-InteractionID = 10
Part-Species1-TRef =273
Part-Species1-DRef = 3.0E-10
Part-Species1-omega=0.24
! =============================================================================== !
! Species2, e
! =============================================================================== !
Part-Species2-SpeciesName = electron
Part-Species2-InteractionID = 4
Part-Species2-TRef =273
Part-Species2-DRef = 2.817920E-15
Part-Species2-omega=0.24

Empty file.
Original file line number Diff line number Diff line change
@@ -0,0 +1,2 @@
MPI=1,6
cmd_suffix=DSMC.ini
47 changes: 47 additions & 0 deletions regressioncheck/NIG_DSMC/Ambipolar_Diffusion_SF/hopr.ini
Original file line number Diff line number Diff line change
@@ -0,0 +1,47 @@
!=============================================================================== !
! MAKEFILE PARAMETER (put a "#" in front, NO blanks!)
!=============================================================================== !
! This is only a dummy parameter needed for the regression check
#MPI=

!=============================================================================== !
! OUTPUT
!=============================================================================== !
ProjectName =cube ! name of the project (used for filenames)
Debugvisu =F ! Write debug mesh to tecplot file
Logging =F ! Write log files

!=============================================================================== !
! MESH
!=============================================================================== !
Mode =1 ! 1 Cartesian 2 gambit file 3 CGNS
nZones =1 ! number of zones
Corner =(/0.,0.,0.,,1E-6,0.,0.,,1E-6,1E-6,0.,,0.,1E-6,0. ,,0.,0.,1E-6,,1E-6,0.,1E-6,,1E-6,1E-6,1E-6,,0.,1E-6,1E-6/) ! [0,1]x[0,1]x[0,0.05]
nElems =(/10,10,10/) ! Anzahl der Elemente in jede Richtung (nfine 4:16 5:32 6:64 7:128)
BCIndex =(/1,1,2,2,3,3/) ! Indices of UserDefinedBoundaries
elemtype =108 ! Elementform (108: Hexaeder)
useCurveds =F ! T if curved boundaries defined
SpaceQuandt =1. ! characteristic length of the mesh
ConformConnect=T
jacobianTolerance = 1e-27
!=============================================================================== !
! BOUNDARY CONDITIONS
!=============================================================================== !
nUserDefinedBoundaries=3
BoundaryName=BC1 ! Outflow: open (absorbing) [for MAXWELL]
BoundaryType=(/4,0,0,0/) ! Outflow: open (absorbing) [for MAXWELL]
BoundaryName=BC2 ! Outflow: open (absorbing) [for MAXWELL]
BoundaryType=(/4,0,0,0/) ! Outflow: open (absorbing) [for MAXWELL]
BoundaryName=BC3 ! Outflow: open (absorbing) [for MAXWELL]
BoundaryType=(/4,0,0,0/) ! Outflow: open (absorbing) [for MAXWELL]

!=============================================================================== !
! BASIS
!=============================================================================== !
NVisu = 7

!=============================================================================== !
! SEARCH
!=============================================================================== !
! nElemsNodeSearch=50
! RefineSideSearch=50
108 changes: 108 additions & 0 deletions regressioncheck/NIG_DSMC/Ambipolar_Diffusion_SF/parameter.ini
Original file line number Diff line number Diff line change
@@ -0,0 +1,108 @@
CFLscale = 0.2
IniExactFunc = 0
N = 1
NAnalyze = 1
! =============================================================================== !
! MESH
! =============================================================================== !
MeshFile = cube_mesh.h5
Logging = F
useCurveds = F
! if boundaries have to be changed (else they are used from Mesh directly):
TrackingMethod = triatracking!, tracing, refmapping
! =============================================================================== !
! OUTPUT / VISUALIZATION
! =============================================================================== !
ProjectName = Reservoir
IterDisplayStep = 1
Part-AnalyzeStep = 1
CalcNumDens = T
CalcTemp = T
! =============================================================================== !
! CALCULATION
! =============================================================================== !
Particles-ManualTimeStep = 1E-12
tend = 1E-10
Analyze_dt = 1E-7 ! Timestep of analyze outputs
! =============================================================================== !
! PARTICLES
! =============================================================================== !
Part-maxParticleNumber=3500000
Part-nSpecies=2
Part-nBounds=3
Part-Boundary1-SourceName=BC1
Part-Boundary1-Condition=open
Part-Boundary2-SourceName=BC2
Part-Boundary2-Condition=open
Part-Boundary3-SourceName=BC3
Part-Boundary3-Condition=open
Part-FIBGMdeltas=(/1E-6,1E-6,1E-6/)
Particles-HaloEpsVelo=5000
! =============================================================================== !
! DSMC
! =============================================================================== !
UseDSMC = T
Particles-DSMC-CollisMode = 1
Particles-DSMC-CalcQualityFactors = F

Particles-DSMC-AmbipolarDiffusion = T

Particles-DSMC-ElectronicModel = 0
! =============================================================================== !
! Weighting Factor
! =============================================================================== !
Part-Species1-MacroParticleFactor = 5
Part-Species2-MacroParticleFactor = 5
! =============================================================================== !
! Species1 | N+
! =============================================================================== !
Part-Species1-ChargeIC = 1.60217653E-19
Part-Species1-MassIC = 2.32600E-26

Part-Species1-nInits = 0
Part-Species1-nSurfacefluxBCs=3
Part-Species1-Surfaceflux1-BC=1
Part-Species1-SurfaceFlux1-velocityDistribution = maxwell_lpn
Part-Species1-SurfaceFlux1-PartDensity = 2E23
Part-Species1-SurfaceFlux1-VeloIC = 1000.
Part-Species1-SurfaceFlux1-VeloVecIC = (/0.,1.,0./)
Part-Species1-SurfaceFlux1-MWTemperatureIC = 1000.
Part-Species1-Surfaceflux2-BC=2
Part-Species1-SurfaceFlux2-velocityDistribution = maxwell_lpn
Part-Species1-SurfaceFlux2-PartDensity = 2E23
Part-Species1-SurfaceFlux2-VeloIC = 1000.
Part-Species1-SurfaceFlux2-VeloVecIC = (/0.,1.,0./)
Part-Species1-SurfaceFlux2-MWTemperatureIC = 1000.
Part-Species1-Surfaceflux3-BC=3
Part-Species1-SurfaceFlux3-velocityDistribution = maxwell_lpn
Part-Species1-SurfaceFlux3-PartDensity = 2E23
Part-Species1-SurfaceFlux3-VeloIC = 1000.
Part-Species1-SurfaceFlux3-VeloVecIC = (/0.,1.,0./)
Part-Species1-SurfaceFlux3-MWTemperatureIC = 1000.
! =============================================================================== !
! Species2 | e
! =============================================================================== !
Part-Species2-ChargeIC = -1.60217653E-19
Part-Species2-MassIC = 9.1093826E-31
Part-Species2-nInits = 0
Part-Species2-nSurfacefluxBCs=3
Part-Species2-Surfaceflux1-BC=1
Part-Species2-SurfaceFlux1-velocityDistribution = maxwell_lpn
Part-Species2-SurfaceFlux1-PartDensity = 2E23
Part-Species2-SurfaceFlux1-VeloIC = 1000.
Part-Species2-SurfaceFlux1-VeloVecIC = (/0.,1.,0./)
Part-Species2-SurfaceFlux1-MWTemperatureIC = 1000.
Part-Species2-Surfaceflux2-BC=2
Part-Species2-SurfaceFlux2-velocityDistribution = maxwell_lpn
Part-Species2-SurfaceFlux2-PartDensity = 2E23
Part-Species2-SurfaceFlux2-VeloIC = 1000.
Part-Species2-SurfaceFlux2-VeloVecIC = (/0.,1.,0./)
Part-Species2-SurfaceFlux2-MWTemperatureIC = 1000.
Part-Species2-Surfaceflux3-BC=3
Part-Species2-SurfaceFlux3-velocityDistribution = maxwell_lpn
Part-Species2-SurfaceFlux3-PartDensity = 2E23
Part-Species2-SurfaceFlux3-VeloIC = 1000.
Part-Species2-SurfaceFlux3-VeloVecIC = (/0.,1.,0./)
Part-Species2-SurfaceFlux3-MWTemperatureIC = 1000.


4 changes: 4 additions & 0 deletions regressioncheck/NIG_DSMC/Ambipolar_Diffusion_SF/readme.md
Original file line number Diff line number Diff line change
@@ -0,0 +1,4 @@
# Surface Flux, Ambipolar Diffusion
* Empty computational Domains with surface flux of ions and electrons applied on all sides
* Electrons are treated with the ambipolar diffusion model (attached to the ion movement, but with their own velocity vector)
* Sanaty check will fail if there are a different number of electrons and ions in an element and leads to an abort
Loading

0 comments on commit e9a0f00

Please sign in to comment.