Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INCLINATION”

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 757 records · Page 42

DISCOVERY OF AN EDGE-ON DEBRIS DISK WITH A DUST RING AND AN OUTER DISK WING-TILT ASYMMETRY

Using Very Large Telescope/SPHERE near-infrared dual-band imaging and integral field spectroscopy, we discovered an edge-on debris disk around the 17 Myr old A-type member of the Scorpius–Centaurus OB association HD 110058. The edge-on disk can be traced to about 0.″6 or 65 AU projected separation. In its northern and southern wings, the disk shows at all wavelengths two prominent, bright, and symmetrically placed knots at 0.″3 or 32 AU from the star. We interpret these knots as a ring of planetesimals whose collisions may produce most of the dust observed in the disk. We find no evidence for a bow in the disk, but we identify a pair of symmetric, hooklike features in both wings. Based on similar features in the Beta Pictoris disk, we propose that this wing-tilt asymmetry traces either an outer planetesimal belt that is inclined with respect to the disk midplane or radiation-pressure-driven dust blown out from a yet unseen inner belt that is inclined with respect to the disk midplane. The misaligned inner or outer disk may be a result of interaction with a yet unseen planet. Overall, the disk geometry resembles the nearby disk around Beta Pictoris, albeit seen at smaller radial scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Astrometric confirmation and preliminary orbital parameters of the young exoplanet 51 Eridani b with the Gemini Planet Imager

We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b that provide further evidence that the companion is physically associated with 51 Eridani. Combining this new astrometric measurement with those reported in the literature, we significantly reduce the posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance alignment with 51 Eridani to 2 × 10 –7 , an order of magnitude lower than previously reported. If 51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo technique, preliminary constraints are placed on the orbital parameters of the system. The current set of astrometric measurements suggest an orbital semimajor axis of ${14}_{-3}^{+7}$ AU, corresponding to a period of ${41}_{-12}^{+35}$ years (assuming a mass of 1.75 M⊙ for the central star), and an inclination of ${138}_{-13}^{+15}$ deg. The remaining orbital elements are only marginally constrained by the current measurements. As a result, these preliminary values suggest an orbit that does not share the same inclination as the orbit of the distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani.

79 ASTRONOMY AND ASTROPHYSICS↗

SDSS-IV MaNGA: drivers of stellar metallicity in nearby galaxies

The distribution of stellar metallicities within and across galaxies is an excellent relic of the chemical evolution across cosmic time. Here we present a detailed analysis of spatially resolved stellar populations based on >2.6 million spatial bins from 7439 nearby galaxies in the Sloan Digital Sky Survey-IV (SDSS-IV) Mapping Nearby Galaxies at APO (MaNGA) survey. To account for accurate inclination corrections, we derive an equation for morphology-dependent determination of galaxy inclinations. Our study goes beyond the well-known global mass–metallicity relation and radial metallicity gradients by providing a statistically sound exploration of local relations between stellar metallicity [Z/H], stellar surface mass density Σ ⋆ , and galactocentric distance in the global mass–morphology plane. We find a significant resolved mass density–metallicity relation rΣ ⋆ ZR for galaxies of all types and masses above 10 9.8 M ⊙ ⁠. Different radial distances make an important contribution to the spread of the relation. Particularly, in low- and intermediate-mass galaxies, we find that at fixed Σ ⋆ metallicity increases with radius independently of morphology. For high masses, this radial dependence is only observed in high Σ ⋆ regions of spiral galaxies. This result calls for a driver of metallicity, in addition to Σ ⋆ that promotes chemical enrichment in the outer parts of galaxies more strongly than in the inner parts. We discuss gas accretion, outflows, recycling, and radial migration as possible scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterizing the magnetic fields of nearby molecular clouds using submillimeter polarization observations

Of all the factors that influence star formation, magnetic fields are perhaps the least well understood. The goal of this paper is to characterize the 3D magnetic field properties of nearby molecular clouds through various methods of statistically analysing maps of polarized dust emission. Our study focuses on nine clouds, with data taken from the Planck Sky Survey as well as data from the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry observations of Vela C. We compare the distributions of polarization fraction (p), dispersion in polarization angles ($\mathcal {S}$), and hydrogen column density (N H ) for each of our targeted clouds. To broaden the scope of our analysis, we compare the distributions of our clouds’ polarization observables with measurements from synthetic polarization maps generated from numerical simulations. We also use the distribution of polarization fraction measurements to estimate the inclination angle of each cloud’s cloud-scale magnetic field. We obtain a range of inclination angles associated with our clouds, varying from 16° to 69°. We establish inverse correlations between p and both $\mathcal {S}$ and NH in almost every cloud, but we are unable to establish a statistically robust $\mathcal {S}$ versus NH trend. By comparing the results of these different statistical analysis techniques, we are able to propose a more comprehensive view of each cloud’s 3D magnetic field properties. These detailed cloud analyses will be useful in the continued studies of cloud-scale magnetic fields and the ways in which they affect star formation within these molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗

TYC 2990-127-1: An Algol-type SB2 binary system of subgiant and red giant with a probable ongoing mass-transfer

ABSTRACT We present a study of the spectroscopic binary TYC 2990-127-1 from the LAMOST survey. We use full-spectrum fitting to derive radial velocities and spectral parameters. The high mass ratio indicates that the system underwent mass transfer in the past. We compute the orbital solution and find that it is a very close sub-giant/red giant pair on circular orbit, slightly inclined to the sky-plane. Fitting of the TESS photometrical data confirms this and suggests an inclination of i ∼ 39.8°. The light curve and spectrum around Hα show signs of irregular variability, which supports ongoing mass transfer. The binary evolution simulations suggest that the binary may experience non-conservative mass transfer with accretion efficiency 0.3, and the binary will enter into common envelope (CE) phase in the subsequent evolution. The remnant product after the ejection of CE may be a detached double helium white dwarf (He WD) or a merger.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring the mass of the black widow PSR J1555-2908

ABSTRACT Accurate measurements of the masses of neutron stars are necessary to test binary evolution models, and to constrain the neutron star equation of state. In pulsar binaries with no measurable post-Keplerian parameters, this requires an accurate estimate of the binary system’s inclination and the radial velocity of the companion star by other means than pulsar timing. In this paper, we present the results of a new method for measuring this radial velocity using the binary synthesis code Icarus. This method relies on constructing a model spectrum of a tidally distorted, irradiated star as viewed for a given binary configuration. This method is applied to optical spectra of the newly discovered black widow PSR J1555–2908. By modeling the optical spectroscopy alongside optical photometry, we find that the radial velocity of the companion star is 397 ± 4 km s−1 (errors quoted at 95 per cent confidence interval), as well as a binary inclination of >75°. Combined with γ-ray pulsation timing information, this gives a neutron star mass of 1.67$^{+0.15}_{-0.09}$ M⊙ and a companion mass of 0.060$^{+0.005}_{-0.003}$ M⊙, placing PSR J1555–2908 at the observed upper limit of what is considered a black widow system.

Kennedy, M. R.↗

CHANG-ES XXV: H i imaging of nearby edge-on galaxies – Data Release 4

