Numerical Solutions of the Classical and Modified Buckley-Leverett Equations Applied to Two-Phase Fluid Flow

Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.

GLOBAL SOLUTIONS TO 1D COMPRESSIBLE NAVIER-STOKES/ALLEN-CAHN SYSTEM WITH DENSITY-DEPENDENT VISCOSITY AND FREE-BOUNDARY

This paper is concerned with the Navier-Stokes/Allen-Cahn system,which is used to model the dynamics of immiscible two-phase flows.We consider a 1D free boundary problem and assume that the viscosity coefficient depends on the density in the form ofη(ρ)=ρ^(α).The existence of unique global H^(2m)-solutions(m∈N)to the free boundary problem is proven for when 0<α<1/4.Furthermore,we obtain the global C^(∞)-solutions if the initial data is smooth.

Design of low-alloying and high-performance solid solution-strengthened copper alloys with element substitution for sustainable development

Solid solution-strengthened copper alloys have the advantages of a simple composition and manufacturing process,high mechanical and electrical comprehensive performances,and low cost;thus,they are widely used in high-speed rail contact wires,electronic component connectors,and other devices.Overcoming the contradiction between low alloying and high performance is an important challenge in the development of solid solution-strengthened copper alloys.Taking the typical solid solution-strengthened alloy Cu-4Zn-1Sn as the research object,we proposed using the element In to replace Zn and Sn to achieve low alloying in this work.Two new alloys,Cu-1.5Zn-1Sn-0.4In and Cu-1.5Zn-0.9Sn-0.6In,were designed and prepared.The total weight percentage content of alloying elements decreased by 43%and 41%,respectively,while the product of ultimate tensile strength(UTS)and electrical conductivity(EC)of the annealed state increased by 14%and 15%.After cold rolling with a 90%reduction,the UTS of the two new alloys reached 576 and 627MPa,respectively,the EC was 44.9%IACS and 42.0%IACS,and the product of UTS and EC(UTS×EC)was 97%and 99%higher than that of the annealed state alloy.The dislocations proliferated greatly in cold-rolled alloys,and the strengthening effects of dislocations reached 332 and 356 MPa,respectively,which is the main reason for the considerable improvement in mechanical properties.

Localized wave solutions and interactions of the (2+1)-dimensional Hirota-Satsuma-Ito equation

This paper studies the(2+1)-dimensional Hirota-Satsuma-Ito equation.Based on an associated Hirota bilinear form,lump-type solution,two types of interaction solutions,and breather wave solution of the(2+1)-dimensional Hirota-Satsuma-Ito equation are obtained,which are all related to the seed solution of the equation.It is interesting that the rogue wave is aroused by the interaction between one-lump soliton and a pair of resonance stripe solitons,and the fusion and fission phenomena are also found in the interaction between lump solitons and one-stripe soliton.Furthermore,the breather wave solution is also obtained by reducing the two-soliton solutions.The trajectory and period of the one-order breather wave are analyzed.The corresponding dynamical characteristics are demonstrated by the graphs.

Evaluation of frictional pressure drop correlations for air-water and air-oil two-phase flow in pipeline-riser system

Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.

PDE Standardization Analysis and Solution of TypicalMechanics Problems

A numerical approach is an effective means of solving boundary value problems(BVPs).This study focuses on physical problems with general partial differential equations(PDEs).It investigates the solution approach through the standard forms of the PDE module in COMSOL.Two typical mechanics problems are exemplified:The deflection of a thin plate,which can be addressed with the dedicated finite element module,and the stress of a pure bending beamthat cannot be tackled.The procedure for the two problems regarding the three standard forms required by the PDE module is detailed.The results were in good agreement with the literature,indicating that the PDE module provides a promising means to solve complex PDEs,especially for those a dedicated finite element module has yet to be developed.

THE LIMITING PROFILE OF SOLUTIONS FOR SEMILINEAR ELLIPTIC SYSTEMS WITH A SHRINKING SELF-FOCUSING CORE

