Numerical Modeling of Multi-Frequency Complex Dielectric Permittivity Dispersion of Sedimentary Rocks

Geophysics ◽  
2021 ◽  
pp. 1-69
Author(s):  
Artur Posenato Garcia ◽  
Zoya Heidari

The dielectric response of rocks results from electric double layer (EDL), Maxwell-Wagner (MW), and dipolar polarizations. The EDL polarization is a function of solid-fluid interfaces, pore water, and pore geometry. MW and dipolar polarizations are functions of charge accumulation at the interface between materials with contrasting impedances and the volumetric concentration of its constituents, respectively. However, conventional interpretation of dielectric measurements only accounts for volumetric concentrations of rock components and their permittivities, not interfacial properties such as wettability. Numerical simulations of dielectric response of rocks provides an ideal framework to quantify the impact of wettability and water saturation ( Sw) on electric polarization mechanisms. Therefore, in this paper we introduce a numerical simulation method to compute pore-scale dielectric dispersion effects in the interval from 100 Hz to 1 GHz including impacts of pore structure, Sw, and wettability on permittivity measurements. We solve the quasi-electrostatic Maxwell's equations in three-dimensional (3D) pore-scale rock images in the frequency domain using the finite volume method. Then, we verify simulation results for a spherical material by comparing with the corresponding analytical solution. Additionally, we introduce a technique to incorporate α-polarization to the simulation and we verify it by comparing pore-scale simulation results to experimental measurements on a Berea sandstone sample. Finally, we quantify the impact of Sw and wettability on broadband dielectric permittivity measurements through pore-scale numerical simulations. The numerical simulation results show that mixed-wet rocks are more sensitive than water-wet rocks to changes in Sw at sub-MHz frequencies. Furthermore, permittivity and conductivity of mixed-wet rocks have weaker and stronger dispersive behaviors, respectively, when compared to water-wet rocks. Finally, numerical simulations indicate that conductivity of mixed-wet rocks can vary by three orders of magnitude from 100 Hz to 1 GHz. Therefore, Archie’s equation calibrated at the wrong frequency could lead to water saturation errors of 73%.

2021 ◽  
Author(s):  
Zulkuf Azizoglu ◽  
◽  
Zoya Heidari ◽  

Broadband relative dielectric dispersion measurements are considered interesting options for assessment of water-filled pore volume. Conventional models such as Complex Refractive Index Model (CRIM) and Maxwell Garnett (MG), often overlook or oversimplify the complexity of pore structure, geometrical distribution of the constituting fluids, and spatial distribution of minerals. This yields to significant errors in assessment of water saturation especially in rocks with complex pore structure. Therefore, it becomes important to quantify the impacts of pore structure and spatial distribution of minerals on broadband relative dielectric dispersion measurements to be able to make decisions about reliability of water saturation estimates from these measurements in a given formation. The objectives of this paper are (a) to quantify the impacts of pore structure and spatial distribution of minerals on relative dielectric permittivity measurements in a wide range of frequencies, (b) to propose a new simple and physically meaningful workflow, which honors pore geometry and spatial distribution of minerals to enhance fluid saturation assessment using relative dielectric permittivity measurements, (c) to verify the reliability of the introduced model in the pore-scale domain. First, we perform numerical simulations of relative dielectric dispersion measurements in the frequency range of 20 MHz to 1 GHz in the pore-scale domain. The input to the numerical simulator includes pore-scale images of actual complex carbonate rock samples. We use a physically meaningful model which honors spatial distribution of the rock constituents for the multi-frequency interpretation of relative dielectric response. To verify the reliability of the model in multiple frequencies, we apply the model to the results of relative dielectric simulations in the pore-scale domain on 3D computed tomography scan (CT-scan) images of carbonate rock samples, which are synthetically saturated to obtain a wide range of water saturation. We successfully verified the reliability of the introduced model in the pore-scale domain using carbonate rock samples with multi-modal pore-size distribution. Estimated water saturations from the results of simulations at 1 GHz resulted in an average relative error of less than 4%. We observed measurable improvements in fluid saturation estimates compared to the cases which CRIM or MG models are used. Results demonstrated that application of conventional models to estimate water saturation from relative dielectric response is not reliable in frequencies below 1 GHz.


2016 ◽  
Vol 30 (02) ◽  
pp. 1550268 ◽  
Author(s):  
Jinwei Shi ◽  
Xingbai Luo ◽  
Jinming Li ◽  
Jianwei Jiang

