A Practical Gas Apparent Permeability Model: Multi-Scale Simulations of Rarefied Gas Flow in Matrix

2021 ◽  
Author(s):  
Di Chai ◽  
Xiaoli Li

Abstract A common practice in gas-shale reservoir simulation, which arbitrarily increases intrinsic matrix permeability to match the production data, has been proven inefficient and unreliable. Alternatively, accurate estimations of gas apparent permeability (AP) in matrix is desired. This work presents an analytical AP model considering rarefaction in nanopores and coupling experimentally confirmed mechanisms in shale matrix for theoretical completeness. Meanwhile, physical terms in AP model are simplified with semi-empirical correlations for the practicability in large-scale field simulation. Compared with other gas transport models in nanopores, the newly-developed analytical model has been successfully validated against molecular dynamic (MD) simulation, direct simulation Monte Carlo (DSMC), Lattice Boltzmann (LB) simulation, and experimental flux results for five types of gases (i.e., methane, nitrogen, helium, argon, and oxygen) with the minimum deviation. It is observed that analytical models excluding Knudsen diffusion mechanism cannot fully characterize rarefaction effect. Next, Knudsen diffusion cannot be explained as the only underlying mechanism of rarefaction because the mass flux is largely underestimated in transition flow regime. However, the weighted superposition of second-order slip boundary and Knudsen diffusion can provide the satisfactory fitting with data. This work provides an analytical model which not only considers non-negligible multi-physics in shale reservoirs (i.e., rarefaction effect, multilayer adsorption, surface diffusion and confinement effect) but also simplifies non-linear physical terms using semi-empirical linear correlations to facilitate AP calculations in core-scale simulations.

2012 ◽  
Vol 710 ◽  
pp. 641-658 ◽  
Author(s):  
Hamed Darabi ◽  
A. Ettehad ◽  
F. Javadpour ◽  
K. Sepehrnoori

AbstractWe study the gas flow processes in ultra-tight porous media in which the matrix pore network is composed of nanometre- to micrometre-size pores. We formulate a pressure-dependent permeability function, referred to as the apparent permeability function (APF), assuming that Knudsen diffusion and slip flow (the Klinkenberg effect) are the main contributors to the overall flow in porous media. The APF predicts that in nanometre-size pores, gas permeability values are as much as 10 times greater than results obtained by continuum hydrodynamics predictions, and with increasing pore size (i.e. of the order of the micrometre), gas permeability converges to continuum hydrodynamics values. In addition, the APF predicts that an increase in the fractal dimension of the pore surface leads to a decrease in Knudsen diffusion. Using the homogenization method, a rigorous analysis is performed to examine whether the APF is preserved throughout the process of upscaling from local scale to large scale. We use the well-known pulse-decay experiment to estimate the main parameter of the APF, which is Darcy permeability. Our newly derived late-transient analytical solution and the late-transient numerical solution consistently match the pressure decay data and yield approximately the same estimated value for Darcy permeability at the typical core-sample initial pressure range and pressure difference. Other parameters of the APF may be determined from independent laboratory experiments; however, a pulse-decay experiment can be used to estimate the unknown parameters of the APF if multiple tests are performed and/or the parameters are strictly constrained by upper and lower bounds.


Aerospace ◽  
2021 ◽  
Vol 8 (6) ◽  
pp. 150
Author(s):  
Yeon-Kyu Park ◽  
Geuk-Nam Kim ◽  
Sang-Young Park

The CANYVAL-C (CubeSat Astronomy by NASA and Yonsei using a virtual telescope alignment for coronagraph) is a space science demonstration mission that involves taking several images of the solar corona with two CubeSats—1U CubeSat (Timon) and 2U CubeSat (Pumbaa)—in formation flying. In this study, we developed and evaluated structural and thermal designs of the CubeSats Timon and Pumbaa through finite element analyses, considering the nonlinearity effects of the nylon wire of the deployable solar panels installed in Pumbaa. On-orbit thermal analyses were performed with an accurate analytical model for a visible camera on Timon and a micro propulsion system on Pumbaa, which has a narrow operating temperature range. Finally, the analytical models were correlated for enhancing the reliability of the numerical analysis. The test results indicated that the CubeSats are structurally safe with respect to the launch environment and can activate each component under the space thermal environment. The natural frequency of the nylon wire for the deployable solar panels was found to increase significantly as the wire was tightened strongly. The conditions of the thermal vacuum and cycling testing were implemented in the thermal analytical model, which reduced the differences between the analysis and testing.


Sensors ◽  
2021 ◽  
Vol 21 (12) ◽  
pp. 4206
Author(s):  
Farhan Nawaz ◽  
Hemant Kumar ◽  
Syed Ali Hassan ◽  
Haejoon Jung

