Impact of climate change and human activities on the spatiotemporal dynamics of surface water area in Gansu Province, China

Understanding the dynamics of surface water area and their drivers is crucial for human survival and ecosystem stability in inland arid and semi-arid areas.This study took Gansu Province,China,a typical area with complex terrain and variable climate,as the research subject.Based on Google Earth Engine,we used Landsat data and the Open-surface Water Detection Method with Enhanced Impurity Control method to monitor the spatiotemporal dynamics of surface water area in Gansu Province from 1985 to 2022,and quantitatively analyzed the main causes of regional differences in surface water area.The findings revealed that surface water area in Gansu Province expanded by 406.88 km2 from 1985 to 2022.Seasonal surface water area exhibited significant fluctuations,while permanent surface water area showed a steady increase.Notably,terrestrial water storage exhibited a trend of first decreasing and then increasing,correlated with the dynamics of surface water area.Climate change and human activities jointly affected surface hydrological processes,with the impact of climate change being slightly higher than that of human activities.Spatially,climate change affected the'source'of surface water to a greater extent,while human activities tended to affect the'destination'of surface water.Challenges of surface water resources faced by inland arid and semi-arid areas like Gansu Province are multifaceted.Therefore,we summarized the surface hydrology patterns typical in inland arid and semi-arid areas and tailored surface water'supply-demand'balance strategies.The study not only sheds light on the dynamics of surface water area in Gansu Province,but also offers valuable insights for ecological protection and surface water resource management in inland arid and semi-arid areas facing water scarcity.

Casson Nanofluid Flow with Cattaneo-Christov Heat Flux and Chemical Reaction Past a Stretching Sheet in the Presence of Porous Medium

In the current work,inclined magnetic field,thermal radiation,and the Cattaneo-Christov heat flux are taken into account as we analyze the impact of chemical reaction on magneto-hydrodynamic Casson nanofluid flow on a stretching sheet.Modified Buongiorno’s nanofluid model has been used to model the flow governing equations.The stretching surface is embedded in a porousmedium.By using similarity transformations,the nonlinear partial differential equations are transformed into a set of dimensionless ordinary differential equations.The numerical solution of transformed dimensionless equations is achieved by applying the shooting procedure together with Rung-Kutta 4th-order method employing MATLAB.The impact of significant parameters on the velocity profile f(ζ),temperature distributionθ(ζ),concentration profileϕ(ζ),skin friction coefficient(Cf),Nusselt number(Nux)and Sherwood number(Shx)are analyzed and displayed in graphical and tabular formats.With an increase in Casson fluid 0.5<β<2,the motion of the Casson fluid decelerates whereas the temperature profile increases.As the thermal relation factor expands 0.1<γ1<0.4,the temperature reduces,and consequently thermal boundary layer shrinks.Additionally,by raising the level of thermal radiation 1<Rd<7,the temperature profile significantly improves,and an abrupt expansion has also been observed in the associated thermal boundary with raise thermal radiation strength.It was observed that higher permeability 0<K<4 hinders the acceleration of Casson fluid.Higher Brownian motion levels 0.2<Nb<0.6 correspond to lower levels of the Casson fluid concentration profile.Moreover,it is observed that chemical reaction 0.2<γ2<0.5 has an inverse relation with the concentration level of Casson fluid.The current model’s significant uses include heat energy enhancement,petroleum recovery,energy devices,food manufacturing processes,and cooling device adjustment,among others.Furthermore,present outcomes have been found in great agreementwith already publishedwork.

Additive Design and Manufacturing of Jet Engine Parts

The additive design （AD） and additive manufacturing （AM） of jet engine parts will revolutionize the traditional aerospace industry. The unique characteristics of AM, such as gradient materials and micro-structures, have opened up a new direction in jet engine design and manufacturing. Engineers have been liberated from many constraints associated with traditional methodologies and technologies. One of the most significant features of the AM process is that it can ensure the consistency of parts because it starts from point（s）, continues to line（s） and layer（s）, and ends with the competed part. Collaboration between design and manufacturing is the key to success in fields including aerodynamics, thermodynamics, structural integration, heat transfer, material development, and machining. Engineers must change the way they design a part, as they shift from the traditional method of ＂subtracting material＂ to the new method of ＂adding material＂ in order to manufacture a part. AD is not the same as designing for AM, A new method and new tools are required to assist with this new way of designing and manufacturing. This paper discusses in detail what is required in AD and AM, and how current problems can be solved.

Spatial artificial neural network model for subgrid-scale stress and heat flux of compressible turbulence

