Initial Anisotropy Research Articles

Investigation on the Influencing Factors of K0 of Granular Materials Using Discrete Element Modelling

Earth pressure coefficient at rest K0 is commonly estimated by empirical equations, which to date has had insufficient accuracy and universality. For better prediction, the investigation on the factors influencing K0 is required. A series of discrete element method (DEM) simulations of oedometer tests are conducted to verify the key factors influencing K0 of granular materials. The influences of initial fabric anisotropy, particle shape, initial void ratio, inter-particle friction angle is investigated. The evolution of microstructure is monitored during the tests to reveal the relationship between the microstructure evolution and K0 values. The results show that the effect of fabric anisotropy exists but is limited. Particle shape, initial void ratio, and inter-particle friction angle all significantly affect the K0 values alone. According to the DEM results, an attempt is made to propose a more reasonable empirical equation in which K0 is a function of relative density, critical state friction angle, and “shape factor”. This new empirical equation has higher accuracy and can consider the effect of particle shape, inspiring the determination of K0 values in practical engineering.

Coupled Geomechanics and Flow Modeling of Fractured Reservoirs considering Matrix Permeability Anisotropy

The effect of geomechanics is crucial in the modeling of fractured reservoirs since fractures can be more stress-sensitive than the rock matrix. In fractured reservoirs, the microscale fractures are often homogenized into the matrix continuum leading to matrix permeability anisotropy. The embedded discrete fracture model (EDFM) is becoming increasingly popular for numerical simulation of fractured reservoirs and has been successfully employed in the simulation of coupled geomechanics and flow systems. However, the anisotropic permeability of matrix cannot be considered in traditional EDFM. An approach of coupled geomechanics and flow modeling is proposed for fractured reservoirs considering anisotropic matrix permeability. In order to calculate the fluid exchange between fractures and rock matrix in an anisotropic formation, an integrally embedded discrete fracture model (IEDFM) is used. The geomechanics model of the proposed approach uses an equivalent continuum approach, which introduces an equivalent material to represent the overall deformation of the fractured rock under normal and shear stresses. The constitutive relations of the equivalent continuum are derived rigorously from stress-strain analysis, where the stress-dependent moduli of natural fractures are included. The coupled geomechanics and flow system is solved using the fixed-stress split iterative coupling strategy with the dynamic hydraulic parameters of matrix and fractures updated separately. Several examples are performed to demonstrate the applicability of the proposed approach for modeling the coupled geomechanics and flow system in fractured reservoirs considering anisotropic permeability. The effect of anisotropy is investigated, which indicates that the dynamic behavior of a fracture is highly orientation-related due to initial stress anisotropy and matrix permeability anisotropy. Simulations also show that the anisotropic matrix permeability affects the compaction in the reservoir domain, which reflects on the performance of production.

Damage Analysis of Third-Generation Advanced High-Strength Steel Based on the Gurson–Tvergaard–Needleman (GTN) Model

The third generation of advanced high-strength steels (AHSS) brought attention to the steel and automotive industries due to its good compromise between formability and production costs. This work evaluated a third-generation AHSS (USS CR980XG3TM) through microstructural and X-ray diffraction (XRD) analyses, uniaxial tensile and plane-strain tension testing, and numerical simulations. The damage behavior of this steel is described with the Gurson–Tvergaard–Needleman (GTN) model using an identification procedure based on the uniaxial tensile and initial microvoids data. The microstructure of the CR980XG3TM steel is composed of ferrite, martensite–austenite islands, and retained austenite with a volume fraction of 12.2%. The global formability of the CR980XG3TM steel, namely the product of the uniaxial tensile strength and total elongation values, is 24.3 GPa%. The Lankford coefficient shows a weak initial plastic anisotropy of the CR980XG3TM steel with the in-plane anisotropy close to zero (−0.079) and the normal anisotropy close to unity (0.917). The identified GTN parameters for the CR980XG3TM steel provided a good forecast for the limit strains defined according to ISO 12004-2 standard from the uniaxial tensile and plane-strain tension data.

