3D Flow Model Research Articles

3D numerical and experimental modelling of multiphase flow through an annular geometry applied for cuttings transport

Accurate estimation of volume fraction and pressure gradient is considered indicating parameters of efficient cuttings transportation. It is vital, in this regard, to consider all parameters that can affect cuttings volume fraction and pressure drop during enrollment of the drilling process. The analysis was conducted based on the turbulent flow of a power-law fluid through an annular domain by employing the Finite Volume Method. In addition, dimensional relationships were developed with the Buckingham-π theorem. Before carrying out simulations, the numerical schemes were validated using actual measurements made with the flow loop system available in Texas A&M University at Qatar. The simulation results demonstrated the followings: (i) using a power-law type drilling fluid with a shear-thinning character would reduce energy consumption for an inclination greater than 45°; (ii) inclination angles from 45° to 60° would be least desirable for an effective cuttings transportation with a turbulent Ostwald-de Waele fluid.

Numerical study of the weir angle on the flow pattern and scour around the submerged weirs

In this study, the flow pattern and scouring around the submerged sharp edge weir, with different slope angles, have been numerically simulated using the Flow-3D model. After validating the results of the numerical model in simulating the flow pattern and scour around the weir with a 90∘ angle, the effect of the weir angles relative to the upstream and downstream on the flow and scour field had been investigated. The results showed upstream and downstream of the weir, vortex flows were formed and by taking an angle of 90∘, the vortex flows were weakened. In addition, the downstream slope of the weir had little effect on the upstream water level and most changes occurred in the downstream flow depth with a 90∘ angle. The downstream scour depth was independent of the upstream slope angle and the upstream scour depth was independent of the downstream slope. When the upstream angle became closer to the horizon, scouring was reduced in this area. The depth of scour downstream was greater than the depth of scouring upstream. Angles of 50∘, 70∘, 110∘ and 130∘ had slight changes in the depth of the downstream scour, but compared to the angle of 90∘, the depth of the downstream scour for these angles decreased significantly.

Numerical Study on the Effect of Water Waves and Depths on Inclined Braces with Respect to the Stability of VLFS Platforms in the Caspian Sea

Abstract Very large floating structures (VLFSs) have various applications, such as recreational applications, port facilities, etc. A surge in the population, the advantages of building floating structures compared to traditional methods of land extraction from the sea, and the development of construction technologies, have led to engineers paying attention to very large floating structures. Bracing systems are capable of controlling and reducing the horizontal responses of a floating platform, but they have no major impact on its vertical responses. In the present study, the semi-floating platform was numerically designed to be least affected by the three factors of wave force, horizontal torsion, and horizontal displacement. In order to optimize the design, the semi-floating platform was simulated and subjected to the three wave directions with collision angles of 40, 45 and 55 degrees in the environmental conditions of the Caspian Sea and by exerting the wave effect in a Flow-3D model. Examination of the platform’s movements has demonstrated that the arrangement of an eight-way restraint system with a 40-degree restraint angle responds better to the impact of waves and is more economical compared to other designs.

Subject-Specific Computational Fluid-Structure Interaction Modeling of Rabbit Vocal Fold Vibration.

A full three-dimensional (3D) fluid-structure interaction (FSI) study of subject-specific vocal fold vibration is carried out based on the previously reconstructed vocal fold models of rabbit larynges. Our primary focuses are the vibration characteristics of the vocal fold, the unsteady 3D flow field, and comparison with a recently developed 1D glottal flow model that incorporates machine learning. The 3D FSI model applies strong coupling between the finite-element model for the vocal fold tissue and the incompressible Navier-Stokes equation for the flow. Five different samples of the rabbit larynx, reconstructed from the magnetic resonance imaging (MRI) scans after the in vivo phonation experiments, are used in the FSI simulation. These samples have distinct geometries and a different inlet pressure measured in the experiment. Furthermore, the material properties of the vocal fold tissue were determined previously for each individual sample. The results demonstrate that the vibration and the intraglottal pressure from the 3D flow simulation agree well with those from the 1D flow model based simulation. Further 3D analyses show that the inferior and supraglottal geometries play significant roles in the FSI process. Similarity of the flow pattern with the human vocal fold is discussed. This study supports the effective usage of rabbit larynges to understand human phonation and will help guide our future computational studies that address vocal fold disorders.