The subgrid-scale(SGS)stress and SGS heat flux are modeled by using an artificial neural network(ANN)for large eddy simulation(LES)of compressible turbulence.The input features of ANN model are based on the first-order and second-order derivatives of filtered velocity and temperature at different spatial locations.The proposed spatial artificial neural network(SANN)model gives much larger correlation coefficients and much smaller relative errors than the gradient model in an a priori analysis.In an a posteriori analysis,the SANN model performs better than the dynamic mixed model(DMM)in the prediction of spectra and statistical properties of velocity and temperature,and the instantaneous flow structures.

Numerical simulation of two-dimensional granular shearing flows and the friction force of a moving slab on the granular media

This paper studies some interesting features of two-dimensional granular shearing flow by using molecular dynamic approach for a specific granular system. The obtained results show that the probability distribution function of velocities of particles is Gaussian at the central part, but diverts from Gaussian distribution nearby the wall. The macroscopic stress along the vertical direction has large fluctuation around a constant value, the non-zero average velocity occurs mainly near the moving wall, which forms a shearing zone.. In the shearing movement, the volume of the granular material behaves in a random manner. The equivalent fl＇iction coefficient between moving slab and granular material correlates with the moving speed at low velocity, and approaches constant as the velocity is large enough.

Disturbance rejection and H_∞ pinning control of linear complex dynamical networks

This paper concerns the disturbance rejection problem of a linear complex dynamical network subject to external disturbances. A dynamical network is said to be robust to disturbance, if the H∞ norm of its transfer function matrix from the disturbance to the performance variable is satisfactorily small. It is shown that the disturbance rejection problem of a dynamical network can be solved by analysing the H∞ control problem of a set of independent systems whose dimensions are equal to that of a single node. A counter-intuitive result is that the disturbance rejection level of the whole network with a diffusive coupling will never be better than that of an isolated node. To improve this, local feedback injections are applied to a small fraction of the nodes in the network. Some criteria for possible performance improvement are derived in terms of linear matrix inequalities. It is further demonstrated via a simulation example that one can indeed improve the disturbance rejection level of the network by pinning the nodes with higher degrees than pinning those with lower degrees.

Cost and effect of pinning control for network synchronization

The problem of pinning control for the synchronization of complex dynamical networks is discussed in this paper. A cost function of the controlled network is defined by the feedback gain and the coupling strength of the network. An interesting result is that a lower cost is achieved by using the control scheme of pinning nodes with smaller degrees. Some strict mathematical analyses are presented for achieving a lower cost in the synchronization of different star-shaped networks. Numerical simulations on some non-regular complex networks generated by the Barabasi-Albert model and various star-shaped networks are performed for verification and illustration.

Thermal buckling analysis of truss-core sandwich plates

Truss-core sandwich plates have received much attention in virtue of the high values of strength-to-weight and stiffness-to-weight as well as the great ability of impulseresistance recently. It is necessary to study the stability of sandwich panels under the influence of the thermal load. However, the sandwich plates are such complex threedimensional （3D） systems that direct analytical solutions do not exist, and the finite element method （FEM） cannot represent the relationship between structural parameters and mechanical properties well. In this paper, an equivalent homogeneous continuous plate is ideMized by obtaining the effective bending and transverse shear stiffness based on the characteristics of periodically distributed unit cells. The first order shear deformation theory for plates is used to derive the stability equation. The buckling temperature of a simply supported sandwich plate is given and verified by the FEM. The effect of related parameters on mechanical properties is investigated. The geometric parameters of the unit cell are optimized to attain the maximum buckling temperature. It is shown that the optimized sandwich plate can improve the resistance to thermal buckling significantly.

Suppression of thermal postbuckling and nonlinear panel flutter motions of variable stiffness composite laminates using piezoelectric actuators

Variable stiffness composite laminates(VSCLs)are promising in aerospace engineering due to their designable material properties through changing fiber angles and stacking sequences.Aiming to control the thermal postbuckling and nonlinear panel flutter motions of VSCLs,a full-order numerical model is developed based on the linear quadratic regulator(LQR)algorithm in control theory,the classical laminate plate theory(CLPT)considering von Kármán geometrical nonlinearity,and the first-order Piston theory.The critical buckling temperature and the critical aerodynamic pressure of VSCLs are parametrically investigated.The location and shape of piezoelectric actuators for optimal control of the dynamic responses of VSCLs are determined through comparing the norms of feedback control gain(NFCG).Numerical simulations show that the temperature field has a great effect on aeroelastic tailoring of VSCLs;the curvilinear fiber path of VSCLs can significantly affect the optimal location and shape of piezoelectric actuator for flutter suppression;the unstable panel flutter and the thermal postbuckling deflection can be suppressed effectively through optimal design of piezoelectric patches.

HIGH PERFORMANCE SPARSE SOLVER FOR UNSYMMETRICAL LINEAR EQUATIONS WITH OUT-OF-CORE STRATEGIES AND ITS APPLICATION ON MESHLESS METHODS

A new direct method for solving unsymmetrical sparse linear systems（USLS） arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin （MLPG） method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.

