Search results

Understanding the fluid flow through rock masses, which commonly consist of rock matrix and fractures, is a fundamental issue in many application areas of rock engineering. As the equivalent porous medium approach is the dominant approach for engineering applications, it is of great significance to estimate the equivalent permeability tensor of rock masses. This study aims to develop a novel numerical approach to estimate the equivalent permeability tensor for fractured porous rock masses.

The radial point interpolation method (RPIM) and finite element method (FEM) are coupled to simulate the seepage flow in fractured porous rock masses. The rock matrix is modeled by the RPIM, and the fractures are modeled explicitly by the FEM. A procedure for numerical experiments is then designed to determinate the equivalent permeability tensor directly on the basis of Darcy’s law.

The coupled RPIM-FEM method is a reliable numerical method to analyze the seepage flow in fractured porous rock masses, which can consider simultaneously the influences of fractures and rock matrix. As the meshes of rock matrix and fracture network are generated separately without considering the topology relationship between them, the mesh generation process can be greatly facilitated. Using the proposed procedure for numerical experiments, which is designed directly on the basis of Darcy’s law, the representative elementary volume and equivalent permeability tensor of fractured porous rock masses can be identified conveniently.

A novel numerical approach to estimate the equivalent permeability tensor for fractured porous rock masses is proposed. In the approach, the RPIM and FEM are coupled to simulate the seepage flow in fractured porous rock masses, and then a numerical experiment procedure directly based on Darcy’s law is introduced to estimate the equivalent permeability tensor.

A study on heat and mass transfer behavior for an anisotropic porous medium embedded in square cavity/annulus is conducted using incompressible smoothed particle hydrodynamics (ISPH) method. In the case of square cavity, the left wall has hot temperature T_h and mass C_h and the right wall have cool temperature T_c and mass C_c and both of the top and bottom walls are adiabatic. While in the case of square annulus, the inner surface wall is considered to have a cool temperature T_c and mass C_c while the outer surface is exposed to a hot temperature T_h and mass C_h. The paper aims to discuss these issues.

The governing partial differential equations are transformed to non-dimensional governing equations and are solved using ISPH method. The results present the influences of the Dufour and Soret effects on the fluid flow and heat and mass transfer.

The effects of various physical parameters such as Darcy parameter, permeability ratio, inclination angle of permeability and Rayleigh numbers on the temperature and concentration profiles together with the local Nusselt and Sherwood numbers are presented graphically. The results from the current ISPH method are well-validated and have favorable comparisons with previously published results and solutions by the finite volume method.

A study on heat and mass transfer behavior on an anisotropic porous medium embedded in square cavity/annulus is conducted using Incompressible Smoothed Particle Hydrodynamics (ISPH) method. In the ISPH algorithm, a semi-implicit velocity correction procedure is utilized, and the pressure is implicitly evaluated by solving pressure Poisson equation (PPE). The evaluated pressure has been improved by relaxing the density invariance condition to formulate a modified PPE.

Perfectly matched layers (PMLs) are used for reflectionless truncation of the problem boundaries in finite element methods applications. In this paper, the method of PMLs is extended to truncate any lossless medium and the method is implemented for the T formulation.

The basic concept behind PMLs is to create an artificial material with a complex and diagonal anisotropic permittivity and permeability. For the A, V formulation PMLs are well known.

It is shown that it is possible to truncate any linear lossless material with PML layers, and if the material has small losses the PML works fine.

In the present paper, an artificial anisotropic lossy material is applied to a 3D edge finite‐element T, formulation to form perfectly matched layers.