Numerical Simulation of Alpine Flash Flood Flow and Sedimentation in Gullies With Large Gradient Variations

Compared with floods occurring over plains, alpine flash floods are formed over scattered locations with complex terrain and data is often lacking regarding the land topography, flow and sedimentation, causing difficulties when developing mathematical models to predict flash floods. The existing flash flood models mainly focus on the influence of water flow, such as the sharp increase in flow discharge caused by convergence at steep slopes, disregarding the sediment load carried by the water flow. However, under the effect of high-intensity sediment transport, the sedimentation in gullies may lead to surges in water level, causing the phenomenon of “great disasters resulting from minor flooding”. In this study, an efficient and accurate water-sediment coupling model was established with a Godunov-type finite volume method based 2D flow model, and a sediment module and OPENMP parallel computing module were added as well. Firstly, a common gully in mountainous area with large gradient variations was used as a generalized model to explore the impact of sedimentation on the flow field of the gully and compared with the physical model. Then, the alpine flash flood incident in Gengdi Village was simulated with the computer model. The calculation results show that the high-intensity sedimentation significantly increased the magnitude of alpine flash floods. Calculated by this mathematical model, the research results verified that this mathematic model can efficiently, accurately and concisely predict the occurrence of flash floods in gullies with large gradient variations. The model also provides a flow and sediment modeling method that incorporates the effect of high-intensity sediment transport into the traditional flash flood flow model. Thus, this model can be a powerful tool for detecting flash floods.

Modelling Groundwater Flow and Contaminant Migration in Heterogeneous Fractured Media at a Municipal Solid Waste Landfill in Nanjing Lishui, China

The migration of groundwater flow and contaminants in fractured medium is complicated owing to the strong heterogeneity and anisotropy of fractured rock mass. Taking the environmental restoration and groundwater protection of the Lishui domestic waste landfill in Nanjing as the background, the groundwater environmental impact assessment and prediction are conducted for the groundwater environmental pollution that may be caused by the leakage of the landfill leachate after the closure of the domestic waste landfill. The strata of the landfill site are clay-cobble gravel, strongly and moderately weathered breccia, with obvious anisotropy and significant differences in rock mass permeability. A 3D numerical model of groundwater flow and contaminant migration in the landfill area is established by integrating the hydrogeological field tests and a conceptual model in the study area. Based on the parametric inversion method, the heterogeneous anisotropic permeability coefficient of the fractured medium is calibrated, and the temporal and spatial migration characteristics of contaminants such as ammonia nitrogen and mercury are predicted using the corrected model under the normal and failure conditions of the antiseepage curtain. The calculated results show that when the antiseepage fails, the maximum migration distances of contaminants in the horizontal direction after 100 days in the old and new landfills are 7.66 m and 15.64 m, respectively, and the maximum migration distances after 20 years are 192.5 m and 113.89 m, respectively. The migration direction and distances of contaminants are consistent with the hydrogeological conditions of the study area. The model calculation results provide a corresponding basis for the antiseepage control of contaminants.

Stormwater management in urban areas using dry gallery infiltration systems

The increase in the frequency of extreme precipitation events due to climate change, together with the continuous development of cities and surface sealing that hinder water infiltration into the subsoil, is accelerating the search for new facilities to manage stormwater. The Canary Islands (Spain) are taking advantage of the knowledge acquired in the construction of water mines to exploit a novel stormwater management facility, which we have defined as a dry gallery. Dry galleries are constituted by a vertical well connected to a horizontal gallery dug into highly permeable volcanic layers of the vadose zone, from where infiltration takes place. However, the lack of scientific knowledge about these facilities prevents them from being properly dimensioned and managed. In this work, we simulate for the first time the infiltration process and the wetting front propagation from dry galleries based on a 3D unsaturated flow model and provide some recommendations for the installation and sizing of these facilities. The fastest advance of the wetting front takes place during the earliest times of infiltration (<2 h), with plausible propagation velocities and infiltration rates higher than 1000 m∙d−1 and 2 m3∙s−1. As time progresses, the propagation velocity and infiltration rate decrease as a consequence of the hydraulic gradient attenuation between the gallery and the aquifer. Therefore, stormwater infiltration is a highly transient process in which a sizing underestimation of 100% may be committed if unsaturated conditions or geological configuration are neglected.

