Search results

The aim of this study is to numerically simulate the density-driven convection in heterogeneous porous media associated with anisotropic permeability field, which is important to the safe and stable long term CO₂ storage in laminar saline aquifers.

The study uses compact finite difference and the pseudospectral method to solve Darcy’s law.

The presence of heterogeneous anisotropy may result in non-monotonic trend of the breakthrough time and quantity of CO₂ dissolved in the porous medium, which are important to the CO₂ underground storage.

The manuscript numerically study the convective phenomena of mixture contained CO₂ and brine. The phenomena are important to the process of CO₂ enhanced oil recovery. Interesting qualitative patterns and quantitative trends are revealed in the manuscript.

Develop a local radial point interpolation method (LRPIM) to analyze the dissipation process of excess pore water pressure in porous media and verify its numerical capability.

Terzaghi's consolidation theory is used to describe the dissipation process. A local residual form is formulated over only a sub‐domain. This form is spatially discretized by radial point interpolation method (RPIM) with basis of multiquadrics (MQ) and thin‐plate spline (TPS), and temporally discretized by finite difference method. One‐dimensional (1D) and two‐dimensional consolidation problems are numerically analyzed.

The LRPIM is suitable, efficient and accurate to simulate this dissipation process. The shape parameters, q=1.03, R=0.1 for MQ and η=4.001 for TPS, are still valid.

The asymmetric system matrix in LRPIM spends more resources in storage and CPU time.

Local residual form requires no background mesh, thus being a truly meshless method. This provides a fast and practical algorithm for engineering computation.

This paper provides a simple, accurate and fast numerical algorithm for the dissipation process of excess pore water pressure, largely simplifies data preparation, shows that the shape parameters from solid mechanics are also suitable for the dissipation process.

In petroleum industry, hydraulic fracturing is essential to enhance oil productivity. The hydraulic fractures are usually generated in the process of hydraulic fracturing. Although some mathematical models were proposed to analyze the well-flow behavior of conventional fracture, there are few models to depict unconventional fracture like reorientation fracture. To figure out the effect of reorientation fracture on production enhancement and guide the further on-site operating, this paper aims to investigate the well-flow behavior of vertical reorientation fracture in horizontal permeability anisotropic reservoir.

Based on the governing equation considering horizontal permeability anisotropy, the mathematical models for reorientation fractures in infinite reservoir are developed by using the principle of superposition. Furthermore, a rectangular closed drainage area is also considered to investigate the well-flow behavior of reorientation fracture, and the mathematical models are developed by using Green’s and source functions.

Computational results indicate that the flux distribution of infinite conductivity fracture is uniform at very early times. After a period, it will stabilize eventually. High permeability anisotropy and small inclination angle of reorientation will cause significant end point effect in the infinite conductivity fracture. The reorientation fractures with small inclination angle in high anisotropic reservoir are capable of improving 1-1.5 times more oil productivity in total.

This paper develops the mathematical methods to study the well-flow behavior for unconventional fracture, especially for reorientation fracture. The results validate the production enhancement effect of reorientation fracture and identify the sensitive parameters of productivity.

Natural convection from a semi‐infinite vertical plate embedded in a fluid saturated porous medium is studied both analytically and numerically. The plate is assumed to be heated isothermally or by a constant heat flux. The porous medium, modeled according to Darcy’s law, is anisotropic in permeability with its principal axes oriented in a direction that is oblique to the gravity vector. In the large Rayleigh number limit, the governing boundary‐layer equations are solved in closed form, using a similarity transformation. Comparisons between the numerical solution of the full equations and analytical solutions are presented for a wide range of the governing parameters. The effects of the anisotropic permeability ratio K^*, of the orientation angle of the principal axes θ, and of the Rayleigh number R_H on the flow and heat transfer are investigated. Results indicate that the anisotropic properties of the porous medium considerably modify the heat transfer, velocity and temperature profiles from that expected under isotropic conditions.

The natural two‐dimensional convection taking place between horizontal concentric cylinders filled with a satured anisotropic porous medium is studied numerically. The anisotropy concerns exclusively the permeability. Isothermal boundary conditions are applied on both inner and outer boundaries, with the outer boundary being warmer. The effects of the anisotropic permeability ratio K^*, of the orientation angle of the principal axes defined by γ, and of the Rayleigh number Ra^* on the flow and heat transfer are investigated. Results indicate that a net circulating flow around the annulus is generated, except for values of γ that preserve the symmetry of flow conditions with respect to the vertical diameter. It is also shown that the anisotropic part of the resistivity tensor is equivalent to a magnetic resistivity tensor.

This paper aims to investigate the impact of combining grain-oriented electrical steel (GOES) grades on specific iron losses and the flux density distribution within a single-phase magnetic core.

This paper presents the results of finite-element method (FEM) simulations investigating the impact of mixing two different GOES grades on losses of a single-phase magnetic core. The authors used different models: a 3D model with a highly detailed geometry including both saturation and anisotropy, as well as a simplified 2D model to save computation time. The behavior of the flux distribution in the mixed magnetic core is analyzed. Finally, the results from the numerical simulations are compared with experimental results.

The specific iron losses of a mixed magnetic core exhibit a nonlinear decrease with respect to the GOES grade with the lowest losses. Analyzing the magnetic core behavior using 2D and 3D FEM shows that the rolling direction of the GOES grades plays a critical role on the nonlinearity variation of the specific losses.

The novelty of this research lies in achieving an optimum trade-off between the manufacturing cost and the core efficiency by combining conventional and high-performance GOES grade in a single-phase magnetic core.

