Linear numerical calculation method for obtaining critical point,pore fluid,and framework parameters of gas-bearing media

Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method.

SOLVERS FOR SYSTEMS OF LARGE SPARSE LINEAR AND NONLINEAR EQUATIONS BASED ON MULTI-GPUS

Numerical treatment of engineering application problems often eventually results in a solution of systems of linear or nonlinear equations.The solution process using digital computational devices usually takes tremendous time due to the extremely large size encountered in most real-world engineering applications.So,practical solvers for systems of linear and nonlinear equations based on multi graphic process units（GPUs）are proposed in order to accelerate the solving process.In the linear and nonlinear solvers,the preconditioned bi-conjugate gradient stable（PBi-CGstab）method and the Inexact Newton method are used to achieve the fast and stable convergence behavior.Multi-GPUs are utilized to obtain more data storage that large size problems need.

Oscillation of Systems of Parabolic Differential Equations with Deviating Arguments

To study a class of boundary value problems of parabolic differential equations with deviating arguments, averaging technique, Green’s formula and symbol function sign(·) are used. The multi dimensional problem was reduced to a one dimensional oscillation problem for ordinary differential equations or inequalities. Two oscillatory criteria of solutions for systems of parabolic differential equations with deviating arguments are obtained.

Iterative Solution for Systems of Binary Non-mixed Monotone Operator Equations

With using the cone and partial ordering t heory and mixed monotone operator theory, the existence and uniqueness for solut ion of systems of non-monotone binary nonliear operator equations are discussed. And the iterative sequences which converge to solution of systems of operator e quations and the error estimates are also given. Some corresponding results for the mixed monotone operations and the unary operator equations are improved and generalized.

Hot deformation behaviors of 35% SiC_p/2024Al metal matrix composites

The hot deformation behaviors of 35%SiCp/2024 aluminum alloy composites were studied by hot compression tests using Gleeble-1500D thermo-mechanical simulator at temperatures ranging from 350 to 500 °C under strain rates of 0.01-10 s-1. The true stress-true strain curves were obtained in the tests. Constitutive equation and processing map were established. The results show that the flow stress decreases with the increase of deformation temperature at a constant strain rate, and increases with the increase of strain rate at constant temperature, indicating that composite is a positive strain rate sensitive material. The flow stress behavior of composite during hot compression deformation can be represented by a Zener-Hollomon parameter in the hyperbolic sine form. Its activation energy for hot deformation Q is 225.4 kJ/mol. To demonstrate the potential workability, the stable zones and the instability zones in the processing map were identified and verified through micrographs. Considering processing map and microstructure, the hot deformation should be carried out at the temperature of 500 °C and the strain rate of 0.1-1 s-1.

Global Attractor of a Reaction-diffusion System with Hamiltonian Structure

In this paper, we consider a reaction diffusion system with Hamitonian structure, we first prove the existence of an invariant region for system and the continuity of the semigroup, then establish the absorbing sets and global attractor.

Blow-Up Rate for a Semilinear Parabolic System

This paper deals with the blow-up rate of positive solution for a semilinearparabolic system coupled in the equations and boundary condition. The upper and lower bounds ofblow-up rates are obtained.

Dynamic spheroidization kinetics behavior of Ti-6.5Al-2Zr-1Mo-1V alloy with lamellar microstructure

Abstract： The dynamic spheroidization kinetics behavior of Ti-6.5Al-2Zr-1Mo-1V alloy with a lamellar initial microstructure was studied by isothermal hot compression tests in the temperature range of 750-950℃ and strain rates of 0.001-10 s^-1. The results show that the spheroidized fraction of alpha lamellae increases with the increase of temperature and the decrease of strain rate. The spheroidization kinetics curves predicted by JMAK equation agree well with experimental ones. The corresponding SEM and TEM observations indicate that the dynamic spheroidization process can be divided into two stages. The primary stage is boundary splitting formed by two competing mechanisms which are dynamic recrystallization and mechanical twin. In the second stage, the penetration of beta phase into the alpha/alpha grain boundaries is dominant, which is controlled in nature by diffusion of the chemical elements such as Al, Mo and V.