Flood characterization based on forensic analysis of bridge collapse using UAV reconnaissance and CFD simulations

Flash floods are common manifestations of extreme weather events and one of the most severe natural hazards. In Europe, they have been responsible for 359 fatalities and an economic loss totalling 67 million USD in the past decade (EM-DAT), while their increasing severity is linked to climate change. Nevertheless, flash floods remain a poorly documented natural phenomenon due to the lack of flow intensity data in many of the affected watersheds. Based on a thorough field investigation, including UAV-based 3D mapping and material characterization with on-site testing, we carry out a numerical study of a notable flood that caused the collapse of bridges and buildings in Central Greece, following a recent Mediterranean hurricane. Focusing on a carefully selected case study, we combine 3D modelling of flow–structure interaction with detailed mechanical modelling of the nonlinear structural response to reproduce the flood-induced fracture of a bridge abutment. Back-analysis of this failure responds to the fundamental problem of estimating the undocumented magnitude of this extreme event. The paper estimates a lower bound value of the flow velocity at the studied location. This can be valuable input for the interpretation of the extensive damage that took place downstream and for the re-assessment of flood risk in a region where similar events are expected to become more frequent because of climate change. The approach, where disaster forensics and engineering analysis are used to fill the gap of missing real-time measurements, can be implemented for the a posteriori estimation of flood intensity in similar events. The well-documented case study of a bridge failure due to extreme flooding can also be used for validation of future numerical and experimental methods and motivate investigations of the mechanisms governing flow–soil–structure interaction in river crossings.

Flow and heat transfer of nitrogen during liquid nitrogen fracturing in coalbed methane reservoirs

Liquid nitrogen (LN2) fracturing has the potential to induce complex fracture networks, avoid formation damage, and eliminate water consumption. However, the flow and heat transfer of nitrogen in fractures and its effect on the fracture generation during LN2 fracturing have not been studied, and the fracturing mechanism remains unclear. In this paper, the nitrogen flow during LN2 fracturing in a straight fracture was simulated. First, a 3D unsteady-state fluid flow and heat transfer model for LN2 fracturing in coalbed methane (CBM) reservoir was developed, which considered the phase transition of nitrogen, thermophysical properties variation of coal and the heat transfer between nitrogen and formation. This model was then validated against published analytical solutions. Subsequently, the model was applied to elucidate the phase distribution of nitrogen and its influence on fracture generation. Finally, the factors that affect the flow and heat transfer of nitrogen were analyzed. The results showed that the nitrogen at the fracture tip was in a supercritical state. Thermal stress had minor effects on the propagation of the main fracture. In addition, fracture aperture, injection velocity, reservoir temperature, injection fluid temperature, fracture propagation pressure and coal cleat porosity could affect the effectiveness of LN2 fracturing in a coal seam. The main findings of this study are the keys to the research of liquid nitrogen fracturing mechanisms in CBM reservoirs. • Nitrogen at the fracture tip is in supercritical state during the fracture propagation. • Thermal stress had minor effects on the propagation process of the main fracture. • LN2 injection parameters and coal cleat porosity could affect the heat transfer efficiency in LN2 fracturing in coal seam.

Numerical Analysis of Velocity Magnitude on Wave Energy Converter System in Perforated Breakwater

Waves are an alternative energy source that can be used for electricity generation. Wave Energy Converter (WEC) system in perforated breakwater is potentially applicable WEC system for coastal area. The magnitude of wave energy generated is determined by the volume of sea water inside the perforated breakwater. This volumetric flow rate is calculated using the flow velocity at perforated holes on the structure slope. Therefore, this research aims to study the velocity magnitude by analyzing the interrelation among wave steepness, wave run-up and relative velocity. The method used consists of applying numeric 3D flow model in the perforated structure of the breakwater with the variation of wave height, wave period and structure slope. The result shows that, the steeper the structure, the bigger is the relative run up (Ru/H). The higher the relative run up, the higher are the relative run-up velocities (V/Vru). As the velocity increase, the volumetric flow rate inside perforated breakwater will be higher, which leads to higher wave energy. Hence, it can be concluded that the higher the velocities (V/Vru), the higher is the wave energy generated.

Numerical modeling of changes in groundwater storage and nitrate load in the unconfined aquifer near a river receiving reclaimed water.

