Ember Source problem

71 views
Skip to first unread message

anuragsu...@gmail.com

unread,
Feb 11, 2021, 7:48:23 AM2/11/21
to Cantera Users' Group
Dear Ember user, 

I plan to use ember to tabulate unsteady flamelets, I start from a pure mixing solution and record the flamelets in time.  This is ember configurations i use 
################################################
 output = 'save/'
 conf = Config(
     Chemistry(mechanismFile='Base_OMEs_PAH.xml',
               transportModel='UnityLewis',
               kineticsModel='standard'),
     Paths(outputDir=output,
           # logFile='ex_diffusion.log'
           ),
     InitialCondition(flameType='diffusion',
                      restartFile='./prof000103.h5',  # Pure mixing solution
                      fuel='NC12H26:1',
                      oxidizer='N2:0.85, O2:0.15',
                      Tfuel=363,
                      pressure=6000000,
                      Toxidizer=900,
                      xLeft=0.00,
                      xRight=0.04,
                      centerWidth=0.002,
                      slopeWidth=0.001),
     StrainParameters(initial=500,
                      final=500),
     General(nThreads=48,chemistryIntegrator='qss',errorStopCount=1000,
             splittingMethod='strang'),
     Grid(vtol=0.13, dvtol=0.14, gridMin=5e-7, gridMax=0.0001),
     Times(globalTimestep=1.0e-8,
           profileStepInterval=300,
           regridStepInterval=600),
     TerminationCondition(tEnd= 0.8,measurement=None),
     OutputFiles(fileExtension = 'h5',firstFileNumber=0))
 if __name__ == '__main__':
     conf.run()

##############################################################

This is what the flamelets look in temperature 
temp_Z.png
and this is what I have as progress variable (combination of CO2,CO,H2 and H2O) vs Z in time.yc_z.png
However when I compute chemical sources from the functionalities provided in Utils.py of ember i get fluctuations and peaks. Z_Yc_OmegaC.png
a better view of the problem can be this 
omegaYc.png
So I wanted to know if I am doing something wrong or there could be some other way to avoid these spikes 

Thanks and Regards 
Reply all
Reply to author
Forward
0 new messages