Jan Mandel/Blog/2011 Feb Mar
February 5: WRF-Fire paper for GMD
Jan Mandel, Jonathan D. Beezley, and Adam K. Kochanski, Coupled atmosphere-wildland fire modeling with WRF-Fire, UCD CCM Report 295, February 2011. Geoscientific Model Development (GMD), submitted. Also available as arXiv:1102.1343.
February 7: Data assimilation seminar: IEEE Prague paper
- Jon leads
- For 6th IEEE International Conference on Intelligent Data Acquisition and Advanced Computing Systems: Technology and Applications (IDAACS'2011), September 15-17, 2011, Prague, Czech Republic
- Abstract submission Feb 27, 2011, From Author Instructions: Prospective authors are invited to submit an extended abstract of 3-4 pages (in English only). The manuscript should define the scope of the paper, emphasize on new advances, theories and/or applications and include an analysis of results and findings. It must be structured such that the program committee will be able to understand the originality and the value of the work. It is not clear if the final camera-ready paper is going to be the same.
- Camera ready paper submission Jun 15, 2011
- Proceedings: IEEE
- Reference: Jan Mandel, Jonathan D. Beezley, and Loren Cobb, Spectral and morphing ensemble Kalman filters, abstract presentation files draft paper directory (will not continue, this paper and the next one supersede the draft)
February 8: Jon's matlab osimorph classes
February 9: Revisit registration and morphing
- Now split off the morphing part pdf files
- See also AMS 2011 EnKF presentation and osimorph notes
- Laroche-1994-VAM Marzban-2010-OFV displacement correction review notes
February 14: CCM Colloquium: Olivier Desjardins
- Olivier Desjardins, A Discontinuous Galerkin Conservative Level Set Method for Simulating Turbulent Atomization
- Olsson and Kreis, 2005 Olsson et al, 2007 Desjardins el at, 2008
February 18: My colloquium at Applied Math, CU Boulder
- Presentation "Coupled atmosphere-wildland ﬁre numerical simulation and data sources" pdf files See also
- Papers: 2011 preprint at arXiv:1102.1343, 2009 paper in IEEE Control Systems Magazine, also available as arXiv:0712.3965.
- Presentations of our group at AGU Fall Meeting, 2010
- A previous presentation at IMA, December 2010
- WRF-Fire software at openwfm.org
- YouTube channels openwfm and wildfireproject
February 21: Data assimilation seminar: Nonparametric and wavelet methods
- Nonparametric statistics (Loren)
- Wavelet estimation of covariance (Jon) photos
February 24: WRF-Fire: Fuel height
- Finishing up vertical interpolation of the wind to different heights for different fuels
- Code is in branch jm2/fuelheight. See the February 24, 2011 commit
- How to get midflame height from wind reduction factors pdf files
February 27: Wavelet EnKF paper submitted
February 28: Data assimilation seminar: Probability measures - Basic concepts and types of convergence of random variables
- Convergence of random variables from Billingsley's book.
- My new notes on random variables and old notes on probability
- Recall our paper On the Convergence of the Ensemble Kalman Filter, CCM Report 278, January 2009. arXiv:0901.2951 Slightly expanded version accepted in Applications of Mathematics
March 4: Abstract for ISF 2011
- The 31st Annual International Symposium on Forecasting, June 26-29, 2011, Prague, The University of Economics
- submitted abstract
March 5: Bedrich's optical flow bibliography search
March 7: Data assimilation seminar: Loren Cobb: Ergodic theory
March 8: Proofs of the GMD paper
March 9: WRF-Fire: GMDD paper published
- GMD paper published
- wind interploation
- Merging branch jm2/vktest2 branch vk/test2 branch master
- Harmanli fire paper
- Janus cluster
March 13: Finishing the Harmanli paper
March 14: Data assimilation seminar: Jon Beezley: Kalman filter in Hilbert space and white noise
March 16: WRF-Fire
March 20: WRF-Fire: merging improved heat flux code
- The improved ignition time interpolation is described on page 7 of the AMS 2011 paper
- Volodymyr's code for ignition time interpolation
- Made first various changes to enforce state consistency: lfn>=0 implies tign>=time_now. Cannot exclude tign=time_now and lfn<0 because of rounding errors. No eps anywhere in the tign computations.
- cosmetic changes to Volodymyr's code completed
- Removes spurious peaks in heat flux on the fire grid
- edited, no substantial changes yet
- fuel_left_method=2 should be exact if lfn and tign are linear (allowing for the fireline). Implemented in subroutine fuel_left_cell_2, my Matlab prototype is standalone/matlab/fuel_burnt.m
fixed: but flux has hotspots:
probably because method 2 is not monotonous in tign. Max difference in fuel fraction between method 1 and 2 cca 11%.
- merged into master (fast forward), but do not use method 2 yet
- Is method 1 is monotonous in tign? It would be enough to show that tign interpolation is monotonous.
- Next: Could there be spurious jumps in flux if the (tign,lfn) -> flux mapping is not smooth?
- This all needs to be resolved before implementing canopy fire, because spurious heat flux hotspots will cause spurious canopy ignition.
- More problems in fuel calculation found: fuel_left was 0.25 on domain boudary after first time step, causing heat released from 75% of the fuel in the first time step along domain boundary. Reason: should be divided to average only where computed. fuel_left can increase due to approximation errors, causing large negative heat flux spikes. Also, fire_area not initialized on boundary when computed, relying on default initialization in the caller instead. Fixed Sun, 27 Mar 2011. Compare with the "After at 10s" above and note different scale, 105 not 106.
- Mar 28, 2011 tested and merged into master (fast forward)
March 22: "WRF-Fire Applied in Bulgaria" paper appeared
- Nina Dobrinkova, Georgi Jordanov, and Jan Mandel, WRF-Fire Applied in Bulgaria, Numerical Methods and Applications, I. Dimov, S. Dimova, and N. Kolkovska, eds., Lecture Notes in Computer Science 6046, Springer, Berlin/Heidelberg, 2011, pp. 133-140. Proceedings of 7th International Conference on Numerical Methods and Applications - NM&A'10, August 20 - 24, 2010, Borovets, Bulgaria. doi:10.1007/978-3-642-18466-6_15. Preprint arXiv:1007.5347, July 2010.
March 24: WRF-Fire: Perimeter ignition
- Input ignition times over whole domain from file input_tign_g, replay fire until fire_perimeter_time, atmosphere follows, then fire simulation takes over
- The ignition times serve also as as the initial level set function
- Test problem: cd WRFV3/test/em_fire/hill, see file README.txt
- merged into master