Order Of Magnitude Difference Research Articles

Research on the accuracy of the generalized Lévy’s solution for deflection of open cylindrical shell structures

In the last century, the plate and shell theory has been developed and matured through the research of many theoretical workers at home and abroad. The generalized Lévy’s solution can be used to obtain the static response of a semi-cylindrical shell with curved ends simply supported under arbitrary load. However, the solution is complex and difficult to calculate, so on this basis, a series of simplified methods under specific working conditions are developed, such as the moment-free theory and semi-moment-free theory of cylindrical shells. In this paper, MATLAB will be used to calculate the generalized Levi’s solution and semi-moment-free theory, and compare them with the finite element results of ANSYS. Finally, it is found that the two theoretical solutions can only roughly predict the distribution form of displacement, with an order of magnitude difference. The engineering significance is limited, and the results still need to be improved.

Adaptive Ranking-Based Constraint Handling for Explicitly Constrained Black-Box Optimization.

We propose a novel constraint-handling technique for the covariance matrix adaptation evolution strategy (CMA-ES). The proposed technique is aimed at solving explicitly constrained black-box continuous optimization problems, in which the explicit constraint is a constraint whereby the computational time for the constraint violation and its (numerical) gradient are negligible compared to that for the objective function. This method is designed to realize two invariance properties: invariance to the affine transformation of the search space, and invariance to the increasing transformation of the objective and constraint functions. The CMA-ES is designed to possess these properties for handling difficulties that appear in black-box optimization problems, such as non-separability, ill-conditioning, ruggedness, and the different orders of magnitude in the objective. The proposed constraint-handling technique (CHT), known as ARCH, modifies the underlying CMA-ES only in terms of the ranking of the candidate solutions. It employs a repair operator and an adaptive ranking aggregation strategy to compute the ranking. We developed test problems to evaluate the effects of the invariance properties, and performed experiments to empirically verify the invariance of the algorithm. We compared the proposed method with other CHTs on the CEC 2006 constrained optimization benchmark suite to demonstrate its efficacy. Empirical studies reveal that ARCH is able to exploit the explicitness of the constraint functions effectively, sometimes even more efficiently than an existing box-constraint handling technique on box-constrained problems, while exhibiting the invariance properties. Moreover, ARCH overwhelmingly outperforms CHTs by not exploiting the explicit constraints in terms of the number of objective function calls.

Comparison of macroscale and microscale mechanical properties of fresh and fixed-frozen porcine colonic tissue

Mechanical changes to the microenvironment of the extracellular matrix (ECM) in tissue have been hypothesised to elicit a pathogenic response in the surrounding cells. Hence, 3D scaffolds are a popular method of studying cellular behaviour under conditions that mimic in vivo microenvironment. To create a 3D biomimetic scaffold that captures the in vivo ECM microenvironment a robust mechanical characterisation of the whole ECM at the microscale is necessary. This study examined the multiscale methods of characterising the ECM microenvironment using porcine colon tissue. To facilitate fresh tissue microscale mechanical characterisation, a protocol for sectioning fresh, unfixed, soft biological tissue was developed. Four experiments examined both the microscale and macroscale mechanics of both fresh (Fr) and fixed-frozen (FF) porcine colonic tissue using microindentation for microscale testing and uniaxial compression testing for macroscale testing. The results obtained in this study show a significant difference in elastic modulus between Fr and FF tissue at both the macroscale and microscale. There was an order of magnitude difference between the Fr and FF tissue at the microscale between each of the three layers of the colon tested i.e. the muscularis propria (MP), the submucosa (SM) and the mucosa (M). Macroscale testing cannot capture these regional differences. The findings in this study suggest that the most appropriate method for mechanically characterising the ECM is fresh microscale mechanical microindentation. These methods can be used on a range of biological tissues to create 3D biomimetic scaffolds that are more representative of the in vivo ECM, allowing for a more in-depth characterisation of the disease process.

Towards an AI-based understanding of the solar wind: A critical data analysis of ACE data

