Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

Codesign of an integrated metal–insulator–semiconductor photocathode for photoelectrochemical reduction of CO 2 to ethylene

Photoelectrochemical carbon-dioxide reduction (PEC CO 2 R) is a potentially attractive means for producing chemicals and fuels using sunlight, water, and carbon dioxide; however, this technology is in its infancy. To date, most studies of PEC CO 2 R have reported products containing one carbon atom (C 1 products) but the production of valuable products containing two or more carbons (C 2+ products), such as ethylene, ethanol, etc., is rarely demonstrated. Metal–semiconductor–insulator (MIS) photocathode/catalyst structures offer a promising approach for this purpose, since they integrate the functions of light absorption, charge separation, and catalysis. In this study, we have investigated a Cu/TiO 2 /p-Si photocathode/catalyst structure with the aim of establishing the effects of semiconductor–insulator interactions on the performance of the photocathode and the influence of the direction of illumination of the MIS structure on the total current density and the distribution of products formed by on the Cu catalyst. Here, we have also examined the influence of ionomer coatings deposited on the Cu surface on the total current density and the distribution of products formed. A major finding is that for a fixed Cu potential the distribution of products formed by PEC CO 2 R are the same, irrespective of the direction of illumination, and are identical to those obtained by electrochemical reduction of CO 2 (EC CO 2 R). Another important finding is that the total current density and the faradaic efficiency to ethylene are enhanced significantly by deposition of a thin bilayer of Sustainion/Nafion onto the surface of the Cu.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Switching current reduction in magnetoresistive random access memories

We present an approach for minimizing the critical current for the magnetization switching in magnetic tunnel junctions by optimizing the spatial distribution of the current density. We show that such a minimization is possible because critical current is determined by the condition of making one of the magnetization eigenstates grow in time. The excitation of the eigenstates is enhanced when the spatial distributions of the eigenstates and current density overlap. Critical current can be viewed as a functional of the current density spatial distribution and it can be minimized by optimizing this distribution. Such an optimization results in a major reduction of the critical current and increase of the switching efficiency, viz. the ratio between the energy barrier and critical current. The minimized critical current increases approximately linearly with the magnetic tunnel junction size, which is much slower than critical current for the case of a uniform current density. The optimized efficiency can be approximately a constant with respect to the magnetic tunnel junction size, which is much higher than the efficiency for the uniform current density. Additional optimization can be achieved by spatially modulating the material parameters, e.g., the saturation magnetization. The presented approach and obtained scaling of the critical current and efficiency offers opportunities for the magnetic tunnel junction optimization.

36 MATERIALS SCIENCE↗

Effect of applied potential on metal surfaces: Surface energy, Wulff shape and charge distribution

Here we use grand canonical density functional theory to predict the surface energies, Wulff shapes, charge distributions and catalytically active sites of different metal surfaces under electrochemical conditions. We propose a method for computing surface energies from grand canonical density functional theory (GC-DFT) calculations of periodic slab models and use it to compute the surface energies of the facets of Pt, Cu, and Ag crystals to predict their Wulff shapes under electrochemical conditions. GC-DFT predicts that, for the pure metals studied, solvation only slightly affects the Wulff shape while applied potentials considerably affect the surface energies and corresponding Wulff shapes. We used Bader charge analysis of GC-DFT computed electron densities to investigate the effect of applied potential on the distribution of electron density over the atoms of the surfaces of Pt, Cu, Ag, and the 75–25 Ag-Pt and Au-Ni alloys. This analysis shows that, under an applied potential, the electron density is unevenly distributed over the surface atoms and that the charges of atoms more exposed to solvent are more sensitive to bias. Our results show that the most sensitive atom to bias can be used to identify the most favorable adsorption site and thus, the active sites of electrochemical reactions, which is computationally less demanding than calculating the adsorption energies on all possible adsorption sites.

36 MATERIALS SCIENCE↗

