Deep learning(DL)has become a crucial technique for predicting the El Niño-Southern Oscillation(ENSO)and evaluating its predictability.While various DL-based models have been developed for ENSO predictions,many f...Deep learning(DL)has become a crucial technique for predicting the El Niño-Southern Oscillation(ENSO)and evaluating its predictability.While various DL-based models have been developed for ENSO predictions,many fail to capture the coherent multivariate evolution within the coupled ocean-atmosphere system of the tropical Pacific.To address this three-dimensional(3D)limitation and represent ENSO-related ocean-atmosphere interactions more accurately,a novel this 3D multivariate prediction model was proposed based on a Transformer architecture,which incorporates a spatiotemporal self-attention mechanism.This model,named 3D-Geoformer,offers several advantages,enabling accurate ENSO predictions up to one and a half years in advance.Furthermore,an integrated gradient method was introduced into the model to identify the sources of predictability for sea surface temperature(SST)variability in the eastern equatorial Pacific.Results reveal that the 3D-Geoformer effectively captures ENSO-related precursors during the evolution of ENSO events,particularly the thermocline feedback processes and ocean temperature anomaly pathways on and off the equator.By extending DL-based ENSO predictions from one-dimensional Niño time series to 3D multivariate fields,the 3D-Geoformer represents a significant advancement in ENSO prediction.This study provides details in the model formulation,analysis procedures,sensitivity experiments,and illustrative examples,offering practical guidance for the application of the model in ENSO research.展开更多
In this article,we develop the Laplace transform(LT)based Chebyshev spectral collocation method(CSCM)to approximate the time fractional advection-diffusion equation,incorporating the Atangana-Baleanu Caputo(ABC)deriva...In this article,we develop the Laplace transform(LT)based Chebyshev spectral collocation method(CSCM)to approximate the time fractional advection-diffusion equation,incorporating the Atangana-Baleanu Caputo(ABC)derivative.The advection-diffusion equation,which governs the transport of mass,heat,or energy through combined advection and diffusion processes,is central to modeling physical systems with nonlocal behavior.Our numerical scheme employs the LT to transform the time-dependent time-fractional PDEs into a time-independent PDE in LT domain,eliminating the need for classical time-stepping methods that often suffer from stability constraints.For spatial discretization,we employ the CSCM,where the solution is approximated using Lagrange interpolation polynomial based on the Chebyshev collocation nodes,achieving exponential convergence that outperforms the algebraic convergence rates of finite difference and finite element methods.Finally,the solution is reverted to the time domain using contour integration technique.We also establish the existence and uniqueness of the solution for the proposed problem.The performance,efficiency,and accuracy of the proposed method are validated through various fractional advection-diffusion problems.The computed results demonstrate that the proposed method has less computational cost and is highly accurate.展开更多
This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the acc...This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the accuracy mismatch between tem-poral low-order finite difference and spatial high-order discre tization,the ir time collocation points must increase dramatically to solve highly oscillatory solutions of structural vibration,which results in a surge in computing time and a decrease in accuracy.To address this problem,we introduced the step-by-step idea in the space-time spectral method.The Chebyshev polynomials and Lagrange's equation were applied to derive discrete spatial goverming equations,and a matrix projection method was used to map the calculation results of prev ious steps as the initial conditions of the subsequent steps.A series of numerical experiments were carried out.The results of the proposed method were compared with those obtained by traditional space-time spectral methods,which showed that higher accuracy could be achieved in a shorter computation time than the latter in highly oscillatory cases.展开更多
In this paper, we proposed a spectral gradient-Newton two phase method for constrained semismooth equations. In the first stage, we use the spectral projected gradient to obtain the global convergence of the algorithm...In this paper, we proposed a spectral gradient-Newton two phase method for constrained semismooth equations. In the first stage, we use the spectral projected gradient to obtain the global convergence of the algorithm, and then use the final point in the first stage as a new initial point to turn to a projected semismooth asymptotically newton method for fast convergence.展开更多
Spectral conjugate gradient method is an algorithm obtained by combination of spectral gradient method and conjugate gradient method,which is characterized with global convergence and simplicity of spectral gradient m...Spectral conjugate gradient method is an algorithm obtained by combination of spectral gradient method and conjugate gradient method,which is characterized with global convergence and simplicity of spectral gradient method,and small storage of conjugate gradient method.Besides,the spectral conjugate gradient method was proved that the search direction at each iteration is a descent direction of objective function even without relying on any line search method.Spectral conjugate gradient method is applied to full waveform inversion for numerical tests on Marmousi model.The authors give a comparison on numerical results obtained by steepest descent method,conjugate gradient method and spectral conjugate gradient method,which shows that the spectral conjugate gradient method is superior to the other two methods.展开更多
As a generalization of the two-term conjugate gradient method(CGM),the spectral CGM is one of the effective methods for solving unconstrained optimization.In this paper,we enhance the JJSL conjugate parameter,initiall...As a generalization of the two-term conjugate gradient method(CGM),the spectral CGM is one of the effective methods for solving unconstrained optimization.In this paper,we enhance the JJSL conjugate parameter,initially proposed by Jiang et al.(Computational and Applied Mathematics,2021,40:174),through the utilization of a convex combination technique.And this improvement allows for an adaptive search direction by integrating a newly constructed spectral gradient-type restart strategy.Then,we develop a new spectral CGM by employing an inexact line search to determine the step size.With the application of the weak Wolfe line search,we establish the sufficient descent property of the proposed search direction.Moreover,under general assumptions,including the employment of the strong Wolfe line search for step size calculation,we demonstrate the global convergence of our new algorithm.Finally,the given unconstrained optimization test results show that the new algorithm is effective.展开更多
This study employed tri-component continuous monitoring data from 10 measurement points on both sides of a base isolation layer in the basement of a large-span high-rise building in Beijing,as well as from a free-fiel...This study employed tri-component continuous monitoring data from 10 measurement points on both sides of a base isolation layer in the basement of a large-span high-rise building in Beijing,as well as from a free-field station and roof frame,during a Mw 5.5 magnitude earthquake in Pingyuan,Shandong,in 2023.The H/V spectral ratio method was used to evaluate the structural dynamic response characteristics of the building and analyze the regulatory effect of the base-isolation layer on seismic waves.The results indicate that during the earthquake,the peak frequency of the free-field and the measurement points below the base-isolation layer was stable at 0.17 Hz,whereas the main frequency of the measurement points above the base-isolation layer increased to 0.75–1.18 Hz,which is 4–6 times greater than that of the points below.The amplitude was suppressed by more than 70%,confirming that the base isolation layer effectively isolated the low-frequency energy from the ground and increased the response frequency of the building.When the building was excited by an earthquake,a three-tier frequency gradient was formed throughout the building:“base-isolation layer(0.17 Hz)-main body(1.18 Hz)-roof frame(3.83 Hz)”,which can effectively avoid resonance of the entire building.In addition,the composite base-isolation device changed the dynamic characteristics of the structure.The resonance period was extended from 0.74 s(theoretical value without base isolation)to 5.9 s(calculated value),and the resonance frequency was reduced from 1.35 to 0.17 Hz.This finding indicates that the base-isolation layer can enhance seismic performance by increasing flexibility and damping.展开更多
This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SE...This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.展开更多
Radiative cooling textiles with spectrally selective surfaces offer a promising energy-efficient approach for sub-ambient cooling of outdoor objects and individuals.However,the spectrally selective mid-infrared emissi...Radiative cooling textiles with spectrally selective surfaces offer a promising energy-efficient approach for sub-ambient cooling of outdoor objects and individuals.However,the spectrally selective mid-infrared emission of these textiles significantly hinders their efficient radiative heat exchange with self-heated objects,thereby posing a significant challenge to their versatile cooling applicability.Herein,we present a bicomponent blow spinning strategy for the production of scalable,ultra-flexible,and healable textiles featuring a tailored dual gradient in both chemical composition and fiber diameter.The gradient in the fiber diameter of this textile introduces a hierarchically porous structure across the sunlight incident area,thereby achieving a competitive solar reflectivity of 98.7%on its outer surface.Additionally,the gradient in the chemical composition of this textile contributes to the formation of Janus infrared-absorbing surfaces:The outer surface demonstrates a high mid-infrared emission,whereas the inner surface shows a broad infrared absorptivity,facilitating radiative heat exchange with underlying self-heated objects.Consequently,this textile demonstrates multi-scenario radiative cooling capabilities,enabling versatile outdoor cooling for unheated objects by 7.8℃ and self-heated objects by 13.6℃,compared to commercial sunshade fabrics.展开更多
A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue prob...A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).展开更多
Fluid dynamic research on rectangular and trapezoidal fins is aimed at increasing heat transfer by means of large surfaces.The trapezoidal cavity form is compared with its thermal and flow performance,and it is reveal...Fluid dynamic research on rectangular and trapezoidal fins is aimed at increasing heat transfer by means of large surfaces.The trapezoidal cavity form is compared with its thermal and flow performance,and it is revealed that trapezoidal fins tend to be more efficient,particularly when material optimization is critical.Motivated by the increasing need for sustainable energy management,this work analyses the thermal performance of inclined trapezoidal and rectangular porous fins utilising a unique hybrid nanofluid.The effectiveness of nanoparticles in a working fluid is primarily determined by their thermophysical properties;hence,optimising these properties can significantly improve overall performance.This study considers the dispersion of Graphene Oxide(GO)and Molybdenum Disulfide in the base fluid,engine oil.Temperature profiles are analysed by altering the radiative,porosity,wet porous,and angle of inclination parameters.Surface and contour plots are constructed by using the Lobatto IIIa Collocation Method with BVP5C solver in MATLAB and Gradient Descent Optimisation to predict the combined heat transfer rate.According to the study,fluid temperature consistently decreases when the angle of inclination,wet porous parameter,porosity parameter,and radiative parameter increase,suggesting significantly improved heat dissipation.The trapezoidal fin consistently exhibits a superior heat transfer mechanism than a rectangular fin.It is found that the trapezoidal fin transmits heat at a rate that is 0.05%higher than that of the rectangular fin.Validation of the present study is done through the comparison of previous studies.This research provides useful design insights for sophisticated engineering uses,including electrical cooling devices,heat exchangers,radiators,and solar heaters.展开更多
We analyzed the infrared 0R）-near infrared （NIR） 2D correlation spectra of drugs perturbed by temperature. By identification of functional groups by IR spectrum and by the correlation analysis of IR-NIR spectrum, w...We analyzed the infrared 0R）-near infrared （NIR） 2D correlation spectra of drugs perturbed by temperature. By identification of functional groups by IR spectrum and by the correlation analysis of IR-NIR spectrum, we identified the characteristic spectral bands that were closely related to the structure of a drug substance of interest. These characteristic spectral bands were relatively less interfered by other ingredients for analysis by the NIR correlation coefficient method. With these characteristic spectral bands, the accuracy of screening illegally added Sildenafil citrate, Tadalafil and Metforrnin hydrochloride in Chinese patent drugs and healthcare products reached about 90%, which met the requirements of rapid screening.展开更多
The conventional gravity gradient method to plot the geologic body location is fuzzy. When the depth is large and the geologic body is small, the Vzz and Vzx derivative errors are also large. We describe that using th...The conventional gravity gradient method to plot the geologic body location is fuzzy. When the depth is large and the geologic body is small, the Vzz and Vzx derivative errors are also large. We describe that using the status distinguishing factor to optimally determine the comer location is more accurate than the conventional higher-order derivative method. Thus, a better small geologic body and fault resolution is obtained by using the gravity gradient method and trial theoretical model calculation. The actual data is better processed, providing a better basis for prospecting and determination of subsurface geologic structure.展开更多
The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical ...The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.展开更多
In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident bou...In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.展开更多
In this paper, we present a new hybrid conjugate gradient algorithm for unconstrained optimization. This method is a convex combination of Liu-Storey conjugate gradient method and Fletcher-Reeves conjugate gradient me...In this paper, we present a new hybrid conjugate gradient algorithm for unconstrained optimization. This method is a convex combination of Liu-Storey conjugate gradient method and Fletcher-Reeves conjugate gradient method. We also prove that the search direction of any hybrid conjugate gradient method, which is a convex combination of two conjugate gradient methods, satisfies the famous D-L conjugacy condition and in the same time accords with the Newton direction with the suitable condition. Furthermore, this property doesn't depend on any line search. Next, we also prove that, moduling the value of the parameter t,the Newton direction condition is equivalent to Dai-Liao conjugacy condition.The strong Wolfe line search conditions are used.The global convergence of this new method is proved.Numerical comparisons show that the present hybrid conjugate gradient algorithm is the efficient one.展开更多
In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Comb...In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Combining the quasi-Newton method with the new method, the former is modified to have global convergence property. Numerical results show that the new algorithm is efficient.展开更多
The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a m...The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.展开更多
We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorou...We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in L^∞ norm and weighted L^2-norm. The numerical examples are given to illustrate the theoretical results.展开更多
A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite differen...A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.展开更多
基金Supported by the Laoshan Laboratory(No.LSKJ202202402)the National Natural Science Foundation of China(No.42030410)+2 种基金the Startup Foundation for Introducing Talent of Nanjing University of Information Science&Technology,and Jiangsu Innovation Research Group(No.JSSCTD 202346)supported by the China National Postdoctoral Program for Innovative Talents(No.BX20240169)the China Postdoctoral Science Foundation(No.2141062400101)。
文摘Deep learning(DL)has become a crucial technique for predicting the El Niño-Southern Oscillation(ENSO)and evaluating its predictability.While various DL-based models have been developed for ENSO predictions,many fail to capture the coherent multivariate evolution within the coupled ocean-atmosphere system of the tropical Pacific.To address this three-dimensional(3D)limitation and represent ENSO-related ocean-atmosphere interactions more accurately,a novel this 3D multivariate prediction model was proposed based on a Transformer architecture,which incorporates a spatiotemporal self-attention mechanism.This model,named 3D-Geoformer,offers several advantages,enabling accurate ENSO predictions up to one and a half years in advance.Furthermore,an integrated gradient method was introduced into the model to identify the sources of predictability for sea surface temperature(SST)variability in the eastern equatorial Pacific.Results reveal that the 3D-Geoformer effectively captures ENSO-related precursors during the evolution of ENSO events,particularly the thermocline feedback processes and ocean temperature anomaly pathways on and off the equator.By extending DL-based ENSO predictions from one-dimensional Niño time series to 3D multivariate fields,the 3D-Geoformer represents a significant advancement in ENSO prediction.This study provides details in the model formulation,analysis procedures,sensitivity experiments,and illustrative examples,offering practical guidance for the application of the model in ENSO research.
基金extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through Large Research Project under grant number RGP2/174/46.
文摘In this article,we develop the Laplace transform(LT)based Chebyshev spectral collocation method(CSCM)to approximate the time fractional advection-diffusion equation,incorporating the Atangana-Baleanu Caputo(ABC)derivative.The advection-diffusion equation,which governs the transport of mass,heat,or energy through combined advection and diffusion processes,is central to modeling physical systems with nonlocal behavior.Our numerical scheme employs the LT to transform the time-dependent time-fractional PDEs into a time-independent PDE in LT domain,eliminating the need for classical time-stepping methods that often suffer from stability constraints.For spatial discretization,we employ the CSCM,where the solution is approximated using Lagrange interpolation polynomial based on the Chebyshev collocation nodes,achieving exponential convergence that outperforms the algebraic convergence rates of finite difference and finite element methods.Finally,the solution is reverted to the time domain using contour integration technique.We also establish the existence and uniqueness of the solution for the proposed problem.The performance,efficiency,and accuracy of the proposed method are validated through various fractional advection-diffusion problems.The computed results demonstrate that the proposed method has less computational cost and is highly accurate.
基金supported by the Advance Research Project of Civil Aerospace Technology(Grant No.D020304)National Nat-ural Science Foundation of China(Grant Nos.52205257 and U22B2083).
文摘This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the accuracy mismatch between tem-poral low-order finite difference and spatial high-order discre tization,the ir time collocation points must increase dramatically to solve highly oscillatory solutions of structural vibration,which results in a surge in computing time and a decrease in accuracy.To address this problem,we introduced the step-by-step idea in the space-time spectral method.The Chebyshev polynomials and Lagrange's equation were applied to derive discrete spatial goverming equations,and a matrix projection method was used to map the calculation results of prev ious steps as the initial conditions of the subsequent steps.A series of numerical experiments were carried out.The results of the proposed method were compared with those obtained by traditional space-time spectral methods,which showed that higher accuracy could be achieved in a shorter computation time than the latter in highly oscillatory cases.
文摘In this paper, we proposed a spectral gradient-Newton two phase method for constrained semismooth equations. In the first stage, we use the spectral projected gradient to obtain the global convergence of the algorithm, and then use the final point in the first stage as a new initial point to turn to a projected semismooth asymptotically newton method for fast convergence.
文摘Spectral conjugate gradient method is an algorithm obtained by combination of spectral gradient method and conjugate gradient method,which is characterized with global convergence and simplicity of spectral gradient method,and small storage of conjugate gradient method.Besides,the spectral conjugate gradient method was proved that the search direction at each iteration is a descent direction of objective function even without relying on any line search method.Spectral conjugate gradient method is applied to full waveform inversion for numerical tests on Marmousi model.The authors give a comparison on numerical results obtained by steepest descent method,conjugate gradient method and spectral conjugate gradient method,which shows that the spectral conjugate gradient method is superior to the other two methods.
基金supported by the National Natural Science Foundation of China(No.72071202)the Key Laboratory of Mathematics and Engineering Applications,Ministry of Education。
文摘As a generalization of the two-term conjugate gradient method(CGM),the spectral CGM is one of the effective methods for solving unconstrained optimization.In this paper,we enhance the JJSL conjugate parameter,initially proposed by Jiang et al.(Computational and Applied Mathematics,2021,40:174),through the utilization of a convex combination technique.And this improvement allows for an adaptive search direction by integrating a newly constructed spectral gradient-type restart strategy.Then,we develop a new spectral CGM by employing an inexact line search to determine the step size.With the application of the weak Wolfe line search,we establish the sufficient descent property of the proposed search direction.Moreover,under general assumptions,including the employment of the strong Wolfe line search for step size calculation,we demonstrate the global convergence of our new algorithm.Finally,the given unconstrained optimization test results show that the new algorithm is effective.
基金supported by the Deep Earth Probe and Mineral Resources Exploration-National Science and Technology Major Project(Grant No.2025ZD005100)by Beijing Geolight Technology Co.,Ltd.under the project“The Impact of Strong Ground Motion on Buildings”(YF-202520).
文摘This study employed tri-component continuous monitoring data from 10 measurement points on both sides of a base isolation layer in the basement of a large-span high-rise building in Beijing,as well as from a free-field station and roof frame,during a Mw 5.5 magnitude earthquake in Pingyuan,Shandong,in 2023.The H/V spectral ratio method was used to evaluate the structural dynamic response characteristics of the building and analyze the regulatory effect of the base-isolation layer on seismic waves.The results indicate that during the earthquake,the peak frequency of the free-field and the measurement points below the base-isolation layer was stable at 0.17 Hz,whereas the main frequency of the measurement points above the base-isolation layer increased to 0.75–1.18 Hz,which is 4–6 times greater than that of the points below.The amplitude was suppressed by more than 70%,confirming that the base isolation layer effectively isolated the low-frequency energy from the ground and increased the response frequency of the building.When the building was excited by an earthquake,a three-tier frequency gradient was formed throughout the building:“base-isolation layer(0.17 Hz)-main body(1.18 Hz)-roof frame(3.83 Hz)”,which can effectively avoid resonance of the entire building.In addition,the composite base-isolation device changed the dynamic characteristics of the structure.The resonance period was extended from 0.74 s(theoretical value without base isolation)to 5.9 s(calculated value),and the resonance frequency was reduced from 1.35 to 0.17 Hz.This finding indicates that the base-isolation layer can enhance seismic performance by increasing flexibility and damping.
基金National Natural Science Foundation of China under Grant Nos.U2139208 and 52278516Key Laboratory of Earthquake Engineering and Engineering Vibration,China Earthquake Administration under Grant No.2024D15Key Laboratory of Soft Soil Characteristic and Engineering Environment,Tianjin Chengjian University under Grant No.2022SCEEKL003。
文摘This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.
基金financial support from the National Natural Science Foundation of China(Grant No.52273067,52233006)the Fundamental Research Funds for the Central Universities(Grant No.2232023A-03)+3 种基金the Shuguang Program of Shanghai Education Development Foundation and Shanghai Municipal Education Commission(Grant No.23SG29)the Natural Science Foundation of Shanghai(Grant No.24ZR1402400)the Shanghai Scientific and Technological Innovation Project(Grant No.24520713000)Innovation Program of Shanghai Municipal Education Commission(Grant No.2021-01-07-00-03-E00108).
文摘Radiative cooling textiles with spectrally selective surfaces offer a promising energy-efficient approach for sub-ambient cooling of outdoor objects and individuals.However,the spectrally selective mid-infrared emission of these textiles significantly hinders their efficient radiative heat exchange with self-heated objects,thereby posing a significant challenge to their versatile cooling applicability.Herein,we present a bicomponent blow spinning strategy for the production of scalable,ultra-flexible,and healable textiles featuring a tailored dual gradient in both chemical composition and fiber diameter.The gradient in the fiber diameter of this textile introduces a hierarchically porous structure across the sunlight incident area,thereby achieving a competitive solar reflectivity of 98.7%on its outer surface.Additionally,the gradient in the chemical composition of this textile contributes to the formation of Janus infrared-absorbing surfaces:The outer surface demonstrates a high mid-infrared emission,whereas the inner surface shows a broad infrared absorptivity,facilitating radiative heat exchange with underlying self-heated objects.Consequently,this textile demonstrates multi-scenario radiative cooling capabilities,enabling versatile outdoor cooling for unheated objects by 7.8℃ and self-heated objects by 13.6℃,compared to commercial sunshade fabrics.
基金supported by the National Natural Science Foundation of China(12401482)the second author was supported by the National Natural Science Foundation of China(12371371,12261160361,11971366)supported by the Open Research Fund of Hubei Key Laboratory of Computational Science,Wuhan University.
文摘A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).
基金supported by the“Regional Innovation System&Education(RISE)”through the Seoul RISE Center,funded by the Ministry of Education(MOE)and the Seoul Metropolitan Government(2025-RISE-01-027-04).
文摘Fluid dynamic research on rectangular and trapezoidal fins is aimed at increasing heat transfer by means of large surfaces.The trapezoidal cavity form is compared with its thermal and flow performance,and it is revealed that trapezoidal fins tend to be more efficient,particularly when material optimization is critical.Motivated by the increasing need for sustainable energy management,this work analyses the thermal performance of inclined trapezoidal and rectangular porous fins utilising a unique hybrid nanofluid.The effectiveness of nanoparticles in a working fluid is primarily determined by their thermophysical properties;hence,optimising these properties can significantly improve overall performance.This study considers the dispersion of Graphene Oxide(GO)and Molybdenum Disulfide in the base fluid,engine oil.Temperature profiles are analysed by altering the radiative,porosity,wet porous,and angle of inclination parameters.Surface and contour plots are constructed by using the Lobatto IIIa Collocation Method with BVP5C solver in MATLAB and Gradient Descent Optimisation to predict the combined heat transfer rate.According to the study,fluid temperature consistently decreases when the angle of inclination,wet porous parameter,porosity parameter,and radiative parameter increase,suggesting significantly improved heat dissipation.The trapezoidal fin consistently exhibits a superior heat transfer mechanism than a rectangular fin.It is found that the trapezoidal fin transmits heat at a rate that is 0.05%higher than that of the rectangular fin.Validation of the present study is done through the comparison of previous studies.This research provides useful design insights for sophisticated engineering uses,including electrical cooling devices,heat exchangers,radiators,and solar heaters.
基金National Key Technology R & D Program-On-site Rapid Identification of Drug Research Project (Grant No. 2008BAI55B06)
文摘We analyzed the infrared 0R）-near infrared （NIR） 2D correlation spectra of drugs perturbed by temperature. By identification of functional groups by IR spectrum and by the correlation analysis of IR-NIR spectrum, we identified the characteristic spectral bands that were closely related to the structure of a drug substance of interest. These characteristic spectral bands were relatively less interfered by other ingredients for analysis by the NIR correlation coefficient method. With these characteristic spectral bands, the accuracy of screening illegally added Sildenafil citrate, Tadalafil and Metforrnin hydrochloride in Chinese patent drugs and healthcare products reached about 90%, which met the requirements of rapid screening.
基金support by the "Eleventh Five-Year" National Science and Technology Support Program (No. 2006BAB01A02)the Pivot Program of the National Natural Science Fund (No. 40930314)
文摘The conventional gravity gradient method to plot the geologic body location is fuzzy. When the depth is large and the geologic body is small, the Vzz and Vzx derivative errors are also large. We describe that using the status distinguishing factor to optimally determine the comer location is more accurate than the conventional higher-order derivative method. Thus, a better small geologic body and fault resolution is obtained by using the gravity gradient method and trial theoretical model calculation. The actual data is better processed, providing a better basis for prospecting and determination of subsurface geologic structure.
基金National Natural Science Foundation of China under Grant No.51078337,No.51108431 and No.91315301
文摘The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51309050 and 51221961)the National Basic Research Program of China(973 Program,Grant Nos.2013CB036101 and 2011CB013703)
文摘In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.
文摘In this paper, we present a new hybrid conjugate gradient algorithm for unconstrained optimization. This method is a convex combination of Liu-Storey conjugate gradient method and Fletcher-Reeves conjugate gradient method. We also prove that the search direction of any hybrid conjugate gradient method, which is a convex combination of two conjugate gradient methods, satisfies the famous D-L conjugacy condition and in the same time accords with the Newton direction with the suitable condition. Furthermore, this property doesn't depend on any line search. Next, we also prove that, moduling the value of the parameter t,the Newton direction condition is equivalent to Dai-Liao conjugacy condition.The strong Wolfe line search conditions are used.The global convergence of this new method is proved.Numerical comparisons show that the present hybrid conjugate gradient algorithm is the efficient one.
文摘In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Combining the quasi-Newton method with the new method, the former is modified to have global convergence property. Numerical results show that the new algorithm is efficient.
基金supported by The National Key Research and Development Program Plane(No.2017YFC0601505)National Natural Science Foundation(No.41672325)Science&Technology Department of Sichuan Province Technology Project(No.2017GZ0393)
文摘The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.
基金supported by NSFC Project(11301446,11271145)China Postdoctoral Science Foundation Grant(2013M531789)+3 种基金Specialized Research Fund for the Doctoral Program of Higher Education(2011440711009)Program for Changjiang Scholars and Innovative Research Team in University(IRT1179)Project of Scientific Research Fund of Hunan Provincial Science and Technology Department(2013RS4057)the Research Foundation of Hunan Provincial Education Department(13B116)
文摘We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in L^∞ norm and weighted L^2-norm. The numerical examples are given to illustrate the theoretical results.
基金the National Natural Science Foundation of China
文摘A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.