All artificial intelligence models today require preprocessed and cleaned data to work properly. This crucial step depends on the quality of the data analysis being done. The Space Weather community increased its use of AI in the past few years, but a thorough data analysis addressing all the potential issues is not always performed beforehand. Here is an analysis of a largely used dataset: Level-2 Advanced Composition Explorer’s SWEPAM and MAG measurements from 1998 to 2021 by the ACE Science Center. This work contains guidelines and highlights issues in the ACE data that are likely to be found in other space weather datasets: missing values, inconsistency in distributions, hidden information in statistics, etc. Amongst all specificities of this data, the following can seriously impact the use of algorithms: Histograms are not uniform distributions at all, but sometime Gaussian or Laplacian. Algorithms will be inconsistent in the learning samples as some rare cases will be underrepresented. Gaussian distributions could be overly brought by Gaussian noise from measurements and the signal-to-noise ratio is difficult to estimate. Models will not be reproducible from year to year due to high changes in histograms over time. This high dependence on the solar cycle suggests that one should have at least 11 consecutive years of data to train the algorithm. Rounding of ion temperatures values to different orders of magnitude throughout the data, (probably due to a fixed number of bits on which measurements are coded) will bias the model by wrongly over-representing or under-representing some values. There is an extensive number of missing values (e.g., 41.59% for ion density) that cannot be implemented without pre-processing. Each possible pre-processing is different and subjective depending on one’s underlying objectives A linear model will not be able to accurately model the data. Our linear analysis (e.g., PCA), struggles to explain the data and their relationships. However, non-linear relationships between data seem to exist. Data seem cyclic: we witness the apparition of the solar cycle and the synodic rotation period of the Sun when looking at autocorrelations.Some suggestions are given to address the issues described to enable usage of the dataset despite these challenges.

Properties of polyamide 6/multiwalled carbon nanotubes composites: The influence of processing methods

AbstractThe objective of this work was to elucidate the influence of shear rates on the properties of polyamide 6/multiwalled carbon nanotube (PA6/CNT) composites which was realized by adopting different types of processing methods that feature different orders of magnitude in shear rates, such as compression molding (CM, ~0 s−1), conventional injection molding (CIM, ~102 s−1) and microinjection molding (μIM, ~105 s−1). Electrical conductivity (σ) results indicated that the prevailing high shearing conditions in injection molding was unfavorable for the formation of intact filler network, thereby resulting in a much lower σ than CM counterparts. Moreover, the σ of PA6/CNT microparts was higher than that of CIM macroparts when the filler content was less than 5 wt%, otherwise the σ of CIM macroparts prevailed over that of μIM counterparts. A better filler distribution was observed when PA6/CNT composites were processed under higher shearing conditions, as corroborated by SEM. In addition, CNTs were preferentially aligned along flow direction and a higher degree of CNT orientation was expected with increasing shear rates, as confirmed by Raman spectral analysis. The tensile strength of injection molded PA6/CNT samples increased with increasing filler concentrations, and the more preferential orientation and better distribution of CNT were considered to be the contributing factors. The comparative study of the properties of PA6/CNT composites that processed using different methods was important for their practical applications in industrial sectors.

The impact of bimodal pore size distribution and wettability on relative permeability and capillary pressure in a microporous limestone with uncertainty quantification

Pore-scale X-ray imaging combined with a steady-state flow experiment was used to study the displacement processes during waterflooding in an altered-wettability carbonate, Ketton limestone, with more than two orders of magnitude difference in pore size between macropores and microporosity. We simultaneously characterized macroscopic and local multiphase flow parameters, including relative permeability, capillary pressure, wettability, and fluid occupancy in pores and throats. An accurate method was applied for porosity and fluid saturation measurements using greyscale based differential imaging without image segmentation. The relative permeability values were corrected by considering the measured saturation profile along the sample length to account for the so-called capillary end effect. The behaviour of relative permeability and capillary pressure was compared to other measurements in the literature to demonstrate the combined effects of wettability and pore structure. Typical oil-wet behaviour in resolvable macropores was measured from contact angle, fluid occupancy and curvature. The capillary pressure was negative while the oil relative permeability dropped quickly as oil was drained to low saturation and flowed through connected oil layers. Brine initially largely flowed through water-wet microporosity, and then filled the centre of large oil-wet pore bodies. Thus, the brine relative permeability remained exceptionally low until brine formed a connected flow path in the macropores leading to a substantial increase in relative permeability. Overall, this work demonstrates that not only wettability but also pore size distribution and microporosity have significant impact on displacement processes.

Self-charging of sprays