Presents a fully coupled numerical model to simulate the slow transient phenomena involving heat and mass transfer in deforming partially saturated porous materials. Makes use of the modified effective stress concept together with the capillary pressure relationship. Examines phase changes (evaporation‐condensation(, heat transfer through conduction and convection, as well as latent heat transfer. The governing equations in terms of gas pressure, capillary pressure, temperature and displacements are coupled non‐linear differential equations and are discretized by the finite element method in space and by finite differences in the time domain. The model is further validated with respect to a documented experiment on partially saturated soil behaviour, and the effects of two‐phase flow, as compared to the one‐phase flow solution, are analysed. Two other examples involving drying of a concrete wall and thermoelastic consolidation of partially saturated clay demonstrate the importance of proper physical modelling and of appropriate choice of the boundary conditions.

The use of gradient-based methods in finite element schemes can be prevented by undefined derivatives, which are encountered when modeling hysteresis in constitutive material laws. This paper aims to present a method to deal with this problem.

Non-smooth Newton methods provide a generalized framework for the treatment of minimization problems with undefined derivatives. Within this paper, a magnetostatic finite element formulation that includes hysteresis is presented. The non-linear equations are solved using a non-smooth Newton method.

The non-smooth Newton method shows promising convergence behavior when applied to a model problem. The numbers of iterations for magnetization curves with and without hysteresis are within the same range.

Mathematical tools like Clarke's generalized Jacobian are applied to magnetostatic field problems with hysteresis. The relation between the non-smooth Newton method and other methods for solving non-linear systems with hysteresis like the M(B)-iteration is established.

The tensor complex permeability of a multi-turn coil with elliptic cross-section is analytically expressed. In field analysis, a multi-turn coil can be modeled by the uniform material that has the present tensor complex permeability. It is shown that the frequency characteristic of the present tensor complex permeability is in good agreement with that evaluated by finite element method applied to a unit cell of the multi-turn coil region.

The authors introduce a new method to evaluate the complex permeability of a multi-turn rectangular coil. To obtain the complex permeability of a rectangular coil in a closed form, it is approximated as an elliptic coil. Because the rectangular coil has different complex permeabilities in the vertical and horizontal directions, the complex permeability have to be defined in a tensor form. It suffices to discretize the coil region into rather coarse finite elements without considering the skin depth in contrast to the conventional finite element method.

The proposed method is shown to give the impedance of multi-turn coils which is in good agreement with results obtained by the conventional finite element (FE) analysis. By extending the proposed approach, the authors can easily perform 3D FE analysis without difficulty in discretization of the coil region with fairly fine finite elements. Moreover, they found that the approximation of rectangular coils as the elliptic coils is valid for analysis of quasi-static fields using this homogenization method.

The novelty of this study is in the approximation of the rectangular coils with elliptic coils, and the complex permeability for them is formulated here in a closed form. The proposed formula includes that for the round coils. Using the present method, the authors analyze the rectangular coils without fine discretization.

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.

The resin transfer molding process for composites manufacturing consists of either of two considerations, namely, the fluid flow analysis through a porous fiber preform where the location of the flow front is of fundamental importance, and the combined flow/heat transfer/cure analysis. In this paper, the continuous sensitivity formulations are developed for the process modeling of composites manufactured by RTM to predict, analyze, and optimize the manufacturing process. Attention is focused here on developments for isothermal flow simulations, and various illustrative examples are presented for sensitivity analysis of practical applications which help serve as a design tool in the process modeling stages.

The purpose of this paper is to present a fully automated numerical tool for computing the effective permeability of porous media from digital images which come from the modern imagery technique.

The permeability is obtained by the homogenization process applied to a periodic rigid solid in which the fluid flow is described by the Stokes equations. The unit cell problem is solved by using the Fast Fourier Transform (FFT) algorithm, well adapted for the microstructures defined by voxels.

Various 3-D examples are considered to show the capacity of the method. First, the case of flow through regular arrays of aligned cylinders or spheres are considered as benchmark problems. Next, the method is applied to some more complex and realistic porous solids obtained with assemblies of overlapping spherical pores having identical or different radii, regularly or randomly distributed within the unit cell.

The use of FFT allows the resolution of high-dimension problems and open various possibilities for computing the permeability of porous microstructures coming from X-ray microtomography.

The purpose of this paper is to study the microwave behaviour of effective magnetic permeability for two‐component ferrite like metamaterial medium in the direction of a biasing magnetic field. The metamaterial medium is presented as an infinite host dielectric material (air) with periodically embedded ferric cylindrical and spherical inclusions saturated with an external dc magnetic field. The study is based on the effective medium theory developed for polycrystalline metaferrites. The simulations show that the presented metamaterial can exhibit the ultra‐low refractive index (ULI) phenomenon and the phenomenon of negative magnetic permeability for the case of microwave propagation in the direction of bias.

The obtained results are based on the wave long approximation of permeability tensor of the presented metamaterial media obtained earlier by the first author. Using the standard approach, the authors apply the above expressions for the microwave propagation in direction of biasing dc magnetic field considering different polarization of the incident microwave.

The considered artificial material media can become either material with a ULI or with negative values in the GHz frequencies.

The paper is concerned with part of the theory of a new generation of artificial ferrites.