ABSTRACT We present the ${\rm H}\, {\small I}$ distribution of galaxies from the Continuum Haloes in Nearby Galaxies – an EVLA Survey (CHANG-ES). Though the observational mode was not optimized for detecting ${\rm H}\, {\small I}$, we successfully produce ${\rm H}\, {\small I}$ cubes for 19 galaxies. The moment-0 maps from this work are available on CHANG-ES data release website (i.e. https://www.queensu.ca/changes). Our sample is dominated by star-forming, ${\rm H}\, {\small I}$-rich galaxies at distances from 6.27 to 34.1 Mpc. ${\rm H}\, {\small I}$ interferometric images on two of these galaxies (NGC 5792 and UGC 10288) are presented here for the first time, while 12 of our remaining sample galaxies now have better ${\rm H}\, {\small I}$ spatial resolutions and/or sensitivities of intensity maps than those in existing publications. We characterize the average scale heights of the ${\rm H}\, {\small I}$ distributions for a subset of most inclined galaxies (inclination > 80 deg), and compare them to the radio continuum intensity scale heights, which have been derived in a similar way. The two types of scale heights are well correlated, with similar dependence on disc radial extension and star formation rate surface density but different dependence on mass surface density. This result indicates that the vertical distribution of the two components may be governed by similar fundamental physics but with subtle differences.

79 ASTRONOMY AND ASTROPHYSICS↗

The LMC geometry and outer stellar populations from early DES data

The Dark Energy Camera has captured a large set of images as part of Science Verification (SV) for the Dark Energy Survey. The SV footprint covers a large portion of the outer Large Magellanic Cloud (LMC), providing photometry 1.5 magnitudes fainter than the main sequence turn-off of the oldest LMC stellar population. We derive geometrical and structural parameters for various stellar populations in the LMC disc. For the distribution of all LMC stars, we find an inclination of i = –38.14°±0.08° (near side in the North) and a position angle for the line of nodes of θ₀ = 129.51°±0.17°. We find that stars younger than ~4 Gyr are more centrally concentrated than older stars. Fitting a projected exponential disc shows that the scale radius of the old populations is R >4Gyr = 1.41 ± 0.01 kpc, while the younger population has R <4Gyr = 0.72 ± 0.01 kpc. However, the spatial distribution of the younger population deviates significantly from the projected exponential disc model. The distribution of old stars suggests a large truncation radius of R t = 13.5 ± 0.8 kpc. If this truncation is dominated by the tidal field of the Galaxy, we find that the LMC is ≃24 +9 –6 times less massive than the encircled Galactic mass. By measuring the Red Clump peak magnitude and comparing with the best-fit LMC disc model, we find that the LMC disc is warped and thicker in the outer regions north of the LMC centre. As a result, our findings may either be interpreted as a warped and flared disc in the LMC outskirts, or as evidence of a spheroidal halo component.

79 ASTRONOMY AND ASTROPHYSICS↗

Birefringence tests of gravity with multimessenger binaries

Extensions to General Relativity (GR) allow the polarization of gravitational waves (GW) from astrophysical sources to suffer from amplitude and velocity birefringence, which respectively induce changes in the ellipticity and orientation of the polarization tensor. Here, we introduce a multi-messenger approach to test this polarization behavior of GWs during their cosmological propagation using binary sources, for which the initial polarization is determined by the inclination and orientation angles of the orbital angular momentum vector with respect to the line of sight. In particular, we use spatially resolved radio imaging of the jet from a binary neutron star (BNS) merger to constrain the orientation angle and hence the emitted polarization orientation of the GW signal at the site of the merger, and compare to that observed on Earth by GW detectors. For GW170817, using past measurements of the inclination angle, we constrain the deviation from GR due to amplitude birefringence to κ A = $-0.12^{+0.60}_ {-0.61}$, while the velocity birefringence parameter κ V remains unconstrained. The inability to constrain κ V is due to the low amplitude of GW170817 in the Virgo detector, and measurements of the polarization orientation require information from a combination of multiple detectors with different alignments. For this reason, we also mock future BNS mergers with resolved afterglow proper motion and project that κ V could be constrained to a precision of 5 rad (corresponding to an angular shift of the GW polarization of δ φV ≈ 0.2 rad for a BNS at 100 Mpc) by a future network of third-generation ground-based GW detectors such as Cosmic Explorer and the radio High Sensitivity Array. Crucially, this velocity birefringence effect cannot be constrained with dark binary mergers as it requires polarization information at the emission time, which can be provided only by electromagnetic emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Microscopic Distributions of Defect Luminescence From Subgrain Boundaries in Multicrystalline Silicon Wafers

We investigate the microscopic distributions of sub-band-gap luminescence emission (the so-called D-lines D1/D2/D3/D4) and the band-to-band luminescence intensity, near recombination-active sub-grain boundaries in multicrystalline silicon wafers for solar cells. We find that the sub-band-gap luminescence from decorating defects/impurities (D1/D2) and from intrinsic dislocations (D3/D4) have distinctly different spatial distributions, and are asymmetric across the sub-grain boundaries. The presence of D1/D2 is correlated with a strong reduction in the band-to-band luminescence, indicating a higher recombination activity. In contrast, D3/D4 emissions are not strongly correlated with the band-to-band intensity. Based on spatially-resolved, synchrotron-based micro-X-ray fluorescence measurements of metal impurities, we confirm that high densities of metal impurities are present at locations with strong D1/D2 emission but low D3/D4 emission. Finally, we show that the observed asymmetry of the sub-band-gap luminescence across the sub-grain boundaries is due to their inclination below the wafer surface. Based on the luminescence asymmetries, the sub-grain boundaries are shown to share a common inclination locally, rather than be orientated randomly.

14 SOLAR ENERGY↗

Shear-based deformation processing of age-hardened aluminum alloy for single-step sheet production

Shear-based deformation processing by hybrid cutting-extrusion and free machining are used to make continuous strip, of thickness up to one millimeter, from low-workability AA6013-T6 in a single deformation step. The intense shear can impose effective strains as large as 2 in the strip without pre-heating of the workpiece. The creation of strip in a single step is facilitated by three factors inherent to the cutting deformation zone: highly confined shear deformation, in situ plastic deformation-induced heating and high hydrostatic pressure. The hybrid cutting-extrusion, which employs a second die located across from the primary cutting tool to constrain the chip geometry, is found to produce strip with smooth surfaces (S a < 0.4 μm) that is similar to cold-rolled strip. The strips show an elongated grain microstructure that is inclined to the strip surfaces – a shear texture – that is quite different from rolled sheet. Furthermore, this shear texture (inclination) angle is determined by the deformation path. Through control of the deformation parameters such as strain and temperature, a range of microstructures and strengths could be achieved in the strip. When the cutting-based deformation was done at room temperature, without workpiece pre-heating, the starting T6 material was further strengthened by as much as 30% in a single step. In elevated-temperature cutting-extrusion, dynamic recrystallization was observed, resulting in a refined grain size in the strip. Implications for deformation processing of age-hardenable Al alloys into sheet form, and microstructure control therein, are discussed.

36 MATERIALS SCIENCE↗

Poster - 23: Dosimetric Characterization and Transferability of an Accessory Mounted Mini-Beam Collimator

Objective: The dosimetric characterization of an accessory-mounted mini-beam collimator across three beam matched linear accelerators. Materials and Methods: Percent depth dose and profiles were measured for the open and mini-beam collimated fields. The average beam quality and peak-to-valley dose ratio (PVDR), the ratio of average peak dose to average valley dose, were obtained from these measurements. The open field relative output and the mini-beam collimator factor, the ratio of the mini-beam dose to open field dose at the beam center, were measured for square fields of side 2, 3, 4, and 5 cm. Mini-beam output as a function of collimator inclination angle relative to the central axis was also investigated. Results and Discussion: Beam quality for both the open and mini-beam collimated fields agreed across all linacs to within ±1.0%. The PVDR was found to vary by up to ±6.6% from the mean. For the 2, 3, and 4 cm fields the average open field relative output with respect to the 5 cm field was 0.874±0.4%, 0.921±0.3%, and 0.962±0.1%. The average collimator factors were 0.450±3.9%, 0.443±3.9%, 0.438±3.9%, and 0.434±3.9%. A decrease in collimator factor greater than 7% was found for an inclination angle change of 0.09°. Conclusion: The mini-beam collimator has revealed a difference between the three linacs not apparent in the open field data, yet transferability can still be attained through thorough dosimetric characterization.

07 ISOTOPE AND RADIATION SOURCES↗

Plans for Liquid Metal Divertor in Tokamak Compass

The COMPASS tokamak (R = 0.56 m, a = 0.2 m, B{sub T} = 1.3 T, I{sub p} ~ 300 kA, pulse duration 0.4 s) operates in ITER-like plasma shape in H-mode with Type-I ELMs. In 2019, we plan to install into the divertor a test target based on capillary porous system filled with liquid lithium/tin. This single target will be inclined toroidally in order to be exposed to ITER-relevant surface heat flux (20 MW/m{sup 2}). Based on precisely measured actual heat fluxes, our simulations predict (for 45° inclination, without accounting for the lithium vapor shielding) the surface temperature rises up to 700°C within 120 ms of the standard ELMy H-mode heat flux with ELM filaments reaching hundreds MW/m{sup 2}. Significant lithium vaporization is expected. The target surface will be observed by spectroscopy, fast visible and infrared cameras. The scientific program will be focused on operational issues (redeposition of the evaporated metal, ejection of droplets, if any) as well as on the effect on the plasma physics (improvement of plasma confinement, L–H power threshold, Z{sub eff}, etc.). After 2024, a closed liquid divertor may be installed into the planned COMPASS Upgrade tokamak (R = 0.84 m, a = 0.3 m, B{sub T} = 5 T, I{sub p} = 2 MA, P{sub in} = 8 MW, pulse duration ~2 s) with ITER-relevant heat fluxes loading the entire toroidal divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Light Water Reactor Sustainability Program Advanced Seismic Soil Structure Modeling

Risk calculations should focus on providing best estimate results, and associated insights, for evaluation and decision-making. Specifically, seismic probabilistic risk assessments (SPRAs) are intended to provide best estimates of the various combinations of structural and equipment failures that can lead to a seismic induced core damage event. However, in some instances the current SPRA approach has large uncertainties, and potentially masks other important events (for instance, it was not the seismic motions that caused the Fukushima core melt events, but the tsunami ingress into the facility). SPRA’s are performed by convolving the seismic hazard (this is the estimate of all likely damaging earthquakes at the site of interest) with the seismic fragility (the conditional probability of failure of a structure, system, or component given the occurrence of earthquake ground motion). In this calculation, there are three main pieces to seismic risk quantification, 1) seismic hazard and nuclear power plants (NPPs) response to the hazard, 2) fragility or capacity of structures, systems and components (SSC), and 3) systems analysis. Two areas where NLSSI effects may be important in SPRA calculations are, 1) when calculating in-structure response at the area of interest, and 2) calculation of seismic fragilities (current fragility calculations assume a lognormal distribution for probability of failure of components). Some important effects when using NLSSI in the SPRA calculation process include, 1) gapping and sliding, 2) inclined seismic waves coupled with gapping and sliding of foundations atop soil, 3) inclined seismic waves coupled with gapping and sliding of deeply embedded structures, 4) soil dilatancy, 5) soil liquefaction, 6) surface waves, 7) buoyancy, 8) concrete cracking and 9) seismic isolation The focus of the research task presented here-in is on implementation of NLSSI into the SPRA calculation process when calculating in-structure response at the area of interest. The specific nonlinear soil behavior included in the NLSSI calculation presented in this report is gapping and sliding. Other NLSSI effects are not included in the calculation. The results presented in this report document initial model runs in the linear and nonlinear analysis process. Final comparisons between traditional and advanced SPRA will be presented in the September 30th deliverable.