The charging of poorly conducting liquids due to flows is a well-known phenomenon, yet the precise charging mechanism is not fully understood. This is especially relevant for sprays, where the spray plume dynamics and maximum distance travelled of a spray dramatically changes for different levels of charging: charging of the droplets makes them repel, thereby preventing drop coalescence and altering the shape of the spray plume. As the charging depends on many factors including the flow and the interactions between the liquid and the nozzle, many models and scaling laws exist in the literature. In this work we focus on perhaps the simplest flow regime, laminar jets created by ultra short channels, and quantify the charging as a function of the different parameters. We present a simple model that collapses all the data for over 4 orders of magnitude difference in streaming currents for various nozzle sizes, flow velocities and surface treatments. We further show that the charging polarity can even be reversed by applying an oppositely charged coating to the nozzle, an important step for any application.

Differences in water and vapor transport through angstrom-scale pores in atomically thin membranes

The transport of water through nanoscale capillaries/pores plays a prominent role in biology, ionic/molecular separations, water treatment and protective applications. However, the mechanisms of water and vapor transport through nanoscale confinements remain to be fully understood. Angstrom-scale pores (~2.8–6.6 Å) introduced into the atomically thin graphene lattice represent ideal model systems to probe water transport at the molecular-length scale with short pores (aspect ratio ~1–1.9) i.e., pore diameters approach the pore length (~3.4 Å) at the theoretical limit of material thickness. Here, we report on orders of magnitude differences (~80×) between transport of water vapor (~44.2–52.4 g m−2 day−1 Pa−1) and liquid water (0.6–2 g m−2 day−1 Pa−1) through nanopores (~2.8–6.6 Å in diameter) in monolayer graphene and rationalize this difference via a flow resistance model in which liquid water permeation occurs near the continuum regime whereas water vapor transport occurs in the free molecular flow regime. We demonstrate centimeter-scale atomically thin graphene membranes with up to an order of magnitude higher water vapor transport rate (~5.4–6.1 × 104 g m−2 day−1) than most commercially available ultra-breathable protective materials while effectively blocking even sub-nanometer (>0.66 nm) model ions/molecules.

Estimation of the space weathering timescale on (25143) Itokawa: Implications on its rejuvenation process

Context. The space weathering timescales of near-Earth S-type asteroids have been investigated using several approaches (i.e., experiments, sample analyses, and theoretical approaches), yet there are orders of magnitude differences. Aims. We aim to examine the space weathering timescale on a near-Earth S-type asteroid, Itokawa, using Hayabusa/AMICA images, and further investigate the evolutional process of the asteroid. Methods. We focused on bright mottles on the boulder surfaces generated via impacts with interplanetary dust particles (IDPs). We compared the bright mottle size distribution with an IDP flux model to determine the space weathering timescale. Results. We found that the space weathering timescale on Itokawa’s boulder surfaces is 103 yr (in the range of 102–104 yr), which is consistent with the timescales of space weathering by light ions from the solar wind. Conclusions. From this result, we conclude that Itokawa’s surface has been weathered rapidly in 103 yr but portions of the surface are exposed due to seismic shaking that was triggered by a recent impact, which created the Kamoi crater.

Quenching transitions for the rovibrational transitions of water: Ortho-H2O in collision with ortho- and para-H2.

We present here the first full computation of the rovibrational quenching of a polyatomic molecule (water) by a rotating molecular projectile (H2). The computation is performed for quenching from the first bending mode of water at ν ≃ 1595cm-1 with a rotation energy of up to ∼400cm-1 in the bending mode. Molecular hydrogen is in its para and ortho modifications; it is rotating with a rotational quantum number of up to 4 and 3, respectively. All computations are performed on a very reliable and fully tested potential water-hydrogen energy surface of full dimensionality. Dynamics is performed in the full coupled channel formalism in the rigid bender approximation with a decoupling of the water rotation and vibration bases. Rate coefficients are converged for a kinetic temperature range 50-500K. The crucial importance of the proper treatment of the projectile rotation is emphasized with orders of magnitude differences between the different channels for the H2 rotation. Sensitivity to the actual rovibrational initial state of water exists but in a weaker manner. Overall quenching rate coefficients are about 10-12cm3s-1, remaining one to three orders of magnitude lower than pure rotational quenching. They should be employed to model denser and warmer astrophysical media, such as high atmospheres or star and planet forming regions, which are to be explored by infrared space telescopes, such as JWST.

The impact of time-dependent stellar activity on exoplanet atmospheres

