Using fracture-based continuum modeling of coupled geomechanical-hydrological processes for numerical simulation of hydraulic fracturing

This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.

Nonlinear wave dispersion in monoatomic chains with lumped and distributed masses:discrete and continuum models

The main objective of this paper is to investigate the influence of inertia of nonlinear springs on the dispersion behavior of discrete monoatomic chains with lumped and distributed masses.The developed model can represent the wave propagation problem in a non-homogeneous material consisting of heavy inclusions embedded in a matrix.The inclusions are idealized by lumped masses,and the matrix between adjacent inclusions is modeled by a nonlinear spring with distributed masses.Additionally,the model is capable of depicting the wave propagation in bi-material bars,wherein the first material is represented by a rigid particle and the second one is represented by a nonlinear spring with distributed masses.The discrete model of the nonlinear monoatomic chain with lumped and distributed masses is first considered,and a closed-form expression of the dispersion relation is obtained by the second-order Lindstedt-Poincare method(LPM).Next,a continuum model for the nonlinear monoatomic chain is derived directly from its discrete lattice model by a suitable continualization technique.The subsequent use of the second-order method of multiple scales(MMS)facilitates the derivation of the corresponding nonlinear dispersion relation in a closed form.The novelties of the present study consist of(i)considering the inertia of nonlinear springs on the dispersion behavior of the discrete mass-spring chains;(ii)developing the second-order LPM for the wave propagation in the discrete chains;and(iii)deriving a continuum model for the nonlinear monoatomic chains with lumped and distributed masses.Finally,a parametric study is conducted to examine the effects of the design parameters and the distributed spring mass on the nonlinear dispersion relations and phase velocities obtained from both the discrete and continuum models.These parameters include the ratio of the spring mass to the lumped mass,the nonlinear stiffness coefficient of the spring,and the wave amplitude.

Static-to-kinematic modeling and experimental validation of tendon-driven quasi continuum manipulators with nonconstant subsegment stiffness

Continuum robots with high flexibility and compliance have the capability to operate in confined and cluttered environments. To enhance the load capacity while maintaining robot dexterity, we propose a novel non-constant subsegment stiffness structure for tendon-driven quasi continuum robots(TDQCRs) comprising rigid-flexible coupling subsegments.Aiming at real-time control applications, we present a novel static-to-kinematic modeling approach to gain a comprehensive understanding of the TDQCR model. The analytical subsegment-based kinematics for the multisection manipulator is derived based on screw theory and product of exponentials formula, and the static model considering gravity loading,actuation loading, and robot constitutive laws is established. Additionally, the effect of tension attenuation caused by routing channel friction is considered in the robot statics, resulting in improved model accuracy. The root-mean-square error between the outputs of the static model and the experimental system is less than 1.63% of the arm length(0.5 m). By employing the proposed static model, a mapping of bending angles between the configuration space and the subsegment space is established. Furthermore, motion control experiments are conducted on our TDQCR system, and the results demonstrate the effectiveness of the static-to-kinematic model.

High-order harmonic generation of ZnO crystals in chirped and static electric fields

High harmonic generation in ZnO crystals under chirped single-color field and static electric field are investigated by solving the semiconductor Bloch equation(SBE). It is found that when the chirp pulse is introduced, the interference structure becomes obvious while the harmonic cutoff is not extended. Furthermore, the harmonic efficiency is improved when the static electric field is included. These phenomena are demonstrated by the classical recollision model in real space affected by the waveform of laser field and inversion symmetry. Specifically, the electron motion in k-space shows that the change of waveform and the destruction of the symmetry of the laser field lead to the incomplete X-structure of the crystal-momentum-resolved(k-resolved) inter-band harmonic spectrum. Furthermore, a pre-acceleration process in the solid four-step model is confirmed.

High-order targeted essentially non-oscillatory scheme for two-fluid plasma model

The weakly ionized plasma flows in aerospace are commonly simulated by the single-fluid model,which cannot describe certain nonequilibrium phenomena by finite collisions of particles,decreasing the fidelity of the solution.Based on an alternative formulation of the targeted essentially non-oscillatory(TENO)scheme,a novel high-order numerical scheme is proposed to simulate the two-fluid plasmas problems.The numerical flux is constructed by the TENO interpolation of the solution and its derivatives,instead of being reconstructed from the physical flux.The present scheme is used to solve the two sets of Euler equations coupled with Maxwell's equations.The numerical methods are verified by several classical plasma problems.The results show that compared with the original TENO scheme,the present scheme can suppress the non-physical oscillations and reduce the numerical dissipation.

Calibration of quantitative rescattering model for simulating vortex high-order harmonic generation driven by Laguerre–Gaussian beam with nonzero orbital angular momentum