Reclaimed water (RW) has been widely used as an alternative water resource to recharge rivers in mega-city Beijing. At the same time, the RW also recharges the ambient aquifers through riverbank filtration and modifies the subsurface hydrodynamic system and hydrochemical characteristics. To assess the impact of RW recharge on the unconfined groundwater system, we conducted a 3D groundwater flow and solute transport model based on 10years of sequenced groundwater monitoring data to analyze the changes of the groundwater table, Cl- loads, and NO3-N loads in the shallow aquifer after RW recharge to the river channel. The results show that the groundwater table around the river channel elevated by about 3-4m quickly after RW recharge from Dec. 2007 to Dec. 2009, and then remained stable due to the continuous RW infiltration. However, the unconfined groundwater storage still declined overall from 2007 to 2014 due to groundwater exploitation. The storage began to recover after groundwater extraction reduction, rising from 3.76 × 108 m3 at the end of 2014 to 3.85 × 108 m3 at the end of 2017. Cl- concentrations varied from 5-75mg/L before RW recharge to 50-130mg/L in 2years (2007-2009), and then remained stable. The zones of the unconfined groundwater quality affected by RW infiltration increased from 11.7 km2 in 2008 to 26.7 km2 in 2017. Cl- loads in the zone increased from 1.8×103t in 2008 to 3.8 × 103 t in 2017, while NO3-N loads decreased from 29.8 t in 2008 to 11.9 t in 2017 annually. We determined the maximum area of the unconfined groundwater quality affected by RW, and groundwater outside this area not affected by RW recharge keeps its original state. The RW recharge to the river channel in the study area is beneficial to increase the groundwater table and unconfined groundwater storage with lesser environmental impacts.

Research on Seepage Control of Jurong Pumped Storage Hydroelectric Power Station

Based on the geological and hydrogeological conditions of the Jurong Pumped Storage Hydroelectric Power Station (JPSHP), a 3D groundwater flow model was developed in the power station area, which took into account the heterogeneity and anisotropy of fractured rocks. A control inversion method for fractured rock structural planes was proposed, where larger-scale fractures were used as water-conducting media and the relatively intact rock matrix was used as water-storage media. A statistical method was used to obtain the geometric parameter values of the structural planes, so as to obtain the hydraulic conductivity tensor of the fractured rocks. Combining the impermeable drainage systems of the upper storage reservoir, underground powerhouse and lower storage reservoir, the 3D groundwater seepage field in the study area was predicted using the calibrated model. The leakage amounts of the upper storage reservoir, powerhouse and lower storage reservoir were 710.48 m3/d, 969.95 m3/d and 1657.55 m3/d, respectively. The leakage changes of the upper storage reservoir, powerhouse and lower storage reservoir were discussed under the partial and full failure of the anti-seepage system. The research results provide a scientific basis for the seepage control of the power station, and it is recommended to strengthen the seepage control of the upper and lower storage reservoirs and the underground powerhouse to avoid excessive leakage and affect the efficiency of the reservoir operation.

LOCAL WELL-POSEDNESS FOR A 3D LIQUID-GAS TWO PHASE MODEL WITH VACUUM

In this paper we prove the local well-posedness of strong solutions to a 3D liquid-gas two-phase flow model with vacuum in a bounded domain without the standard compatibility conditions.

Wavelet Operational Matrices and Lagrange Interpolation Differential Quadrature‐Based Numerical Algorithms for Simulation of Nanofluid in Porous Channel

This work analyses the features of nanofluid flow and thermal transmission (NFTT) in a rectangular channel which is asymmetric by developing two numerical algorithms based on scale‐2 Haar wavelets (S2HWs), Lagrange’s interpolation differential quadrature technique (LIDQT), and quasilinearization process (QP). In the simulation procedure, first of all, using similarity transformation (ST), the governing unsteady 2D flow model is changed into two highly non‐linear ODEs. After that, QP is applied to linearize the non‐linear ODEs, and finally S2HWs and LIDQT are used to simulate the non‐linear system of ODEs. In results and discussion section, the parameters Reynolds number (R), expansion ratio and Nusselt (Nu), and nanoparticle volume fraction (φ) are analysed with respect to velocity and temperature profiles. The proposed techniques are easy to implement for fluid and heat transfer (FHT) problems.

Examples on Simulation Model