Probing the alignment-dependent mechanical behaviors and time-evolutional aligning process of collagen scaffolds.

Efficiently manipulating and reproducing collagen (COL) alignment in vitro remains challenging because many of the fundamental mechanisms underlying and guiding the alignment process are not known. We reconcile experiments and coarse-grained molecular dynamics simulations to investigate the mechanical behaviors of a growing COL scaffold and assay how changes in fiber alignment and various cross-linking densities impact their alignment dynamics under shear flow. We find higher cross-link densities and alignment levels significantly enhance the apparent tensile/shear moduli and strength of a bulk COL system, suggesting potential measures to facilitate the design of stronger COL based materials. Since fibril alignment plays a key factor in scaffold mechanics, we next investigate the molecular mechanism behind fibril alignment with Couette flow by computationally investigating the effects of COL's structural properties such as chain lengths, number of chains, tethering conditions, and initial COL conformations on the COL's final alignment level. Our computations suggest that longer chain lengths, more chains, greater amounts of tethering, and initial anisotropic COL conformations benefit the final alignment, but the effect of chain lengths may be more dominant over other factors. These results provide important parameters for consideration in manufacturing COL-based scaffolds where alignment and cross-linking are necessary for regulating performance.

The behavior of elastic moduli with fluid content in the VTI media

For addressing the behavior of a reservoir with different fluid types, the Biot-Gassmann equation often is the base of the practical simulation. Despite the prevalence application of this equation with the isotropy conditions, the unforeseen errors always expose in simulation results because of the anisotropy state in reality. We investigated the anisotropy model with the integration of two analytical strategies using a three-component VSP data set. We obtained an initial anisotropy model in the region of acquired walkaway VSP using slowness polarization inversion, and updated the anisotropy model with the application of anisotropy ray-tracing and tomography. We applied a layer-stripping approach to the anisotropy model during raytracing to optimize the inversion. Given a computed geomechanical model and extracted rock properties of a carbonate reservoir, we developed the anisotropy Biot-Gassmann model, for finding the elastic moduli. We used the substitution strategy to generate the dynamic model of elastic moduli. We showed how the compressional modulus and rigidity change with the anisotropy model in different fluid content. We found that integration of slowness polarization and raytracing tomography increases the maneuverability to control the predicted anisotropy model and intensifies the convergence rate of the inverse problem. We observed that the isotropy assumption in modeling the elastic parameters makes around 8–10 % drift value in compressional modulus relevant to the reality, whereas rigidity showed reluctant behavior to fluid.

DEM models using direct and indirect shape descriptions for Toyoura sand along monotonous loading paths

Two different DEM models are proposed for quantitatively simulating Toyoura sand macroscopic response along various monotonous loading paths and for a wide range of initial densities. The first model adopts spherical particles and compensates for the irregular shapes of Toyoura sand grains by adding an additional rolling resistance stiffness to the classical linear contact model. The second model follows a different strategy whereby rolling stiffness is abandoned in favor of more complex shapes in the form of a few different 3D polyhedrons defined from a 2D micrograph of Toyoura particles. After a preliminary analysis of the number of particles for optimal REV simulations, the two different modeling approaches are calibrated using triaxial compression in so-called drained conditions, adopting a common contact friction angle for the two models. Similar predictive abilities are then obtained along so-called undrained (constant volume) triaxial compression and extension paths. Although it leads to 9-times longer simulations, the polyhedral approach is easier to calibrate regarding the contact parameters. It also enables a more precise description of the microstructure in terms of particle shapes and initial fabric anisotropy, whose crucial role is evidenced in a parametric analysis.

Evolution towards energy equipartition in star clusters: effects of the tidal field, primordial binaries, and internal velocity anisotropy

