Modeling of Solute Transport in a 3D Rough-Walled Fracture–Matrix System

Fluid flow and solute transport in a 3D rough-walled fracture–matrix system were simulated by directly solving the Navier–Stokes equations for fracture flow and solving the transport equation for the whole domain of fracture and matrix with considering matrix diffusion. The rough-walled fracture–matrix model was built from laser-scanned surface tomography of a real rock sample, by considering realistic features of surfaces roughness and asperity contacts. The numerical modeling results were compared with both analytical solutions based on simplified fracture surface geometry and numerical results by particle tracking based on the Reynolds equation. The aim is to investigate impacts of surface roughness on solute transport in natural fracture–matrix systems and to quantify the uncertainties in application of simplified models. The results show that fracture surface roughness significantly increases heterogeneity of velocity field in the rough-walled fractures, which consequently cause complex transport behavior, especially the dispersive distributions of solute concentration in the fracture and complex concentration profiles in the matrix. Such complex transport behaviors caused by surface roughness are important sources of uncertainty that needs to be considered for modeling of solute transport processes in fractured rocks. The presented direct numerical simulations of fluid flow and solute transport serve as efficient numerical experiments that provide reliable results for the analysis of effective transmissivity as well as effective dispersion coefficient in rough-walled fracture–matrix systems. Such analysis is helpful in model verifications, uncertainty quantifications and design of laboratorial experiments.     

Behaviors of Solutions for a Singular Prey–Predator Model and its Shadow System

We study the asymptotic behaviors and quenching of the solutions for a two-component system of reaction–diffusion equations modeling prey–predator interactions in an insular environment. First, we give a global existence result for the solutions to the corresponding shadow system. Then, by constructing some suitable Lyapunov functionals, we characterize the asymptotic behaviors of global solutions to the shadow system. Also, we give a finite time quenching result for the shadow system. Finally, some global existence results for the original reaction–diffusion system are given.     

Analytical study of the load transfer in fibre-reinforced 2D composite materials

The phenomenon of the load diffusion from a fibre to a surrounding matrix is analysed for the 2D case. We use an approximate analytical approach based on the asymptotic reduction of the governing biharmonic problem into two harmonic problems. The comparison of the obtained solutions with known results of other authors shows an acceptable accuracy of the proposed asymptotic simplifications. All solutions are obtained in closed analytical form.The case of perfect bonding between fibre and matrix for a single fibre and for a periodic system of fibres is firstly considered. Then we study the influence of the interface stiffness. In the case when only a single fibre is loaded, the influence of all other fibres is predicted by means of a three-phase model. The proposed approach gives a possibility to solve the problems for a broken fibre and for a broken matrix, as well as for partly debonded fibres. The important problem of infinite matrix cracks is also solved in the present paper.The obtained results can be used for the calculation of pull-out and push-out tests, as well as for the investigation of the fracture of composite materials.     

一类三维随机动力系统的最大Lyapunov指数

基于一维扩散过程的奇异边界理论,使用摄动方法研究了白噪声参激的一类余维二分岔系统的最大Lyapunov指数渐近表达式和数值解,主要讨论了一维相扩散过程同时存在两类奇异边界以及FPK方程存在平稳解的一般性条件. 通过对参激噪声作用项系数矩阵的分析,给出了不变测度的解析解及其相应的Monte Carlo数值仿真结果,并导出了一维相扩散过程P分岔点的确定方法. 对于一类特殊情形,给出了最大Lyapunov指数的渐近表达式;对于参激噪声作用项系数矩阵的一般情形,则给出了系统最大Lyapunov指数的数值结果.      

Antiplane interaction of line crack with an arbitrarily located elliptical inclusion

The antiplane problem of the interaction between a main crack and an arbitrarily located elastic elliptical inclusion near its tip is addressed in the current study. The analysis is based on the use of the complex potentials for the antiplane problem, Laurent series expansion method and an appropriate superposition scheme. The stress intensity factor at the main crack is obtained in a general series form. Explicit asymptotic solutions are also derived by using a perturbation technique and retaining the leading order terms in series expansion. The present solutions are shown to coincide with the Taylor expansion of exact solutions for special cases available in the literature. Discussed are changes in the crack tip stress intensity which can be enhanced or suppressed depending on the location of the elliptical inclusion. The explicit solutions provided herein are well suited for the further quantitative analysis of toughening mechanisms in ceramic composite materials.     

Attractivity in non-autonomous systems

We discuss some sufficient conditions for attractivity of a closed set by using Liapunov functions and show how to apply our results to the study of asymptotic behaviors of solutions in nonautonomous systems. We also discuss the convergence of solutions by the same idea. These results are extendable to functional differential equations.     

Time and Space Fractional Diffusion in Finite Systems