Abstract: A new methodology was developed Further real-time determination gate control operations of a river-reservoir system to minimize flooding conditions. The methodology is based upon an optimization-simulation model approach interfacing the genetic algorithm within simulation software for short-term rainfall forecasting, rainfall–runoff modeling (HEC-HMS), and a one-dimensional (1D), two-dimensional (2D), and combined 1D and 2D combined unsteady flow models (HEC-RAS). Both realtime rainfall data from next-generation radar (NEXRAD) and gaging stations, and forecasted rainfall are needed to make gate control decisions (reservoir releases) in real-time so that at timet, rainfall is known and rainfall over the future timeperiod(∆t)totimet+ ∆t can be forecasted. This new model can be used to manage reservoir release schedules (optimal gate operations) before, during, and after a rainfall event. Through real-time observations and optimal gate controls, downstream water surface elevations are controlled to avoid exceedance of threshold flood levels at target locations throughout a riverreservoir system to minimize the damage. In an example application, an actual river reach with a hypothetical upstream flood control reservoir is modeled in real-time to test the optimization-simulation portion of the overall model. Keywords: Simulation – Random numbers- Steps for simulation – Problems.

Slip length measurement in rectangular graphene nanochannels with a 3D flow analysis

Although many molecular dynamics simulations have been conducted on slip flow on graphene, experimental efforts remain very limited and our understanding of the flow friction on graphene remains far from sufficient. Here, to accurately measure the slip length in rectangular nanochannels, we develop a 3D capillary flow model that fully considers the nonuniform cross-section velocity profile, slip boundary conditions, and the dynamic contact angle. We show that the 3D analysis is necessary even for a channel with a width/height ratio of 100. We fabricated graphene nanochannels with 45-nm depth and 5-μm width, and measured slip lengths of about 30–40 nm using this 3D flow model. We also reevaluated the slip-length data for graphene obtained from capillary filling experiments in the literature: 30 nm instead of originally claimed 45 nm for a 25-nm-deep channel, and 47 nm instead of 60 nm for an 8.5-nm-deep channel. We discover a smaller slip length than existing experimental measurements due to our full 3D flow analysis considered in our method. This work presents a rigorous analysis approach while also providing a better understanding of slip flow in graphene nanochannels, which will benefit further innovation in nanofluidic applications, including electronics cooling and biomedical chips.

Using depth specific electrical conductivity estimates to improve hydrological simulations in a heterogeneous tile-drained field

In agricultural fields, tile drains represent potential pathways for the migration of solutes, such as nitrates, in groundwater and surface water bodies. Tile drain flow is controlled by the temporal and spatial dynamics of the shallow groundwater table, which results from complex interactions between climate, topography and soil heterogeneity. Studies on the effect of topsoil heterogeneity on shallow water and drainage dynamics by fully 3D surface water and groundwater flow modeling are limited. The objective of our study is to demonstrate the use of depth specific electrical conductivity (EC) estimates to improve hydrological simulations in a tile-drained field. The model was applied to a field site in Denmark where times series of drainage discharge and water table elevations are available. Clay-rich soil zones were identified in a tile-drained field using depth specific electrical conductivity estimates generated by the inversion of apparent electrical conductivity data measured using an electromagnetic induction instrument. One model that included the low-permeability clayey zones in the soil layers down to a depth of 1.2 m was compared to a simpler model that assumed homogeneous soil layers. Both models simulate drainage discharge that compares well to the observations. However, including the clayey zones improves the simulation of hydraulic heads, and water table fluctuations, and generates flooded areas that are more representative of those observed during the wet seasons. Our results suggest that the simulation of water table fluctuations can be improved when the soil heterogeneity determined from depth specific EC estimates is included in integrated hydrological models. A better representation of the subsurface flow dynamics will also improve subsequent simulations of the transport and fate of agrochemical substances leaching from fields such as nitrate, which may deteriorate the quality of groundwater and surface water bodies.

Numerical Modeling of Failure Mechanisms in Articulated Concrete Block Mattress as a Sustainable Coastal Protection Structure