Enabled by the fifth-generation (5G) and beyond 5G communications, large-scale deployments of Internet-of-Things (IoT) networks are expected in various application fields to handle massive machine-type communication (mMTC) services. Device-to-device (D2D) communications can be an effective solution in massive IoT networks to overcome the inherent hardware limitations of small devices. In such D2D scenarios, given that a receiver can benefit from the signal-to-noise-ratio (SNR) advantage through diversity and array gains, cooperative transmission (CT) can be employed, so that multiple IoT nodes can create a virtual antenna array. In particular, Opportunistic Large Array (OLA), which is one type of CT technique, is known to provide fast, energy-efficient, and reliable broadcasting and unicasting without prior coordination, which can be exploited in future mMTC applications. However, OLA-based protocol design and operation are subject to network models to characterize the propagation behavior and evaluate the performance. Further, it has been shown through some experimental studies that the most widely-used model in prior studies on OLA is not accurate for networks with networks with low node density. Therefore, stochastic models using quasi-stationary Markov chain are introduced, which are more complex but more exact to estimate the key performance metrics of the OLA transmissions in practice. Considering the fact that such propagation models should be selected carefully depending on system parameters such as network topology and channel environments, we provide a comprehensive survey on the analytical models and framework of the OLA propagation in the literature, which is not available in the existing survey papers on OLA protocols. In addition, we introduce energy-efficient OLA techniques, which are of paramount importance in energy-limited IoT networks. Furthermore, we discuss future research directions to combine OLA with emerging technologies.


2010 ◽  
Vol 88 (8) ◽  
pp. 815-830 ◽  
Author(s):  
Lesley R. Rutledge ◽  
Stacey D. Wetmore

The present work uses 129 nucleobase – amino acid CCSD(T)/CBS stacking and T-shaped interaction energies as reference data to test the ability of various density functionals with double-zeta quality basis sets, as well as some semi-empirical and molecular mechanics methods, to accurately describe noncovalent DNA–protein π–π and π+–π interactions. The goal of this work is to identify methods that can be used in hybrid approaches (QM/MM, ONIOM) for large-scale modeling of enzymatic systems involving active-site (substrate) π–π contacts. Our results indicate that AMBER is a more appropriate choice for the lower-level method in hybrid techniques than popular semi-empirical methods (AM1, PM3), and suggest that AMBER accurately describes the π–π interactions found throughout DNA–protein complexes. The M06–2X and PBE-D density functionals were found to provide very promising descriptions of the 129 nucleobase – amino acid interaction energies, which suggests that these may be the most suitable methods for describing high-level regions. Therefore, M06–2X and PBE-D with both the 6–31G(d) and 6–31+G(d,p) basis sets were further examined through potential-energy surface scans to better understand how these techniques describe DNA–protein π–π interactions in both minimum and nonminimum regions of the potential-energy surfaces, which is critical information when modeling enzymatic reaction pathways. Our results suggest that studies of stacked nucleobase – amino acid systems should implement the PBE-D/6–31+G(d,p) method. However, if T-shaped contacts are involved and (or) smaller basis sets must be considered due to limitations in computational resources, then M06–2X/6–31G(d) provides an overall excellent description of both nucleobase – amino acid stacking and T-shaped interactions for a range of DNA–protein π–π and π+–π interactions.


Author(s):  
Yancheng Li ◽  
Jianchun Li

This paper presents a recent research breakthrough on the development of a novel adaptive seismic isolation system as the quest for seismic protection for civil structures, utilizing the field-dependent property of the magnetorheological elastomer (MRE). A highly-adjustable MRE base isolator was developed as the key element to form smart seismic isolation system. The novel isolator contains unique laminated structure of steel and MRE layers, which enable its large-scale civil engineering applications, and a solenoid to provide sufficient and uniform magnetic field for energizing the field-dependent property of MR elastomers. With the controllable shear modulus/damping of the MR elastomer, the developed adaptive base isolator possesses a controllable lateral stiffness while maintaining adequate vertical loading capacity. Experimental results show that the prototypical MRE base isolator provides amazing increase of lateral stiffness up to 1630%. Such range of increase of the controllable stiffness of the base isolator makes it highly practical for developing new adaptive base isolation system utilizing either semi-active or smart passive controls. To facilitate the structural control development using the adaptive MRE base isolator, an analytical model was developed to stimulate its behaviors. Comparison between the analytical model and experimental data proves the effectiveness of such model in reproducing the behavior of MRE base isolator, including the observed strain stiffening effect.


Author(s):  
C-M Chen ◽  
R-F Fung

The dynamic equations of a micro-positioning Scott—Russell (SR) mechanism associated with two flexible hinges and an offset are developed to calculate output responses. Both rigid and flexible hinges are considered to explore the results. The main features in the kinematics of the SR mechanism are its displacement amplification and straight-line motion, which are widely needed in practical industries. The manufacturing inaccuracy of the SR mechanism definitely causes geometric offsets of flexure hinges, and affects displacement amplification and straight-line output motion. Analytical models based on kinematics and Hamilton's principle are derived to explore the variation of linearity ratio, magnification factor, and deviation factor due to various offsets and link lengths. From numerical simulations for the SR mechanism with various offsets of flexible hinges in the conditions of different link lengths, it is found that offsets of flexure hinges obviously affect the amplifying factor and linearity ratio, and appear to dominate the changes of magnification factors. Moreover, an analytical model is also used to predict magnification factors due to various offsets. Finally, some conclusions concerning the effects of offset on the performance of the SR mechanism are drawn.