Spatiotemporal nonlocal diffusion in a bounded system is addressed by considering fractional diffusion in a linear, composite system. By considering limiting conditions, solutions for combinations of Neumann and Dirichlet boundary conditions (either zero or nonzero) at the ends of a finite system are derived in terms of Mittag–Leffler functions by the Laplace transformation. Computational viability is demonstrated by inverting the solutions numerically and comparing resulting calculations with asymptotic solutions. Time and space fractional derivatives, defined by variables \(\alpha \) and \(\beta \), respectively, are employed in the Caputo sense; a single-sided, asymmetric space derivative is used. Inspection of the asymptotic solutions leads to insights on the structure of the solutions that may not be available otherwise; the resulting deductions are verified through the numerical inversions. For pure superdiffusion, characteristics of some of the solutions presented here are very similar to those of classical diffusion but combined effects for the corresponding situation result in power-law behaviors. Incidentally, to our knowledge, the pressure distribution for space fractional diffusion at long enough times in a finite system is derived based on first principles for the first time.     

Transient dynamic fracture analysis by an extended meshfree method with different crack-tip enrichments

Investigation of transient dynamic stress intensity factors (DSIFs) of two-dimensional fracture problems of isotropic solids and orthotropic composites by an extended meshfree method is described. We adopt the recently developed extended meshfree radial point interpolation method (X-RPIM), which combines either the standard branch functions or the new linear ramp function associated with Heaviside functions to capture crack-tip behaviors. It is the first time the linear ramp function integrating into meshfree X-RPIM has been presented in a dynamical fracture context. We are particularly interested in exploring insight into the behaviors of DSIFs under dynamic impact loadings (e.g., step, blast and sine loading types) using our meshfree method. For some of these problems numerical examples have been performed using the new ramp functions, and the obtained results of DSIFs have also been compared with those using the standard enrichment functions under which the two schemes have the same setting. In each case it is found that the numerical solutions delivered using the X-RPIM associated with the ramp enrichments are in good agreement with those with the standard functions. The paper first describes formulations and then provides verification of our developed approach through a series of numerical examples in transient dynamic fracture for both solids and orthotropic composites. Illustration of scattered elastic stress waves propagating in the cracked body is depicted to take an insight look at the behavior of dynamic response.     

Équation énergétique locale pour la prevision des niveaux vibratoires et acoustiques instationnaires

Insights and additionnal clarifications concerning the energy diffusion model [4] are in concern herein. We showed in former times under steady state conditions, that the diffusion model is well suited for an uncorrelated plane wave dynamics [7], energy behaviour modelling. We will prove, in this short note, that the later is not enough when transient dynamics is in concern. We thus propose a new equation alternative to the diffusion equation.     

Conduction-radiation effect on transient natural convection with thermophoresis

This paper discusses the effect of thermophoretic particle deposition on the transient natural convection laminar flow along a vertical flat surface,which is immersed in an optically dense gray fluid in the presence of thermal radiation.In the analysis,the radiative heat flux term is expressed by adopting the Rosseland diffusion approximation.The governing equations are reduced to a set of parabolic partial differential equations.Then,these equations are solved numerically with a finite-difference scheme in the entire time regime.The asymptotic solutions are also obtained for sufficiently small and large time.The obtained asymptotic solutions are then compared with the numerical solutions,and they are found in excellent agreement.Moreover,the effects of different physical parameters,i.e.,the thermal radiation parameter,the surface temperature parameter,and the thermophoretic parameter,on the transient surface shear stress,the rate of surface heat transfer,and the rate of species concentration,as well as the transient velocity,temperature,and concentration profiles are shown graphically for a fluid(i.e.,air) with the Prandtl number of 0.7 at 20 C and 1.013 × 10 5 Pa.     

Exponential Propagation for Fractional Reaction–Diffusion Cooperative Systems with Fast Decaying Initial Conditions

We study the time asymptotic propagation of solutions to the reaction–diffusion cooperative systems with fractional diffusion. We prove that the propagation speed is exponential in time, and we find the precise exponent of propagation. This exponent depends on the smallest index of the fractional laplacians and on the principal eigenvalue of the matrix DF(0) where F is the reaction term. We also note that this speed does not depend on the space direction.     

An asymptotic solution for fluid production from an elliptical hydraulic fracture at early-times

The early-time transient flow during the start-up of fluid production from a porous medium by a well intersected by a vertical elliptical hydraulic fracture is studied using an asymptotic analysis. The analysis is focused on the situation of practical interest where the fracture conductivity is high so that production from the fracture dominates. The first three terms in a short-time asymptotic expansion for the production rate during constant-pressure production, and for the well-pressure during constant-rate production, are obtained. It is shown that the fracture tip starts to influence the production rate only when the dimensionless time is increased to the square of the reciprocal of the dimensionless fracture conductivity. The asymptotic results also show that geometric factors of an elliptical fracture introduce non-negligible corrections to the so-called bilinear flow in the early times, which were previously erroneously associated with the effect of the fracture tip.     

Analysis of fractional anomalous diffusion caused by an instantaneous point source in disordered fractal media

