Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

Mobility and dynamic erosion process of granular flow:insights from numerical investigation using material point method

In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility and dynamic erosion process of granular flow subjected to the complex settings,i.e.,the aspect ratio,granular mass,friction and dilatancy resistance,gravity and presence of water.A set of power scaling laws were proposed to describe the final deposit characteristics of granular flow by the relations of the normalized run-out distance and the normalized final height of granular flow against the aspect ratio,being greatly affected by the complex geological settings,e.g.,granular mass,the friction and dilatancy resistance of granular soil,and presence of water in granular flow.An index of the coefficient of friction of granular soil was defined as a ratio of the target coefficient of friction over the initial coefficient of friction to quantify the scaling extent of friction change(i.e.,friction strengthening or weakening).There is a characteristic aspect ratio of granular column corresponding to the maximum mobility of granular flow with the minimum index of the apparent coefficient of friction.The index of the repose coefficient of friction of granular flow decreased gradually with the increase in aspect ratio because higher potential energy of granular column at a larger aspect ratio causes a larger kinetic energy of granular soil to weaken the friction of granular soil as a kind of velocity-related friction weakening.An increase in granular mass reduces gradually the indexes of the apparent and repose coefficients of friction of granular soil to enhance the mobility of granular flow.The mobility of granular flow increases gradually with the decrease in friction angle or increase in dilatancy angle of granular soil.However,the increase of gravity accelerates granular flow but showing the same final deposit profile without any dependence on gravity.The mobility of granular flow increases gradually by lowering the indexes of the apparent and repose coefficients of friction of granular flow while changing the surroundings,in turn,the dry soil,submerged soil and saturated soil,implying a gradually increased excessive mobility of granular flow with the friction weakening of granular soil.Presence of water in granular flow may be a potential catalyzer to yield a long run-out granular flow,as revealed in comparison of water-absent and water-present granular flows.In addition,the dynamic erosion and entrainment of based soil induced by granular flow subjected to the complex geological settings,i.e.,the aspect ratio,granular mass,gravity,friction and dilatancy resistance,and presence of water,were comprehensively investigated as well.

Experimental investigation of closed-loop active control to modulate coherent structures by mu-level method

The experimental research on zero-net-mass-flux jet closed-loop active control was conducted in the wind tunnel.The mu-level method successfully detected burst events of the coherent structures. The streamwise velocity signals in the turbulent boundary layer were measured by HWA. The drag reduction rate of 16.7% is obtained comparable to that of the open-loop control and saves 75% of the input energy at the asynchronous 100 V/160 Hz control case, which reflects the advantages of the closed-loop control. The experimental findings indicate that the intensity increases in the near-wall region.The perturbation of the PZT vibrators on the skewness factor is concentrated in the region y+< 60. The generation of highspeed fluids is depressed and the downward effect of high-speed fluids weakens. The alteration of energy distribution and the discernible impact of modulation between structures of varying scales are observed. The correlation coefficient exhibits a strong positive correlation, which indicates that the large-scale structures produce modulation effect on small-scale ones.The occurrence of burst events is effectively suppressed. The disturbance has the characteristics of stable periodicity,positive and negative symmetry, low intermittency, and high pulsation strength. The conditional phase waveform shows that the fluctuation amplitude increases, indicating amplitude modulation effects on coherent structures.

A Numerical Investigation Based on Exponential Collocation Method for Nonlinear SITR Model of COVID-19

In this work,the exponential approximation is used for the numerical simulation of a nonlinear SITR model as a system of differential equations that shows the dynamics of the new coronavirus(COVID-19).The SITR mathematical model is divided into four classes using fractal parameters for COVID-19 dynamics,namely,susceptible(S),infected(I),treatment(T),and recovered(R).The main idea of the presented method is based on the matrix representations of the exponential functions and their derivatives using collocation points.To indicate the usefulness of this method,we employ it in some cases.For error analysis of the method,the residual of the solutions is reviewed.The reported examples show that the method is reasonably efficient and accurate.