54 ENVIRONMENTAL SCIENCES↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗

Influence of the Fluctuating Velocity Field on the Surface Pressures in a Jet/Fin Interaction

The mechanism by which aerodynamic effects of jet/fin interaction arise from the flow structure of a jet in crossflow is explored using particle image velocimetry measurements of the crossplane velocity field as it impinges on a downstream fin instrumented with high-frequency pressure sensors. A Mach 3.7 jet issues into a Mach 0.8 crossflow from either a normal or inclined nozzle, and three lateral fin locations are tested. Conditional ensemble-averaged velocity fields are generated based upon the simultaneous pressure condition. Additional analysis relates instantaneous velocity vectors to pressure fluctuations. The pressure differential across the fin is driven by variations in the spanwise velocity component, which substitutes for the induced angle of attack on the fin. Pressure changes at the fin tip are strongly related to fluctuations in the streamwise velocity deficit, wherein lower pressure is associated with higher velocity and vice versa. The normal nozzle produces a counter-rotating vortex pair that passes above the fin, and pressure fluctuations are principally driven by the wall horseshoe vortex and the jet wake deficit. In conclusion, the inclined nozzle produces a vortex pair that impinges the fin and yields stronger pressure fluctuations driven more directly by turbulence originating from the jet mixing.

42 ENGINEERING↗