To analyze the process of jet penetration in water medium quantitatively, the properties of jet penetration spaced target with water interlayer were studied through test and numerical simulation. Two theoretical models of jet penetration in water were proposed. The theoretical model 1 was established considering the impact of the shock wave, combined with the shock equation Rankine–Hugoniot and the virtual origin calculation method. The theoretical model 2 was obtained by fitting theoretical analysis and numerical simulation results. The effectiveness and universality of the two theoretical models were compared through the numerical simulation results. Both the models can reflect the relationship between the penetration velocity and the penetration distance in water well, and both the deviation and stability of theoretical model 1 are better than 2, the lower penetration velocity, and the larger deviation of the theoretical model 2. Therefore, the theoretical model 1 can reflect the properties of jet penetration in water effectively, and provide the reference of model simulation and theoretical research.


Polymers ◽  
2021 ◽  
Vol 13 (7) ◽  
pp. 1173
Author(s):  
Ilze Beverte ◽  
Ugis Cabulis ◽  
Sergejs Gaidukovs

As a non-metallic composite material, widely applied in industry, rigid polyurethane (PUR) foams require knowledge of their dielectric properties. In experimental determination of PUR foams’ dielectric properties protection of one-side capacitive sensor’s active area from adverse effects caused by the PUR foams’ test objects has to be ensured. In the given study, the impact of polytetrafluoroethylene (PTFE) films, thickness 0.20 mm and 0.04 mm, in covering or simulated coating the active area of one-side access capacitive sensor’ electrodes on the experimentally determined true dielectric permittivity spectra of rigid PUR foams is estimated. Penetration depth of the low frequency excitation field into PTFE and PUR foams is determined experimentally. Experiments are made in order to evaluate the difference between measurements on single PUR foams’ samples and on complex samples “PUR foams + PTFE film” with two calibration modes. A modification factor and a small modification criterion are defined and values of modifications are estimated in numerical calculations. Conclusions about possible practical applications of PTFE films in dielectric permittivity measurements of rigid PUR foams with one-side access capacitive sensor are made.


2021 ◽  
pp. 204141962110377
Author(s):  
Yaniv Vayig ◽  
Zvi Rosenberg

A large number of 3D numerical simulations were performed in order to follow the trajectory changes of rigid CRH3 ogive-nosed projectiles, impacting semi-infinite metallic targets at various obliquities. These trajectory changes are shown to be related to the threshold ricochet angles of the projectile/target pairs. These threshold angles are the impact obliquities where the projectiles end up moving in a path parallel to the target’s face. They were found to depend on a non-dimensional entity which is equal to the ratio between the target’s resistance to penetration and the dynamic pressure exerted by the projectile upon impact. Good agreement was obtained by comparing simulation results for these trajectory changes with experimental data from several published works. In addition, numerically-based relations were derived for the penetration depths of these ogive-nosed projectiles at oblique impacts, which are shown to agree with the simulation results.


2016 ◽  
Vol 10 (11) ◽  
pp. 203
Author(s):  
Mohd Zaid Othman ◽  
Qasim H. Shah ◽  
Muhammad Akram Muhammad Khan ◽  
Tan Kean Sheng ◽  
M. A. Yahaya ◽  
...  

A series of numerical simulations utilizing LS-DYNA was performed to determine the mid-point deformations of V-shaped plates due to blast loading. The numerical simulation results were then compared with experimental results from published literature. The V-shaped plate is made of DOMEX 700 and is used underneath an armour personal carrier vehicle as an anti-tank mine to mitigate the effects of explosion from landmines in a battlefield. The performed numerical simulations of blast loading of V-shaped plates consisted of various angles i.e. 60°, 90°, 120°, 150° and 180°; variable mass of explosives located at the central mid-point of the V-shaped vertex with various stand-off distances. It could be seen that the numerical simulations produced good agreement with the experimental results where the average difference was about 26.6%.


2014 ◽  
Vol 1042 ◽  
pp. 188-193 ◽  
Author(s):  
Xing Jun Hu ◽  
Jing Chang

