forked from bernardparent/CFDWARP
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathCHANGELOG
20044 lines (16221 loc) · 777 KB
/
CHANGELOG
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
_____________________________________________________________________________________________________
CFDWARP CHANGELOG
_____________________________________________________________________________________________________
2022 MAR 16a
* in model/transport/_dixonlewis/transport.c and model/transport/doc/report.tex
- in _mue_from_Nn_Ni_Te() and find_dmue_from_Nn_Ni_Te(), changed the
way the mobility is calculated from muei=1.061e16*pow(Te,1.5)/Ni to muei=1.9e16*pow(Te,1.5)/Ni
2022 MAR 14a
* in model/transport/.config
- changed the gnoffo option to experimental
2022 MAR 07a
* in /model/chem/_Air_Plasma_11s/farbar2013.c
- fixed a bug in reaction 3
* in /model/chem/_Air_Plasma_11s/kim2021.c
- changed the controlling temperatures of the forward and backward reactions and removed the
radiative recombination reactions
- clipped controlling temperatures of backward reactions to the range 300 - 32000 K using a
minmax function
- changed reaction 11 and 12 to match standard Park model
- fixed a typo in reaction 26
* in /model/chem/_Air_Plasma_11s/doc/report.tex
- updated the Kim 11-species chemical model with the correct controlling temperatures
- changed reaction 11 and 12 to match standard Park model and removed the radiative
recombination reactions
- added references for all the reactions of the Kim 11-species chemical model
* in /model/chem/_Air_Plasma_11s/doc/all.bib
- updated references
2022 FEB 28b
* in model/transport/_gnoffo/transport.c
- took away "_gnoffo" at the end of each function name
2022 FEB 28a
* in lib/soap/doc/report.tex
- added log10 function to table
* in lib/soap/soap.c
- in functions_builtin(), added log10 function
* in model/fluid/_*/fluid.c, model/fluid/_*/fluid_conv.c, model/fluid/_*/fluid_post.c
- added find_rhok()
- in the calls to find_nuk_eta_kappa(), _eta(), _kappan_from_rhok_T_Te(), _etan_from_rhok_T_Te()
changed w for rhok
* in model/fluid/_*plasma/fluid.c
- in find_prim_fluid_mem(), fixed a bug related to the determination of rhok (rhok was mistakenly
set to rhok/rho not to rhok)
* in model/share/fluid_share.c
- changed the parameters to find_nuk_from_rhok_T_Te() from w to rhok
* in model/thermo/_generic/thermo.c
- removed all the functions related to the transport coefficients (mobility, diffusion coefficient,
thermal conductivity, viscosity)
- removed the functions related to finding the enthalpy using the NASA Glenn polynomials
* in model/share
- added transport_share.c
* in model/thermo/_generic/
- added enthalpy.c
* in model
- added transport
* in src/test.c
- in test_eta() and test_kappa() and test_Pr(), added rhok instead of w as a parameter to
find_nuk_eta_kappa();
2022 FEB 19a
* in model/thermo/_generic/thermo.c
- set INCLUDE_ELECTRONS_IN_ETA to TRUE
- set INCLUDE_ELECTRONS_IN_KAPPA to TRUE
2022 FEB 18a
* in model/fluid/_Navier-Stokes/fluid_post.c
- in find_post_variable_value_fluid(), fixed a bug related to the parameters passed to _eta and _kappa
2022 FEB 11a
* in model/thermo/_generic/thermo.c
- in _muk_from_rhok_Tk_Ek() and find_dmuk_from_rhok_Tk_Ek(), removed the mobilities for Oplus, Nplus
- in _muk_from_rhok_Tk_Ek() and find_dmuk_from_rhok_Tk_Ek(), modified the default mobility to
match the one outlined in the thermo document
- in _kappac_from_rhok_Tk_Ek(), now use _cpk_from_T_equilibrium instead 2.5*kB/m or 3.5*kB/m
- in _etac_from_rhok_Tk_Ek(), now use _cpk_from_T_equilibrium instead 2.5*kB/m or 3.5*kB/m
* in model/thermo/doc/report.tex
- added section on transport coefficients for partially-ionized plasmas
* in model/fluid/_Navier-Stokes_plasma/fluid_plasma.c
- in find_Dstarmat(), find_Dstarmat_test(), find_Dstar(), find_Dstar_interface()
replaced _hk_from_T() with _hk_from_T_equilibrium()
* in model/fluid/_Favre-Reynolds_plasma/fluid_plasma.c
- in find_Dstar(), replaced _hk_from_T() with _hk_from_T_equilibrium()
2022 FEB 06a
* in lib/soap/eg/calc.c
- added constants eps0, epsilon0, and NA
* in model/chem/_Air_Plasma_11s/chem.hh
- added parent2022a citation
* in model/emfield/_Epotential-Bfixed/emfield.c
- uncommented the function _Ekmag()
* in model/fluid/_*/fluid.c,fluid.h,fluid_post.c
- changed the parameters to _kappa(), _eta(), _ktilde(), _omega()
* in model/fluid/_*_plasma/fluid_source.c
- changed the call from find_dmue_from_rhok_Te() to find_dmuk_from_rhok_Tk_Ek()
- changed the call from _mue_from_rhok_Te() to _muk_from_rhok_Tk_Ek()
* in model/fluid/_*_plasma/fluid.c
- added the function _kappan()
- changed the parameters to the function _kappai()
- in _kappai(), now calls _kappac_from_rhok_Tk_Ek()
- in _kappa(), added the charged species kappa contributions using _kappac_from_rhok_Tk_Ek()
- in _eta(), added the charged species eta contributions using _etac_from_rhok_Tk_Ek()
- in find_Kstar_interface(), now determines the thermal conductivity as kappan+kappai instead
of kappa+kappai
- in find_prim_fluid_mem(), etamem is now determined as the sum of etan and etai found through
_etac_from_rhok_Tk_Ek()
- in find_prim_fluid_mem(), kappamem is now determined as the sum of kappan and kappai found through
_kappac_from_rhok_Tk_Ek()
* in model/fluid/doc/Transport_Equations/report.tex
- Rearranged the various boundary conditions within one section
* in model/share/fluid_share.c
- changed the parameters to functions _eta() and _kappa() and _ktilde()
* in model/thermo/_generic/thermo.c
- modified the charged species viscosities in _etak_from_T()
- added the functions _etan_from_w_T(), _etac_from_rhok_Tk_Ek(), _kappac_from_rhok_Tk_Ek()
- added the function _mui_from_Nn_Ni_Ti() which includes the contribution to the ion mobility
coming from ion-ion collisions
- added the function find_dmui_from_Nn_Ni_Ti()
* in model/thermo/doc/report.tex
- Added section on Transport Coefficients of Fully-Ionized Plasmas
* in src/control.c
- added parent2022a citation
2022 JAN 22a
* in model/fluid/doc/Transport_Equations/report.tex
- improved the derivation of the thermal conductivity of the electrons in the presence
of a magnetic field (now function of the tensor mobility mutilde, not the mobility mu)
2022 JAN 16a
* in src/control.c
- added kim2021a citation
2022 JAN 14a
* in model/chem/_Air_Plasma_11s/
- added kim2021.c and kim2021.h
* in model/chem/_Air_Plasma_11s/chem.c and model/chem/_Air_Plasma_11s/chem.hh
- incorporated the kim2021 model
* in model/chem/_Air_Plasma_11s/doc/report.tex
- added the kim2021 model documentation
* in src/control.c
- changed the copyright years for Bernard Parent to include 2022
2021 DEC 07a
* in model/fluid/doc/Transport_Equations/report.tex
- added a section on the Thermionic Emission Boundary Condition
* in model/metrics/_generic/metrics.c
- in _distance_between_near_bdry_node_and_boundary_plane(), added i,j,k position of problematic
node before the fatal_error
2021 NOV 17a
* in model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- in update_bdry_wall(), fixed a bug related to CATALYTIC and INJECTION wall conditions:
np[lA].numbdryparam-2 should be np[lA].numbdryparam-1
2021 NOV 12a
* in model/fluid/_Navier-Stokes_plasma/fluid.c
- in write_model_fluid_template() and read_model_fluid_actions(), took away reading in
gamma (secondary electron emission coefficient at the surface) as a global parameter
* in model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- in write_bdry_fluid_template() and update_bdry_wall() made changes so that gamma is now read
as an additional argument to the boundary condition call and not as a global parameter
- in update_bdry_wall(), took away secondary electron emission in the case of a thermionic
boundary condition
2021 OCT 30a
* in model/fluid/_Navier-Stokes/fluid_bdry.c and model/fluid/_Favre-Reynolds_Tv/fluid_bdry.c
- in update_bdry_wall(), now sets charged densities to zero even when INJECTION is TRUE
2021 OCT 29a
* in model/metrics/_generic/metrics.c
- in find_unit_vector_normal_to_boundary_plane(), got rid of the fatal_error and replaced
with a return(FALSE) in case of not being able to find a boundary plane
* in model/share/fluid_share.c
- in _w_product_at_catalytic_wall(), added a if condition so that find_Vk is called only if
specR is an ion or electron
- in update_w_at_catalytic_wall(), now finds wwall[specp] using _w_product_at_catalytic_wall()
* in model/fluid/doc/Transport_Equations/report.tex
- in the wall catalytic BC, fixed a typo in the equation determining the boundary condition
for the product species
* in model/fluid/_Favre-Reynolds_Tv/fluid_bdry.c, model/fluid/_Navier-Stokes/fluid_bdry.c,
and model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- in update_bdry_wall(), fixed a bug: when adjusting w at the wall so that no mass is lost
or created by the catalicity of the electrons and ions, the contribution is now added to
(instead of replacing) the previous contributions
* in ./model/share/fluid_share.c, ./model/fluid/_Navier-Stokes/fluid_bdry.c
and ./model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- changed the parameters to update_w_at_catalytic_wall()
2021 OCT 28b
* in model/fluid/_Favre-Reynolds_Tv/fluid_bdry.c
- added BDRY_WALLTFIXEDINJECTION1
- in update_bdry_wall(), added code to make sure that surface catalyticity of electrons with
ions leads to production of neutrals
* in model/fluid/_Navier-Stokes/fluid_bdry.c
- in write_bdry_fluid_template(), changed the default template of bdry conditions: need to make
sure a wall node is never in a corner
2021 OCT 28a
* in model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- in _Vk_drift_perp_to_bdry(), changed the algorithm so that only nodes perpendicular to
the surface are used to determie Vk_drift
- added _mu_effective_perpendicular_boundary_plane()
- in update_bdry_wall(), changed the way the thermionic electrode BC is applied: now
when the electrode is an anode, the electron flux out of the surface is set to diffusion
when the electrode is a cathode, the electron flux out of the surface is set to convection
* in bin/MHD_earth_entry2.wrp
- made various changes to the grid and cycle code.
2021 OCT 14a
* in model/chem/_Air_Plasma_11s/doc/report.tex
- in the Parent-Park model, fixed a typo related to the third body lists
* in model/share/fluid_share.c
- in _w_product_at_catalytic_wall(), fixed a variable definition bug: dim_t VionRB
should be EXM_vec3D_t VionRB
* in src/common.c
- added _smooth3()
2021 OCT 06a
* in model/emfield/_Epotential-Bfixed/emfield.c and model/emfield/_Epotential/emfield.c
- added find_DeltaVk_drift_from_mu()
* in model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- renamed _Vk_perp_to_bdry() to _Vk_drift_perp_to_bdry()
- in _Vk_drift_perp_to_bdry(), now calls find_Vk_drift() and
find_unit_vector_normal_to_boundary_plane()
* in model/share/emfield_share.c
- added find_Vk_drift()
* in bin/MHD_earth_entry2.wrp
- made improvements
2021 OCT 03a
* in model/fluid/_Navier-Stokes_plasma/fluid_post.c
- got rid of the last 4 post processing variables related to cesium
* in model/fluid/doc/Transport_Equations/report.tex
- in model/fluid/doc/Transport_Equations/report.tex, cleaned up the text about the
surface injection boundary condition
* in model/share/fluid_share.c
- in write_disc_fluid_template(), set default zetaD to 0.0 instead of 1.0
2021 OCT 01a
* in cycle/res/_generic/res.c
- changed YSTARDH_SECONDORDER to FALSE
- changed TVDLIMITER to LIMITER_MINMOD
* in model/chem/_Cs_Air_Plasma_15s/chem.c
- in add_dW_dx_Cesium(), fixed a bug: when CESIUMIONIZATIONMODEL was not set to
CESIUMIONIZATIONMODEL_LENARD1964, the linearization of W was not correct
- in add_dW_dx_Cesium(), fixed a bug: when CESIUMRECOMBINATIONMODEL was not set to
CESIUMRECOMBINATIONMODEL_LENARD1964, the linearization of W was not correct
2021 SEP 28a
* in model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- in update_bdry_wall(), rewrote the code related to wwall[specneutral]; this was done
for clarity and should not change the model
2021 SEP 26a
* in model/fluid/_Navier-Stokes/fluid_bdry.c
- in update_bdry_wall(), made w at wall node equal to w at near wall node no matter
the order of accuracy specified
- in update_bdry_wall() now does not set w of the charged species to zero at injection surface
- in update_bdry_wall(), at non injection surfaces added correction to the mass fraction
gradient of the neutral species which is the product of a positive ion recombining
with electrons through the call to _w_product_at_catalytic_wall()
* in model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- in update_bdry_wall(), made w at wall node equal to w at near wall node no matter
the order of accuracy specified
- in update_bdry_wall(), at both injection and non injection surfaces added correction
to the mass fraction gradient of the neutral species which is the product of a positive
ion recombining with electrons through the call to _w_product_at_catalytic_wall()
* in model/fluid/_Navier-Stokes_plasma/fluid_post.c
- in find_post_variable_value_fluid(), added 4 new post variables related to cesium
(NOTE: THESE MUST BE DELETED BECAUSE THEY ARE USED FOR TESTING CESIUM INJECTION ONLY)
* in model/fluid/doc/Transport_Equations/report.tex
- in the section about surface mass injection, started documentation of how to specify
the boundary conditions for a non-neutral solver
* in model/metrics/_generic/metrics.c
- in find_unit_vector_normal_to_boundary_plane(), removed a fatal_error and rather return
a FALSE if n couldn't be found
- in _distance_between_near_bdry_node_and_boundary_plane(), added fatal_error in case n
couldn't be found
* in model/share/fluid_share.c
- added _w_product_at_catalytic_wall()
- in update_w_V_at_injection_wall(), added fatal_error in case n couldn't be found
* in model/thermo/_generic/thermo.c
- added find_neutral_spec_from_ion_spec()
- in find_nuk_from_w_rho_T_Te(), fixed a bug related to the determination of nu for the electrons
in the case of a neutral plasma: when the mass fraction of all positive ions was zero,
nu for the electrons couldn't be determined and this lead to NaN.
2021 SEP 16a
* in model/chem/_Cs_Air_Plasma_15s/chem.c
- added parentpark2021 option as a air chemical solver
- added CESIUMIONIZATIONMODEL option
- added CESIUMRECOMBINATIONMODEL option
* in model/chem/_Cs_Air_Plasma_15s/doc/report.tex
- added park1993 and parentpark2021 models
- added various cesium ionization and recombination models
* in model/metrics/_generic/metrics.c
- in find_unit_vector_normal_to_boundary_plane(), fixed a bug related to inner node B
being misaligned with the boundary plane
2021 SEP 15b
* in bin
- added rampinjector.wrp
* in model/fluid/doc/Transport_Equations/report.tex
- fixed some typesetting mistakes in the section on Wall Injection Boundary Conditions
* in model/metrics/_generic/metrics.c
- in find_unit_vector_normal_to_boundary_plane(), fixed a bug for the 3D version of the algorithm
- in find_unit_vector_normal_to_boundary_plane(), added a fatal_error should lB not be
perpendicular to the surface along the generalized coordinates
2021 SEP 15a
* in model/fluid/_Navier-Stokes/fluid_bdry.c and model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- added BDRY_WALLTFIXEDINJECTION1
- changed the call to update_w_at_catalytic_wall() to include paramstart and paramend
* in model/fluid/doc/Transport_Equations/report.tex
- added section for wall boundary condition with injection of species
* in model/metrics/_generic/metrics.c
- added find_unit_vector_normal_to_boundary_plane()
- added _distance_between_near_bdry_node_and_boundary_plane()
* in model/share/fluid_share.c
- in update_w_at_catalytic_wall(), now calls _distance_between_near_bdry_node_and_boundary_plane()
instead of calculating the non-normal distance between the near wall node and the wall node
- in update_w_at_catalytic_wall(), added the function parameters paramstart and paramend and
used the latter when calling _bdry_param()
- added update_w_V_at_injection_wall()
- made sure that update_w_at_catalytic_wall() is not compiled unless the fluid module has
multiple species
* in src/control.c
- in read_control(), added the constants specO2, specN2, specH2, etc.
2021 AUG 13a
* moved parent2021.c to parentpark2021.c
* in model/chem/_Air_Plasma_11s
- removed parent2021.h
* in model/chem/_Air_Plasma_11s/doc/report.tex
- deleted parent2021 model and updated the parentpark2021 documentation to reflect
the new parentpart2021.c code
2021 AUG 11a
* in model/chem/_Air_Plasma_11s/parent2021.c
- wrote the reactions as function of EXPON_TTV, EXPON_TTE, and EXPON_TVTE which are defined
at the beginning of the file.
- The constants EXPON_* are defined such that the reactions are function of only Te or T but
this can be changed easily
2021 AUG 08a
* in model/chem/_Air_Plasma_11s
- added the module parent2021.c
* in model/chem/_Air_Plasma_11s/parentpark2021.c
- in reactions 10, 11, 12, made the reverse controlling temperature equal to sqrt(Tv*Te)
which follows the logic of the Park model
* in model/chem/_Air_Plasma_11s/chem.c
- added module parent2021
* in model/chem/_Air_Plasma_11s/doc/report.tex
- changed the Parent-Park solver to reflect changes to the controlling temperature of the
reverse reactions 10, 11, and 12
- added the Parent 2021 solver
2021 AUG 06a
* in model/chem/_Air_Plasma_11s/doc/report.tex
- in the parentdunn table, changed the forward reactions 7, 13, and 16 to those in Boyd2007
* in model/chem/_Air_Plasma_11s/parentdunn2021.c
- changed reactions 7, 13, and 16 (forward rates) to those in Boyd2007
* in model/chem
- added _Air_Plasma_12s
- moved _Air_Plasma_8s_28r_Macheret to old
2021 AUG 04a
* in model/chem/_Air_Plasma_11s/doc/report.tex
- updated the Parent-Dunn model to include the two-body dissociative recombination rates
from Ajjay
* in model/chem/_Air_Plasma_11s
- added parentdunn2021.c and parentdunn2021.h
* in model/chem/_Air_Plasma_11s/chem.c
- added calls to parentdunn2021
- removed calls to parent2021
* in model/chem/_Air_Plasma_11s
- removed parent2021.c
2021 AUG 03a
* in model/chem/_Air_Plasma_11s/doc/report.tex
- finalized the documentation of the Parent-Park model with Ajjay's curve fits based on BOLSIG+
* in model/chem/_Air_Plasma_11s
- added parentpark2021.c and parentpark2021.h
* in model/chem/_Air_Plasma_11s/chem.c
- added the parentpark2021 model
* in model/chem/_Air_Plasma_11s/dunnkang1973.c
- changed TEMAX_TOWNSEND from 60000 to 1e99 K
* in model/chem/_Air_Plasma_11s/parent2021.c
- changed TEMAX_TOWNSEND from 100e6 to 1e99 K
2021 JUL 29a
* in model/chem/_Air_Plasma_11s/doc/report.tex
- updated documentation for Dunn-Kang chemical model to indicate whether the translational or
electron temperature is used in the reaction rates.
2021 JUL 23a
* in model/chem/_Air_Plasma_8s/rajendran2022.c
- fixed a bug related to the units of the pre-exponential factor of the Townsend ionization reactions
2021 JUL 17a
* in model/chem/
- added the Air Plasma 8 species model and updated the .config and Makefile to include the model
- removed the code associated with _Air_Plasma_8s_28r_Macheret since it has been implemented as
one of the modules in the Air Plasma 8 species model
* in model/share/chem_share.c
- Fixed a bug in add_to_dW_3r3p() and add_to_dW_3r2p() functions: multiplied the number density
of the third reactant for dWdTe, dWdTv and dWdT elements
* in src/control.c
- updated the reference for the Air Plasma 8 species Macheret (2007) chemical model
* in model/chem/_Air_Plasma_11s/doc/report.tex
- added documentation for Park-Parent chemical model
- updated documentation for Park-Dunn-Kang chemical model
2021 JUL 09b
* in model/share/fluid_share.c
- in update_w_at_catalytic_wall(), added error messages when gamma is not within physically
valid bounds
2021 JUL 09a
* in model/fluid/_Navier-Stokes/fluid_bdry.c
- implemented the BDRY_WALLTFIXEDCATALYTIC1 boundary condition
* in model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- implemented the BDRY_WALLTFIXEDCATALYTIC1 boundary condition
* in model/share/fluid_share.c
- added the function update_w_at_catalytic_wall()
* in src/control.c
- in read_control_functions(), added the treatment of the function _spec()
* in model/fluid/doc/Transport_Equations/report.tex
- added the section "Catalytic Surface Boundary Conditions"
2021 JUN 16a
* in model/fluid/_Navier-Stokes_plasma/fluid_bdry.c
- in write_bdry_fluid_template(), added boundary condition BDRY_WALLTFIXEDTHERMIONIC1
- added function _Vk_perp_to_bdry()
- in update_bdry_wall(), added parameter THERMIONIC
- in update_bdry_wall(), made modifications for thermionic BC
2021 JUN 12a
* in model/chem/_Air_Plasma_11s/dunnkang1973.c
- reaction 10 had the wrong rate: this has been fixed
* in model/chem/_Air_Plasma_11s
- added parent2021.c
* in model/chem/_Air_Plasma_11s/doc/report.tex
- added parent2021 model
- fixed rate of reaction 10 of Dunn Kang model
2021 JUN 09a
* in model/emfield/doc/Electrodynamics/report.tex
- added Lorenz condition and electrodynamic potential
2021 APR 15b
* updated config/removeproprietary.sh through make config
2021 APR 15a
* in model/chem/.config
- made the Air_Plasma_8s_28r_Macheret option non proprietary
* in model/chem/_Air_Plasma_8s_28r_Macheret/chem.c
- changed the license to BSD 2 clause
2021 APR 13a
* in src/control.c
- added reference farbar2013a
* in model/chem/_Air_Plasma_11s/chem.hh
- added reference farbar2013a
* in model/chem/_Air_Plasma_11s/chem.c
- added modifications for the Farbar model
* in model/chem/_Air_Plasma_11s/
- added farbar2013.c and farbar2013.h and made relevant modifications to Makefile
* in model/chem/_Air_Plasma_11s/doc/
- added the Farbar model to report.tex and corresponding reference to all.bib
2021 APR 06b
* in src/post.c
- in read_post_actions() function, fixed incorrect condition for updating postfilename due
to incorrect reference to filename pointer address instead of filename pointer
2021 APR 06a
* in model/share/emfield_share.c
- in _E_dot_J(), changed find_E(np[l],gl,E); to find_E_smoothed(np[l],gl,E);
* in src/test.c
- added test_dmuk()
* in model/fluid/_*_plasma/fluid_plasma.c
- added _mu_test()
- in find_H1() and find_H2(), added call to _mu_test instead of _mu when TEST is defined
- in find_dH1_dUstar() and find_dH2_dUstar(), updated the Jacobian terms present only
when TEST is defined
* in model/fluid/_*_plasma/fluid_source.c
- in find_Semfield(), now calls _E_dot_J() instead of the minmod between _E_dot_J() and
_E_dot_J_recast()
2021 APR 05a
* in cycle/share/cycle_share.c
- added header <src/post.h>
- in runtime_actions(), added code to enable Post() to be called from Cycle()
* in cycle/share/cycle_share.h
- added function prototype for runtime_actions()
* in src/post.c
- added header <cycle/share/cycle_share.h>
- in read_post_actions(), added condition to read postfilename only when -opm flag is used
- in read_post(), removed line which prevented Post() to be executed with more than one MPI process
- in read_post(), added condition to prevent newline to be printed when Post() is called within Cycle()
2021 MAR 27a
* in model/chem/_Air_Plasma_11s/chem.c
- added QEISOURCETERMS within the variables to read in Model(AirPlasma11s(););
* in find_Qei() and find_dQei_dx()
- set Qei to 0 if QEISOURCETERMS is set to FALSE
2021 MAR 23b
* in model/chem/_Air_Plasma_11s/boyd2007.c, park1993.c, lenard1964.c
- added find_Qei_*() and find_dQei_dx_*()
* in model/chem/_Air_Plasma_11s/dunnkang1973.c
- in find_Qei_DunnKang1973() and find_dQei_dx_DunnKang1973(), fixed a typo: 3.6e31/sqr(calA)
should be 3.6e31/calA and 1.1e32/sqr(calA) should be 1.1e32/calA
* in model/chem/_Air_Plasma_11s/chem.c
- renamed TOWNSENDIONIZATION to ADDITIONALREACTION
- changed the rates for ADDITIONALREACTION[5] and ADDITIONALREACTION[6] using BOLSIG+
* in model/chem/_Air_Plasma_11s/doc/report.tex
- changed the rates for the 5th and 6th additional reaction using BOLSIG+
2021 MAR 23a
* in read_data_file(), fixed a few bugs related to read_data_file_mpi when BDF is the chosen
time stepping method
2021 MAR 21a
* in src/common.h
- added dim_t Esmoothedmem to npbs_t
* in model/thermo/_generic/thermo.c
- in _ionizationpotential(), added the ionization energies of species O and N
* in model/share/emfield_share.c
- in _E_dot_Je_over_mueNe(), now calls find_E_smoothed() instead of find_E()
* in model/emfield/_*/emfield.c
- added find_E_smoothed() which returns Esmoothedmem
- in find_prim_emfield_mem_1(), added code to find Esmoothedmem
* in model/chem/_Air_Plasma_11s/dunnkang1973.c
- implemented clipping of Te so that it doesn't not exceed 60000 K when determining the
electron impact ionization reactions
- added find_Qei_DunnKang1973() and find_dQei_dx_DunnKang1973()
* in model/chem/_Air_Plasma_11s/chem.c
- added TOWNSENDREACTION[5] and TOWNSENDREACTION[6] for Townsend ionization of N and O function
of the reduced electric field
- in find_Qei(), added call to find_Qei_DunnKang1973()
- in find_dQei_dx(), added call to find_dQei_dx_DunnKang1973()
2021 MAR 19a
* in model/chem/_CO2_Plasma_23s_Park/chem.c
- in find_dW_dx(), fixed a bug in reactions 40 and 41 where dWdTv should be dWdTe
2021 MAR 17a
* in model/chem/_Air_Plasma_11s/chem.c
- in add_W_Townsend() and add_dW_dx_Townsend(), changed the chemical reaction for
NO+ + 2e- -> NO + e- to the one specified in Dunn-Kang (1973)
* in model/chem/_Air_Plasma_11s/doc/report.tex
- in the table "Additional Townsend Ionization and recombination rates for NO, O2.."
changed the chemical reaction for NO+ + 2e- -> NO + e- to the one specified in
Dunn-Kang (1973)
2021 MAR 15a
* in model/chem/_Air_Plasma_11s/chem.c
- added a fourth reaction to the Townsend reactions e- + NO -> 2 e- + NO+
* in model/chem/_Air_Plasma_11s/doc/report.tex
- added a fourth reaction to the Townsend reactions e- + NO -> 2 e- + NO+
2021 MAR 11a
* in model/chem/_H2_He_CH4_Plasma_13s
- added doc directory
* in model/chem/_H2_He_CH4_Plasma_13s/chem.c
- implemented Chris Johnston's Neptune atmosphere model
* in model/share/fluid_share.c
- in find_Saxi(), added the effective pressure to the total energy flux axisymmetric source
as it should be
2021 MAR 07a
* in model/chem/_CO2_Plasma_23s_Park/chem.c
- made modifications to the reactions so that the backward and forward temperatures
* in model/chem/_CO2_Plasma_23s_Park/doc/report.tex
- made modifications to the reactions so that the backward and forward temperatures
2021 FEB 27c
* in config/config.c
- in read_config_file(), made changes so that the options that are both EXPERIMENTAL
and PROPRIETARY are not hidden when doing make config or make configp
2021 FEB 27b
* in bin
- updated Starikovskii.wrp
2021 FEB 27a
* in model/fluid/_Favre-Reynolds_plasma/fluid.h
- added _kappae()
2021 FEB 22a
* in model/thermo/_generic/thermo.c
- in _kappa_from_w_T() and _eta_from_w_T(), fixed a bug introduced in 21feb21a related
to including the electron species
2021 FEB 21c
* in model/chem/_CO2_Plasma_23s_Park/chem.c
- fixed a bug: reaction 13 used N instead of X as a parameter
2021 FEB 21b
* in src/common.h
- in gl_t, added boolean variable OUTPUTBINARYMPI
* in src/data.c
- added header unistd.h
- added functions _ai_mpidatafile(),read_data_file_mpi(),write_data_file_mpi() to implement
parallel I/O file format WARPMPBFORMAT010
- modified functions write_data_file(),read_data_file() to implement parallel I/O file format
WARPMPBFORMAT010
* in src/data.h
- in input_t, added boolean variables BINARYMPI,M1,MM1,M2,MM2,M3,MM3
* in src/main.c
- renamed -im1,-im2,-im3 flags to -i1,-i2,-i3
- added code to implement parallel I/O fileformat to be used via -im and -om flags
- repurposed -im1,-im2,-im3 flags to read time level minus 1, time level minus 2 and time
level minus 3 data of type BINARYMPI
2021 FEB 21a
* in model/fluid/_*_plasma/fluid.c
- in _kappa() and find_prim_fluid_mem(), now calls _kappah() to include only the heavy species
- added _kappai()
- in find_Kstar_interface(), added kappai to kappa for Kstar[fluxet][fluxet]
* in model/fluid/_*_plasma/fluid_plasma.c
- in find_Dstar(), took away the charged species thermal diffusion from the velocity
(now added to find_Kstar_interface)
* in model/fluid/_Favre-Reynolds_plasma/fluid_post.c
- added kappae and kappav to the post variables
* in model/fluid/_Navier-Stokes_plasma/fluid_post.c
- added kappai to the post variables
* in model/fluid/doc/Transport_Equations/report.tex
- in the charged species energy transport equation section, made changes and took away thermal
diffusion (already included within pressure gradients)
* in model/share/emfield_share.c
- added _dTkdXj()
* in model/thermo/_generic/thermo.c
- added _etak_from_T()
- added _kappak_from_T()
- added _kappan_from_w_T(), _kappa_from_w_T(), _kappaw_from_w_T()
- added _eta_from_w_T()
- added find_nuk_from_w_rho_T_Te()
- rewrote find_nuk_eta_kappa() so that it calls _kappa_from_w_T(), _eta_from_w_T(),
find_nuk_from_w_rho_T_Te()
- electrons are now excluded when finding kappa, eta, and nu for the mixture
2021 FEB 18a
* in model/thermo/_generic/thermo.c
- in find_nuk_eta_kappa(), took away contributions of charged species to thermal conductivity
kappa because such are already accounted for through the enthalpies
2021 FEB 17a
* renamed
- find_Dstarplus() -> find_Dstarmatplus()
- find_Dstarminus() -> find_Dstarmatminus()
- find_Dstar() -> find_Dstarmat()
* in cycle/res/_generic/res.c
- renamed add_dDstarU_residual_FVS() to add_dDstarmatU_residual_FVS()
- added new add_dDstarU_residual_FVS() where Dstar is a vector not a matrix and with metricsint
given as a parameter to find_Dstar
* in cycle/share/ts_share.c
- renamed find_DstarUstar_jacobians_FVS() to find_DstarmatUstar_jacobians_FVS()
- added new find_DstarUstar_jacobians_FVS() where Dstar is a vector not a matrix and with metricsint
given as a parameter to find_Dstar
* in model/fluid/_Favre-Reynolds_plasma/fluid_plasma.c
- implemented new find_Dstar*() functions and erased previous ones
* in model/fluid/_Navier-Stokes_plasma/fluid_plasma.c
- added new find_Dstar(), find_Dstarplus(), and find_Dstarminus() which accept metricsint as an
argument
* in model/fluid/_drift-diffusion/fluid.c
- added new find_Dstar(), find_Dstarplus(), and find_Dstarminus() which accept metricsint as an
argument
* in model/share/emfield_share.c
- in _dTedXj(), now calls _Tk() instead of _Te()
2021 FEB 15a
* in cycle/res/_generic/res.c
- in add_dDstarU_residual_FDS(), changed TVD limiter to FIRSTORDER
- in add_dDstarU_residual(), now calls add_dDstarU_residual_FVS
* in cycle/share/ts_share.c
- in find_DstarUstar_jacobians(), now calls find_DstarUstar_jacobians_FVS()
* renamed
- _dTedXj_interface -> _dTedXj_interface_upwind
- _dTedXj -> _dTedXj_interface
* in model/fluid/_*_plasma/fluid_plasma.c
- in find_Dstar(), added code to take away thermal diffusion from electron velocity
* in model/fluid/_Navier-Stokes_plasma/fluid_plasma.c
- in find_Dstar_interface(), made minor modifications to averaging of mu and Ve at interface
* in model/fluid/_drift-diffusion/fluid.c
- added find_Dstar_interface()
* in model/share/emfield_share.c
- added _dTedXj()
2021 FEB 12b
* in model/fluid/_*_plasma/fluid_plasma.c
- in _DeltaVkstar_interface(), changed arith average to minmod
2021 FEB 12a
* in cycle/res/_generic/res.c
- renamed add_dDstarU_residual() to add_dDstarU_residual_FVS()
- added add_dDstarU_residual_FDS()
- added add_dDstarU_residual() which calls add_dDstarU_residual_FDS()
* in cycle/share/res_share.c
- in _limiter_TVD(), added LIMITER_FIRSTORDER
* in cycle/share/ts_share.c
- renamed find_Dstar_interface() to find_DstarUstar_jacobians_FVS()
- added find_DstarUstar_jacobians_FDS()
- added find_DstarUstar_jacobians()
- in add_Dstar_to_TDMA_check() and find_TDMA_jacobians_conservative(), now calls
find_DstarUstar_jacobians() instead of find_Dstar_interface()
* in model/emfield/_Epotential-Te/emfield.c
- in add_dUstar_to_U_emfield(), added the addition of Temin_EMF and Temax_EMF to the clipped
variables
- in find_Ve_drift_at_interface(), took away thermal diffusion contribution to electron velocity
- in find_linearization_coefficients_inner_node_emfield_OHMBASED(), set the e thermal diffusion
coefficients to 2 times their normal value instead of 1.25 times
- in find_Fstar_interface_emfield(), added to Fstar the full electron thermal diffusion terms
- in find_Fstar_interface_emfield(), changed the limiter from VANLEER to MINMOD in determining
TeL and TeR
* in model/fluid/_*_plasma/fluid.c
- in find_Kstar_interface(), changed the electron thermal diffusion coefficients to 2 times
the normal value for CYCLELEVEL_TS (up from 1.25) and to the normal value for CYCLELEVEL_RES
(up from 0.25)
* in model/fluid/_*_plasma/fluid_plasma.c
- added _DeltaVkstar_interface()
- added find_Dstar_interface()
* in model/fluid/_Favre-Reynolds_plasma/fluid_plasma.c
- in find_Dstar(), changed find_Ve_from_Je(np, gl, l, Ve); to find_Ve_from_J(np, gl, l, Ve);
2021 FEB 05a
* in model/emfield/doc/Electrodynamics/report.tex
- improved document
* in model/emfield/_*/emfield.c
- in read_model_emfield_actions(), now read in AXISYMMETRIC only is _2D is defined
* in model/emfield/_*/emfield.hh
- in gl_model_emfield_t, AXISYMMETRIC is defined only if _2D is defined
2021 FEB 04a
* in model/emfield/doc/
- added directory Electrodynamics
2021 FEB 02a
* in lib/soap/soap.c
- in BA_for_parallel(), fixed a bug that caused the incorrect update of variables when using
the for_parallel built-in action
2021 JAN 29a
* in bin
- updated Starikovskii.wrp
* in model/emfield/_*/emfield.c
- in write_model_emfield_template(), added a newline after AXISYMMETRIC=NO;
* in model/fluid/_Navier-Stokes_perfect/fluid.c
- in read_model_fluid_actions(), added new variables ETAMODEL, eta_ref, eta_ref_T, eta_C, Pr
and removed kappa
- in find_Kstar_interface(), now takes average of eta and kappa at interface
- added function _kappa(), _eta_from_T(), _eta()
* in model/fluid/_Navier-Stokes_perfect/fluid_bdry.c
- in write_bdry_fluid_template(), added BDRY_INFLOWFLUCTUATING
- added find_V_P_T_bdry_fluctuating(), update_bdry_inflow_fluctuating()
* in model/fluid/_Navier-Stokes_perfect/fluid_conv.c
- in find_jacvars(), now finds jacvars->eta from _eta(np,gl) and not from gl->model.fluid.eta
- in find_jacvars_from_musclvars() and find_jacvars_from_U(), now finds jacvars->eta from
_eta_from_T(T,gl) and not from gl->model.fluid.eta
* in model/fluid/_Navier-Stokes_perfect/fluid_post.c
- in find_post_variable_name_fluid() and find_post_variable_value_fluid(), added kappa and eta
to post variables
2021 JAN 22a
* in model/chem/_Air_Plasma_11s/dunnkang1973.c and model/chem/_Air_Plasma_11s/doc/report.tex
- modified the Dunn-Kang model: backward reactions are now explicitly written as in the
original paper instead of using the equilibrium constant
* in model/fluid/_Navier-Stokes/fluid_source.c
- added missing #include fluid_share.h
* in model/fluid/_Navier-Stokes_perfect/fluid_source.c
- removed functions find_Saxi() and find_dSaxi_dU()
- added missing #include fluid_share.h
* in model/fluid/_drift-diffusion/fluid.c
- in write_model_fluid_template() and read_model_fluid_actions(), added code to read
AXISYMMETRIC from control file
- in find_Sstar() and find_dSstar_dUstar(), added calls to find_Saxi and find_dSaxi_dU
2021 JAN 21a
* in model/emfield/_Epotential*/emfield.c
- read_model_emfield_actions(), write_model_emfield_template(), added AXISYMMETRIC flag
- in find_Sstar_emfield(), added axisymmetric source terms
* in model/fluid/_Navier-Stokes/fluid_source.c and model/fluid/_Navier-Stokes_plasma/fluid_source.c
- deleted find_Saxi() and find_dSaxi_dU()
* in model/fluid/_Navier-Stokes_plasma/fluid_source.c
- in find_Vk(), added an assert(spec<ncs)
* in model/share/fluid_share.c
- in _Qk(), turned off the axisymmetric contributions by Jason Etele
- deleted _Qk_old()
- renamed find_Saxi() to find_Saxi_FavreReynolds()
- added find_Saxi() and find_dSaxi_dU() to be used by all fluid modules
2021 JAN 19b
* in src/control.c
- in read_control(), added "EMFIELD" variable to control file variables
* in model/fluid/_*_plasma/fluid.c and model/fluid/_drift-diffusion/fluid.c
- in write_model_fluid_template(), changed Temin=Tmin; to Temin=300.0;
* in bin/RAM-C-II.wrp
- made modifications so that SteadyState() can be read either by emfield or non-emfield
solvers
2021 JAN 19a
* in src/control.c
- in find_metrics_on_all_nodes(), deleted the variable newzone (not needed anymore)
* in model/chem/_Air_Plasma_11s/dunnkang1973.c
- made reactions with electrons in reactants function of Te not of T
2021 JAN 18a
* in config/config.c
- in output_makefileheader() added -rdynamic flag when using debugging symbols and gcc or mpicc
* in lib/exm/exm.c
- added output_backtrace() when linux is defined
* in lib/gridg/*.h and src/common.h
- removed the include assert.h
* in lib/soap/soap.c
- removed the include assert.h
- added the include exm.h
* in src/common.c
- removed backtrace() function
* in src/control.c
- in find_metrics_on_all_nodes(), fixed a bug: the metrics at the interface were not properly
stored in RAM for all interfaces with MPI, and this lead to the current density not being
determined properly on some nodes when using MPI. Note that this bug only occurred when using
MPI and using a emfield module other than none.
2021 JAN 16a
* in cycle/share/cycle_share.c
- in update_linked_nodes(), fixed a bug: _al(gl,i,j,k) should be _ai(gl,i,j,k)
* in cycle/tsemf/_AF/tsemf.c
- replaced in a few locations _al(gl,l,dim,+0) by l
- added assert(dim<nd) in a few locations
* in model/chem/_Cs_Air_Plasma_11s/chem.c
- added option CESIUMIONIZATIONMODEL within write_model_chem_template(),
read_model_chem_actions(), add_dW_dx_Cesium(), add_W_Cesium(): now the ionization due to
electron impact (and its backward reaction) can be set to either Lenard (1964)
or Andrienko (2021)
* in model/chem/_Cs_Air_Plasma_11s/doc/report.tex
- added the Andrienko (2021) cesium ionization and recombination reactions
* in model/fluid/_Navier-Stokes_plasma/fluid.c
- in read_model_fluid_actions(), changed "1" for TEMODEL_CONSTANT, "2" for TEMODEL_TEQUILIBRIUM,
etc
* in src/common.c, added the function output_backtrace() but left it commented so that it does
not lead to compilation problems on non-Linux systems
2021 JAN 07a
* in lib/soap/soap.c
- in functions_builtin(), got rid of the spline function
- in functions_builtin(), renamed "Interpolate" to "interpolate"
* in lib/soap/doc/report.tex
- deleted spline function and renamed Interpolate to interpolate
2021 JAN 06a
* in lib/soap/doc/report.tex
- added description of function Interpolate
2021 JAN 04a
* in src/test.c
- in test_s(), test_cp(), test_hmolar(), test_h(), initialized the species mass fractions
to 1e-10 instead of 0
- added test_eta(), test_kappa(), test_Pr()
* in src/control.c
- in write_license(), changed the copyright year to 1998-2021
* in model/thermo/_generic/thermo.c
- got rid of MAX_T_FOR_NUK_ETA_KAPPA_POLYNOMIALS=6000.0 and replaced with TSTAR_MAX=100.0
- in _Omega11() and _Omega22(), limited Tstar to TSTAR_MAX
- in find_nuk_eta_kappa(), got rid of T limit to MAX_T_FOR_NUK_ETA_KAPPA_POLYNOMIALS
- in find_nuk_eta_kappa() commented line "if (w[spec]<1.0E-12) w[spec]=1.0E-9;"
* in model/emfield/_Epotential*/emfield.c
- commented out _Tk_smooth()
- in _mu_mem(), now calls _smooth2(0 instead of _Tk_smooth
* in src/common.c
- in _smooth(), added a fatal_error in case node l is not an inner node
- added function _smooth2()
2020 DEC 29a
* in model/emfield/_Epotential-Bfixed/emfield.c
- changed maxiterEk from 1000 to 2000
- took away constant EFIELD_SPECIES_FRAME and incorporated it within the control file
- commented function _Ekmag() because it is not used any more
- in find_Ek_muk_mem(), implemented a faster and more robust algorithm to find the local
electric field which tries various relaxation factors (instead of only one as before) before
quitting
- in find_Ek_muk_mem_nonlocal(), fixed a mistake:
Tk=max(Tk,_T(np[l],gl)+gl->model.emfield.Te_buffer) should be
Tk=max(Tk,gl->model.emfield.Te_buffer)
* in model/share/emfield_share.c
- in find_Ve_from_J(), added a call to fatal_error in case no inner node can be found near a
boundary node
* in src/common.c
- in _al(), added an assert to verify that theta is less than nd
- in _all(), fixed a bug: ltmp should have a type of long not double
- added function _alll()
- in find_l_of_nearest_inner_node(), added a fourth check to find an inner node in case the bdry
node is in a 3D corner
* in model/chem/_CO2_Plasma_23s_Park/chem.hh
- renamed #define _CHEM_CSAIRPLASMA to #define _CHEM_CO2PLASMA23S
2020 DEC 27a
* in model/chem/_Cs_Air_Plasma_*s/chem.c
- in find_Qei() and find_dQei_dx(), added missing definition and initialization of R
* in model/share/fluid_share.c
- in write_disc_fluid_template(), changed betag[spec]=-0.5; to betag[spec]=-0.01;
2020 DEC 26a
* in model/chem/_Cs_Air_Plasma_11s/chem.c and model/chem/_Cs_Air_Plasma_15s/chem.c
- in find_Qei() and find_dQei_dx() added contributions from CESIUMREACTION[1]
* in model/chem/_H2_Air_Plasma_*/chem.c and model/chem/_H2_Air_Plasma_*/doc/report.tex