In this paper,we consider the semilinear elliptic equation systems{△u+u=αQ_(n)(x)|u|^(α-2)|v|^(β)u in R^(N),-△v+v=βQ(x)|u|^(α)|v|^(β-2)v in R^(N),where N≥3,α,β>1,α+β<2^(*),2^(*)=2N/N-2 and Q_(n) are bounded given functions whose self-focusing cores{x∈R^(N)|Q_(n)(x)>0} shrink to a set with finitely many points as n→∞.Motivated by the work of Fang and Wang[13],we use variational methods to study the limiting profile of ground state solutions which are concentrated at one point of the set with finitely many points,and we build the localized concentrated bound state solutions for the above equation systems.

Solution approximations for a mathematical model of relativistic electrons with beta derivative

The main aim of this paper is to obtain the exact and semi-analytical solutions of the nonlinear Klein-Fock-Gordon(KFG)equation which is a model of relativistic electrons arising in the laser thermonuclear fusion with beta derivative.For this purpose,both the modified extended tanh-function(mETF)method and the homotopy analysis method(HAM)are used.While applying the mETF the chain rule for beta derivative and complex wave transform are used for obtaining the exact solution.The advantage of this procedure is that discretization or normalization is not required.By applying the mETF,the exact solutions are obtained.Also,by applying the HAM semi-analytical results for the considered equation are acquired.In HAM?curve gives us a chance to find the suitable value of the for the convergence of the solution series.Also,comparative graphical representations are given to show the effectiveness,reliability of the methods.The results show that the m ETF and HAM are reliable and applicable tools for obtaining the solutions of non-linear fractional partial differential equations that involve beta derivative.This study can bring a new perspective for studies on fractional differential equations.On the other hand,it can be said that scientists can apply the considered methods for different mathematical models arising in physics,chemistry,engineering,social sciences and etc.which involves fractional differentiation.Briefly the results may cause a new insight who studies on relativistic electron modelling.

Bifurcations, Analytical and Non-Analytical Traveling Wave Solutions of (2 + 1)-Dimensional Nonlinear Dispersive Boussinesq Equation

For the (2 + 1)-dimensional nonlinear dispersive Boussinesq equation, by using the bifurcation theory of planar dynamical systems to study its corresponding traveling wave system, the bifurcations and phase portraits of the regular system are obtained. Under different parametric conditions, various sufficient conditions to guarantee the existence of analytical and non-analytical solutions of the singular system are given by using singular traveling wave theory. For certain special cases, some explicit and exact parametric representations of traveling wave solutions are derived such as analytical periodic waves and non-analytical periodic cusp waves. Further, two-dimensional wave plots of analytical periodic solutions and non-analytical periodic cusp wave solutions are drawn to visualize the dynamics of the equation.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Analytic Solutions for Hindered Rotation Quantum Problem

Normalizable analytic solutions of the quantum rotor problem with divergent potential are presented here as solution of the Schrödinger equation. These solutions, unknown to the literature, represent a mathematical advance in the description of physical phenomena described by the second derivative operator associated with a divergent interaction potential and, being analytical, guarantee the optimal interpretation of such phenomena.

Numerical Simulations of Wave Impact Forces on the Open-Type Sea Access Road Using A Two-Phase SPH Model

A numerical study based on a two-dimensional two-phase SPH(Smoothed Particle Hydrodynamics)model to analyze the action of water waves on open-type sea access roads is presented.The study is a continuation of the analyses presented by Chen et al.(2022),in which the sea access roads are semi-immersed.In this new configuration,the sea access roads are placed above the still water level,therefore the presence of the air phase becomes a relevant issue in the determination of the wave forces acting on the structures.Indeed,the comparison of wave forces on the open-type sea access roads obtained from the single and two-phase SPH models with the experimental results shows that the latter are in much better agreement.So in the numerical simulations,a two-phaseδ-SPH model is adopted to investigate the dynamical problems.Based on the numerical results,the maximum horizontal and uplifting wave forces acting on the sea access roads are analyzed by considering different wave conditions and geometries of the structures.In particular,the presence of the girder is analyzed and the differences in the wave forces due to the air cushion effects which are created below the structure are highlighted.

Existence of Entropy Solution for Degenerate Parabolic-Hyperbolic Problem Involving p(x)-Laplacian with Neumann Boundary Condition

We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.

Malargüe Site Remediation:A Successful Solution for Uranium Mill Tailings

Comisión Nacional de Energía Atómica (CNEA) has the responsibility for restoring uranium mining facilities once the operations have finished.CNEA,within its Environmental Program and in compliance with its legal responsibilities,decides to implement a restoration project for all sites related to the mining and processing of uranium ores.The Malargüe Site is located within the Province of Mendoza in the city of Malargüe.It is the first site to successfully complete its remediation.The activities consist of relocation of tailings to an engineering repository.The tailings management(encapsulation) and rehabilitation of the area was finished in June 2017.The remediation alternative for the ore tailings was selected after conducting comparative studies and submitted the project to the society for consideration.The objective of the encapsulation of the mineral tails is to isolate them from the environment,also proceeding with the decontamination and rehabilitation of the area (landscaping,post-closure monitoring and 20 years monitoring period).Encapsulation consisted of the construction of a containment cell for the mine tailings,to isolate them and prevent pollutants from entering the environment through the transfer routes.To clean the impacted areas,the soil was removed,it was incorporated into the encapsulation,and the filling was carried out with natural soils from the area.Remediation prevents radon transfer to the environment,as ^(222)Ra is an alpha emitter with a half-life of four days,which produces its own radioactive progeny.Radon progeny are solids,and when a ^(222)Ra nucleus emits an alpha particle into the air,the resulting ^(218)Po nucleus,momentarily electrically charged,adheres to any dust particle.Remediation prevents the discharge into the air containing radon and also containing dust particles charged with intensely radioactive radon progeny.The tasks mentioned make it possible to decrease radon emanation,reduce radiological risks to the public and prevent the entry of rainwater into the system.In addition,the containment system prevents the discharge of contaminated liquids into the environment,avoiding contamination of the groundwater.All these activities are according to the concepts of sustainability.

GLOBAL WEAK SOLUTIONS FOR AN ATTRACTION-REPULSION CHEMOTAXIS SYSTEM WITH p-LAPLACIAN DIFFUSION AND LOGISTIC SOURCE

This paper is concerned with the following attraction-repulsion chemotaxis system with p-Laplacian diffusion and logistic source:■The system here is under a homogenous Neumann boundary condition in a bounded domainΩ ■ R^(n)(n≥2),with χ,ξ,α,β,γ,δ,k_(1),k_(2)> 0,p> 2.In addition,the function f is smooth and satisfies that f(s)≤κ-μs~l for all s≥0,with κ ∈ R,μ> 0,l> 1.It is shown that(ⅰ)if l> max{2k_(1),(2k_(1)n)/(2+n)+1/(p-1)},then system possesses a global bounded weak solution and(ⅱ)if k_(2)> max{2k_(1)-1,(2k_(1)n)/(2+n)+(2-p)/(p-1)} with l> 2,then system possesses a global bounded weak solution.

Dancing bubble sonoluminescence in phosphoric acid solution

Sonoluminescence is more distinctly observed in phosphoric and sulfuric acid,which exhibit high viscosity and lower vapor pressures relative to water.Within an 85-wt%phosphoric acid solution saturated with argon(Ar),variations in the light-emitting regimes of bubbles were noted to correspond with increments in the driving acoustic intensity.Specifically,the bubbles were observed to perform a dance-like motion 2 cm below the multi-bubble sonoluminescence(MBSL)cluster,traversing a 25-mm^(2) grid during the camera exposure period.Spectral analysis conducted at the beginning of the experiment showed a gradual attenuation of CN(B^(2)S–X^(2)S)emission concurrent with a strengthening of Ar(4p–4s)atom emission lines.The application of a theoretical temperature model to the spectral data revealed that the internal temperature of the bubbles escalates swiftly upon their implosion.This study is instrumental in advancing the comprehension of the underlying mechanisms of sonoluminescence and in the formulation of a dynamic model for the behavior of the bubbles.

Hamiltonian system for the inhomogeneous plane elasticity of dodecagonal quasicrystal plates and its analytical solutions

A Hamiltonian system is derived for the plane elasticity problem of two-dimensional dodecagonal quasicrystals by introducing the simple state function. By using symplectic elasticity approach, the analytic solutions of the phonon and phason displacements are obtained further for the quasicrystal plates. In addition, the effectiveness of the approach is verified by comparison with the data of the finite integral transformation method.

On entire solutions of some Fermat type differential-difference equations

On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].

Simulation of Two-Phase Flowback Phenomena in Shale Gas Wells

The gas-water two-phaseflow occurring as a result of fracturingfluidflowback phenomena is known to impact significantly the productivity of shale gas well.In this work,this two-phaseflow has been simulated in the framework of a hybrid approach partially relying on the embedded discrete fracture model(EDFM).This model assumes the region outside the stimulated reservoir volume(SRV)as a single-medium while the SRV region itself is described using a double-medium strategy which can account for thefluid exchange between the matrix and the micro-fractures.The shale gas adsorption,desorption,diffusion,gas slippage effect,fracture stress sensitivity,and capillary imbibition have been considered.The shale gas production,pore pressure distribution and water saturation distribution in the reservoir have been simulated.The influences of hydraulic fracture geometry and nonorthogonal hydraulic fractures on gas production have been determined and discussed accordingly.The simulation results show that the daily gas production has an upward and downward trend due to the presence of a large amount of fracturingfluid in the reservoir around the hydraulic fracture.The smaller the angle between the hydraulic fracture and the wellbore,the faster the daily production of shale gas wells decreases,and the lower the cumulative production.Nonplanar fractures can increase the control volume of hydraulic fractures and improve the production of shale gas wells.