Investigation of the epidemiology,pathogenicity and immunogenicity of Bordetella bronchiseptica isolated from cats and dogs in China from 2021 to 2023

Bordetella bronchiseptica(Bb)is recognized as a leading cause of respiratory diseases in dogs and cats.However,epidemiological data on Bb in dogs and cats in China are still limited,and there is no commercially available vaccine.Live vaccines containing Bb that are widely used abroad are generally efective but can establish latency and potentially reactivate to cause illness in some immunodefcient vaccinated recipients,raising safety concerns.In this study,34 canine-derived and two feline-derived Bb strains were isolated from 1809 canine and 113 feline nasopharyngeal swab samples collected from eight provinces in China from 2021 to 2023.The PCR results showed that the percentage of positive Bb was 22.94%(441/1922),and more than 90%of the Bb isolates had four virulence factor-encoding genes(VFGs),namely,fhaB,prn,betA and dnt.All the isolated strains displayed a multidrug-resistant phenotype.The virulence of 10 Bb strains isolated from dogs with respiratory symptoms was tested in mice,and we found that eight isolates were highly virulent.Furthermore,the eight Bb isolates with high virulence were inactivated and intramuscularly injected into mice,and three Bb strains(WH1218,WH1203 and WH1224)with the best protective efcacy were selected.Dogs immunized with these three strains exhibited strong protection against challenge with the Bb feld strain WH1218.Ultimately,the WH1218 strain with the greatest protection in dogs was selected as the vaccine candidate.Dogs and cats that received a vaccine containing 109 CFU of the inactivated WH1218 strain showed complete protection against challenge with the Bb feld strain WH1218.This study revealed that Bb is an important pathogen that causes respiratory diseases in domestic dogs and cats in China,and all the isolates exhibited multidrug resistance.The present work contributes to the current understanding of the prevalence,antimicrobial resistance,and virulence genes of Bb in domestic dogs and cats.Additionally,our results suggest that the WH1218 strain is a promising candidate safe and efcacious inactivated Bb vaccine.

A Review on Sources,Extractions and Analysis Methods of a Sustainable Biomaterial:Tannins

Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.

Drilling-based measuring method for the c-φ parameter of rock and its field application

The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.

A speckle noise suppression method based on surface waves investigation and monitoring data

The internal energy distribution of waves can be described using ocean-wave spectra.In many ways,obtaining wave spectra on a global scale is critical.Surface waves investigation and monitoring onboard the Chinese-French oceanography satellite is the first space-borne instrument for detecting wave spectra specially,which was launched on October 29,2018.It can avoid the shortage of synthetic aperture radar detection results while still having some problems,especially with the effects of speckle noise.In this study,a method to suppress the speckle noise is proposed.First,the empirical formula for background speckle noise is established.Second,many spatio-temporal representative fluctuation spectra are classified and averaged.Third,rational transfer function filtering is used to obtain speckle noise close to the along-track direction.Finally,a signal-to-noise ratio threshold is used to suppress the abnormal speckle noise.This method solves the problems existing in previous denoising methods,such as excessive denoising in the along-track direction and the inability of some abnormal noises to be denoised in the two-dimensional directional wave spectra.

A comparison study on structure-function relationship of polysaccharides obtained from sea buckthorn berries using different methods:antioxidant and bile acid-binding capacity

In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

Investigation of the block toppling evolution of a layered model slope by centrifuge test and discrete element modeling

Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.

Material point method simulation of hydro-mechanical behaviour in twophase porous geomaterials: A state-of-the-art review

The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.

Sparse Modal Decomposition Method Addressing Underdetermined Vortex-Induced Vibration Reconstruction Problem for Marine Risers