Out-of-distribution detection with non-parametric density estimation for models predicting processing history of uranium ore concentrates

The rapid advancement in machine learning (ML) and computer vision (CV) coincides with the growth of interest in deploying these ML/CV models in numerous fields from medicine to social science. Similar to those areas, we have witnessed a great number of works in materials science employing ML/CV models – neural networks in particular – in their studies in recent years. These models have proven to obtain accurate performance in various tasks. However, these models struggle to attain a similar performance when encountering test samples coming from a distribution that is different from the training set. More importantly, they fail without providing any warning to the users. Therefore, we propose a framework for detecting out-of-distribution (OOD) samples to alert users when a human intervention might be necessary in this work. Specifically, we explore the use of a non-parametric density estimation method to detect OOD samples. Here, we assess OOD detection capability of the proposed framework on ML models developed for categorizing precipitation routes of U 3 O 8 when encountering OOD datasets that contain samples (1) undergone different imaging acquisition process, (2) undergone different material synthesis process, and (3) different materials than ID set. Through those experiments, we achieve an average area under the receiver operating characteristic (AUROC) of at least 91% on average in detecting OOD samples. With minimal overhead cost and superior performance, the proposed framework enables a reliable and safe system when deploying in real-world scenarios.

Convolutional neural networks↗

Particle Size and Density Measurement

This project has developed a miniature gravity settling classifier (flow cell) and flow system that separates slurry particles by settling rate into six different Cuts. The settling rates for particles in each Cut were determined from COMSOL computational fluid dynamics (CFD) modeling of the flow in the flow cell together with particle tracking for particles of differing settling rates. For particles collected in each Cut, the particle concentration and size distribution were measured. The effective particle densities from each Cut are calculated from their respective CFD modeled settling rates using a relationship between the unhindered settling rate, particle size distribution (PSD), and effective density. The particle size and density distribution (PSDD) of the slurry sample is then calculated by combining the size and density distributions determined for the particles collected in each of the settling rate Cuts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A predictive discrete-continuum multiscale model of plasticity with quantified uncertainty

Multiscale models of materials, consisting of upscaling discrete simulations to continuum models, are unique in their capability to simulate complex materials behavior. The fundamental limitation in multiscale models is the presence of uncertainty in the computational predictions delivered by them. In this work, a sequential multiscale model has been developed, incorporating discrete dislocation dynamics (DDD) simulations and a strain gradient plasticity (SGP) model to predict the size effect in plastic deformations of metallic micro-pillars. The DDD simulations include uniaxial compression of micro-pillars with different sizes and over a wide range of initial dislocation densities and spatial distributions of dislocations. An SGP model is employed at the continuum level that accounts for the size-dependency of flow stress and hardening rate. Sequences of uncertainty analyses have been performed to assess the predictive capability of the multiscale model. The variance-based global sensitivity analysis determines the effect of parameter uncertainty on the SGP model prediction. The multiscale model is then constructed by calibrating the continuum model using the data furnished by the DDD simulations. A Bayesian calibration method is implemented to quantify the uncertainty due to microstructural randomness in discrete dislocation simulations (density and spatial distribution of dislocations) on the macroscopic continuum model prediction (size effect in plastic deformation). Here, the outcomes of this study indicate that the discrete-continuum multiscale model can accurately simulate the plastic deformation of micro-pillars, despite the significant uncertainty in the DDD results. Additionally, depending on the macroscopic features represented by the DDD simulations, the SGP model can reliably predict the size effect in plasticity responses of the micropillars with below 10% of error.

36 MATERIALS SCIENCE↗

Geostatistical interpolation of streambed hydrologic attributes with addition of left censored data and anisotropy

