Implicit Finite Difference Method Research Articles

Finite difference simulation of natural convection of two-phase hybrid nanofluid along a vertical heated wavy surface

This study employs finite difference modelling to investigate the natural convection of a two-phase hybrid nanofluid consisting of Al 2 O 3 and Cu nanoparticles dispersed in water along a vertically heated wavy surface. The hybrid nanofluid has unique features that influence convective heat transfer. Numerical solutions analyze fluid flow patterns and temperature distributions, providing insights into heat exchange and thermal management. The research highlights the importance of hybrid nanofluids, notably those containing Al 2 O 3 and Cu in a water base, in improving overall heat transfer efficiency. The mathematical model accounts for laminar and incompressible fluid flow with a Prandtl number of Pr = 6.2, Lewis number Le = 10, and maximum 10 % concentration of hybrid nanoparticles. After transforming the governing equations into a non-dimensional form, the implicit finite difference method is used to solve them. The solution employs the in-house FORTRAN 90 code and compares its outcomes with the benchmark results. Various parameters, such as the Schmidt number (Sc = 1 to 10), volume fraction of nanoparticles ( ϕ = 0.0 to 0.1), wavy amplitude (A = 0.0 to 0.3), and N BT = (0.05 to 0.2), are investigated regarding temperature, velocity, local skin friction coefficient ( C f ) , local Nusselt number (Nu), streamlines, and isotherms. Elevated volume fraction generally decreases skin friction, increases temperature, and may reduce velocity. Meanwhile, higher skin friction, temperature, velocity, and local Nusselt numbers often correlate with elevated Schmidt numbers. Changes in amplitude and N BT affect skin friction, temperature, velocity, and local Nusselt numbers, providing information on the dynamics of fluid systems. For example, when the volume fraction (ϕ) increases from 0 to 0.1 at Y = 2, the temperature rises by 75 % , while the velocity gradually decreases by 11.11 % . This is in contrast to findings according to when the N BT rises from 0.05 to 0.2 at Y = 1. In this case, the temperature drops by 18.75 % , followed by a 9.68 % fall in fluid velocity. Understanding these interactions enhances comprehension of heat transfer and fluid dynamics. By investigating these interactions, we not only improve our understanding of heat transfer and fluid dynamics, but also open the door to new strategies for improving thermal systems and manufacturing procedures.

Finite element modeling of thermal‐hydro‐mechanical coupled processes in unsaturated freezing soils considering air‐water capillary pressure and cryosuction

AbstractThis paper presents a comprehensive computational model for analyzing thermo‐hydro‐mechanical coupled processes in unsaturated porous media under frost actions. The model employs the finite element method to simulate multiphase fluid flows, heat transfer, phase change, and deformation behaviors. A new soil freezing characteristic curve model is proposed to consider the suctions from air‐water capillary pressure and water‐ice cryosuction. A total pore pressure with components from liquid water pressure, air pressure, and ice pressure is used in the effective stress law. Vapor and dry air are considered miscible gases, utilizing the ideal gas law and Dalton's law. The governing equations encompass the linear momentum balance equation, the energy balance equation, and mass conservation equations for water species (ice, liquid, and vapor) and dry air. Weak forms are formulated based on primary variables of displacement, water pressure, air pressure, and temperature. The spatial discretization is achieved through the finite element method, while temporal discretization employs the fully implicit finite difference method, resulting in a system of fully coupled nonlinear equations. To verify the proposed computational model, a numerical implementation is developed and validated against a set of experimental data from the literature. The successful verification demonstrates the robustness of the model. A detailed discussion of the contributions from phase change strain and different sources of pore pressure is also addressed.

Finite difference analyses of bioconvective Williamson nanofluids stream over irregular plates with effects of nonlinear radiation and viscous dissipation

The purpose of this work is to provide a comprehensive analysis of the bioconvective Williamson nanofluid stream while considering the effects of radiation flux and viscous dissipation. The two-phase nanofluids model is applied where the impacts of the Brownian motion and thermophores are significant. This study starts with introducing the governing system in case of the irregular plane then it transforms to the regular plane. The solution methodology is based on the fully implicit finite difference method. The major results disclosed that maximizing the amplitude parameter causes a diminishing in the velocity and skin-friction coefficients. Also, lower behaviors are noted for the velocity in the case of Williamson nanofluids compared to a Newtonian nanofluid case (We = 0). Furthermore, the considered range of Biot number (1 ≤ Bi ≤ 5)causes an improvement in the surface temperature up to 71.59%.