Determination of pollution point source in parabolic system model

This paper considers an inverse problem for a partial differential equation to identify a pollution point source in a watershed. The mathematical model of the problem is a weakly coupled system of two linear parabolic equations for the concentrations u（x, t） and v（x, t） with an unknown point source F（x, t） = A（ t）δ（x- s） related to the concentration u（x, t）, where s is the point source location and A（t） is the amplitude of the pollution point source. Assuming that source F becomes inactive after time T＊, it is proved that it can be uniquely determined by the indirect measurements { v（0, t）, v（ a, t）, v（ b, t）, v（ l, t）, 0 〈 t ≤ T, T＊ 〈 T}, and, thus, the local Lipschitz stability for this inverse source problem is obtained. Based on the proof of its uniqueness, an inversion scheme is presented to determine the point source. Finally, two numerical examples are given to show the feasibility of the inversion scheme.

Global Classical Solutions for One Dimensional Navier Stokes Equations

In this article the author considers Cauchy problem for one dimensional Navier Stokes equations and the global smooth resolvablity for classical solutions is obtained.

Existence of Positive Solutions for a Class of Degenerate and Nondegenerate Stable Diffusion Models

This paper is concerned with a class of degenerate and nondegenerate stable diffusion models.By using the upper and lower solution method and Schauder fixed point principle,the author studies the existence of positive solutions for these stable_diffusion models under some conditions.

Minor Self-conjugate and Skewpositive Semidefinite Solutions to a System of Matrix Equations over Skew Fields

Minor self conjugate (msc) and skewpositive semidefinite (ssd) solutions to the system of matrix equations over skew fields [A mn X nn =A mn ,B sn X nn =O sn ] are considered. Necessary and sufficient conditions for the existence of and the expressions for the msc solutions and the ssd solutions are obtained for the system.

Advancing Hierarchical Equations of Motion for Efficient Evaluation of Coherent Two-dimensional Spectroscopy

To advance hierarchical equations of motion as a standard theory for quantum dissipative dynamics, we put forward a mixed Heisenberg-SchrSdinger scheme with block-matrix implementation on efficient evaluation of nonlinear optical response function. The new approach is also integrated with optimized hierarchical theory and numerical filtering algorithm. Different configurations of coherent two-dimensional spectroscopy of model excitonic dimer systems are investigated, with focusing on the effects of intermolecular transfer coupling and bi-exciton interaction.

Global existence and blow up of a degenerate parabolic system

This paper deals with positive solutions of a degenerate parabolic system: u t= Δ u m+ v p ln α(h+u), v t= Δ v n+u q ln β(h+v) with homogeneous Dirichlet boundary conditions and positive initial conditions. This system describes the processes of diffusion of heat and burning in two component continuous media with nonlinear conductivity and volume energy release. We obtain the global existence and blow up results of the solution relying on comparison with carefully constructed upper solutions and lower solutions.

Variable damping constraint tomography and its application in VSP Data

First arrival travel time tomography has achieved wide application. However, tomographic resolution is insufficient because geometry constraints cause rays to be unevenly distributed in the velocity model. The variable damping constraint method adopts uneven priori information to match uneven data distribution which can lessen the correlation between velocity correction values and ray coverage density. In this paper, we combine the variable damping constraint with a smoothness constraint which is added into the regularization equations in velocity inversion to avoid instability caused by only using the variable damping constraint method. The alpha-trimmed-mean filter is used to smooth and denoise intermediate results in the velocity inversion process. We use the LSQR algorithm to enhance the convergence rate and suppress error propagation in solving linear equations. In this paper, we apply the proposed tomographic method to perform velocity inversion using VSP data. The application in recovery test of the checkerboard model and velocity inversion of real VSP data show that the variable damping constraint method can improve tomographic quality because it can solve the effects of uneven ray coverage. In addition, the examples show that the tomographic result near geophones is much more reliable than other areas in the velocity model.