Spatial geostatistical interpolation of point measurements of streambed attributes in the hyporheic zone may be constrained by the streambed anisotropy, and data density and spatial distribution may significantly impact the results. Spatial clustering and low spatial data density can be caused by bedrock outcropping at the streambed limiting installation of in-stream piezometers. This study examines parameter error variability of the geostatistical interpolation using anisotropic interpolation methods and increasing the data density by adding left censored values (i.e., data below measurement limit) to locations where measurements were limited by exposed bedrock lining the streambed. The reduction in relative standard error of the interpolation was determined for the spatial distributions of streambed attributes including hydraulic conductivity, seepage flux, and mercury solute flux measured in two different years along a study reach in East Fork Poplar Creek, Tennessee, USA. Here, two methods to impute the left censored values were compared including the conventional half the detection limit substitution method, and the Stochastic Approximation of Expectation-Maximization (SAEM) algorithm, which both had comparable results. Imputing left censored data increased the data density to recommended ranges, reduced data clustering, increased the spatial dependence for some attributes, and reduced the standard error for each of the three attributes. For the reach considered herein, addition of the left censored values resulted in a larger error reduction than the consideration of anisotropy within the interpolation, which confirms the benefit of data addition to increase data density within data-limited river corridors.

58 GEOSCIENCES↗

Cuspy dark matter density profiles in massive dwarf galaxies

ABSTRACT Rotation curves of galaxies probe their total mass distributions, including dark matter. Dwarf galaxies are excellent systems to investigate the dark matter density distribution, as they tend to have larger fractions of dark matter compared to higher mass systems. The core-cusp problem describes the discrepancy found in the slope of the dark matter density profile in the centres of galaxies (β*) between observations of dwarf galaxies (shallower cores) and dark matter-only simulations (steeper cusps). We investigate β* in six nearby spiral dwarf galaxies for which high-resolution CO J = 1–0 data were obtained with ALMA (Atacama Large Millimeter/submillimeter Array). We derive rotation curves and decompose the mass profile of the dark matter using our CO rotation curves as a tracer of the total potential and 4.5 $\mu$m photometry to define the stellar mass distribution. We find 〈β*〉 = 0.6 with a standard deviation of ±0.1 among the galaxies in this sample, in agreement with previous measurements in this mass range. The galaxies studied are on the high stellar mass end of dwarf galaxies and have cuspier profiles than lower mass dwarfs, in agreement with other observations. When the same definition of the slope is used, we observe steeper slopes than predicted by the FIRE and NIHAO simulations. This may signal that these relatively massive dwarfs underwent stronger gas inflows towards their centres than predicted by these simulations, that these simulations overpredict the frequency of accretion or feedback events, or that a combination of these or other effects are at work.

Cooke, Lauren H. (ORCID:0000000284337553)↗

Spatial Control of Substitutional Dopants in Hexagonal Monolayer WS 2 : The Effect of Edge Termination

The ability to control the density and spatial distribution of substitutional dopants in semiconductors is crucial for achieving desired physicochemical properties. Substitutional doping with adjustable doping levels has been previously demonstrated in 2D transition metal dichalcogenides (TMDs); however, the spatial control of dopant distribution remains an open field. In this work, edge termination is demonstrated as an important characteristic of 2D TMD monocrystals that affects the distribution of substitutional dopants. Particularly, in chemical vapor deposition (CVD)-grown monolayer WS 2 , it is found that a higher density of transition metal dopants is always incorporated in sulfur-terminated domains when compared to tungsten-terminated domains. Here two representative examples demonstrate this spatial distribution control, including hexagonal iron- and vanadium-doped WS 2 monolayers. Density functional theory (DFT) calculations are further performed, indicating that the edge-dependent dopant distribution is due to a strong binding of tungsten atoms at tungsten-zigzag edges, resulting in the formation of open sites at sulfur-zigzag edges that enable preferential dopant incorporation. Based on these results, it is envisioned that edge termination in crystalline TMD monolayers can be utilized as a novel and effective knob for engineering the spatial distribution of substitutional dopants, leading to in-plane hetero-/multi-junctions that display fascinating electronic, optoelectronic, and magnetic properties.

2D transition metal dichalcogenides↗

Freeze-in from preheating