When investigating the vortex-induced vibration(VIV)of marine risers,extrapolating the dynamic response on the entire length based on limited sensor measurements is a crucial step in both laboratory experiments and fatigue monitoring of real risers.The problem is conventionally solved using the modal decomposition method,based on the principle that the response can be approximated by a weighted sum of limited vibration modes.However,the method is not valid when the problem is underdetermined,i.e.,the number of unknown mode weights is more than the number of known measurements.This study proposed a sparse modal decomposition method based on the compressed sensing theory and the Compressive Sampling Matching Pursuit(Co Sa MP)algorithm,exploiting the sparsity of VIV in the modal space.In the validation study based on high-order VIV experiment data,the proposed method successfully reconstructed the response using only seven acceleration measurements when the conventional methods failed.A primary advantage of the proposed method is that it offers a completely data-driven approach for the underdetermined VIV reconstruction problem,which is more favorable than existing model-dependent solutions for many practical applications such as riser structural health monitoring.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

Semi-analytical investigation of heat transfer in a porous convective radiative moving longitudinal fin exposed to magnetic field in the presence of a shape-dependent trihybrid nanofluid

The thermal examination of a non-integer-ordered mobile fin with a magnetism in the presence of a trihybrid nanofluid(Fe_3O_4-Au-Zn-blood) is carried out. Three types of nanoparticles, each having a different shape, are considered. These shapes include spherical(Fe_3O_4), cylindrical(Au), and platelet(Zn) configurations. The combination approach is utilized to evaluate the physical and thermal characteristics of the trihybrid and hybrid nanofluids, excluding the thermal conductivity and dynamic viscosity. These two properties are inferred by means of the interpolation method based on the volume fraction of nanoparticles. The governing equation is transformed into a dimensionless form, and the Adomian decomposition Sumudu transform method(ADSTM) is adopted to solve the conundrum of a moving fin immersed in a trihybrid nanofluid. The obtained results agree well with those numerical simulation results, indicating that this research is reliable. The influence of diverse factors on the thermal overview for varying noninteger values of γ is analyzed and presented in graphical representations. Furthermore, the fluctuations in the heat transfer concerning the pertinent parameters are studied. The results show that the heat flux in the presence of the combination of spherical, cylindrical, and platelet nanoparticles is higher than that in the presence of the combination of only spherical and cylindrical nanoparticles. The temperature at the fin tip increases by 0.705 759% when the value of the Peclet number increases by 400%, while decreases by 11.825 13% when the value of the Hartman number increases by 400%.

Experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering boundary effect

The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.

Robustness Study and Superior Method Development and Validation for Analytical Assay Method of Atropine Sulfate in Pharmaceutical Ophthalmic Solution

Background: The robustness is a measurement of an analytical chemical method and its ability to contain unaffected by little with deliberate variation of analytical chemical method parameters. The analytical chemical method variation parameters are based on pH variability of buffer solution of mobile phase, organic ratio composition changes, stationary phase (column) manufacture, brand name and lot number variation;flow rate variation and temperature variation of chromatographic system. The analytical chemical method for assay of Atropine Sulfate conducted for robustness evaluation. The typical variation considered for mobile phase organic ratio change, change of pH, change of temperature, change of flow rate, change of column etc. Purpose: The aim of this study is to develop a cost effective, short run time and robust analytical chemical method for the assay quantification of Atropine in Pharmaceutical Ophthalmic Solution. This will help to make analytical decisions quickly for research and development scientists as well as will help with quality control product release for patient consumption. This analytical method will help to meet the market demand through quick quality control test of Atropine Ophthalmic Solution and it is very easy for maintaining (GDP) good documentation practices within the shortest period of time. Method: HPLC method has been selected for developing superior method to Compendial method. Both the compendial HPLC method and developed HPLC method was run into the same HPLC system to prove the superiority of developed method. Sensitivity, precision, reproducibility, accuracy parameters were considered for superiority of method. Mobile phase ratio change, pH of buffer solution, change of stationary phase temperature, change of flow rate and change of column were taken into consideration for robustness study of the developed method. Results: The limit of quantitation (LOQ) of developed method was much low than the compendial method. The % RSD for the six sample assay of developed method was 0.4% where the % RSD of the compendial method was 1.2%. The reproducibility between two analysts was 100.4% for developed method on the contrary the compendial method was 98.4%.

Cascading multi-segment rupture process of the 2023 Turkish earthquake doublet on a complex fault system revealed by teleseismic P wave back projection method

In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.