ABSTRACT This paper is the second in a series investigating the evolution of star clusters towards energy equipartition (EEP). Here, we focus on the effects of the external tidal field of the host galaxy, initial anisotropy in the velocity distribution, and primordial binary star population. The results of our N-body simulations show that regardless of the strength of the tidal field or the fraction of primordial binaries: (i) the evolution towards EEP in the intermediate and outer regions of initially anisotropic systems is more rapid than for isotropic systems; (ii) this evolution also proceeds at different rates for the tangential and radial components of the velocity dispersion; and (iii) the outer regions of the initially isotropic systems show a tendency to evolve towards a state of ‘inverted’ EEP in which low-mass stars have smaller velocity dispersion than high-mass stars. We also find that the clusters with primordial binaries stay even farther from EEP than systems containing only single stars. Finally, we show that all these results also hold when the degree of EEP is calculated using quantities measured in projection as it is done in observational studies, and that our findings could be tested with current and upcoming observational data.

The Structuring and Empowerment Policy of Losari Beach Street Vendors Based on The Concept of Smart City

Mining activities in open pit and underground mines will always be associated with rock breaking or stripping activities (both mechanical and blasting), so that this can affect the structure and strength of rocks. The strength of the rock is strongly influenced by the presence of initial cracks (pre-existing cracks) and rock anisotropy conditions associated with discontinuous plane conditions. Fracture mechanics is a science that illustrates how a fracture can occur and propagate during applied stress on material. The main parameter in fracture mechanics is called fracture toughness which shows the resistance of the material to propagate the crack. There are several mode in determining type I fracture toughness, one of which is type I fracture toughness Flattened Brazilian Disc (FBD) mode. Type I fracture toughness test is carried out using a compression machine in a laboratory and is conducted on concrete samples consisting of 3 (three) various samples, with a ratio of cement and sand composition of 1:1, 1:2, and 2:1. This test also uses different loading rate values, namely 2.50 mm/min, 2.70 mm/min, and 2.83 mm/min. The results of the type I fracture toughness value from each loading rate will be compared to determine the effect of the loading rate on the value of type I fracture toughness. The obtained fracture toughness value is also related to the physical and mechanical properties of the samples. Based on the results of tests, it can be seen that the loading rate affects the value of fracture toughness, the increase in fracture toughness value is followed by the higher loading rate. In addition, it can be seen that the fracture toughness value is directly proportional to the uniaxial compressive strength value and the indirect tensile strength value. The average correlation value obtained is R2 = 0.9884 (indicating a strong relationship).

Design of Groundwater Sustain Irrigation Network System in Panca Lautang Area, Sidrap Regency, South Sulawesi

Mining activities in open pit and underground mines will always be associated with rock breaking or stripping activities (both mechanical and blasting), so that this can affect the structure and strength of rocks. The strength of the rock is strongly influenced by the presence of initial cracks (pre-existing cracks) and rock anisotropy conditions associated with discontinuous plane conditions. Fracture mechanics is a science that illustrates how a fracture can occur and propagate during applied stress on material. The main parameter in fracture mechanics is called fracture toughness which shows the resistance of the material to propagate the crack. There are several mode in determining type I fracture toughness, one of which is type I fracture toughness Flattened Brazilian Disc (FBD) mode. Type I fracture toughness test is carried out using a compression machine in a laboratory and is conducted on concrete samples consisting of 3 (three) various samples, with a ratio of cement and sand composition of 1:1, 1:2, and 2:1. This test also uses different loading rate values, namely 2.50 mm/min, 2.70 mm/min, and 2.83 mm/min. The results of the type I fracture toughness value from each loading rate will be compared to determine the effect of the loading rate on the value of type I fracture toughness. The obtained fracture toughness value is also related to the physical and mechanical properties of the samples. Based on the results of tests, it can be seen that the loading rate affects the value of fracture toughness, the increase in fracture toughness value is followed by the higher loading rate. In addition, it can be seen that the fracture toughness value is directly proportional to the uniaxial compressive strength value and the indirect tensile strength value. The average correlation value obtained is R2 = 0.9884 (indicating a strong relationship).