We consider the production of dark matter during the process of reheating after inflation. The relic density of dark matter from freeze-in depends on both the energy density and energy distribution of the inflaton scattering or decay products composing the radiation bath. Here, we compare the perturbative and non-perturbative calculations of the energy density in radiation. We also consider the (likely) possibility that the final state scalar products are unstable. Assuming either thermal or non-thermal energy distribution functions, we compare the resulting relic density based on these different approaches. We show that the present-day cold dark matter density can be obtained through freeze-in from preheating for a large range of dark matter masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evolved Gas Analysis–Mass Spectrometry Exposes Polymer Network Structures

Polymer network structures in epoxy thermosets play an important role in the final thermoset material properties. However, analytical characterization of these network structures is difficult due to their amorphous nature. In this work, the application of evolved gas analysis–mass spectrometry (EGA-MS) to characterize the polymer network structures of bisphenol A (BPA)-based thermosets is demonstrated. Analytical characterization of the polymer network structures is accomplished by monitoring the Product-Specific Kinetics (PSK) of BPA monomer formation during thermal degradation investigations. We relate observed differences in the activation energy (E a ) of BPA monomer formation to the local packing environment around the BPA monomer units within the polymer network. Variations in the local environment related to the polymer networks manifest qualitatively as broadening in the thermal profile of the BPA monomer evolution and quantitatively as changes in the activation energy (E a ). Three BPA thermoset formulations were investigated; two amine-cured thermoset with 4,4′-diaminodiphenylmethane (DDM) or poly(propylene glycol) bis(2-amino-propyl ether) (PPG400) and a homopolymerized thermoset via curing with Epikure 3253 catalyst (3253). Results revealed that the 3253 thermoset contained two distinct packing densities in the polymer network, while DDM and PPG400 thermosets had uniform distributions of packing densities. Results from the DDM thermoset revealed a gradually decreasing E a , while the apparent E a of PPG400 was consistent over the entire degradation. Furthermore, these differences in E a were concluded to stem from the flexibility of the corresponding polymer networks and the ability of the network components to rearrange and occupy formed voids. Due to the minimal sample required for analysis (100–200 μg), this EGA-MS technique has great potential for postproduction evaluation of composite parts to identify changes in the polymer networks from use and aging, which could signal compromised performance.

Degradation↗

Runaway electron beam dynamics at low plasma density in DIII-D: energy distribution, current profile, and internal instability

