VFLUX: Vertical Fluid Heat Transfer Solver (VFlu[H]X Solver)

VFLUX is a program that calculates one-dimensional vertical fluid flow (seepage flux) through saturated porous media, using heat transport equations. It uses temperature time series data measured by multiple temperature sensors in a vertical profile in order to calculate flux at specific times and depths. VFLUX is written as a MATLAB toolbox, a set of functions that run in the MATLAB environment.

Please cite as: (Gordon et al., 2012; Irvine et al., 2015)
- Irvine, DJ, LK Lautz, MA Briggs, RP Gordon, JM McKenzie. 2015. Experimental evaluation of the applicability of phase, amplitude, and combined methods to determine water flux and thermal diffusivity from temperature time series using VFLUX 2. Journal of Hydrology, 531(3):728-737. doi:10.1016/j.jhydrol.2015.10.054.
- Gordon, RP, LK Lautz, MA Briggs, JM McKenzie. 2012. Automated calculation of vertical pore-water flux from field temperature time series using the VFLUX method and computer program. Journal of Hydrology, 420-421:142-158. doi: 10.1016/j.jhydrol.2011.11.053.
- (Two minor typographical errors were recently discovered in the publication Gordon et al. (2012). The definition for H in Equation 12 (p. 147) should be corrected to read “H = Cw/λo”. Also, Equation 9 (p. 147) should read “κe = (λo/C) + β|v|”, where v is the thermal front velocity. Please note that the code of the VFLUX program has always contained the correct forms of these two equations. Thanks go to Chengpeng Lu of Hohai University and Dylan Irvine of Monash University for bringing these typographical errors to our attention.)

VFLUX 1.2.5: The original VFLUX code, with amplitude- and phase-based solutions to determine flux from temperature time series data, as well as signal processing methods, data visualization, sensitivity analysis, and Monte Carlo error analysis modules.

VFLUX 2.0.0: Includes all functionality of the original VFLUX, with the addition of solutions for the “combined” amplitude ratio and phase lag methods (Luce et al., 2013; McCallum et al., 2012). More information can be found in Irvine et al. (2015).

VFLUX requires the Captain Toolbox and the ‘resample’ function from the MATLAB Signal Processing Toolbox. The Captain Toolbox is available for free at http://www.lancs.ac.uk/staff/taylorcj/tdc/download.php. The Signal Processing toolbox is available from MathWorks at http://www.mathworks.com/products/signal/. If the Signal Processing Toolbox is not available to you, then an alternate function may be substituted (for more information, contact the authors).

VFLUX Add-ons
(1) vflux_qar_opt is an add on program for VFLUX2 (v 2.0.0 and greater). vflux_qar_opt has two main applications: 1) to fine tune flux estimates using the benefits of two analytical solutions, and 2) to provide a workflow to assess potential streambed scour. vflux_qar_opt can be run after the main vflux functions, where the user can optimize flux estimates by refining the thermal diffusivity (κe) value that best reproduces the known sensor spacing (Δz). More information can be found in the vflux_qar_opt documentation, which is included with the m-file in the download, and associated manuscript (Irvine et al. 2017).
Irvine, DJ, MA Briggs, I Cartwright, CR Scruggs, LK Lautz. 2017. Improved vertical streambed flux estimation using multiple diurnal temperature methods in series. Groundwater, 55(1): 73-80. doi: 10.1111/gwat.12436

VFLUX是一款基于热传输方程计算饱和多孔介质中一维垂直流体流动（渗透通量）的程序。该程序通过分析垂直剖面中多个温度传感器所测量的温度时间序列数据，以计算特定时间和深度的通量。VFLUX以MATLAB工具箱的形式编写，是一组在MATLAB环境中运行的函数集。

引用格式：[Gordon et al., 2012; Irvine et al., 2015]

- Irvine, DJ, LK Lautz, MA Briggs, RP Gordon, JM McKenzie. 2015. 相位、振幅及综合方法在确定水通量和热扩散率中的应用实验评估——以VFLUX 2.0为例. 水文，531(3):728-737. doi:10.1016/j.jhydrol.2015.10.054.
- Gordon, RP, LK Lautz, MA Briggs, JM McKenzie. 2012. 基于VFLUX方法和计算机程序自动计算现场温度时间序列的垂直孔隙水通量. 水文，420-421:142-158. doi: 10.1016/j.jhydrol.2011.11.053.
（近期在Gordon et al. (2012)的出版物中发现了两项细微的印刷错误。方程12（第147页）中H的定义应更正为“H = Cw/λo”。方程9（第147页）应读作“κe = (λo/C) + β|v|”，其中v为热前沿速度。请注意，VFLUX程序的代码始终包含这两个方程的正确形式。感谢河海大学的Chengpeng Lu和莫纳什大学的Dylan Irvine指出这些印刷错误。）

VFLUX 1.2.5：原始的VFLUX代码，包括基于振幅和相位的解法来确定温度时间序列数据的通量，以及信号处理方法、数据可视化、敏感性分析和蒙特卡洛误差分析模块。

VFLUX 2.0.0：包含原始VFLUX的所有功能，并新增了“综合”振幅比和相位滞后方法的解法（Luce et al., 2013; McCallum et al., 2012）。更多详细信息可参考Irvine et al. (2015)。

VFLUX需要Captain Toolbox和MATLAB信号处理工具箱中的‘resample’函数。Captain Toolbox可免费从http://www.lancs.ac.uk/staff/taylorcj/tdc/download.php下载。信号处理工具箱可从MathWorks（http://www.mathworks.com/products/signal/）获取。如果您无法访问信号处理工具箱，则可以替换为其他函数（更多信息请联系作者）。

VFLUX附加组件
(1) vflux_qar_opt是VFLUX2（v 2.0.0及更高版本）的附加程序。vflux_qar_opt具有两个主要应用：1）利用两种解析解的优势，微调通量估计；2）提供一项工作流程，以评估潜在的河床冲刷。vflux_qar_opt可以在主vflux函数运行之后运行，用户可以通过优化热扩散率（κe）值，该值最佳地再现已知的传感器间距（Δz），以细化通量估计。更多详细信息可在vflux_qar_opt文档中找到，该文档包含在下载的m文件中，并附有相关论文（Irvine et al. 2017）。
Irvine, DJ, MA Briggs, I Cartwright, CR Scruggs, LK Lautz. 2017. 利用串联的多个日温度方法改进垂直河床通量估计. 地下水，55(1): 73-80. doi: 10.1111/gwat.12436