Nonlinear Convection Flow of a Micropolar Nanofluid Past a Stretching Sphere with Convective Heat Transfer

An incompressible, steady combined nonlinear convective transport system on a micropolar nanofluid through a stretching sphere with convective heat transfer was investigated. The conservation equations corresponding to momentum, microrotation, thermal energy, and concentration particles have been formulated with suitable boundary constraints. By using the required non-dimensional variables, the conservation equations have been turned into a set of high-order standard differential equations. Then, an implicit finite difference method, also known as the Keller-Box Method (KBM), was used to numerically solve the flow problem. The obtained outcomes are displayed through graphs and tables to explain the impact of various governing variables over velocity, temperature, concentration, number of skin friction, wall coupled stress, Nusselt number, and Sherwood number. The findings demonstrate that increasing the convective heat parameter Bi enhances the factor of skin friction, local Nusselt number, Sherwood number, velocity field, and temperature profile while lowering the wall-coupled stress. It is observed that for high values of the material parameter β, the fluid velocity and the spin of the micro-elements both increase, which causes the dynamic viscosity and microrotation velocity to decrease. In addition, as the rates of magnetic constant Ma, thermophoresis Nt and Brownian movement Nb rise, the thermal distribution and its thickness of boundary layer increase. However, it decline along the enlarging quantities of nonlinear convection parameter λ, Prandtl number Pr, material parameter β, and solutal Grashof number Gm, which agrees to increase fluid density. When the range of thermophoresis Nt surges, it causes an increment in the nanoparticle species, but the opposite behavior have seen in the case of Brownian number Nb, and Lewis number Le. The comparison made with the related published paper achieves a significant agreement. The numerical result is generated through the implementation of the computational software MATLAB R2023a.

The Generalized Fractional-Order Fisher Equation: Stability and Numerical Simulation

This study examines the stability and numerical simulation of the generalized fractional-order Fisher equation. The equation serves as a mathematical model describing population dynamics under the influence of factors such as natural selection and migration. We propose an implicit exponential finite difference method to solve this equation, considering the conformable fractional derivative. Furthermore, we analyze the stability of the method through theoretical considerations. The method involves transforming the problem into systems of nonlinear equations at each time since our method is an implicit method, which is then solved by converting them into linear equations systems using the Newton method. To test the accuracy of the method, we compare the results obtained with exact solutions and with those available in the literature. Additionally, we examine the symmetry of the graphs obtained from the solution to examine the results. The findings of our numerical simulations demonstrate the effectiveness and reliability of the proposed approach in solving the generalized fractional-order Fisher equation.

Transient magnetohydrodynamic heat and mass transfer analysis of chemically reacting fluid flow over oscillatory permeable media

This study focuses on the numerical observation of the convective motion of chemically active magnetohydrodynamic (MHD) fluid through a vertically oriented permeable medium, incorporating variations in mass and heat transfer. The fluid type is assumed to be incompressible, chemically strongly ionized and viscous with some mass infusibility. The model associated with this problem is solved by a highly stable Implicit Finite Difference Method (IFDM). The method is used for small and large deflection of the physical parameters, which results in a noticeable fluid flow behavior. Numerical configuration is graphically depicted to scrutinize the fluid behavior. The momentum, energy, concentration diffusion, skin friction, Nusselt number and Sherwood number are investigated for numerous factors such as magnetic field, permeability and chemical reaction rate. The current study unveils significant findings, demonstrating that a heightened rate of chemical reaction in the presence of magnetic effects, coupled with specific porosity, diminishes ionization energy, resulting in a concurrent decrease in the concentration and momentum profiles of the fluid flow. The rise in the viscous diffusion rate is attributed to escalating values of the Schmidt number, causing an augmentation in dynamic viscosity and consequently resulting in an overall reduction in the momentum of the fluid flow.

A novel finite difference scheme for numerical solution of fractional order population growth model

In this paper, we propose a new scheme based on the implicit finite difference method for solving the fractional population growth model (FPGM). We use the well-known L1 finite difference method to approximate the Caputo fractional derivative of order 0 < α ≤ 1, and the linear interpolation to approximate the integral part. We provide a study on the stability and convergence of the scheme. We present the numerical solution of the proposed method and compare it with three other methods to demonstrate its validity and efficiency.

Analyzing thermophoresis phenomenon in a transient free convective MHD chemically radioactive fluid over oscillatory permeable media