We calibrate the macroscopic vortex high-order harmonic generation(HHG)obtained by the quantitative rescattering(QRS)model to compute single-atom induced dipoles against that by solving the time-dependent Schr?dinger equation(TDSE).We show that the QRS perfectly agrees with the TDSE under the favorable phase-matching condition,and the QRS can accurately predict the main features in the spatial profiles of vortex HHG if the phase-matching condition is not good.We uncover that harmonic emissions from short and long trajectories are adjusted by the phase-matching condition through the time-frequency analysis and the QRS can simulate the vortex HHG accurately only when the interference between two trajectories is absent.This work confirms that it is an efficient way to employ the QRS model in the single-atom response for precisely simulating the macroscopic vortex HHG.

General Kinetostatic Modeling and Deformation Analysis of a Two-Module Rod-Driven Continuum Robot with Friction Considered

Continuum robots actuated by flexible rods have large potential applications,such as detection and operation tasks in confined environments,since the push and pull actuation of flexible rods withstand tension and compressive force,and increase the structure's rigidity.In this paper,a generalized kinetostatics model for multi-module and multi-segment continuum robots considering the effect of friction based on the Cosserat rod theory is established.Then,the model is applied to a two-module rod-driven continuum robot with winding ropes to analyze its deformation and load characteristics.Four different in-plane configurations under the external load term as S1,S2,C1,and C2 are defined.Taking a bending plane as an example,the tip deformation along thex-axis of these shapes is simulated and compared,which shows that the load capacity of C1 and C2 is generally larger than that of S1 and S2.Furthermore,the deformation experiments and simulations show that the maximum error ratio without external loads relative to the total length is no more than 3%,and it is no more than 4.7%under the external load.The established kinetostatics model is proven sufficient to accurately analyze the rod-driven continuum robot with the consideration of internal friction.

On iterative learning control design for tracking iteration-varying trajectories with high-order internal model

In this paper, iterative learning control （ILC） design is studied for an iteration-varying tracking problem in which reference trajectories are generated by high-order internal models （HOLM）. An HOlM formulated as a polynomial operator between consecutive iterations describes the changes of desired trajectories in the iteration domain and makes the iterative learning problem become iteration varying. The classical ILC for tracking iteration-invariant reference trajectories, on the other hand, is a special case of HOlM where the polynomial renders to a unity coefficient or a special first-order internal model. By inserting the HOlM into P-type ILC, the tracking performance along the iteration axis is investigated for a class of continuous-time nonlinear systems. Time-weighted norm method is utilized to guarantee validity of proposed algorithm in a sense of data-driven control.

Numerical modeling of wave equation by a truncated high-order finite-difference method

Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.

Variable Curvature Modeling Method of Soft Continuum Robots with Constraints

The inherent compliance of continuum robots holds great promise in the fields of soft manipulation and safe human–robot interaction.This compliance reduces the risk of damage to the manipulated object and its surroundings.However,continuum robots possess theoretically infinite degrees of freedom,and this high flexibility usually leads to complex deformations when subjected to external forces and positional constraints.Describing these complex deformations is the main challenge in modeling continuum robots.In this study,we investigated a novel variable curvature modeling method for continuum robots,considering external forces and positional constraints.The robot configuration curve is described using the developed mechanical model,and then the robot is fitted to the curve.A ten-section continuum robot prototype with a length of 1 m was developed in order to validate the model.The feasibility and accuracy of the model were verified by the ability of the robot to reach target points and track complex trajectories with a load.This work was able to serve as a new perspective for the design analysis and motion control of continuum robots.

A High-Order Conservative Semi-Lagrangian Solver for 3D Free Surface Flows with Sediment Transport on Voronoi Meshes

In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons.

Multi-factor high-order intuitionistic fuzzy timeseries forecasting model

Fuzzy sets theory cannot describe the neutrality degreeof data, which has largely limited the objectivity of fuzzy time seriesin uncertain data forecasting. With this regard, a multi-factor highorderintuitionistic fuzzy time series forecasting model is built. Inthe new model, a fuzzy clustering algorithm is used to get unequalintervals, and a more objective technique for ascertaining membershipand non-membership functions of the intuitionistic fuzzy setis proposed. On these bases, forecast rules based on multidimensionalintuitionistic fuzzy modus ponens inference are established.Finally, contrast experiments on the daily mean temperature ofBeijing are carried out, which show that the novel model has aclear advantage of improving the forecast accuracy.

A High-Order Conservative Numerical Method for Gross-Pitaevskii Equation with Time-Varying Coefficients in Modeling BEC

We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation （BEC）. This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.

The density wave in a new anisotropic continuum model