Hot compressive deformation behavior and microstructure evolution of HIPed FGH96 superalloy

Hot deformation behavior and microstructure evolution of hot isostatically pressed FGH96 P/M superalloy were studied using isothermal compression tests. The tests were performed on a Gleeble-1500 simulator in a temperature range of 1000-1150 °C and strain rate of 0.001-1.0 s-1, respectively. By regression analysis of the stress—strain data, the constitutive equation for FGH96 superalloy was developed in the form of hyperbolic sine function with hot activation energy of 693.21 kJ/mol. By investigating the deformation microstructure, it is found that partial and full dynamical recrystallization occurs in specimens deformed below and above 1100 °C, respectively, and dynamical recrystallization （DRX） happens more readily with decreasing strain rate and increasing deformation temperature. Finally, equations representing the kinetics of DRX and grain size evolution were established.

The stability problem of reverse time migration for viscoacoustic VTI media

In real strata anisotropy and viscosity extensively exists. They degraded waveforms in amplitude, resulting in which reducing of image resolution. To obtain high-precision imaging of deep reservoirs, we extended the separated viscous and anisotropic reverse time migration （RTM） to a stable viscoacoustic anisotropic RTM for vertical transverse isotropic （VTI） media, based on single generalized standard and linear solid （GSLS） media theory.. We used a pseudo-spectral method to develop the numerical simulation. By introducing a regularization operator to eliminate the high-frequency instability problem, we built a stable inverse propagator and achieved viscoacoustic VTI media RTM. High-resolution imaging results were obtained after correcting for the effects of anisotropy and viscosity. Synthetic tests verify the validity and accuracy of algorithm.

Exact Solutions of a Generalized Multi-Fractional Nonlinear Diffusion Equation in Radical Symmetry

This paper is devoted to investigating exact solutions of a generalized fractional nonlinear anomalousdiffusion equation in radical symmetry.The presence of external force and absorption is also considered.We firstinvestigate the nonlinear anomalous diffusion equations with one-fractional derivative and then multi-fractional ones.Inboth situations,we obtain the corresponding exact solutions,and the solutions found here can have a compact behavioror a long tailed behavior.

Existence of multi-bump solutions for coupled Schr dinger systems

The Schrodinger equation -△u＋λ2u=｜u｜2q-2u has a unique positive radial solution Uλ, which decays exponentially at infinity. Hence it is reasonable that the Schrolinger system -△u1＋u1=｜u1｜2q-1u1-εb（x）｜u2｜1｜u1｜q-1u1,-△u2＋u2=｜u2｜2q-2u2-εb（x）｜u1｜1｜u2｜q-1u2 has multiple-bump solutions which behave like Uλ in the neighborhood of some points. For u=（u1,u2）∈H1（R3）×H1（R3）, a nonlinear functional Iε（u）=I1（u1）＋I2（u2）-ε/q∫R3b（x）｜u1｜q｜u2｜qdx,is defined,where I1（u1）=1/2｜｜u1｜｜2-1/2q∫R3｜u1｜2qdx and I2（u2）=1/2｜｜u2｜｜2ω-1/2q∫R3｜u2｜2qdx. It is proved that the solutions of the system are the critical points of I,. Let Z be the smooth solution manifold of the unperturbed problem and TzZ is the tangent space. The critical point of I is rewritten as the form of z ＋ w, where w ∈ （TzZ）⊥. Using some properties of Iε, it is proved that there exists a critical point of I, close to the form which is a multi-bump solution.

Lienard Equation and Exact Solutions for Some Soliton-Producing NonlinearEquations

In this paper, we first consider exact solutions for Lienard equation with nonlinear terms of any order. Then,explicit exact bell and kink profile solitary-wave solutions for many nonlinear evolution equations are obtained by means of results of the Lienard equation and proper deductions, which transform original partial differential equations into the Lienard one. These nonlinear equations include compound KdV, compound KdV-Burgers, generalized Boussinesq,generalized KP and Ginzburg-Landau equation. Some new solitary-wave solutions are found.