A numerical observation of chemically radioactive free convective magneto hydrodynamics fluid passing through a vertical permeable medium with fluctuating mass and heat transfer with thermophoresis is studied here. The fluid type is assumed to be a suspended mixture of particles that are incompressible, chemically strong, ionized, and viscous with some mass infusibility. The model pertaining to this problem is resolved through the application of a highly stable implicit finite difference method. The method is used for small and large deflections of the physical parameters, which result in a noticeable fluid flow behavior depicted graphically. The momentum, energy, concentration, skin friction, Nusselt and Sherwood number are investigated for various factors such as the thermophoresis effect, magnetic field, permeability, and chemical reaction rate. The substantial findings of the present investigation indicate that the impact of thermophoresis is more pronounced in the diffusion of concentration as compared to the rate of chemical reactions. The concentration is considerably influenced by the permeability parameter, surpassing the influence of the chemical reaction rate. For higher radiation parameters and Prandtl number, the molecular interaction of liquids increases, which in turn constrains the extent of radiative heat exchange.

Two layer film flow on an unsteady stretching cylinder

Most of the earlier research works related to the thin film flow over a stretching cylinder are primarily limited to the boundary layer flow of a single liquid. Whereas, multilayer film coating is commonly used to coat microelectronic chips, circuits, etc. This scenario motivated us to investigate the thin double-layer liquid film development over a stretching cylinder. In the present work, the above limitations of earlier models have been removed to investigate the two-layer film flow over a stretching cylinder by considering a full set of momentum equations. The governing set of equations is transformed into a system of nonlinear partial differential equations (PDEs) by applying appropriate similarity transformation. The analytical solutions are obtained by using the asymptotic analysis method for smaller values of Reynolds number (i.e. Re1<<1). Whereas numerical solutions are obtained by implicit finite difference method for smaller as well as moderate values of Re1. It is found that the film thinning rate enhanced for both the liquid layers with raising values of the cylinder’s radius (A) and density ratio (m). It is further seen that both the liquid layers thin faster with a decreasing initial thickness ratio of both liquid layers (δ). It is also observed that the variation of film height for the first layer is very small in comparison to the second layer for increasing values of viscosity ratio parameter (n). The film height for both the liquid layers decreases very fast for a smaller Reynolds number (Re1) at the initial time but a reverse trend is noticed at large time.

Numerical modeling of Carreau–Yasuda fluid with induced magnetic field in a heated curved channel having ciliary walls

Cilia motion is commonly located in mammalian cells of living organisms which continually grow and feature conspicuously. The ciliary apparatus is associated with cell cycle movement and proliferation, and cilia play a dynamic part in human and animal growth and everyday life. The motivation of these applications, this article is to present a metachronal wave analysis for non-Newtonian Carreau–Yasuda fluid inside a curved channel with ciliated walls. Hall effect is encountered in the momentum equations to measure the magnetic field while the viscous dissipation and thermal radiation terms are considered to address the energy analysis. The curvilinear coordinates are used due to the curved nature of flow geometry in the derivation of flow equations. The computations for axial velocity, pressure rise, temperature field, mass concentration, and stream function are carried out under low Reynolds number and long wavelength approximation in the wave frame of reference by utilizing appropriate numerical implicit finite difference method (FDM). The implementation of numerical procedure and graphical representation of the computations are accomplished using MATLAB language. Furthermore, the skin friction and Nusselt number at the channel walls are determined for a variety of critical parameter assessments.

Wave oscillations in thermal boundary layer of Darcy-Forchheimer nanofluid flow along buoyancy-driven porous plate under solar radiation region

The significant contribution of current mechanism is to explore the periodical and fluctuating results of heat and mass flux of Darcy-Forchheimer nanofluid flow along buoyancy-driven porous plate under solar radiation region. Flow through a Darcy-medium has a wide range of applications such as the use of oil in various hydrothermal transfer control, radioactive nuclear disposal systems, water improvement, and filtration of water. The dimensional model is transformed into non-dimension for scaling factors. The primitive-based transformation is applied on steady and oscillatory parts for smooth algorithm in FORTRAN language machine. The convenient model is again reduced into algebraic pattern by using implicit finite difference method. The numerical and graphical results of velocity, temperature and concentration are executed by Gaussian elimination method. To enhance the frequency and wavelength, the impact of solar radiations, thermophoresis, Brownian movement, Schmidt factor, Prandtl factor, porosity and buoyancy pressure is applied on periodic nanoparticles with Darcy Forchheimer relation. The novelty of this proposal to explore the wave oscillations, amplitude and phase angle of thermal and concentration boundary layer of Darcy-Forchheimer nanofluid flow along buoyancy-driven porous plate under solar radiation region. It is noticed that the prominent wavelength and frequency in thermal and concentration boundary layers is generated under porous and solar radiation region. The significance of temperature variation increases as solar radiation, Brownian motion and thermophoresis increases. It is noticed that heat and mass transport enhances as Darcy-Forchheimer and modified buoyancy force increases. In thermodynamic systems, the Darcy-porous materials is prominent assumption that is employed in many domains, especially clinical science, petrochemical and structural engineering, geography, biochemistry and biological physics, and material science. For energy storage and energy conservation, the Darcy-porous material is required for batteries, fuel cells and super capacitors.