The Effect of Loading Rate on Fracture Toughness of Concrete using Type I Fracture Toughness Test Flattened Brazilian Disc (FBD) Mode

Mining activities in open pit and underground mines will always be associated with rock breaking or stripping activities (both mechanical and blasting), so that this can affect the structure and strength of rocks. The strength of the rock is strongly influenced by the presence of initial cracks (pre-existing cracks) and rock anisotropy conditions associated with discontinuous plane conditions. Fracture mechanics is a science that illustrates how a fracture can occur and propagate during applied stress on material. The main parameter in fracture mechanics is called fracture toughness which shows the resistance of the material to propagate the crack. There are several mode in determining type I fracture toughness, one of which is type I fracture toughness Flattened Brazilian Disc (FBD) mode. Type I fracture toughness test is carried out using a compression machine in a laboratory and is conducted on concrete samples consisting of 3 (three) various samples, with a ratio of cement and sand composition of 1:1, 1:2, and 2:1. This test also uses different loading rate values, namely 2.50 mm/min, 2.70 mm/min, and 2.83 mm/min. The results of the type I fracture toughness value from each loading rate will be compared to determine the effect of the loading rate on the value of type I fracture toughness. The obtained fracture toughness value is also related to the physical and mechanical properties of the samples. Based on the results of tests, it can be seen that the loading rate affects the value of fracture toughness, the increase in fracture toughness value is followed by the higher loading rate. In addition, it can be seen that the fracture toughness value is directly proportional to the uniaxial compressive strength value and the indirect tensile strength value. The average correlation value obtained is R2 = 0.9884 (indicating a strong relationship).

Non-flow effects in correlation between harmonic flow and transverse momentum in nuclear collisions

A large anti-correlation signal between elliptic flow v2 and average transverse momentum [pT] was recently measured in small collision systems, consistent with a final-state hydrodynamic response to the initial geometry. This negative v2–[pT] correlation was predicted to change to positive correlation for events with very small charged particle multiplicity Nch due to initial-state momentum anisotropies of the gluon saturation effects. However, the role of non-flow correlations is expected to be important in these systems, which is not yet studied. We estimate the non-flow effects in pp, pPb and peripheral PbPb collisions using Pythia and Hijing models, and compare them with the experimental data. We show that the non-flow effects are largely suppressed using the rapidity-separated subevent cumulant method (details of the cumulant framework are also provided). The magnitude of the residual non-flow is much less than the experimental observation in the higher Nch region, supporting the final-state response interpretation. In the very low Nch region, however, the sign and magnitude of the residual non-flow depend on the model details. Therefore, it is unclear at this moment whether the sign change of v2–[pT] can serve as evidence for initial state momentum anisotropies predicted by the gluon saturation.

Theoretical analysis of long-term setup of driven piles considering soil anisotropy, destructuration and creep effects

The soil adjacent to the pile suffers a severe destructuration during pile installation and then regains its strength during the setup stage. This paper presents a theoretical approach to evaluate the long-term setup of driven piles in soft sensitive clays. An advanced elastic-viscoplastic model that accounts for fabric anisotropy, particle bonding and viscosity, is used in both stages of pile installation and setup. A semi-analytical solution for the undrained cylindrical and spherical cavity expansion in soft sensitive clays is developed to capture the pile installation effects. The changes in the stress state and fabric anisotropy of the soil adjacent to the pile caused by pore pressure dissipation and sustained service loads are analyzed. The evolution of soil mechanical behavior is incorporated into the load-transfer models of shaft friction and base resistance. Comparisons between measured and predicted results of field tests at three different sites demonstrate the applicability of the proposed approach. The calculated results illustrate that the ultimate bearing capacity increases steadily with time during the long-term setup stage. An extensive study is carried out to evaluate the effects of soil's initial anisotropy and bonding on the load-settlement behavior of driven piles.

Numerical simulation of underground excavations in an indurated clay using non-local regularisation. Part 2: sensitivity analysis

