Numerical simulation of three-dimensional fracturing fracture propagation in radial wells

A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.

Molecular dynamics study of thermal stress and heat propagation in tungsten under thermal shock

Using molecular dynamics （MD） simulation, we study the thermal shock behavior of tungsten （W）, which has been used for the plasma facing material （PFM） of tokamaks. The thermo-elastic stress wave, corresponding to the collective displacement of atoms, is analyzed with the Lagrangian atomic stress method, of which the reliability is also analyzed. The stress wave velocity corresponds to the speed of sound in the material, which is not dependent on the thermal shock energy. The peak pressure of a normal stress wave increases with the increase of thermal shock energy. We analyze the temperature evolution of the thermal shock region according to the Fourier transformation. It can be seen that the “obvious” velocity of heat propagation is less than the velocity of the stress wave; further, that the thermo-elastic stress wave may contribute little to the transport of kinetic energy. The heat propagation can be described properly by the heat conduction equation. These results may be useful for understanding the process of the thermal shock of tungsten.

Parametric Instabilities of Parallel Propagating Circularly Polarized Alfven Waves:One-Dimensional Hybrid Simulations

By performing one-dimensional （l-D） hybrid simulations, we analyze in detail the parametric instabilities of the Alfv^n waves with a spectrum in a low beta plasma. The parametric instabilities experience two stages. In the first stage, the density modes are excited and immediately couple with the pump Alfv6n waves. In the second stage, each pump Alfv^n wave decays into a density mode and a daughter Alfv6n mode similar to the monochromatic cases. Ftlrthermore, the proton velocity beam will also be formed after the saturation of the parametric instabilities. When the plasma beta is high, the parametric decay in the second stage will be strongly suppressed.

Numerical simulation of nonlinear propagation of sound waves in a finite horn

Nonlinear acoustic propagation generated by a piston in a finite horn is numerically studied. A quasi-one-dimensional nonlinear model with varying cross-section uses high-order low-dispersion numerical schemes to solve the governing equation. Because of the nonlinear wave distortion and reflected sound waves at the mouth, broadband time-domain impedance boundary conditions are employed. The impedance approximation can be optimized to identify the complex-conjugate pole-residue pairs of the impedance functions, which can be calculated by fast and efficient recursive convolution. The numerical results agree very well with experi- mental data in the situations of weak nonlinear wave propagation in an exponential horn, it is shown that the model can describe the broadband characteristics caused by nonlinear distortion. Moreover, finite-amplitude acoustic propagation in types of horns is simulated, including hyperbolic, conical, exponential and sinusoidal horns. It is found that sound pressure levels at the horn mouth are strongly affected by the horn profiles, the driving velocity and frequency of the piston. The paper also discusses the influence of the horn geometry on nonlinear waveform distortion.

Numerical Investigation of a Refractive Index SPR D-Type Optical Fiber Sensor Using COMSOL Multiphysics

Recently, many programs have been developed for simulation or analysis of the different parameters of light propagation in optical fibers, either for sensing or for communication purposes. In this paper, it is shown the COMSOL Multiphysics as a fairly robust and simple program, due to the existence of a graphical environment, to perform simulations with good accuracy. Results are compared with other simulation analysis, focusing on the surface plasmon resonance （SPR） phenomena for refractive index sensing in a D-type optical fiber, where the characteristics of the material layers, in terms of the type and thickness, and the residual fiber cladding thickness are optimized.

A parameterization scheme for upward lightning in the cloud model and a discussion of the initial favorable environmental characteristics in the cloud

The upward lightning(UL) initiated from the top of tall buildings(at least above 100 m) is a type of atmospheric discharge. Currently, we understand the nature of the UL from ground observations, but the corresponding theoretical research is lacking. Based on an existing bidirectional leader stochastic model, a stochastic parameterization scheme for the UL has been built and embedded in an existing two-dimensional thundercloud charge/discharge model. The ULs simulated from the experiments with two-dimensional high resolution agree generally with the observation results. By analyzing the charge structure of thunderstorm clouds, we determined the in-cloud environmental characteristics that favor the initiation of conventional cloud-to-ground(CG) flashes and analyzed the differences and similarities of some characteristics of the positive and the negative UL. Simulation results indicate that the positive ULs are typically other-lightning-triggered ULs(OLTUL) and are usually a discharge phenomenon between the ground and the lower positive charge region appearing below the main middle negative charge region. The effect of the previous in-cloud lightning(IC) process of space electrical field provides favorable conditions for the initiation of a positive UL. Its entire discharge process is limited, and the branches of the leader are fewer in number as its discharge is not sufficient. A negative UL is generally a discharge phenomenon of the dipole charge structure between the ground and the main negative charge region. The lower temperature stratification and the sinking of the hydrometeors typically initiate a negative UL. Negative ULs develop strongly and have more branches. The OLTUL is initiated mainly during the development stage of a thunderstorm, while the self-triggered UL(STUL) is initiated mainly during the dissipation stage of a thunderstorm.