Thermal radiation effects on oscillating frequency of heat transfer of Darcy–Forchheimer nanofluid with chemical reaction and applications in machining operations

The advancement of cutting tool components and design is presently promoting innovative developments in numerous different machining-related industries. The characteristics of nanofluid are important for machining activities such as the drilling process, grinding, rotating, milling, and cutting. Various machining procedures require distinct lubricating oils and nanofluids for cutting-edge innovations. The significant contribution of the current mechanism is to explore the fluctuating heat and mass flux of Darcy–Forchheimer chemically reactive nanofluid along a buoyancy-driven porous plate under solar radiation region. Flow through a Darcy medium has a wide range of applications such as the use of oil in various hydrothermal transfer control, radioactive nuclear disposal systems, water improvement, and filtration of water. The dimensional model is transformed into non-dimension for scaling factors. The primitive-based transformation is applied on steady and oscillatory parts for smooth algorithm in FORTRAN language machine by using an implicit finite difference method. The numerical and graphical results of velocity, temperature, and concentration are executed by the Gaussian elimination method. To enhance the frequency and wavelength, the impact of solar radiations is applied on periodic nanoparticles with Darcy–Forchheimer relation. The novelty of this proposal is to explore the wave oscillations, amplitude, and phase angle of thermal and concentration boundary layer of Darcy–Forchheimer nanofluid flow under chemical reaction and solar radiation region. It is noticed that the prominent wavelength and frequency in thermal and concentration boundary layers is generated under porous and solar radiation region. The significance of temperature variation increases as solar radiation, chemical reaction, Brownian motion, and thermophoresis increase. It is found that minimum oscillation in heat transport is observed as Pr decreases but maximum oscillation in heat transfer is sketched as Pr enhances.

Application of finite difference method to 2D heat conduction in young hydrating mass concretes structures.

Finite Difference Method is a numerical approximation scheme usually applied in solving differential equations derived from the knowledge of the behaviour of engineering systems. Due to the unsteady state nature of cement hydration and heat conduction in concretes, Crank Nicholson implicit approach was adopted. Concrete block of size 1.10 m x 1.10 m x 1.10 m, cast with mix ratio 1:3:6, probed thermocouples and digital thermometer were used to verify the numerical computational analysis. The temperature profile depicted a relatively hot core, with peak temperature values occurring at 24 hours of placement. Implicit finite difference method was utilized in determining time dependent temperature profile within mass concretes at early ages. The paper could provide a guide towards proactive decision making in terms of controlling thermal cracks.

Implicit finite difference simulation of non-similar conduction–convection interaction of magnetohydrodynamics water-Cu nanofluid flow along a vertical surface

Heat-exchanging devices equipped with conventional thermal fluids are failing to meet the standard in terms of thermal efficiency. The inclusion of nanoparticles within a base fluid to improve thermal efficiency can be a great alternative but a proper understanding on the interactions during the heat transfer phenomena is required. This work examines the interactions between conduction and convection in a nanofluid of water and copper (Cu) particles. The investigation focuses on the behavior of the nanofluid, considering the impact of magnetohydrodynamics (MHDs) flow over a vertical surface. The governing equations were non-dimensionalized by the suitable transformation and then were solved numerically by the implicit finite difference method (FDM). The key findings suggested that increasing magnetic parameter plummeted skin friction coefficient and accelerated the heat dissipation up to by augmenting the peaks of fluid velocity as well as the temperature under a constant volume fraction of nanoparticles. In addition, the thermal boundary layer thickness also increased concurrently. Alternatively, when a constant non-zero magnetic parameter was considered, and volume fraction of nanoparticles was allowed to vary, the convective heat transfer got enhanced. The findings from this study will assist the relevant industry to identify the suitable fraction of Cu nanoparticles in the presence of magnetic field in order to achieve the best thermal efficacy of a heat-exchanging unit.