Analysis on high-temperature oxidation and growth stress of iron-based alloy using phase field method

High-temperature oxidation is an important property to evaluate thermal protection materials. However, since oxidation is a complex process involving microstructure evolution, its quantitative analysis has always been a challenge. In this work, a phase field method （PFM） based on the thermodynamics theory is developed to simulate the oxidation behavior and oxidation induced growth stress. It involves microstructure evolution and solves the problem of quantitatively computational analysis for the oxidation behavior and growth stress. Employing this method, the diffusion process, oxidation performance, and stress evolution axe predicted for Fe-Cr-A1-Y alloys. The numerical results agree well with the experimental data. The linear relationship between the maximum growth stress and the environment oxygen concentration is found. PFM provides a powerful tool to investigate high-temperature oxidation in complex environments.

Microstructure-based three-dimensional characterization of chip formation and surface generation in the machining of particulate-reinforced metal matrix composites

Particulate-reinforced metal matrix composites(PRMMCs)are difficult to machine due to the inclusion of hard,brittle reinforcing particles.Existing experimental investigations rarely reveal the complex material removal mechanisms(MRMs)involved in the machining of PRMMCs.This paper develops a three-dimensional(3D)microstructure-based model for investigating the MRM and surface integrity of machined PRMMCs.To accurately mimic the actual microstructure of a PRMMC,polyhedrons were randomly distributed inside the matrix to represent irregular SiC particles.Particle fracture and matrix deformation and failure were taken into account.For the model’s capability comparison,a two-dimensional(2D)analysis was also conducted.Relevant cutting experiments showed that the established 3D model accurately predicted the material removal,chip morphology,machined surface finish,and cutting forces.It was found that the matrix-particle-tool interactions led to particle fractures,mainly in the primary shear and secondary deformation zones along the cutting path and beneath the machined surface.Particle fracture and dilodegment greatly influences the quality of a machined surface.It was also found that although a 2D model can reflect certain material removal features,its ability to predict microstructural variation is limited.

Effects of Prandtl number in two-dimensional turbulent convection

We report a numerical study of the Prandtl-number(Pr)effects in two-dimensional turbulent Rayleigh-Bénard convection.The simulations were conducted in a square box over the Pr range from 0.25 to 100 and over the Rayleigh number(Ra)range from 10^(7) to 10^(10).We find that both the strength and the stability of the large-scale flow decrease with the increasing of Pr,and the flow pattern becomes plume-dominated at high Pr.The evolution in flow pattern is quantified by the Reynolds number(Re),with the Ra and the Pr scaling exponents varying from 0.54 to 0.67 and-0.87 to-0.93,respectively.It is further found that the non-dimensional heat flux at small Ra diverges strongly for different Pr,but their difference becomes marginal as Ra increases.For the thermal boundary layer,the spatially averaged thicknesses for all the Pr numbers can be described byδθ~Ra^(-0.30) approximately,but the local values vary a lot for different Pr,which become more uniform with Pr increasing.

Effect of turbulence intensity on airfoil flow:numerical simulations and experimental measurements

The effect of the turbulence intensity of the oncoming stream on the aerodynamic characteristics of the NACA-0012 airfoil is investigated by a direct numerical simulation. The numerical results are found to be consistent with the experimental measurements. Based on the finite spectral QUICK scheme, the simulation gets the high accuracy results. Both the simulation and the experiment reveal that the airfoil stall does not exist for the low turbulence intensity, however, occurs when the turbulence intensity increases sufficiently. Besides, the turbulence intensity has a significant effect on both the airfoil boundary layer and the separated shear layer.

Preface: symposium on turbulence structures and aerodynamic heat/force(STSAHF2018)——scientific significance of turbulence research

Most fluid flows in nature and engineering applications are in the state of turbulence.Turbulent motions usually exhibit a wide range of spatial and temporal scales,such as the flow of natural gas and oil in pipelines,the wakes of cars and submarines,the boundary layer of an aircraft,the current in the ocean surface,the atmospheric boundary layer,the interstellar gas clouds(gaseous stars),and the Earth’s wake in the solar wind.Turbulence can greatly improve the heat and mass transfer efficiency of macroscopic flow.For example,chemical engineers use turbulence to mix up and homogenize fluid components and to increase chemical reaction rates in liquids or gases.However,turbulence can also lead to increases in drag,aerodynamic heat,and hydrodynamic and aerodynamic noise.For instance,the aerodynamic loading of high-speed aircraft can be significantly increased due to turbulence.

Explicit formulations and performance of LSFD method on Cartesian mesh

