Skip to content

Commit

Permalink
Merge branch '220-bugfix-dielectric-mortar' into 'master.dev'
Browse files Browse the repository at this point in the history
Resolve "Bugfix dielectric + mortar"

Closes #220

See merge request piclas/piclas!873
  • Loading branch information
scopplestone committed Dec 18, 2023
2 parents 3d1ff06 + d92c431 commit be76595
Show file tree
Hide file tree
Showing 8 changed files with 53 additions and 306 deletions.
Original file line number Diff line number Diff line change
@@ -1,7 +1,2 @@
! Deactivate reggie for now
! Fix by implementing dielectric <-> vacuum mortar interfaces correctly
! Delete the following line afterwards
PICLAS_EQNSYSNAME = poisson

! run only with particles
PICLAS_PARTICLES=OFF
Original file line number Diff line number Diff line change
Expand Up @@ -6,16 +6,11 @@ maxIterCG = 12000

crossProductMetrics = F ! bad for hyperbolic equations (free stream preservation), maybe good for HDG??


MeshFile = MESH3_sphere_center_004_002_NGeo2_mesh.h5, MESH3_sphere_center_008_004_NGeo2_mesh.h5 ! --> NGeo=2



useCurveds = T

!MeshFile = sphere_004_NGeo4_mesh.h5
DielectricRadiusValue = 1.0

! =============================================================================== !
! EQUATION (linearscalaradvection)
! =============================================================================== !
Expand All @@ -26,56 +21,24 @@ TrackingMethod = refmapping
! OUTPUT / VISUALIZATION
! =============================================================================== !
ProjectName = sphere
!NVisu = 4 ! Number of visualization points
!NodeType_visu = VISU ! VISU... Equidistant points
!VarName = ElectricFieldX
!VarName = ElectricFieldY
!VarName = ElectricFieldZ
!Visu3D = T

Logging = F
WriteErrorFiles= F
printRandomSeeds = F
DoCalcErrorNorms = T
! =============================================================================== !
! ANALYZE
! =============================================================================== !
CalcKineticEnergy = T
CalcPotentialEnergy = T
CalcNumSpec = T
CalcInternalEnergy = T

CalcTemp = T
CalcTransTemp = T ! noch nicht im trunk

CalcPartBalance = T! T -> defekt

CalcVelos = F! T
VelocityDirections = (/1,1,1,1/) ! x,y,z,abs

CalcCollRates = T
CalcReacRates = T
Particles-DSMC-CalcQualityFactors = T

Part-WriteMacroValues = T
Part-IterationForMacroVal = 100
Part-WriteFieldsToVTK = T
! =============================================================================== !
! HDF5-Output
! =============================================================================== !
PIC-OutputSource = T

! =============================================================================== !
! CALCULATION
! =============================================================================== !
tend = 0.1!100.0E-12
Analyze_dt = 0.1!100.0E-12
tend = 0.1
Analyze_dt = 0.1
Part-DelayTime = 1 ! delay time for particle push
!ManualTimeStep=1.12E-13

ManualTimestep=1!e-13 ! 1E21: 1.12E-13
! 1E23: 1E-14
Particles-HaloEpsVelo=1!300E6
ManualTimestep=1
Particles-HaloEpsVelo=1

CFLscale = 0.1 ! Scaling of theoretical CFL number
! =============================================================================== !
Expand All @@ -86,7 +49,6 @@ DielectricMuR = 1
xyzDielectricMinMax = (/-1.0,1.0,-1.0,1.0,-1.0,1.0/) !
DielectricCheckRadius = T
Dielectric_E_0 = -1.

! =============================================================================== !
! Field Boundaries
! =============================================================================== !
Expand All @@ -97,123 +59,32 @@ BoundaryType=(/2,200/)
! =============================================================================== !
! PARTICLES
! =============================================================================== !
Part-maxParticleNumber=5!7000000
Part-nSpecies=1!7
PIC-externalField=(/0.,0.,0.,0.,0.,0./)
Part-maxParticleNumber=5
Part-nSpecies=1

Part-FIBGMdeltas=(/0.1,0.1,0.1/)

PIC-Deposition-Type = shape_function
PIC-shapefunction-radius = 0.3
PIC-Depo-Periodic = FALSE

Part-nBounds = 1
Part-Boundary1-SourceName = BC_outer
Part-Boundary1-Condition = open

! =============================================================================== !
! DSMC
! =============================================================================== !
UseDSMC=F
Particles-DSMCElectronicDatabase=DSMCSpecies_electronic_state_full_Data.h5 ! when supplied: doQK=true
Particles-NumberForDSMCOutputs=0
Part-TimeFracForSampling=0.0
Particles-DSMC-CollisMode=3 !(1:elast coll, 2: elast + rela, 3:chem)
Part-NumberOfRandomSeeds =2
Particles-RandomSeed1= 1
Particles-RandomSeed2= 2