On the Fidelity of Computational Models for the Flow of Milled Loblolly Pine: A Benchmark Study on Continuum-Mechanics Models and Discrete-Particle Models

The upstream of bioenergy industry has suffered from unreliable operations of granular biomass feedstocks in handling equipment. Computational modeling, including continuum-mechanics models and discrete-particle models, offers insightful understandings and predictive capabilities on the flow of milled biomass and can assist equipment design and optimization. This paper presents a benchmark study on the fidelity of the continuum and discrete modeling approaches for predicting granular biomass flow. We first introduce the constitutive law of the continuum-mechanics model and the contact law of the coarse-grained discrete-particle model, with model parameters calibrated against laboratory characterization tests of the milled loblolly pine. Three classical granular material flow systems (i.e., a lab-scale rotating drum, a pilot-scale hopper, and a full-scale inclined plane) are then simulated using the two models with the same initial and boundary conditions as the physical experiments. The close agreement of the numerical predictions with the experimental measurements on the hopper mass flow rate, the hopper critical outlet width, the material stopping thickness on the inclined plane, and the dynamic angle of repose, clearly indicates that the two methods can capture the critical flow behavior of granular biomass. The qualitative comparison shows that the continuum-mechanics model outperforms in parameterization of materials and wall friction, and large-scale systems, while the discrete-particle model is more preferred for discontinuous flow systems at smaller scales. Industry stakeholders can use these findings as guidance for choosing appropriate numerical tools to model biomass material flow in part of the optimization of material handling equipment in biorefineries.

Jin, Wencheng↗