Contaminate transport in vadose is a huge concern since the vadose zone is the main passage way for ground water recharge. Understanding this process is crucial in order to prevent contamination, protect and rehabilitate ground water resources. Reactive transport models are instrumental for such purposes and there are numerous solute transport simulation programs for both ground water and vadose zone but most of this models are limited to simple Linear, Langmuir and Freundlich sorption models and first order decay and fail to simulate more complex geochemical reactions that are common in the vadose zone such as cation exchange, surface complexation, redox reaction and biodegradation. So it is necessary to enhance capabilities of solute transport models by incorporating well tested hydrogeochemical models like PHREEQC in to them to be able closely approximate the geochemical transport process in the subsurface. In this PhD research a new reactive transport model called VS2DRT was created by coupling existing public domain solute and heat transport models VS2DT, VS2DH with hydro-chemical model PHREEQC using non-iterative operator splitting technique. VS2DRT was compiled using MinGW compiler using tools like autotools and automake. A graphical user interface was also created using QT creator and Argus ONE numerical development tools. The new model was tested for one dimensional conservative Cl transport, surface complexation, cation exchange, dissolution of calcite and gypsum, heat and solute transport as well as for two dimensional cation exchange cases. Their results were compared with VS2DT, VS2DH, HP1 and HP2 models and the results are in good agreement.

Contaminate transport in vadose is a huge concern since the vadose zone is the main passage way for ground water recharge. Understanding this process is crucial in order to prevent contamination, protect and rehabilitate ground water resources. Reactive transport models are instrumental for such purposes and there are numerous solute transport simulation programs for both ground water and vadose zone but most of this models are limited to simple Linear, Langmuir and Freundlich sorption models and first order decay and fail to simulate more complex geochemical reactions that are common in the vadose zone such as cation exchange, surface complexation, redox reaction and biodegradation. So it is necessary to enhance capabilities of solute transport models by incorporating well tested hydrogeochemical models like PHREEQC in to them to be able closely approximate the geochemical transport process in the subsurface. In this PhD research a new reactive transport model called VS2DRT was created by coupling existing public domain solute and heat transport models VS2DT, VS2DH with hydro-chemical model PHREEQC using non-iterative operator splitting technique. VS2DRT was compiled using MinGW compiler using tools like autotools and automake. A graphical user interface was also created using QT creator and Argus ONE numerical development tools. The new model was tested for one dimensional conservative Cl transport, surface complexation, cation exchange, dissolution of calcite and gypsum, heat and solute transport as well as for two dimensional cation exchange cases. Their results were compared with VS2DT, VS2DH, HP1 and HP2 models and the results are in good agreement.

Non-isothermal compositional two-phase flow is considered to be one of the fundamental physical processes in the field of water resources research. The strong non-linearity and discontinuity emerging from phase transition phenomena pose a serious challenge for numerical modeling. Recently, Lauser et al.[1] has proposed a numerical scheme, namely the Nonlinear Complementary Problem (NCP), to handle this strong non-linearity. In this work, the NCP is implemented at both local and global levels of a Finite element algorithm. In the former case, the NCP is integrated into the local thermodynamic equilibrium calculation. While in the latter one, it is formulated as one of the governing equations. The two different formulations have been investigated through several well established benchmarks and analyzed for their efficiency and robustness. In the second part of the thesis, the presented numerical formulations are applied for application and process studies in the context of nuclear waste disposal in Switzerland. Application studies comprehend the coupling between multiphase transport model and complex bio-geo-chemical process to investigate the degradation of concrete material due to two major reactions: carbonation and Aggregate Silica Reaction(ASR). The chemical processes are simplified into a lookup table and cast into the transport model via source and sink term. The efficiency and robustness of the look-up table are further compared with a fully reactive transport model.