A theoretical analysis of fractional anomalous diffusion caused by an instantaneous point source in disordered fractal media is studied. Using the method of symmetry group of scaling transformations and the H-function, the analytical solutions of concentration distribution are given. At the same time we derive the expressions of scattering function spectrum.The result shows that the scattering function spectra still have the properties of scaling function. The scattering functions of point source, line source and area source in regular Euclidean space can be regarded as particular cases of this paper and are included in this paper. At the end of the paper we discuss the asymptotic behaviors of the solution in detail. The results of this paper can be taken to be the fundamental solutions for every kind of boundary value problems of fractional anomalous diffusion in disordered fractal media.     

Solutions to and Validation of Matrix-Diffusion Models

A model transport system is considered in which a pulse of tracer molecules is advected along a flow channel with porous walls. The advected tracer thus undergoes diffusion, matrix-diffusion, inside the walls, which affects the breakthrough curve of the tracer. Analytical solutions in the form of series expansions are derived for a number of situations which include a finite depth of the porous matrix, varying aperture of the flow channel, and longitudinal diffusion and Taylor dispersion of the tracer in the flow channel. Novel expansions for the Laplace transforms of the concentration in the channel facilitated closed-form expressions for the solutions. A rigorous result is also derived for the asymptotic form of the breakthrough curve for a finite depth of the porous matrix, which is very different from that for an infinite matrix. A specific experimental system was created for validation of matrix-diffusion modeling for a matrix of finite depth. A previously reported fracture-column experiment was also modeled. In both cases model solutions gave excellent fits to the measured breakthrough curves with very consistent values for the diffusion coefficients used as the fitting parameters. The matrix-diffusion modeling performed could thereby be validated.     

Capillary effects in steady-state flow in heterogeneous cores

Effects of capillary heterogeneity at the macroscopic scale have previously been analyzed for static conditions or in the context of outflow end-effects. This paper presents a systematic study for the case of one-dimensional, steady-state flow, that complements recent work on transient displacement. We consider the saturation response to various forms of heterogeneity. Included are analytical results for certain model cases, some general results, and numerical solutions for variously correlated spatial variations. The sensitivity to process parameters, such as rate, heterogeneity length scale and correlation, is studied. Physical interpretations are offered and potential applications in the estimation of heterogeneity are discussed.     

Concentration distribution of fractional anomalous diffusion caused by an instantaneous point source

The Fox function expression and the analytic expression for the concentration distribution of fractional anomalous diffusion caused by an instantaneous point source in n-dimensional space (n= 1, 2 or 3) are derived by means of the condition of mass conservation , the time-space similarity of the solution , Mellin transform and the properties of the Fox function . And the asymptotic behaviors for the solutions are also given .     

A variational approach to the fracture of brittle thin films subject to out-of-plane loading

We address the problem of fracture in homogenous linear elastic thin films using a variational model. We restrict our attention to quasi-static problems assuming that kinetic effects are minimal. We focus on out-of-plane displacement of the film and investigate the effect of bending on fracture. Our analysis is based on a two-dimensional model where the thickness of the film does not need to be resolved. We derive this model through a formal asymptotic analysis. We present numerical simulations in a highly idealized setting for the purpose of verification, as well as more realistic micro-indentation experiments.     

Application of Fractional Differential Equations for Modeling the Anomalous Diffusion of Contaminant from Fracture into Porous Rock Matrix with Bordering Alteration Zone

Solute diffusion from a fracture into a porous rock with an altered zone bordering the fracture is modeled by a system of two diffusion equations (one for the altered zone and another for the intact porous matrix) with different coefficients of effective diffusivity. Since experimental studies of diffusion into rock samples with altered zones indicate that mathematical models of diffusion based on Fick’s law do not adequately describe the concentration field in a sample, fractional order diffusion equations are chosen in this study for modeling the anomalous mass transport in the rocks. In the case of significantly higher porosity of the altered zone (e.g., this is typical for carbonates) the effective diffusivity here can be much higher than the effective diffusivity of non-altered rocks. By introducing a small parameter that is the ratio of effective diffusivities in the non-altered and altered regions and applying the technique of perturbations, approximate analytical solutions for concentrations in the altered zone bordering the fracture and in the intact surrounding rocks are obtained. Based on these solutions, different regimes of diffusion into the rocks with different physical properties are modeled and analyzed. It is shown that, using experimentally obtained data, the orders of the fractional derivatives in the differential equations can be readily calibrated for the every specific rock.     

Homogenization for wave propagation in periodic fibre-reinforced media with complex microstructure. I—Theory

The effective elastic properties of periodic fibre-reinforced media with complex microstructure are determined by the method of asymptotic homogenization via a novel solution to the cell problem. The solution scheme is ideally suited to materials with many fibres in the periodic cell. In this first part of the paper we discuss the theory for the most general situation—N arbitrarily anisotropic fibres within the periodic cell. For ease of exposition we then restrict attention to isotropic phases which results in a monoclinic composite material with 13 effective moduli and expressions for each of these are determined. In the second part of this paper we shall discuss results for a variety of specific microstructures.