A sensitivity study is presented to evaluate the influence of different parameters on the simulation of an underground excavation in the Callovo-Oxfordian (COx) argillaceous formation performed in the Meuse/Haute-Marne underground research laboratory. An elasto-viscoplastic constitutive law representing the characteristic behaviour of indurated mudrocks and stiff clays has been employed. It incorporates anisotropy, strain-softening, creep deformations and dependence of permeability on damage. In addition, a non-local formulation, able to simulate localised deformations objectively, has been incorporated in the analyses. The following features affecting the excavation have been studied: initial stress, strength and stiffness anisotropy, strength parameters, hydraulic and hydromechanical parameters, and scale effects. A simulation reported in a companion paper provides the base case for benchmarking. The results are compared in terms of extent and configuration of the excavation fractured zone, vertical and horizontal tunnel convergences, and the development and evolution of pore pressures in the rock. From the comparisons, an enhanced understanding of the hydromechanical mechanisms associated with underground excavations in COx claystone, and other similar argillaceous materials, has been achieved.

Influences of initial anisotropy and principal stress rotation on the undrained monotonic behavior of a loose silica sand

Triaxial compression and extension tests have been conducted under different initial anisotropy conditions to investigate the undrained response of a crushed silica sand. The loose to medium specimens were prepared using the moist tamping method. Five stress paths with different stress ratios (q/p′) were employed to prepare anisotropically consolidated specimens. Several specimens were consolidated under a specific condition in which a stress rotation occurred under undrained monotonic shearing similar to a reversed cyclic shear stress loading during an earthquake. The effects of initial induced anisotropy at consolidation on the onset of liquefaction, phase transformation, and critical state are investigated within the framework of Anisotropic Critical State Soil Mechanics (ACSSM). In addition, fabric evolution during shearing towards the critical state is evaluated using bidirectional bender element tests. The results illustrate the fact that there is a unique anisotropic critical state representing anisotropic fabric, irrespective of initial anisotropy, and the states of stress. Similar to the critical state line, the phase transformation line has the same loci for different initial anisotropies.

An elastoplastic solution for cylindrical cavity expansion under constant water content conditions in anisotropic unsaturated soils

In this paper, an elastoplastic solution is developed for cylindrical cavity expansion in anisotropic unsaturated soils under constant water content conditions. The most prominent features of the solution can be attributed to 1) the stress–strain behaviour of unsaturated soil is subtly coupled with the saturation-suction behaviour; 2) both the mechanical and hydraulic behaviours are modelled as elastoplastic processes; and 3) an anisotropic critical state soil model for saturated soil is extended to an unsaturated situation to consider the effects of the initial stress anisotropy and stress-induced anisotropy of the natural unsaturated soils. The problem is formulated as a system of first-order differential equations and solved as an initial value problem. A numerical code that can judge the hydraulic yield condition of soil particles during the cavity expansion process is developed based on the Runge-Kutta algorithm to solve the governing equations. Parametric analyses, which cover a wide range of suctions, overconsolidation ratios, degrees of initial stress anisotropy, are conducted to investigate the effects of these important factors on the expansion responses in unsaturated soils. The proposed solution could potentially be applied to interpret pressuremeter tests and model pile installation effects in anisotropic unsaturated soils.

Anisotropic flow decorrelation in heavy-ion collisions with event-by-event viscous hydrodynamics

Decorrelation of the elliptic flow in rapidity is calculated within a hybrid approach which includes event-by-event viscous fluid dynamics and final state hadronic cascade model. The simulations are performed for Au+Au collisions at center-of-mass collision energies of 27 and 200 GeV per nucleon pair, as well as various colliding systems at 72 GeV per nucleon pair. Initial conditions determined by an extended Monte Carlo Glauber model show better agreement with experimental data than initial conditions from the UrQMD transport model. We show how the observed decorrelation is connected with the decorrelation of initial state spatial anisotropies. We also study how the effect is increased by the final state hadronic cascade.