Implicit Finite Difference Simulation of Hybrid Nanofluid along a Vertical Thin Cylinder with Sinusoidal Wall Heat Flux under the Effects of Magnetic Field

A numerical analysis of magnetohydrodynamic natural convection along a thin vertical cylinder with a sinusoidal heat flux at the wall immersed in copper (Cu) and aluminum-oxide (Al2O3) hybrid nanofluids has been studied. A 2D vertical thin cylinder shape geometry has been considered with a radius of R. The fluid flow is considered laminar and incompressible with the Prandtl number of Pr = 6.2 and 10% concentration of hybrid nanoparticles. The nondimensional governing equations have been solved numerically by using the implicit finite difference method. An in-house FORTRAN 90 code is used for solving this problem and the code is validated with the available benchmark results. Numerical simulations have been performed for a wide range of governing parameters, Hartmann number from Ha = 0 to Ha = 4, nanoparticles volume fractions ϕ = 0.0 to ϕ = 0.1, and the amplitude of the wall heat flux ε = 0.0–0.3. The findings have been illustrated in terms of streamlines, isotherms, local skin friction coefficients, local Nusselt numbers, velocity, and temperature distributions. The flow field and temperature distribution within the boundary layer are deceased by the effects of the wall heat flux amplitudes. It is also noted that the rate of heat transfer increases with particle volume fraction and the amplitude of the wall heat flux. According to the findings, Nu increases by 24.72% as ϕ increases from 0 to 0.1 while ε = 0.3, and 27.66% while ε increases from 0.0 to 0.3 at 5% hybrid nanoparticles. The local skin frictions and Nusselt number diminish with the increment of the Hartman number due to the effects of the Lorenz force. The findings of this study can lead to a better understanding of the fundamental principles regarding the behavior of hybrid nanofluids under complex conditions, such as a vertical thin cylinder with a sinusoidal wall heat flux. Understanding the behavior of hybrid nanofluids in the presence of a magnetic field and a nonuniform wall heat flow can also lead to the development of innovative heat transfer enhancement strategies.

Dynamic analysis of wave scenarios based on enhanced numerical models for the good Boussinesq equation

The good Boussinesq equation, a modification of the Boussinesq equation, aims to enhance predictions about shallow water wave behavior. This paper introduces two finite difference schemes for solving the good Boussinesq equation including linear and nonlinear implicit finite difference methods. Both schemes utilize the pseudo-compact difference approach, delivering second-order precision with an additional term to boost numerical simulation accuracy while maintaining the grid points of the standard scheme. These schemes rigorously preserve the critical physical characteristics of the good Boussinesq equation, ensuring more precise representation. We establish the existence of solutions with discrete differences and demonstrate, through the discrete energy method, their uniqueness, stability, and second-order convergence in the maximum norm. Furthermore, we propose an iterative algorithm tailored for the nonlinear implicit finite difference scheme, resulting in significant reductions in computational costs compared to the linear scheme. The results of our numerical experiments demonstrate that our methods are competitive and efficient when compared to difference schemes and previously used methods, while maintaining crucial physical qualities. Furthermore, we run relevant numerical simulations to demonstrate the accuracy of the current methods using evidence from the solitary wave interaction with the initial amplitudes of the wave. It is also suggested that the issue has a critical initial wave amplitude for the interaction of two solitary waves, where blow up occurs in a finite amount of time for initial wave amplitudes greater than the new blow-up criteria value.

Analytical and Numerical Models of a New Hybrid System of Earth Air Tunnel Integrated Solar Air Heater With Sensible Storage

Abstract A novel combination of an earth air tunnel (EAT) and a sensible thermal storage-aided solar air heater has been proposed when the sensible storage medium is taken to be rocks filled within the air flow passage. Both time-dependent and steady-state models are reported. The former describes the exact thermal performance of the system for a two-month wintertime for the city of Baghdad, Iraq. The latter describes average performance of the system. The implicit finite difference method and method of separation of variables have been used to solve the pertinent equations in the respective models. Acceptable level of accuracy was obtained between the steady-state numerical and analytical solutions as well as between the numerical and the published data. It is revealed that preheating by a sufficiently long EAT module improves the effective power as well as the output temperature from the rock bed solar air heater by about 35% and 9% respectively for the present set of parameters considered. It is also observed that for the present set of parameters considered, the temperature gradients in the direction normal to air flow in the solar air heater are insignificant and may be ignored. A parametric study is also carried out that assesses the impact of system parameters on the quality as well as quantity of the energy extracted from the system. This work is hence the first to couple an EAT with a sensible thermal storage equipped solar air heater and may pave the way for future studies.