This paper aims to use a history-dependent vector stop hysteresis model incorporated into a two dimensional finite elements (FE) simulation environment to solve the magnetic field problems in electrical machines. The vector stop hysteresis model is valid for representing the anisotropic magnetization characteristics of electrical steel sheets. Comparisons of the simulated results with measurements show that the model is well appropriate for the simulation of electrical machines with alternating, rotating and harmonic magnetic flux densities.

The anisotropy of the permeability of an electrical steel sheet can be represented by integrating anhysteretic surfaces into the elastic element of a vector hysteresis stop model. The parameters of the vector stop hysteresis model were identified by minimizing the errors between the simulated results and measurements. In this paper, a damped Newton method is applied to solve the nonlinear problem, which ensures a robust convergence of the finite elements simulation with vector stop hysteresis model.

Analyzing the measurements of the electrical steel sheets sample obtained from a rotational single sheet tester shows the importance to consider the anisotropic and saturation behavior of the material. Comparing the calculated and measured data corroborates the hypothesis that the presented energy-based vector stop hysteresis model is able to represent these magnetic properties appropriately. To ensure a unique way of hysteresis loops during finite elements simulation, the memory of the vector stop hysteresis model from last time step is kept unchanged during the Newton iterations.

The results of this work demonstrates that the presented vector hysteresis stop model allows simulation of vector hysteresis effects of electrical steel sheets in electrical machines with a limited amount of measurements. The essential properties of the electrical steel sheets, such as phase shifts, the anisotropy of magnetizations and the magnetization characteristics by alternating, rotating, harmonic magnetization types, can be accurately represented.

For gravity dams built on foundations with directional joint sets, the seepage in the foundation possesses anisotropic characteristics and may have adverse effects on the foundation stability. A methodology for system reliability analysis of gravity dam foundations considering anisotropic seepage and multiple sliding surfaces is proposed in this paper.

Anisotropic seepages in dam foundations are simulated using finite element method (FEM) with the equivalent continuum model (ECM), and their effect on dam foundation stability is involved by uplift pressures acting on the potential sliding surfaces. The system failure probability of the dam foundation is efficiently estimated using Monte Carlo method (MCM) combined with response surface method (RSM).

The case study shows that it is necessary to consider the possibly adverse effect of anisotropic seepage on foundation stability of gravity dams and the deterministic analysis of the foundation stability may be misleading. The system reliability analysis of the dam foundation is justified, as the uncertainties in shear strength parameters of the foundation rocks and joint sets as well as aperture, connectivity and spacing of the joint sets are quantified and the system effect of the multiple potential sliding surfaces on the foundation reliability is reasonably considered.

(1) A methodology is proposed for efficient system reliability analysis of foundation stability of gravity dams considering anisotropic seepage and multiple sliding surfaces (2) The influence of anisotropic seepage on the stability of gravity dam foundation  is revealed (3) The influence of estimation errors of RSMs on the system reliability assessment of dam foundation is investigated.

The ordinary vector hysteresis stop model with constant threshold values is not able to prohibit the hysteretic property after the saturation correctly. This paper aims to develop an improved vector hysteresis stop model with threshold surfaces. This advanced anisotropic vector hysteresis stop model can represent the magnetic saturation properties and the hysteresis losses under alternating and rotating magnetizations.

By integrating anhysteretic surfaces into the elastic element of a vector hysteresis stop model, the anisotropy of the permeability of an electrical steel sheet can be represented. Instead of the commonly used constant threshold value for plastic elements of the hysteresis model, threshold surfaces are applied to the stop hysterons. The threshold surfaces can be derived directly from measured alternating major loops of the material sample. By saturated polarization, the constructed threshold surfaces are vanishing. In this way, the reversible magnetic flux density is in the same direction of the applied magnetic flux density. Thus, the saturation properties are satisfied.

Analyzing the measurements of the electrical steel sheets sample obtained from a rotational single sheet tester shows that the clockwise (CW) and counter-CW (CCW) rotational hysteresis losses decrease by saturated flux density. At this state, instead of the domain wall motion, the magnetization rotation is dominant in the material. As a result, the hysteresis losses, which are related to the domain wall motion, are vanished near the saturation. In one stop operator, the plastic element represents the hysteresis part of the model. Integrating threshold surface into the plastic element, the hysteresis part can be modified to zero near the saturation to represent the saturation properties.

The results of this work demonstrate that the presented vector hysteresis stop model allows simulation of anisotropic hysteresis effects, alternating and rotating hysteresis losses. The parameters of the hysteresis model are determined by comparing the measured and modeled minor loops in different alternating magnetization directions. With the identified parameters, the proposed model is excited with rotated excitations in CW and CCW directions. The rotated hysteresis losses, derived from the model, are then compared with those experimentally measured. The modified vector stop model can significantly improve the accuracy of representing hysteresis saturations and losses.

To analyze effectively magnetic shielding effects by shields with fine structure.

Simplification of the fine structure makes it possible to analyze them efficiently. The authors have introduced a homogenization method to estimate effective permeability of magnetic composite structure for the static field. The homogenization method is applied to the analysis of magnetic shields composed of steel plates and steel rods against DC power lines to test its feasibility.

The properties of the magnetic shielding are analyzed by using the homogenization method. The errors of the magnetic fields increase in case of very few layers.

The simplification of the magnetic shields with fine structure by using the homogenization method makes it possible to analyze efficiently magnetic shielding effects, although the accuracy becomes worse in case of very few layers.