Parameters of the post-disruption runaway electron (RE) beam in the collisionless background plasma achieved after deuterium injection are investigated in DIII-D. The spatially resolved RE energy distribution function is measured for the first time during the RE plateau stage by inverting hard X-ray bremsstrahlung spectra. It has maximum energy up to 20 MeV and a non-monotonous feature at 5{6 MeV observed only in the core of the beam supporting the possibility of kinetic instabilities. The RE current profile is reconstructed for the first time using the spatially resolved RE energy distribution. It is found to be more peaked than the pre-disruption plasma current, with higher internal inductance, suggesting preferential formation of REs in the core plasma or potentially a radially inward motion of the REs. The accessed relatively low-current (180 kA) RE beam is found stable, likely due to its elevated q profile. From this base stable equilibrium, instability is accessed by ramping the solenoid current. Under a large applied accelerating loop voltage an internal kink mode is observed. The kink mode leads to sawtooth-like relaxation of the RE current profile, but drives no RE loss. Under conditions of a large decelerating voltage, chirping Alfvénic instabilities driven by REs during the RE plateau stage are observed for the first time. The low-frequency (0.1{7 MHz) Alfvénic instabilities, supposedly compressional Alfvén eigenmodes, are found to be correlated with a RE loss signal. Finally, this provides novel plasma conditions to study these instabilities and investigate opportunities to utilize them for RE control. Bifurcated electron cyclotron emission (ECE) spectrum is measured during the collisionless RE plateau, with a break point at ≈100 GHz suggesting resonant absorption of the ECE at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A catalogue of cataclysmic variables from 20 yr of the Sloan Digital Sky Survey with new classifications, periods, trends, and oddities

ABSTRACT We present a catalogue of 507 cataclysmic variables (CVs) observed in SDSS I to IV including 70 new classifications collated from multiple archival data sets. This represents the largest sample of CVs with high-quality and homogeneous optical spectroscopy. We have used this sample to derive unbiased space densities and period distributions for the major sub-types of CVs. We also report on some peculiar CVs, period bouncers and also CVs exhibiting large changes in accretion rates. We report 70 new CVs, 59 new periods, 178 unpublished spectra, and 262 new or updated classifications. From the SDSS spectroscopy, we also identified 18 systems incorrectly identified as CVs in the literature. We discuss the observed properties of 13 peculiar CVS, and we identify a small set of eight CVs that defy the standard classification scheme. We use this sample to investigate the distribution of different CV sub-types, and we estimate their individual space densities, as well as that of the entire CV population. The SDSS I to IV sample includes 14 period bounce CVs or candidates. We discuss the variability of CVs across the Hertzsprung–Russell diagram, highlighting selection biases of variability-based CV detection. Finally, we searched for, and found eight tertiary companions to the SDSS CVs. We anticipate that this catalogue and the extensive material included in the Supplementary Data will be useful for a range of observational population studies of CVs.

79 ASTRONOMY AND ASTROPHYSICS↗

Explainable AI classification for parton density theory

Quantitatively connecting properties of parton distribution functions (PDFs, or parton densities) to the theoretical assumptions made within the QCD analyses which produce them has been a longstanding problem in HEP phenomenology. To confront this challenge, we introduce an ML-based explainability framework, XAI4PDF, to classify PDFs by parton flavor or underlying theoretical model using ResNet-like neural networks (NNs). By leveraging the differentiable nature of ResNet models, this approach deploys guided backpropagation to dissect relevant features of fitted PDFs, identifying x-dependent signatures of PDFs important to the ML model classifications. By applying our framework, we are able to sort PDFs according to the analysis which produced them while constructing quantitative, human-readable maps locating the x regions most affected by the internal theory assumptions going into each analysis. This technique expands the toolkit available to PDF analysis and adjacent particle phenomenology while pointing to promising generalizations.

Artificial Intelligence↗

Adaptive Sampling for In Situ Cloud Probe (Final Report)

Clouds play a leading role in the Earth's global energy and solar radiation balance and hydrological cycle. Improving cloud models requires detailed information on the cloud microphysical properties, such as droplet size distribution and number density, liquid water content and cloud composition (droplets, ice particles), which can only be provided by aerial in situ measurements. However, for many atmospheric measurement instruments, the lack of flexibility in selecting the operational mode during operation can lead to uncertainties in sampling and measurement characteristics under continuously varying atmospheric conditions. This SBIR project is developing an advanced, compact optical imaging technology for in situ characterization of cloud hydrometeors. The development involves a deep modification of the existing Mesa Photonics’ Cloud Droplet Measurement System (CDMS) in order to implement real-time automatic adaptive sampling based on the acquired in situ data and environmental parameters. The new system, CDMS-2, implements two measurement modes: side-scatter imaging for smaller hydrometeors and direct bright-field-illumination imaging for larger hydrometeors in a significantly larger sample volume. The system measures the droplet size distribution (DSD) and number density with an added capability of discriminating between liquid water and ice hydrometeors (based on polarization-resolved side-scatter imaging). The instrument will implement automatic switching or alternating between the regular side-scatter imaging mode and sparse/large hydrometeor mode (based on the acquired data). Other adaptive sampling capabilities include variable sample volume and dynamic range (based on the measured DSD). The preferred deployment platforms are uncrewed aircraft systems (UAS) and tethered balloon/kite systems (TBS). The Phase I project achieved (or exceeded) the goals listed in the Work Plan. A CDMS-2 laboratory prototype implementing the polarization-resolved side-scatter imaging mode and direct bright-field-illumination imaging mode was designed and built. Additional capabilities included the variable illumination pulse energy and sample volume. The smallest detectable droplet diameter was improved to 3–4 μm (from the nominal 10 μm value specified for the original CDMS). Discrimination between water droplets and ice particles was experimentally demonstrated. The Phase I prototype was extensively tested and calibrated in the laboratory and also tested in the Pi Cloud Chamber at Michigan Technological University (MTU). The two intensive experimental campaigns at MTU provided unique opportunities of testing the CDMS-2 laboratory prototype under realistic warm and mixed-phase cloud conditions (stable for long periods of time), testing different sampling modes and intercomparing the CDMS-2 prototype to other co-located cloud characterization instruments. The Phase I project successfully demonstrated the feasibility of the proposed technology and identified the engineering challenges of designing a field deployable prototype instrument in Phase II. The Phase I study provides a solid basis for development, characterization and field-testing of the proposed advanced cloud probe with adaptive sampling in Phase II followed by commercialization of the technology in Phase III.

47 OTHER INSTRUMENTATION↗

Atomistic Insight into Cation Effects on Binding Energies in Cu-Catalyzed Carbon Dioxide Reduction

Electrolyte properties in general and electrolyte ions in particular have been shown to have significant effects on the kinetics of electrochemical CO 2 and CO reduction at metal electrodes, but these effects have not yet been fully understood. Here, we investigate the effects of cations at the aqueous electrolyte–Cu(211) interface on adsorbate binding energies and the electrolyte structure using density functional theory (DFT). Charging the interface via explicit Na + has systematic effects on adsorbate–electrolyte interactions and conformations. We describe specific local adsorbate–ion interactions, including direct alkali ion–adsorbate coordination and hydrogen bonding via ion-coordinated water molecules. The relative importance of these specific interactions and purely electrostatic field–adsorbate interactions is investigated by comparing the DFT-calculated ion effects to those predicted by purely electrostatic models of the interface and the adsorbates. We find that the trend in ion effects among different adsorbates at constant surface charge density is well explained by a purely electrostatic interaction model. The binding energy of OCCHO is found to depend strongly on the surface charge density as well as the spatial distribution of ions at constant surface charge density. These effects are also explained by a purely electrostatic local electrostatic field–adsorbate model. This indicates that alkali ion effects can be mainly attributed to purely electrostatic field interactions and that the local field at ion-stabilized active sites can depend significantly on both the overall charge density and on the spatial distribution of ions at constant charge density. This work provides new insight on alkali ion effects on a variety of adsorbates relevant to the CO 2 and CO reduction reactions by describing specific local ion–adsorbate interactions, systematic changes in adsorbate–electrolyte interactions induced by explicit surface charge, and comparisons of the calculated effects on binding energies to simple electrostatic field–adsorbate models to explain trends among various adsorbates and at varying surface charge densities and supercell configurations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental characterization of discharge plasma dynamics in a square capillary for prospective applications in laser wakefield acceleration

Research on plasma generated by electrical discharge in gas-filled capillaries plays an important role in advancing laser wakefield accelerators. We present a comprehensive study of the temporal and spatial distribution of plasma density within a short, 1 cm long, square-shaped capillary filled with hydrogen gas. Two transverse capillary sizes, 500 μm and 300 μm, were investigated achieving peak plasma densities of 0.9 × 10 18 cm −3 and 2.5 × 10 18 cm −3 , respectively, at the capillary center. In addition, we explore how these distributions depend on discharge parameters, specifically discharge current and gas flow through the capillary. Plasma density was determined by analyzing the Stark broadening of the hydrogen Balmer-alpha line. The results obtained were compared with the theoretical models and simulations. The comparison between the model predictions and the experimental data at the transient ionization stage of the discharge reveals a discrepancy of a factor of ∼1.3–2.1, depending on the capillary size, which is thoroughly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