Shoreline protection remains a global priority. Typically, coastal areas are protected by armoring them with hard, non-native, and non-sustainable materials such as limestone. To increase the execution speed and environmental friendliness and reduce the weight of individual concrete blocks and reinforcements, concrete blocks can be designed and implemented as Articulated Concrete Block Mattress (ACB Mat). These structures act as an integral part and can be used as a revetment on the breakwater body or shoreline protection. Physical models are one of the key tools for estimating and investigating the phenomena in coastal structures. However, it does have limitations and obstacles; consequently, in this study, numerical modeling of waves on these structures has been utilized to simulate wave propagation on the breakwater, via Flow-3D software with VOF. Among the factors affecting the instability of ACB Mat are breaking waves as well as the shaking of the revetment and the displacement of the armor due to the uplift force resulting from the failure. The most important purpose of the present study is to investigate the ability of numerical Flow-3D model to simulate hydrodynamic parameters in coastal revetment. The run-up values of the waves on the concrete block armoring will multiply with increasing break parameter (0.5&lt;ξm−1,0&lt;3.3) due to the existence of plunging waves until it (Ru2%Hm0=1.6) reaches maximum. Hence, by increasing the breaker parameter and changing breaking waves (ξm−1,0&gt;3.3) type to collapsing waves/surging waves, the trend of relative wave run-up changes on concrete block revetment increases gradually. By increasing the breaker index (surf similarity parameter) in the case of plunging waves (0.5&lt;ξm−1,0&lt;3.3), the low values on the relative wave run-down are greatly reduced. Additionally, in the transition region, the change of breaking waves from plunging waves to collapsing/surging (3.3&lt;ξm−1,0&lt;5.0), the relative run-down process occurs with less intensity.

Heat extraction study of a novel hydrothermal open-loop geothermal system in a multi-lateral horizontal well

The efficient development of hydrothermal resource is an important focus in the global geothermal industry. A novel hydrothermal open-loop geothermal system in a multi-lateral horizontal well is presented to achieve the simultaneous exploitation of multi-layer reservoirs and accelerate the large-scale development of the hydrothermal resource. To investigate the detailed performance of the novel geothermal system, a novel 3D fluid flow and heat transfer numerical model is established based on the local thermal non-equilibrium and dual porosity theories. The flow and temperature fields are analyzed, while the pressure and temperature interferences between the lateral wellbores of the multi-lateral horizontal well system are comprehensively investigated. The influences of the injection flow rate, lateral-wellbore number, spacing between the injection segment and the main wellbore, fracture aperture on the heat extraction are studied. Furthermore, the thermal characteristics of the novel geothermal system and traditional doublet well system are compared for the first time. The investigation indicates that the flow impedance of the novel geothermal system can decrease by 75.38%, and its energy efficiency reaches 8 times that of the doublet well system. The findings can provide a key reference for scientific research and application of the novel and traditional geothermal systems.

Mapping and Monitoring of DNAPL Source Zones With Combined Direct Current Resistivity and Induced Polarization: A Field‐Scale Numerical Investigation

AbstractDirect current (DC) resistivity has been widely investigated for non‐invasive mapping of dense non‐aqueous phase liquids (DNAPLs); however, due to its difficulty in distinguishing DNAPLs from adjacent soils, the DC method is limited for static detection of DNAPLs and more often employed for monitoring DNAPL changes over time. Time‐domain induced polarization (IP) can provide complementary information to better discriminate between DNAPL, water and surrounding soils. Since highly resistive DNAPLs tend to laterally spread and pool on polarizable (chargeable) clay lenses, combined DC and IP (DCIP) has the potential to enhance static detection, and also monitoring of DNAPL source zones (SZs). The objective of this study is to assess DCIP for characterizing and monitoring DNAPL SZs at the field‐scale. A new DNAPL‐DCIP numerical model was first developed that couples a 3D multiphase flow model, which simulates DNAPL release and remediation scenarios, with a 3D DCIP model, which calculates the corresponding resistivity and chargeability response. The sensitivity of the DCIP response to key DNAPL and soil properties was then analyzed at a single subsurface location, closely matching previous theoretical and experimental observations. Finally, a field‐scale simulation of DNAPL release and remediation was conducted, with simultaneous mapping by DCIP surveys. Results demonstrate that chargeability can provide enhanced understanding of the lithological distribution that controls the variability in the DNAPL SZ, with time‐lapse resistivity being used to monitor DNAPL mass changes during SZ remediation. This numerical study suggests that combined DCIP can be valuable for site characterization and time‐lapse monitoring performance at DNAPL‐impacted sites.