Author(s):  
Chao Liu ◽  
Yan He ◽  
Yufeng Li ◽  
Yulin Wang ◽  
Shilong Wang ◽  
...  

Abstract The residual stresses could affect the ability of components to bear loading conditions and also the performance. The researchers considered workpiece surface as a plane and ignored the effect of surface topography induced by the intermittent cutting process when modeling residual stresses. The aim of this research develops an analytical model to predict workpiece residual stresses during intermittent machining by correlating the effect of surface topography. The relative motions of tool and workpiece are analyzed for modeling thermal-mechanical and surface topography. The influence of dynamic cutting force and thermal on different positions of surface topography is also considered in analytical model. Then the residual stresses model with the surface topography effect can be developed in intermittent cutting. The analytical models of dynamic cutting force, surface topography and residual stresses are verified by the experiments. The variation trend of evaluated values of the residual stress of workpiece is basically consistent with that of measured values. The compressive residual stress of workpiece surface in highest point of the surface topography are higher than that in the lowest point.


2019 ◽  
Vol 141 (5) ◽  
Author(s):  
Georg Geiser ◽  
Jens Wellner ◽  
Edmund Kügeler ◽  
Anton Weber ◽  
Anselm Moors

A nonlinear full-wheel time-domain simulation of a two-stage low pressure turbine is presented, analyzed, and compared with the available experimental data. Recent improvements to the computational fluid dynamics (CFD) solver TRACE that lead to significantly reduced wall-clock times for such large scale simulations are described in brief. Since the configuration is characterized by significant unsteady turbulence and transition effects, it is well suited for the validation and benchmarking of frequency-domain methods. Transition, flow separation and wall pressure fluctuations on the stator blades of the second stage are analyzed in detail. A strong azimuthal π-periodicity is observed, manifesting in a significantly varying stability of the midspan trailing edge flow with a quasi-steady closed separation bubble on certain blades and highly dynamic partially open separation bubbles with recurring transition and turbulent reattachment on other blades. The energy spectrum of fluctuating wall quantities in that regime shows a high bandwidth and considerable disharmonic content, which is challenging for frequency-domain-based simulation methods.


2018 ◽  
Vol 611 ◽  
pp. A89 ◽  
Author(s):  
Yueh-Ning Lee ◽  
Patrick Hennebelle

Context. Understanding the origin of the initial mass function (IMF) of stars is a major problem for the star formation process and beyond. Aim. We investigate the dependence of the peak of the IMF on the physics of the so-called first Larson core, which corresponds to the point where the dust becomes opaque to its own radiation. Methods. We performed numerical simulations of collapsing clouds of 1000 M⊙ for various gas equations of state (eos), paying great attention to the numerical resolution and convergence. The initial conditions of these numerical experiments are varied in the companion paper. We also develop analytical models that we compare to our numerical results. Results. When an isothermal eos is used, we show that the peak of the IMF shifts to lower masses with improved numerical resolution. When an adiabatic eos is employed, numerical convergence is obtained. The peak position varies with the eos, and using an analytical model to infer the mass of the first Larson core, we find that the peak position is about ten times its value. By analyzing the stability of nonlinear density fluctuations in the vicinity of a point mass and then summing over a reasonable density distribution, we find that tidal forces exert a strong stabilizing effect and likely lead to a preferential mass several times higher than that of the first Larson core. Conclusions. We propose that in a sufficiently massive and cold cloud, the peak of the IMF is determined by the thermodynamics of the high-density adiabatic gas as well as the stabilizing influence of tidal forces. The resulting characteristic mass is about ten times the mass of the first Larson core, which altogether leads to a few tenths of solar masses. Since these processes are not related to the large-scale physical conditions and to the environment, our results suggest a possible explanation for the apparent universality of the peak of the IMF.


2009 ◽  
Vol 2009 ◽  
pp. 1-15 ◽  
Author(s):  
Bernard Girau ◽  
César Torres-Huitzil ◽  
Nikolaos Vlassopoulos ◽  
José Hugo Barrón-Zambrano

We consider here the feasibility of gathering multiple computational resources by means of decentralized and simple local rules. We study such decentralized gathering by means of a stochastic model inspired from biology: the aggregation of theDictyostelium discoideumcellular slime mold. The environment transmits information according to a reaction-diffusion mechanism and the agents move by following excitation fronts. Despite its simplicity this model exhibits interesting properties of self-organization and robustness to obstacles. We first describe the FPGA implementation of the environment alone, to perform large scale and rapid simulations of the complex dynamics of this reaction-diffusion model. Then we describe the FPGA implementation of the environment together with the agents, to study the major challenges that must be solved when designing a fast embedded implementation of the decentralized gathering model. We analyze the results according to the different goals of these hardware implementations.


Sign in / Sign up

Export Citation Format

Share Document