In order to analyze the impact of engine cabin parts on aerodynamic characteristics, the related parts are divided into three categories except the engine cooling components: front thin plates (average thickness of 2mm), bottom-suspension and interior panels. The aerodynamic drag coefficient (Cd) were obtained upon the combination schemes consisting of the three types of parts by numerical simulation. Results show that Cd by simulation is closer to the test value gained by the wind tunnel experiment when front thin plates were simplified to the two-dimensional interface with zero thickness. The error is only 5.23%. Meanwhile this scheme reduces grid numbers, thus decreasing the calculating time. As the front thin plates can guide the flow, there is no difference on the Cd values gained from the model with or without bottom-suspension or interior panels when the engine cabin contains the front thin plates; while only both bottom-suspension and interior panels are removed, the Cd value can be reduced when the cabin doesn’t contain the front thin plates.


2013 ◽  
Vol 710 ◽  
pp. 320-324
Author(s):  
Ying Zi Jiang ◽  
Wei Li Wang ◽  
Xue Feng Huang ◽  
Lei Fu ◽  
Zhuang Qing Fan

The numerical simulation of shelled Comp.B explosive was studied following the Lee-Tarver ignition and growth model when it was impacted respectively by 4340 Steel, OFHC and 93#W projectile with the same mass; the influences on explosive detonation of the initiation process, the material of projectile and the L/D ratio of projectile were analyzed; the critical initiation speeds of the projectiles of three different materials with different L/D ratio were gained. In order to verify the simulation results, the questions were calculated by the theoretical simplified model, the results of the theoretical calculation and the numerical simulation accorded well based on critical energy criterion. The results show that the capability of igniting explosive, the first is 93#W, the second is OFHC, the last is 4340 Steel; The initiation point were not on the interface of shell and explosive, and the faster of the impacting velocity, the initiation point closer the interface; the bigger of the L/D ratio of projectile, the higher of the critical initiation speed.


2021 ◽  
Vol 2021 ◽  
pp. 1-9
Author(s):  
Qing Dong ◽  
Zheng-hua Zhou ◽  
Su Jie ◽  
Bing Hao ◽  
Yuan-dong Li

At engineering practice, the theoretical basis for the cross-over method, used to obtain shear wave arrival time in the downhole method of the wave velocity test by surface forward and backward strike, is that the polarity of P-wave keeps the same, while the polarity of S-wave transforms when the direction of strike inverted. However, the characteristics of signals recorded in tests are often found to conflict with this theoretical basis for the cross-over method, namely, the polarity of the P-wave also transforms under the action of surface forward and backward strike. Therefore, 3D finite element numerical simulations were conducted to study the validity of the theoretical basis for the cross-over method. The results show that both shear and compression waves are observed to be in 180° phase difference between horizontal signal traces, consistent with the direction of excitation generated by reversed impulse. Furthermore, numerical simulation results prove to be reliable by the analytic solution; it shows that the theoretical basis for the cross-over method applied to the downhole wave velocity test is improper. In meanwhile, numerical simulations reveal the factors (inclining excitation, geophone deflection, inclination, and background noise) that may cause the polarity of the P-wave not to reverse under surface forward and backward strike. Then, as to reduce the influence factors, we propose a method for the downhole wave velocity test under surface strike, the time difference of arrival is based between source peak and response peak, and numerical simulation results show that the S-wave velocity by this method is close to the theoretical S-wave velocity of soil.


2004 ◽  
Vol 108 (1079) ◽  
pp. 15-26 ◽  
Author(s):  
G. R. Whitehouse ◽  
R. E. Brown

In recent years, various strategies for the concurrent operation of fixed-and rotary-wing aircraft have been proposed as a means of increasing airport capacity. Some of these strategies will increase the likelihood of encounters with the wakes of aircraft operating nearby. Several studies now exist where numerical simulations have been used to assess the impact of encounters with the wakes of large transport aircraft on the safety of helicopter operations under such conditions. This paper contrasts the predictions of several commonly-used numerical simulation techniques when each is used to model the dynamics of a helicopter rotor during the same idealised wake encounter. In most previous studies the mutually-induced distortion of the wakes of the rotor and the interacting aircraft has been neglected, yielding the so-called ‘frozen vortex’ assumption. This assumption is shown to be valid only when the helicopter encounters the aircraft wake at high forward speed. At the low forward speeds most relevant to near-airfield operations, however, injudicious use of the frozen vortex assumption may lead to significant errors in predicting the severity of a helicopter’s response to a wake encounter.


Sign in / Sign up

Export Citation Format

Share Document