PT dependence of the correlation between initial spatial anisotropy and final momentum anisotropies in relativistic heavy ion collisions

The particle momentum anisotropy (vn) produced in relativistic nuclear collisions is considered to be a response of the initial geometry or the spatial anisotropy ϵn of the system formed in these collisions. The linear correlation between ϵn and vn quantifies the efficiency at which the initial spatial eccentricity is converted to final momentum anisotropy in heavy ion collisions. We study the transverse momentum and collision centrality dependence of this correlation for charged particles using a hydrodynamical model framework at LHC. The (ϵn−vn) correlation is found to be stronger for central collisions and also for n=2 compared to that for n=3 as expected. However, the transverse momentum (pT) dependent correlation coefficient shows interesting features which strongly depends on the mass as well as pT of the emitted particle. The correlation strength is found to be larger for lighter particles in the lower pT region. We see that the relative fluctuation in anisotropic flow depends strongly on the value of η/s specially in the region pT<1 GeV unlike the correlation coefficient which does not show significant dependence on η/s.

Anisotropic microstructural evolution and coarsening in free sintering and constrained sintering of metal film by using FIB-SEM tomography

The anisotropic microstructure develops during powder processing and constrained sintering. FIB-SEM tomography was used to investigate the microstructural evolutions of Au submicron particles during free sintering and constrained sintering. The decrease in total surface area and the coarsening, i.e. the increase of mean intercept length of solid phase, were observed during the densification process. The anisotropic packing structure of spherical particles induced by casting process was characterized by area weighted fabric tensor, which represented the bond orientation and the anisotropy of contact area. In free sintering, the initial anisotropy in packing structure decreased with densification. In constrained sintering, the initial anisotropy observed by the mean intercept length of pore was reversed in the later stage. The microstructure evolved so as to form the elongated pores along the thickness direction. An open pore structure model was presented to explain the microstructural evolution in the intermediate stage.

Analysis of stress influence and plastic strain on magnetic properties during the forming process

The aim of this paper is to analyze the relation between magnetic and mechanical properties during and after forming processes. For this purpose, several tensile tests were carried out on sheet metal samples up to a defined plastic strain. The specimens were left in the clamping device in order to relieve the force in several steps until the specimen was completely relieved. As a consequence, the gradual relief leads to a reduction of internal stress states. During the forming process, the initial magnetic relative permeability and magnetic anisotropy of the sample were measured several times. Both properties are related to the mechanical states in the material through the effects of magnetic embrittlement and magneto-elasticity. The plastic strain of the specimens was determined by optical measurements and the stresses in the measurement range during the tensile test was determined with the help of a subsequent numerical simulation. This made it possible for the first time to measure the magnetic properties of samples with different plastic strain and different stress states. The evaluation shows that there is a strong correlation between permeability and plastic strain as well as anisotropy and stress. Based on these findings, it has been confirmed, that the determination of the plastic strain by a soft sensor is possible.

Probes of the quark-gluon plasma and plasma instabilities

Penetrating probes in heavy-ion collisions, like jets and photons, are sensitive to the transport coefficients of the produced quark-gluon plasma, such as shear and bulk viscosity. Quantifying this sensitivity requires a detailed understanding of photon emission and jet-medium interaction in a nonequilibrium plasma during the hydrodynamic stages of heavy-ion collisions. Up to now, such an understanding has been hindered by plasma instabilities. These instabilities arise out of equilibrium and lead to spurious divergences when evaluating the rate of interaction of hard probes with the plasma. In this paper, we show that taking into account the time evolution of an unstable plasma cures these divergences. Specifically, we calculate the time evolution of gluon two-point correlators in a setup with a small initial momentum anisotropy and show that the gluon occupation density at first grows exponentially. Based on this calculation, we argue for a phenomenological prescription where instability poles are subtracted. Finally, we show that in the Abelian case instability fields do not affect medium-induced photon emission to our order of approximation.