ABSTRACTM-dwarfs are thought to be hostile environments for exoplanets. Stellar events are very common on such stars. These events might cause the atmospheres of exoplanets to change significantly over time. It is not only the major stellar flare events that contribute to this disequilibrium, but the smaller flares might also affect the atmospheres in an accumulating manner. In this study, we aim to investigate the effects of time-dependent stellar activity on the atmospheres of known exoplanets. We simulate the chemistry of GJ 876c, GJ 581c, and GJ 832c that go from H2-dominated to N2-dominated atmospheres using observed stellar spectra from the MUSCLES collaboration. We make use of the chemical kinetics code vulcan and implement a flaring routine that stochastically generates synthetic flares based on observed flare statistics. Using the radiative transfer code petitradtrans, we also simulate the evolution of emission and transmission spectra. We investigate the effect of recurring flares for a total of 11 d covering 515 flares. Results show a significant change in abundance for some relevant species such as H, OH, and CH4, with factors going up to three orders of magnitude difference with respect to the pre-flare abundances. We find a maximum change of ∼12 ppm for CH4 in transmission spectra on GJ 876c. These changes in the spectra remain too small to observe. We also find that the change in abundance and spectra of the planets accumulates throughout time, causing permanent changes in the chemistry. We conclude that this small but gradual change in chemistry arises due to the recurring flares.

A Closed-Form Equation for Capillary Pressure in Porous Media for All Wettabilities

A saturation–capillary pressure relationship is proposed that is applicable for all wettabilities, including mixed-wet and oil-wet or hydrophobic media. This formulation is more flexible than existing correlations that only match water-wet data, while also allowing saturation to be written as a closed-form function of capillary pressure: we can determine capillary pressure explicitly from saturation, and vice versa. We propose Pc=A+Btanπ2-πSeCfor0≤Se≤1,\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$P_{{\ ext{c}}} = A + B\ an \\left( {\\frac{\\pi }{2} - \\pi S_{e}^{C} } \\right)\\,{\ ext{for}}\\,0 \\le S_{{\ ext{e}}} \\le 1,$$\\end{document}where S_{{text{e}}} is the normalized saturation. A indicates the wettability: A>0 is a water-wet medium, A<0 is hydrophobic while small A suggests mixed wettability. B represents the average curvature and pore-size distribution which can be much lower in mixed-wet compared to water-wet media with the same pore structure if the menisci are approximately minimal surfaces. C is an exponent that controls the inflection point in the capillary pressure and the asymptotic behaviour near end points. We match the model accurately to 29 datasets in the literature for water-wet, mixed-wet and hydrophobic media, including rocks, soils, bead and sand packs and fibrous materials with over four orders of magnitude difference in permeability and porosities from 20% to nearly 90%. We apply Leverett J-function scaling to make the expression for capillary pressure dimensionless and discuss the behaviour of analytical solutions for spontaneous imbibition.

Creep Behavior of Fe–25Ni–15Cr‐Based Alumina‐Forming Austenitic Steel

The creep behavior of Fe–25Ni–15Cr‐based aluminum‐forming austenitic stainless steel is investigated at 700 °C, under the selected stress of 100, 140, 170, and 180 MPa, respectively. Creep properties are analyzed in terms of stress exponent (n) and creep damage tolerance factor (λ). Within the tested stress range, there is an order of magnitude difference between the minimum creep rate and steady creep rate. The evolution of the microstructures after creep test with the stresses of 140, 170, and 180 MPa is characterized by scanning electron microscopy and electron backscatter diffraction detectors, respectively. With the decreased creep stress, the creep void size shows an increasing trend. Higher stress would accelerate the process of recrystallization during the creep process at 700 °C. The proportion of high‐angle grain boundaries increases to a certain extent under high‐stress conditions. A large amount of Laves and NiAl phases are observed after creep test at 140, 170, and 180 MPa, respectively.

Determining baseline radiation levels in marine biota – A comparison of SE Queensland commercial species

