|
From: Steve P. <sj...@gm...> - 2018-11-12 15:56:19
|
Here is your earlier Q, also CCing Michael. Steve On Thu, Nov 8, 2018 at 12:05 PM Runco, John J Capt USAF AETC AFIT/ENY < Joh...@af...> wrote: > Chonglin/Steve, > > > I am currently still getting electron densities higher than the > experimental data indicates (about 3 orders of magnitude). I think I am > coding the fix ambipolar correctly, but was hoping to get a second pair of > eyes on it to see if I was missing anything. Attached is my input dec. > Highlighted in yellow are the lines that have to do with the fix ambipolar > commands. I took the manual's suggestions as well and made 2 different > mixtures, one with the electron and one without. I used a fix in/emit for > the mixture without the electrons, and subsequently a collide command for a > mixture with all the species. > > > Mistakes that I think are possible: > > 1. incorrect coding for the creation of two different groups (ion and > electron) within the mixture containing all species...as the Sparta manual > suggests?: > > mixture air_wo_ID NO+ e group self > > > I am currently trying to run a job using two separate lines to distinguish > the electron and ion (NO+) within the group that contains all the species: > > mixture air_wo_ID NO+ group ion_ID > mixture air_wo_ID e group elec_ID > > > 2. not varying the initial electron velocity used in the mixture command > seen below? (don't think this is the issue due to the air.species and > air.vss and air.tce files, but it came to mind). > species air.species N2 O2 N O NO NO+ e > mixture air_wo_ID N2 O2 N O NO NO+ e vstream 7650. 0. 0. temp > 196.7 > > Appreciate the help, > John > Sparta User > > > > ________________________________ > From: Chonglin Zhang <zha...@gm...> > Sent: Thursday, October 11, 2018 2:05 PM > To: Runco, John J Capt USAF AETC AFIT/ENY > Cc: spa...@li... > Subject: Re: [sparta-users] Fix Ambipolar Question > > Yes, if you are modeling electron density, then you need to use the "fix > ambipolar" command, together with "collide_modify" command. > > You will need to do something like this: > > collide_modify vibrate smooth ambipolar yes # active the > ambipolar approximation for ionization reactions. > fix ID ambipolar e NO+ N+ O+ N2+ O2+ # set which species > is electron, which species are ions; ID is the fix id number. > > > Thanks! > Chonglin > > On Thu, Oct 11, 2018 at 1:42 PM Runco, John J Capt USAF AETC AFIT/ENY < > Joh...@af...<mailto:Joh...@af...>> wrote: > If I am looking to verify electron densities, would accuracy also depend > on incorporating the fix ambipolar command? > > > Thanks, > > John Runco > > > > > _______________________________________________ > sparta-users mailing list > spa...@li...<mailto: > spa...@li...> > https://lists.sourceforge.net/lists/listinfo/sparta-users > |