! for chem reac + HDG
Particles-DSMC-UseOctree=T
Particles-DSMC-UseNearestNeighbour = T
Particles-OctreePartNumNode=80
Particles-OctreePartNumNodeMin=50

EpsMergeElectronicState = 1.e-2 ! merge QK levels when difference falls below eps
Particles-DSMC-BackwardReacRate = T ! automatic reaction backward rate calculation (when foward rate is defined)
Particles-DSMC-PartitionMaxTemp = 1e9! maximale temperatur für allocation von kb=f(T), sonst abort
Particles-DSMC-PartitionInterval = 50 ! dT for Intervall generation of kb

! =============================================================================== !
! Species1 | Al
! =============================================================================== !
Part-Species1-ChargeIC = 1
Part-Species1-MassIC = 100000!4.47908080018E-26
Part-Species1-MassIC = 100000
Part-Species1-MacroParticleFactor = 88
Part-Species1-nInits = 1

Part-Species1-Init1-SpaceIC = point
Part-Species1-Init1-velocityDistribution = maxwell
Part-Species1-Init1-MWTemperatureIC = 400
Part-Species1-Init1-ParticleNumber = 1
Part-Species1-Init1-RadiusIC = 0.!20E-6
Part-Species1-Init1-BasePointIC = (/0.0,0.0,0.0/)
Part-Species1-Init1-BaseVector1IC = (/0.0,0.0,0.0/)
Part-Species1-Init1-BaseVector2IC = (/0.0,0.0,0.0/)
Part-Species1-Init1-CylinderHeightIC = 0.!0E-6
Part-Species1-Init1-VeloIC = 0
Part-Species1-Init1-VeloVecIC = (/0.,0.,1./)
Part-Species1-Init1-Tempelec = 400
! =============================================================================== !
! Species2, e-
! =============================================================================== !
Part-Species2-ChargeIC = -1.60217653E-19
Part-Species2-MassIC = 9.1093826E-31
Part-Species2-MacroParticleFactor = 88
Part-Species2-nInits = 1

Part-Species2-Init1-SpaceIC = cylinder
Part-Species2-Init1-velocityDistribution = maxwell
Part-Species2-Init1-MWTemperatureIC = 1.160e7
Part-Species2-Init1-ParticleNumber = 2000000
Part-Species2-Init1-RadiusIC = 20E-6
Part-Species2-Init1-BasePointIC = (/0.0,0.0,0.0/)
Part-Species2-Init1-BaseVector1IC = (/1.0,0.0,0.0/)
Part-Species2-Init1-BaseVector2IC = (/0.0,1.0,0.0/)
Part-Species2-Init1-CylinderHeightIC = 70E-6
Part-Species2-Init1-VeloIC = 0
Part-Species2-Init1-VeloVecIC = (/0.,0.,1./)
! =============================================================================== !
! Species3, Al+
! =============================================================================== !
Part-Species3-ChargeIC=1.60217653E-19
Part-Species3-MassIC=4.4789897064E-26
Part-Species3-MacroParticleFactor=88
! =============================================================================== !
! Species4, Al+2
! =============================================================================== !
Part-Species4-ChargeIC = 3.20435306000000E-19
Part-Species4-MassIC = 4.4788986125E-26
Part-Species4-MacroParticleFactor = 88
Part-Species4-nInits = 1

Part-Species4-Init1-SpaceIC = cylinder
Part-Species4-Init1-velocityDistribution = maxwell
Part-Species4-Init1-MWTemperatureIC = 1.160e5
Part-Species4-Init1-ParticleNumber = 1000000
Part-Species4-Init1-RadiusIC = 20E-6
Part-Species4-Init1-BasePointIC = (/0.0,0.0,0.0/)
Part-Species4-Init1-BaseVector1IC = (/1.0,0.0,0.0/)
Part-Species4-Init1-BaseVector2IC = (/0.0,1.0,0.0/)
Part-Species4-Init1-CylinderHeightIC = 70E-6
Part-Species4-Init1-VeloIC = 0
Part-Species4-Init1-VeloVecIC = (/0.,0.,1./)
Part-Species4-Init1-Tempelec = 1.160e5
! =============================================================================== !
! Species5, Al+3
! =============================================================================== !
Part-Species5-ChargeIC=4.80652959000000E-19
Part-Species5-MassIC=4.4788075187E-26
Part-Species5-MacroParticleFactor=88
! =============================================================================== !
! Species6, Al+4
! =============================================================================== !
Part-Species6-ChargeIC= 6.40870612000000E-19
Part-Species6-MassIC=4.4787164249E-26
Part-Species6-MacroParticleFactor=88
! =============================================================================== !
! Species7, Al+5
! =============================================================================== !
Part-Species7-ChargeIC= 8.010882650000E-19
Part-Species7-MassIC=4.4786253311E-26
Part-Species7-MacroParticleFactor=88
Binary file not shown.
Binary file not shown.
Binary file not shown.
Loading

0 comments on commit be76595

Please sign in to comment.