MHD flow of third-grade fluid through a vertical micro-channel filled with porous media using semi implicit finite difference method

This study examines the complex dynamics of an incompressible, electrically conducting third-grade fluid (TGF) flow through a vertical micro-channel that is filled with porous media. The research utilizes asymmetrical convective cooling, a transverse magnetic field, and exothermic reactions to introduce a comprehensive model. The viscosity of the fluid, which is determined by the Naham-type law, changes significantly with temperature. Additionally, the convective heat transfer rates, which govern the asymmetric convective cooling at the surfaces of the micro-channel, are modeled using Newton's law of cooling, and Modified Darcy law is used to address the porous medium effect. To solve the complex system of nonlinear partial differential equations, we employ computational solutions derived from reliable and efficient finite difference methods (FDM). The method is tested for different time-steps and mesh sizes and it is found that the algorithms reliably produce the same results. Our investigation includes an evaluation of both spatial and temporal convergence of the FDM scheme. The study present qualitative discussions of graphical results, that highlight the impact of various flow parameters integrated into the system. The velocity and temperature profiles increases for higher values of thermal Grashoff number Gr and Reynolds number Re while opposite effect is noticed for increasing values of Hartman number M.

Experimental and Hybrid Numerical Analytical Approach of an Earth-to-Air Heat Exchanger

Abstract The primary goal of this paper is to assess the thermal behaviour of Earth-to-Air Heat Exchanger (EAHE) under continuous operation mode during the cooling period in warm climatic conditions. An experimental setup was conducted in Biskra University (Algeria) to test the ability to use EAHE as an alternative solution to conventional air conditioning systems. Besides, a transient numerical approach for the flowing air within the EAHE was developed using energy balance equations, discretised by the implicit finite differences method and solved by Thomas method. Furthermore, the transient temperature profile of the soil around the EAHE was investigated in order to determine the adiabatic radius of the soil as a function of the duration of operation. In this context, a transient semi-analytical model was integrated to the air model to estimate the transient EAHE’s performance deterioration. This deterioration is produced by the thermal saturation of the soil, which is influenced by the heat released from the air inside the pipe. Results showed that continuous operation for a duration of four days has minimal impact on the outlet air temperatures due to the implementation of night purging. Furthermore, the design and thermal efficiency of the EAHE are impacted by factors such as higher air velocity and reduced soil thermal conductivity.

Gas composition tracking feasibility using transient finite difference [formula omitted]-scheme model for binary gas mixtures

Implicit finite difference method is popular numerical method for simulation of transient pipeline dynamics, governed by the system of Euler partial differential equations. If advection equation is included in the mathematical formulation, gas composition tracking is also possible. If implicit scheme is used for gas composition tracking, simulated concentration profiles are distorted due to the occurrence of numerical diffusion, a non-physical phenomenon, caused by the 1st order accurate numerical implicit scheme. Intensity of the numerical diffusion can be decreased by use of θ-scheme with arbitrarily chosen implicitness. In order to evaluate accuracy of such numerical scheme and adequacy for usage in actual application, numerical model must be evaluated on actual gas pipeline system topology with actual measured boundary conditions and compared to actual measured results. Aim of this paper is to assess whether hydrogen concentration tracking is feasible using such approach considering complex grid topologies. Paper highlights the peculiarities of the proposed θ-scheme and finite difference method and its impact on the calculated concentration profile for binary gas mixtures. Non-physical flow reversal issue and its impacts on concentration tracking in looped pipeline system is researched and explained. Simulation results shows high accuracy for simulated values of hydraulic variables compared to the measured values, in some cases fit is perfect. Discrepancy between measured and the simulated values of the pressures at nodes on the main transport route is negligible, while the discrepancy increases with the distance from the main transport routes. The largest discrepancy occurred on the most distant node and was between 1.7006 bar (4.85%) and 1.7119 bar (4.88%), depending on the θ value. On the other hand, simulated concentration profile lacks finer details due to the numerical diffusion effects while the simulated concentration profile lags the measured profile, partly by the impacts of the 1st order accurate scheme and partly by effects of non-physical reversed flows. Discrepancy between simulated and measured amplitude of the tracked concentration profile was negligible, while discrepancy due to the time delay of the simulated and measured profile was more expressed and ranged between 2h to 9h. Largest time delay was caused by the non-physical flow reversal.