Search

In this paper, the counter – current imbibition phenomena in a heterogeneous porous media is studied with the consideration of two types of porous materials like volcanic and fine sand and Adomian decomposition method is applied to find the saturation of wetting phase and the recovery rate of the reservoir. A simulation result is developed here to study the effect of heterogeneity, capillarity and relative permeability on saturation rate and to obtain an optimum recovery rate of the reservoir with the choices of some interesting parametric value.

The goal of this work is to analyze a model for the rate-independent evolution of sets with finite perimeter. The evolution of the admissible sets is driven by that of a given time-dependent set, which has to include the admissible sets and hence is to be understood as an external loading. The process is driven by the competition between perimeter minimization and minimization of volume changes. \par In the mathematical modeling of this process, we distinguish the adhesive case, in which the constraint that the (complement of) the `external load' contains the evolving sets is penalized by a term contributing to the driving energy functional, from the brittle case, enforcing this constraint. The existence of Energetic solutions for the adhesive system is proved by passing to the limit in the associated time-incremental minimization scheme. In the brittle case, this time-discretization procedure gives rise to evolving sets satisfying the stability condition, but it remains an open problem to additionally deduce energy-dissipation balance in the time-continuous limit. This can be obtained under some suitable quantification of data. The properties of the brittle evolution law are illustrated by numerical examples in two space dimensions., Comment: Dedicated to Alexander Mielke on the occasion of his 60th birthday

In this paper, the counter - current imbibition phenomenon is discussed in an inclined heterogeneous porous media with the consideration of two types of porous materials like volcanic sand and fine sand. Adomian decomposition method is applied to find the saturation of wetting phase and the recovery rate of the reservoir. Finally, a simulation result is developed to study the saturation of wetting phase and the optimum recovery rate of reservoir with the choices of some interesting parametric values. This problem has a great importance in the field of oil recovery process.

This contribution presents recent results in the modeling and the analysis of delamination problems. It addresses adhesive contact, brittle, and cohesive zone models both in a quasistatic and a viscous, dynamic setting for the bulk part. Also different evolution laws for the delaminating surface are discussed.

We consider the weak solution of the Laplace equation in a planar domain with a straight crack, prescribing a homogeneous Neumann condition on the crack and a nonhomogeneous Dirichlet condition on the rest of the boundary. For every k we express the k-th derivative of the energy with respect to the crack length in terms of a finite number of coefficients of the asymptotic expansion of the solution near the crack tip and of a finite number of other parameters, which only depend on the shape of the domain.

Phase-field models have already been proven to predict complex fracture patterns in two and three dimensions for brittle fracture at small strains. In this paper we discuss a model for phase-field fracture at finite deformations in more detail. Among the identification of crack location and projection of crack growth the numerical stability is one of the main challenges in solid mechanics. We here present a phase-field model at finite strains, which takes into account the anisotropy of damage by applying an anisotropic split and the modified invariants of the right Cauchy-Green strain tensor. We introduce a suitable weak notion of solution that also allows for a spatial and temporal discretization of the model. In this framework we study the existence of solutions %Second the mathematical background of the approach is examined and and we show that the time-discrete solutions converge in a weak sense to a solution of the time-continuous formulation of the model. Numerical examples in two and three space dimensions are carried out in the range of validity of the analytical results.

The focus of this note lies on the numerical analysis of models describing the propagation of a single crack in a linearly elastic material. The evolution of the crack is modeled as a rate-independent process based on the Griffith criterion. We follow two different approaches for setting up mathematically well defined models: the global energetic approach and an approach based on a viscous regularization. &nbsp We prove the convergence of solutions of fully discretized models (i.e. with respect to time and space) and derive relations between the discretization parameters (mesh size, time step size, viscosity parameter, crack increment) which guarantee the convergence of the schemes. Further, convergence rates are provided for the approximation of energy release rates by certain discrete energy release rates. Thereby we discuss both, models with self-contact conditions on the crack faces as well as models with pure Neumann conditions on the crack faces. The convergence proofs rely on regularity estimates for the elastic fields close to the crack tip and local and global finite element error estimates. Finally the theoretical results are illustrated with some numerical calculations.

Brittle Griffith-type delamination of compounds is deduced by means of Gamma-convergence from partial, isotropic damage of three-specimen-sandwich-structures by flattening the middle component to the thickness 0. The models used here allow for nonlinearly elastic materials at small strains and consider the processes to be unidirectional and rate-independent. The limit passage is performed via a double limit: first, we gain a delamination model involving the gradient of the delamination variable, which is essential to overcome the lack of a uniform coercivity arising from the passage from partial damage to delamination. Second, the delamination-gradient is supressed. Noninterpenetration- and transmission-conditions along the interface are obtained.

We discuss a model for crack propagation in an elastic body, where the crack path is described a-priori. In particular, we develop in the framework of finite-strain elasticity a rate-independent model for crack evolution which is based on the Griffith fracture criterion. Due to the nonuniqueness of minimizing deformations, the energy-release rate is no longer continuous with respect to time and the position of the crack tip. Thus, the model is formulated in terms of the Clarke differential of the energy, generalizing the classical crack evolution models for elasticity with strictly convex energies. We prove the existence of solutions for our model and also the existence of special solutions, where only certain extremal points of the Clarke differential are allowed.

A unilateral contact problem between elastic bodies at small strains glued by a brittle adhesive is addressed in the quasistatic rate-independent setting. The delamination process is modeled as governed by stresses rather than by energies. This leads to a specific scaling of an approximating elastic adhesive contact problem, discretized by a semi-implicit scheme and regularized by a BV-type gradient term. An analytical zero-dimensional example motivates the model and a specific local-solution concept. Two-dimensional numerical simulations performed on an engineering benchmark problem of debonding a fiber in an elastic matrix further illustrate the validity of the model, convergence, and algorithmical efficiency even for very rigid adhesives with high elastic moduli.

We consider the nonlocal formulation of continuum mechanics described by peridynamics. We provide a link between peridynamic evolution and brittle fracture evolution for a broad class of peridynamic potentials associated with unstable peridynamic constitutive laws. Distinguished limits of peridynamic evolutions are identified that correspond to vanishing peridynamic horizon. The limit evolution is associated with dynamic brittle fracture and satisfies a dynamic energy inequality expressed in terms of the kinetic energy of the motion together with a bulk elastic energy and a Griffith surface energy. It corresponds to the simultaneous evolution of elastic displacement and brittle fracture with displacement fields satisfying the wave equation inside the cracking domain. The wave equation provides the dynamic coupling between elastic waves and the evolving fracture path inside the media. The elastic moduli, wave speed and energy release rate for the evolution are explicitly determined by moments of the peridynamic influence function and the peridynamic potential energy., Comment: 26 pages, 2 figures, typos and references added

The nonlinear evolution problem for a crack with a kink in elastic body is considered. This nonlinear formulation accounts the condition of mutual non-penetration between the crack faces. The kinking crack is presented with the help of two unknown shape parameters of the kink angle and of the crack length, which minimize an energy due to the Griffith hypothesis. Based on the obtained results of the shape sensitivity analysis, solvability of the evolutionary minimization problem is proved, and the necessary conditions for the optimal crack are derived.

We study the evolution of a single crack in an elastic body and assume that the crack path is known in advance. The motion of the crack tip is modeled as a rate-independent process on the basis of Griffith's local energy release rate criterion. According to this criterion, the system may stay in a local minimum before it performs a jump. The goal of this paper is to prove the existence of such an evolution and to shed light on the discrepancy between the local energy release rate criterion and models which are based on a global stability criterion (as for example the Francfort/Marigo model). We construct solutions to the local model via the vanishing viscosity method and compare different notions of weak, local and global solutions.

By a combination of geometrical and configurational analysis we study the properties of absolute minimal and equilibrium states of general Mumford-Shah functionals, with applications to models of quasistatic brittle fracture propagation. The main results concern the mathematical relations between physical quantities as energy release rate and energy concentration for 3D cracks with complex shapes, seen as outer measures living on the crack edge., Comment: larger version

Griffith's fracture criterion describes in a quasistatic setting whether or not a pre-existing crack in an elastic body is stationary for given external forces. In terms of the energy release rate (ERR), which is the derivative of the deformation energy of the body with respect to a virtual crack extension, this criterion reads: if the ERR is less than a specific constant, then the crack is stationary, otherwise it will grow. In this paper, we consider geometrically nonlinear elastic models with polyconvex energy densities and prove that the ERR is well defined. Moreover, without making any assumption on the smoothness of minimizers, we rigorously derive the well-known Griffith formula and the J-integral, from which the ERR can be calculated. The proofs are based on a weak convergence result for Eshelby tensors. Copyright © 2007 John Wiley & Sons, Ltd.

The purpose of this paper is to fill the gap between the classical treatment of brittle fracture mechanics and the new idea of considering the crack evolution as a free discontinuity problem. Griffith and Irwin criterions of crack propagation are studied and transformed in order to be no longer dependent on any prescription of the geometry of the crack during its evolution. The inequality contained in theorem 6.1. represents the link between generalized Irwin and Griffith criterions of brittle crack propagation. The physical meaning of this inequality is explained in the last section. Further information available at http://irmi.epfl.ch/cag/buliga_bfrac.html ., Comment: unpublished paper from 1997, circulated in manuscript