In this paper the new continuum traffic flow model proposed by Jiang et al is developed based on an improved car-following model, in which the speed gradient term replaces the density gradient term in the equation of motion. It overcomes the wrong-way travel which exists in many high-order continuum models. Based on the continuum version of car-following model, the condition for stable traffic flow is derived. Nonlinear analysis shows that the density fluctuation in traffic flow induces a variety of density waves. Near the onset of instability, a small disturbance could lead to solitons determined by the Korteweg-de-Vries （KdV） equation, and the soliton solution is derived.

Motivation in Foreign Language Learning: Applying the Time Continuum Model of Motivation in FLL

The research is: by using Wdolkowski's Time Continuum Model throughout a lesson plan enables the teacher to increase students'motivation and help them move closer to success in a learning environment. This research supports the theory that instruction is a network of interactions between the teacher and learner that promotes a successful learning experience. It identifies a three-part learning sequence-a beginning, middle and an end. Each part has two of six key motivational factors that when applied correctly by the teacher will maximize the success and continued motivation of the learner.

A review on the application of modified continuum models in modeling and simulation of nanostructures

Analysis of the mechanical behavior of nanos- tructures has been very challenging. Surface energy and non- local elasticity of materials have been incorporated into the traditional continuum analysis to create modified continuum mechanics models. This paper reviews recent advancements in the applications of such modified continuum models in nanostructures such as nanotubes, nanowires, nanobeams, graphenes, and nanoplates. A variety of models for these nanostructures under static and dynamic loadings are men- tioned and reviewed. Applications of surface energy and nonlocal elasticity in analysis of piezoelectric nanomateri- als are also mentioned. This paper provides a comprehensive introduction of the development of this area and inspires fur- ther applications of modified continuum models in modeling nanomaterials and nanostructures.

Aero-engine Blade Fatigue Analysis Based on Nonlinear Continuum Damage Model Using Neural Networks

Fatigue life and reliability of aero-engine blade are always of important significance to flight safety.The establishment of damage model is one of the key factors in blade fatigue research.Conventional linear Miner＇s sum method is not suitable for aero-engine because of its low accuracy.A back propagation neutral network（BPNN） based on the combination of Levenberg-Marquardt（LM） and finite element method（FEM） is used to describe process of nonlinear damage accumulation behavior in material and predict fatigue life of the blade.Fatigue tests of standard specimen made from TC4 are carried out to obtain material fatigue parameters and S-N curve.A nonlinear continuum damage model（CDM）,based on the BPNN with one hidden layer and ten neurons,is built to investigate the nonlinear damage accumulation behavior,in which the results from the tests are used as training set.Comparing with linear models and previous nonlinear models,BPNN has the lowest calculation error in full load range.It has significant accuracy when the load is below 500 MPa.Especially,when the load is 350 MPa,the calculation error of the BPNN is only 0.4%.The accurate model of the blade is built by using 3D coordinate measurement technology.The loading cycle in fatigue analysis is defined from takeoff to cruise in 10 min,and the load history is obtained from finite element analysis（FEA）.Then the fatigue life of the compressor blade is predicted by using the BPNN model.The final fatigue life of the aero-engine blade is 6.55 104 cycles（10 916 h） based on the BPNN model,which is effective for the virtual design of aero-engine blade.

Seepage simulation of high concrete-faced rockfill dams based on generalized equivalent continuum model

This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.

A high-order splitting scheme for the advection diffusion equation

A high-order splitting scheme for the advection-diffusion equation of pollutants is proposed in this paper. The multidimensional advection-diffusion equation is splitted into several one-dimensional equations that are solved by the scheme. Only three spatial grid points are needed in each direction and the scheme has fourth-order spatial accuracy. Several typically pure advection and advection-diffusion problems are simulated. Numerical results show that the accuracy of the scheme is much higher than that of the classical schemes and the scheme can he efficiently solved with little programming effort.

Verification of a laboratory-based dilation model for in situ conditions using continuum models

With respect to constitutive models for continuum modeling applications, the post-yield domain remainsthe area of greatest uncertainty. Recent studies based on laboratory testing have led to thedevelopment of a number of models for brittle rock dilation, which account for both the plastic shearstrain and confining stress dependencies of this phenomenon. Although these models are useful inproviding an improved understanding of how dilatancy evolves during a compression test, there hasbeen relatively little work performed examining their validity for modeling brittle rock yield in situ. Inthis study, different constitutive models for rock dilation are reviewed and then tested, in the context of anumber of case studies, using a continuum finite-difference approach （FLAC）. The uncertainty associatedwith the modeling of brittle fracture localization is addressed, and the overall ability of mobilizeddilation models to replicate in situ deformation measurements and yield patterns is evaluated. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.