Limited data exists on the contribution to radiological dose to members of the public from ingestion of radioactivity in seafood in the Australian diet. There is also a lack of data available to assess the radiological dose to marine fauna in Australian waters. Natural and anthropogenic radionuclides and trace metals were measured in the edible flesh of prawns to determine the radiological dose to humans. The remaining tissues were combined and analysed to enable uptake and environmental radiological doses to be assessed. Although in international studies, the edible flesh is generally measured to determine radiation dose or ingestion of trace metals, the effects of preparation and cooking techniques are rarely assessed. In this study, cooking and preparation techniques that may influence the radiological dose to humans were assessed. Eggs were also removed from a selected number of samples to assess the potential dose in sensitive developmental tissues and possible implications for environmental effects. Order of magnitude differences in 210Po activity concentrations and Cd concentrations were observed between whole animals and the edible flesh of Australian caught King (Melicertus spp.) and Tiger (Penaeus spp.) prawns, with the hepatopancreas primarily responsible for this difference. Most 210Po was unsupported by 210Pb and activity concentrations of all other radionuclides measured (137Cs, 210Pb, 226Ra, 232Th, and 238U) were very low. The major contribution to radiation dose via the ingestion pathway and to the organism itself was from 210Po. Cooking techniques that may lead to leaching of 210Po from the hepatopancreas could substantially increase the radiological dose from ingestion of this isotope. Organism dose estimates using different input assumptions in the radiological assessment tool “ERICA”, including site-specific tissue activity concentrations with site- or region-specific media concentrations, were compared with ERICA default distribution coefficients (Kd) and concentration ratios.

Distribution and carbon isotopic composition of diploptene from epiphytic bryophytes in Wuhan, central China

Diploptene is a ubiquitous hopanoid in the geosphere, synthesized by all hopanoid-containing bacteria. Variations in the concentration and stable carbon isotopic composition (δ13C) of diploptene in ancient peats and lignite can be used to reconstruct certain aspects of the wetland methane cycle in the past. However, the sources and mechanisms that control diploptene δ13C values in wetlands are not fully constrained. To address this, here we determined the distribution and δ13C values of diploptene, as well as n-alkanes, obtained from five genera of epiphytic bryophytes (non-vascular plants such as mosses) that occupy-three different habitats: soil, rock, and tree bark. Our data show that the concentrations of diploptene are highly variable with two order of magnitude differences between the various species. Mosses collected from the soil habitat had higher concentrations compared to those from rock and tree habitats. This suggests that the input from some habitats might dominate the sedimentary signal. The δ13C values of diploptene (δ13Cdip) also vary between species with values ranging between –39.2‰ and –31.2‰. Generally, the δ13C values of diploptene and long chain n-alkanes (i.e., C29 and C31) are similar (±2‰) in most of the bryophyte species. This may suggest that diploptene is produced by heterotrophic bacteria that live in symbiosis with the mosses. However, for some bryophytes the δ13Cdip values are much more 13C depleted (>–2‰) compared to long chain n-alkanes, implying that for some mosses bacterial methanotrophs or methylotrophs may contribute to the diploptene pool. Our findings expand our understanding of the biological sources of diploptene in terrestrial epiphytic bryophytes, which will allow for a more detailed interpretation of the long chain n-alkanes and diploptene (δ13C values) in past environmental and paleoclimatic reconstructions.

Solar cycle and seasonal variability of H in the upper atmosphere of Mars

The complex seasonal, dust, and atmospheric circulation patterns at Mars affect the abundance and spatial extent of water and its photo-dissociated atomic byproduct H. Recent observations by the Mars Atmosphere and Volatile Evolution (MAVEN) mission have been used to constrain the H abundance at the exobase (200 km altitude) as well as its thermal and global escape rate. The MAVEN observations span the last seven years and encompass the effects of the declining phase of Solar Cycle 24 and activating phase of Solar Cycle 25 from mid-perihelion of Mars Year (MY) 32 through most of aphelion of MY36. Results show that there are seasonal as well as solar cycle variations in H exobase temperatures, densities, thermal escape rates, and global escape rates. During perihelion of MYs 32–35, the H exobase density fluctuates with a cadence that corresponds to the approach to perihelion as well as the onset and peak of regional dust storms. Additionally, the abundance of H atoms did not show exceptional variability during the MY34 Global Dust Storm, when compared with trends from adjacent years. Seasonal averages of H exobase densities, thermal escape rates, and global escape rates indicated an order of magnitude difference between perihelion and aphelion, and suggest a water global equivalent layer range of 0.02–13 m lost over 4 billon years. H properties are also found to vary with solar cycle to a lesser extent than with seasons across the MAVEN timeline. Trends show that the exobase density varies inversely with solar activity, while both the thermal and global escape rates vary directly with solar activity. It remains to be seen if the effects of a stronger than recent solar cycle on the upper atmospheric H properties would become more pronounced.