Performance of the LSFD method is compared with conventional FD schemes. Generally, 9-point stencils for 2D cases and 27-point stencils for 3D cases are used for the approximation of the first and second order derivatives obtained with conventional central difference schemes. When the same stencils are used, explicit LSFD formulations for approximation of the first and second order derivatives are presented. The LSFD formulations are actually a combination of conventional central difference schemes along relevant mesh lines. It has been found that LSFD formulations need much less iteration steps than the conventional FD schemes to converge, and the ratio of mesh spacing in the x and y directions is an important parameter in the LSFD application, with a great impact on stability of LSFD computation.

Lattice Boltzmann simulation of flows in bifurcate channel at rotating inflow boundary conditions and resulted different outflow fluxes

The Lattice Boltzmann method （LBM） is used to simulate the flow field in a bifurcate channel which is a simplified model of the draft tube of hydraulic turbine machine. According to the simulation results, some qualitative conclusions can be deduced. The reason of uneven flux in different branches of draft tube is given. Not only the vortex rope itself, but also the attenuation of the rotation strength is important in bringing on the uneven flux. The later leads to adverse pressure gradient, and changes the velocity profile. If the outlet contains more than one exit, the one that contains the vortex rope will lose flux because of this adverse pressure gradient. Several possible methods can be used to minimize the adverse pressure gradient domain in order to improve the efficiency of turbine machine.

Sandwich structures with tapered tubes as core:A quasi-static investigation

In this article,the experimental and finite element analysis is utilized to investigate the quasi-static compression features of sandwich constructions built with tapered tubes.3D printing technology was utilized to create the hollow centers of the tapering tubes,with and without corrugations.The results demonstrate that the energy absorption(EA)and specific energy absorption(SEA)of the single corrugated tapered tube sandwich are 51.6% and 19.8% higher,respectively,than those of the conical tube sandwich.Furthermore,the results demonstrate that energy absorbers can benefit from corrugation in order to increase their efficiency.Additionally,the tapered corrugated tubes'resistance to oblique impacts was studied.Compared to a straight tube,the tapered tube is more resistant to oblique loads and has a lower initial peak crushing force(PCF),according to numerical simulations.After conducting a parametric study,it was discovered that the energy absorption performance of the sandwich construction is significantly affected by the amplitude,number of corrugations,and wall thickness.EA and SEA of DTS with corrugation number of 8 increased by 17.4%and 29.6%,respectively,while PCF decreased by 9.2% compared to DTS with corrugation number of 10.

Bioconvective Hybrid Flow with Microorganisms Migration and Buongiorno’s Model under Convective Condition

Heat transfer improves significantly when the working fluid has high thermal conductivity.Heat transfer can be found in fields such as food processing,solar through collectors,and drug delivery.Considering this notable fact,this work is focused on investigating the bio-convection-enhanced heat transfer in the existence of convective boundary conditions in the flow of hybrid nanofluid across a stretching surface.Buongiorno fluid model with hybrid nanoparticles has been employed along the swimming microorganisms to investigate the mixture base working fluid.The developed nonlinear flow governing equations have been tackled numerically with the help of the bvp4c.The effects of relevant parameters on the flowdynamic have been portrayed in a graphical representation.The velocity profile decreases by raising the levels of buoyancy ratio and mixed convection in the range of 0.1<λ≤0.3.It has been discovered thatwhen bioconvection levels rise,motile microbemigration abruptly slows,which results in a decrease in fluid acceleration.The concentration of fluid flow declined for the Lewis number,but the opposite trend has been observed for the elastic parameter,thermophoresis parameter,and buoyancy ratio.With rising values of Brownian motion and thermophoretic diffusion,the surface drag and Nusselt number decrease significantly.Whereas,the opposite trend has been observed when the values of the thermal Biot number,Prandtl number and buoyancy ratio are enhanced.Additionally,data from this study have been validated by comparison with those that have previously been published,and an appropriate rate of agreement has been observed.

TIME ASYMPTOTIC BEHAVIOR OF THE BIPOLAR NAVIER-STOKES-POISSON SYSTEM

The bipolar Navier-Stokes-Poisson system （BNSP） has been used to simulate the transport of charged particles （ions and electrons for instance） under the influence of electrostatic force governed by the self-consistent Poisson equation. The optimal L^2 time convergence rate for the global classical solution is obtained for a small initial perturbation of the constant equilibrium state. It is shown that due to the electric field, the difference of the charge densities tend to the equilibrium states at the optimal rate （1 ＋ t）^-3/4 in L^2-norm, while the individual momentum of the charged particles converges at the optimal rate （1 ＋ t）^-1/4 which is slower than the rate （1 ＋ t）^-3/4 for the compressible Navier-Stokes equations （NS）. In addition, a new phenomenon on the charge transport is observed regarding the interplay between the two carriers that almost counteracts the influence of the electric field so that the total density and momentum of the two carriers converges at a faster rate （1 ＋ t）^-3/4＋ε for any small constant ε 〉 0. The above estimates reveal the essential difference between the unipolar and the bipolar Navier-Stokes-Poisson systems.