Revisiting Shikhmurzaev’s Approach to the Contact Line Problem

In this paper, we revisit a model for the contact line problem which has been proposed by Shikhmurzaev (Int. J. Multiph. Flow 19(4):589–610, 1993). In the first part, in addition to rederiving the model, we study in detail the assumptions required to obtain the isothermal limit of the model. We also derive in this paper several lubrication approximation models, based on Shikhmurzaev’s approach. The first two lubrication models describe thin film flow of incompressible fluids on solid substrates, based on different orders of magnitude of the slip length parameter. The third lubrication model describes a meniscus formation where a wedge-shaped solid immerses in a thin film of fluid.

A splitting uniformly convergent method for one-dimensional parabolic singularly perturbed convection-diffusion systems

In this paper we deal with solving robustly and efficiently one-dimensional linear parabolic singularly perturbed systems of convection-diffusion type, where the diffusion parameters can be different at each equation and even they can have different orders of magnitude. The numerical algorithm combines the classical upwind finite difference scheme to discretize in space and the fractional implicit Euler method together with an appropriate splitting by components to discretize in time. We prove that if the spatial discretization is defined on an adequate piecewise uniform Shishkin mesh, the fully discrete scheme is uniformly convergent of first order in time and of almost first order in space. The technique used to discretize in time produces only tridiagonal linear systems to be solved at each time level; thus, from the computational cost point of view, the method we propose is more efficient than other numerical algorithms which have been used for these problems. Numerical results for several test problems are shown, which corroborate in practice both the uniform convergence and the efficiency of the algorithm.

Complete characterization of sink-strengths for 1D to 3D mobilities of defect clusters: Bridging between limiting cases with effective sink-strengths calculations.

In a companion paper, new analytical expressions of the cluster sink-strength (CSS) for two one-dimensional diffusers were proposed. These expressions are indispensable to any complete parameterisation of rate-equations cluster dynamics intending to model mobile interstitial clusters such as dislocation loops. Indeed, simulating the long-term microstructural evolution in systems with fast diffusing species such as self-interstitial atom (SIA) clusters relies on establishing the CSS for rate-equation cluster dynamics which critically depends on the complex mobilities of clusters. In this second paper, we will estimate effective CSSs by object kinetic Monte-Carlo (OKMC) on wide ranges of radii, rotation energies, diffusion coefficients, and concentrations of both reaction partners. The symmetric roles of some parameters are used to infer a generic form of a semi-analytical expression of the CSS depending on all these interaction parameters. It is composed of, on one hand, the various analytical expressions established as limiting cases and, on the other hand, of fitted transition functions that allow a gradual switching between them. The analysis of the residuals shows that the overall agreement is reasonably good: it is only in the transition zones that some significant discrepancies are located. This semi-analytical expression answers our practical need for CSSs evaluation in our typical range of conditions, but further extending this range to much smaller diffusion coefficients ratios, it is quite striking to see that the domain for 1D-1D mobility is very extended. It is only when approaching a ratio of 10−6 that the 1D-0 CSSs start becoming more relevant than the adapted 1D-1D CSSs which have completely different orders of magnitude and kinetic orders. Finally, the semi-analytical CSS expression is validated by comparing its implementation in rate-equation cluster dynamics to OKMC simulations of interstitial clusters aggregation. Companion paper: [1]

Spatial sputtering of fused silica after a laser-induced exploding caused by a 355 nm nd: YAG laser

To alleviate laser-induced threats in the high-power laser systems, it is necessary to understand the debris propagation process and distribution characteristics, then take protection and cleaning strategies to sustain an ultra-clean environment. This work presents an experimental investigation of the sputtering behaviors of the laser-induced fused silica debris to clarify their longitudinal propagation and bottom distribution on the millimeter scale. Two types of polished surfaces with more than three orders of magnitude differences in surface roughness were adopted. The result reveals the centralized/decentralized distribution and the ability of long-range movement of debris propelled by the melting and mechanical failure. The maximum largest dispersion angle (LDA) appears at 18J/cm2 for debris within 10–20 μm from the Two Sides polished surface; debris &amp;lt;30 μm dominate the major quantity in the longitudinal propagation, and their centralization ability profoundly relies on the laser fluence and the surface roughness of the target; the diameter of the debris with long-range movement generally below 20 μm. This study could guide for establishing the overall dynamic cleaning strategy inside the high-power laser systems.