Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “software differences”

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 595 records · Page 33

Equivalent Point Sources for Simulating Nuclear/Radiological Search Operations

Modeling and simulation can be used to estimate the detection of different types and strengths of radioactive sources. Simulations of shielded gamma-ray sources in complex environments, such as searches in urban areas using vehicle-mounted detectors, can demonstrate what source strengths and shielding combinations should be able to be detected. For faster calculations, gamma-emitting sources that could have different shielding layers can be simulated ahead of time, resulting in equivalent point source emission distributions that can be used in larger simulation models. Larger-scale models will then run faster using equivalent point sources because the shielding layers around the sources have already been accounted for. Accompanying this paper are 260 equivalent point source emission files for 72 sources with different shielding layers, formatted for three different Monte Carlo codes. Users wanting to request versions of the SCALE and/or MCNP software can find details here: https://rsicc.ornl.gov/

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Design and Development of an Integrated Workstation Automation Hub

Miscellaneous Electronic Loads (MELs) account for one third of all electricity consumption in U.S. commercial buildings, and are drivers for a significant energy use in India. Many of the MEL-specific plug-load devices are concentrated at workstations in offices. The use of intelligence, and integrated controls and communications at the workstation for an Office Automation Hub – offers the opportunity to improve both energy efficiency and occupant comfort, along with services for Smart Grid operations. Software and hardware solutions are available from a wide array of vendors for the different components, but an integrated system with interoperable communications is yet to be developed and deployed. In this study, we propose system- and component-level specifications for the Office Automation Hub, their functions, and a prioritized list for the design of a proof-of-concept system. Leveraging the strength of both the U.S. and India technology sectors, this specification serves as a guide for researchers and industry in both countries to support the development, testing, and evaluation of a prototype product. Further evaluation of such integrated technologies for performance and cost is necessary to identify the potential to reduce energy consumptions in MELs and to improve occupant comfort.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Accelerating Science with the NERSC Burst Buffer Early User Program

NVRAM-based Burst Buffers are an important part of the emerging HPC storage landscape. The National Energy Research Scientific Computing Center (NERSC) at Lawrence Berkeley National Laboratory recently installed one of the first Burst Buffer systems as part of its new Cori supercomputer, collaborating with Cray on the development of the DataWarp software. NERSC has a diverse user base comprised of over 6500 users in 700 different projects spanning a wide variety of scientific computing applications. The use-cases of the Burst Buffer at NERSC are therefore also considerable and diverse. We describe here performance measurements and lessons learned from the Burst Buffer Early User Program at NERSC, which selected a number of research projects to gain early access to the Burst Buffer and exercise its capability to enable new scientific advancements. To the best of our knowledge this is the first time a Burst Buffer has been stressed at scale by diverse, real user workloads and therefore these lessons will be of considerable benefit to shaping the developing use of Burst Buffers at HPC centers.

97 MATHEMATICS AND COMPUTING↗

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗

Cross-instrument Analysis Correlation Software

This program has been designed to assist with the tracking of a sample from one analytical instrument to another such as SEM, microscopes, micro x-ray diffraction and other instruments where particular positions/locations on the sample are examined, photographed, etc. The software is designed to easily enter the position of fiducials and locations of interest such that in a future session in the same of different instrument the positions of interest can be re-found through using the known location fiducials in the current and reference session to transform the point into the current sessions coordinate system. The software is dialog box driven guiding the user through the necessary data entry and program choices. Information is stored in a series of text based extensible markup language (XML) files.

McJunkin, TimothyR.↗

Cross-instrument Analysis Correlation Software

This program has been designed to assist with the tracking of a sample from one analytical instrument to another such as SEM, microscopes, micro x-ray diffraction and other instruments where particular positions/locations on the sample are examined, photographed, etc. The software is designed to easily enter the position of fiducials and locations of interest such that in a future session in the same of different instrument the positions of interest can be re-found through using the known location fiducials in the current and reference session to transform the point into the current sessions coordinate system. The software is dialog box driven guiding the user through the necessary data entry and program choices. Information is stored in a series of text based extensible markup language (XML) files.

McJunkin, Timothy R.↗

Techno-Economic Evaluation of Cellulosic Ethanol Production Based on Pilot Biorefinery Data: a Case Study of Sweet Sorghum Bagasse Processed via L+SScF

Replacing fossil fuels with renewable fuels derived from lignocellulosic biomass can contribute to the mitigation of global warming and the economic development of rural communities. This will require lignocellulosic biofuels to become price competitive with fossil fuels. Techno-economic analyses can provide insights into which parts of the biofuel production process need to be optimized to reduce cost or energy use. Here, we used data obtained from a pilot biorefinery to model a commercial-scale biorefinery that processes lignocellulosic biomass to ethanol, with a focus on the minimum ethanol selling price (MESP). The process utilizes a phosphoric acid-catalyzed pre-treatment of sweet sorghum bagasse followed by liquefaction and simultaneous saccharification and co-fermentation (L+SScF) of hexose and pentose sugars by an engineered Escherichia coli strain. After validating a techno-economic model developed with the SuperPro Designer software for the conversion of sugarcane bagasse to ethanol by comparing it to a published Aspen Plus model, six different scenarios were modeled for sweet sorghum bagasse Under the most optimistic scenario, the ethanol can be produced at a cost close to the energy-equivalent price of gasoline. Aside from an increase in the price of gasoline, the gap between ethanol and gasoline prices could also be bridged by either a decrease in the cost of cellulolytic enzymes or development of value-added products from lignin.

Biofuel↗

Fitness cost of vancomycin-resistant Enterococcus faecium plasmids associated with hospital infection outbreaks

Abstract Background Vancomycin resistance is mostly associated with Enterococcus faecium due to Tn1546-vanA located on narrow- and broad-host plasmids of various families. This study’s aim was to analyse the effects of acquiring Tn1546-carrying plasmids with proven epidemicity in different bacterial host backgrounds. Methods Widespread Tn1546-carrying plasmids of different families RepA_N (n = 5), Inc18 (n = 4) and/or pHTβ (n = 1), and prototype plasmids RepA_N (pRUM) and Inc18 (pRE25, pIP501) were analysed. Plasmid transferability and fitness cost were assessed using E. faecium (GE1, 64/3) and Enterococcus faecalis (JH2-2/FA202/UV202) recipient strains. Growth curves (Bioscreen C) and Relative Growth Rates were obtained in the presence/absence of vancomycin. Plasmid stability was analysed (300 generations). WGS (Illumina-MiSeq) of non-evolved and evolved strains (GE1/64/3 transconjugants, n = 49) was performed. SNP calling (Breseq software) of non-evolved strains was used for comparison. Results All plasmids were successfully transferred to different E. faecium clonal backgrounds. Most Tn1546-carrying plasmids and Inc18 and RepA_N prototypes reduced host fitness (–2% to 18%) while the cost of Tn1546 expression varied according to the Tn1546-variant and the recipient strain (9%–49%). Stability of Tn1546-carrying plasmids was documented in all cases, often with loss of phenotypic resistance and/or partial plasmid deletions. SNPs and/or indels associated with essential bacterial functions were observed on the chromosome of evolved strains, some of them linked to increased fitness. Conclusions The stability of E. faecium Tn1546-carrying plasmids in the absence of selective pressure and the high intra-species conjugation rates might explain the persistence of vancomycin resistance in E. faecium populations despite the significant burden they might impose on bacterial host strains.

Tedim, Ana P.↗

SU-E-T-580: On the Significance of Model Based Dosimetry for Breast and Head and Neck 192Ir HDR Brachytherapy

Purpose: To assess the dosimetric and radiobiological differences between TG43-based and model-based dosimetry in the treatment planning of {sup 192}Ir HDR brachytherapy for breast and head and neck cancer. Methods: Two cohorts of 57 Accelerated Partial Breast Irradiation (APBI) and 22 head and neck (H&N) patients with oral cavity carcinoma were studied. Dosimetry for the treatment plans was performed using the TG43 algorithm of the Oncentra Brachy v4.4 treatment planning system (TPS). Corresponding Monte Carlo (MC) simulations were performed using MCNP6 with input files automatically prepared by the BrachyGuide software tool from DICOM RT plan data. TG43 and MC data were compared in terms of % dose differences, Dose Volume Histograms (DVHs) and related indices of clinical interest for the Planning Target Volume (PTV) and the Organs-At-Risk (OARs). A radiobiological analysis was also performed using the Equivalent Uniform Dose (EUD), mean survival fraction (S) and Tumor Control Probability (TCP) for the PTV, and the Normal Tissue Control Probability (N TCP) and the generalized EUD (gEUD) for the OARs. Significance testing of the observed differences performed using the Wilcoxon paired sample test. Results: Differences between TG43 and MC DVH indices, associated with the increased corresponding local % dose differences observed, were statistically significant. This is mainly attributed to their consistency however, since TG43 agrees closely with MC for the majority of DVH and radiobiological parameters in both patient cohorts. Differences varied considerably among patients only for the ipsilateral lung and ribs in the APBI cohort, with a strong correlation to target location. Conclusion: While the consistency and magnitude of differences in the majority of clinically relevant DVH indices imply that no change is needed in the treatment planning practice, individualized dosimetry improves accuracy and addresses instances of inter-patient variability observed. Research co-financed by the ESF and Greek funds through the Operational Program Education and Lifelong Learning Investing in Knowledge Society of the NSRF. Research Funding Program Aristeia. Nucletron, an Elekta company (Veenendaal, The Netherlands) is gratefully acknowledged for providing Oncentra Brachy v4.4 for research purposes.

61 RADIATION PROTECTION AND DOSIMETRY↗

EXAGRAPH: Graph and combinatorial methods for enabling exascale applications

Combinatorial algorithms in general and graph algorithms in particular play a critical enabling role in numerous scientific applications. However, the irregular memory access nature of these algorithms makes them one of the hardest algorithmic kernels to implement on parallel systems. With tens of billions of hardware threads and deep memory hierarchies, the exascale computing systems in particular pose extreme challenges in scaling graph algorithms. The codesign center on combinatorial algorithms, ExaGraph, was established to design and develop methods and techniques for efficient implementation of key combinatorial (graph) algorithms chosen from a diverse set of exascale applications. Algebraic and combinatorial methods have a complementary role in the advancement of computational science and engineering, including playing an enabling role on each other. In this paper, we survey the algorithmic and software development activities performed under the auspices of ExaGraph from both a combinatorial and an algebraic perspective. In particular, we detail our recent efforts in porting the algorithms to manycore accelerator (GPU) architectures. We also provide a brief survey of the applications that have benefited from the scalable implementations of different combinatorial algorithms to enable scientific discovery at scale. We believe that several applications will benefit from the algorithmic and software tools developed by the ExaGraph team.

97 MATHEMATICS AND COMPUTING↗

Methods, software and datasets to verify DVH calculations against analytical values: Twenty years late(r)

Purpose: The authors designed data, methods, and metrics that can serve as a standard, independent of any software package, to evaluate dose-volume histogram (DVH) calculation accuracy and detect limitations. The authors use simple geometrical objects at different orientations combined with dose grids of varying spatial resolution with linear 1D dose gradients; when combined, ground truth DVH curves can be calculated analytically in closed form to serve as the absolute standards. Methods: DICOM RT structure sets containing a small sphere, cylinder, and cone were created programmatically with axial plane spacing varying from 0.2 to 3 mm. Cylinders and cones were modeled in two different orientations with respect to the IEC 1217 Y axis. The contours were designed to stringently but methodically test voxelation methods required for DVH. Synthetic RT dose files were generated with 1D linear dose gradient and with grid resolution varying from 0.4 to 3 mm. Two commercial DVH algorithms—PINNACLE (Philips Radiation Oncology Systems) and PlanIQ (Sun Nuclear Corp.)—were tested against analytical values using custom, noncommercial analysis software. In Test 1, axial contour spacing was constant at 0.2 mm while dose grid resolution varied. In Tests 2 and 3, the dose grid resolution was matched to varying subsampled axial contours with spacing of 1, 2, and 3 mm, and difference analysis and metrics were employed: (1) histograms of the accuracy of various DVH parameters (total volume, D{sub max}, D{sub min}, and doses to % volume: D99, D95, D5, D1, D0.03 cm{sup 3}) and (2) volume errors extracted along the DVH curves were generated and summarized in tabular and graphical forms. Results: In Test 1, PINNACLE produced 52 deviations (15%) while PlanIQ produced 5 (1.5%). In Test 2, PINNACLE and PlanIQ differed from analytical by >3% in 93 (36%) and 18 (7%) times, respectively. Excluding D{sub min} and D{sub max} as least clinically relevant would result in 32 (15%) vs 5 (2%) scored deviations for PINNACLE vs PlanIQ in Test 1, while Test 2 would yield 53 (25%) vs 17 (8%). In Test 3, statistical analyses of volume errors extracted continuously along the curves show PINNACLE to have more errors and higher variability (relative to PlanIQ), primarily due to PINNACLE’s lack of sufficient 3D grid supersampling. Another major driver for PINNACLE errors is an inconsistency in implementation of the “end-capping”; the additional volume resulting from expanding superior and inferior contours halfway to the next slice is included in the total volume calculation, but dose voxels in this expanded volume are excluded from the DVH. PlanIQ had fewer deviations, and most were associated with a rotated cylinder modeled by rectangular axial contours; for coarser axial spacing, the limited number of cross-sectional rectangles hinders the ability to render the true structure volume. Conclusions: The method is applicable to any DVH-calculating software capable of importing DICOM RT structure set and dose objects (the authors’ examples are available for download). It includes a collection of tests that probe the design of the DVH algorithm, measure its accuracy, and identify failure modes. Merits and applicability of each test are discussed.

60 APPLIED LIFE SCIENCES↗

SU-E-T-464: Implementation and Validation of 4D Acuros XB Dose Calculations

Purpose: In this abstract we implement and validate a 4D VMAT Acuros XB dose calculation using Gafchromic film. Special attention is paid to the physical material assignment in the CT dataset and to reported dose to water and dose to medium. Methods: A QUASAR phantom with a 3 cm sinusoidal tumor motion and 5 second period was scanned using 4D computed tomography. A CT was also obtained of the static QUASAR phantom with the tumor at the central position. A VMAT plan was created on the average CT dataset and was delivered on a Varian TrueBeam linear accelerator. The trajectory log file from this treatment was acquired and used to create 10 VMAT subplans (one for each portion of the breathing cycle). Motion for each subplan was simulated by moving the beam isocentre in the superior/inferior direction in the Treatment Planning System on the static CT scan. The 10 plans were calculated (both dose to medium and dose to water) and summed for 1) the original HU values from the static CT scan and 2) the correct physical material assignment in the CT dataset. To acquire a breathing phase synchronized film measurements the trajectory log was used to create a VMAT delivery plan which includes dynamic couch motion using the Developer Mode. Three different treatment start phases were investigated (mid inhalation, full inhalation and full exhalation). Results: For each scenario the coronal dose distributions were measured using Gafchromic film and compared to the corresponding calculation with Film QA Pro Software using a Gamma test with a 3%/3mm distance to agreement criteria. Good agreement was found between calculation and measurement. No statistically significant difference in agreement was found between calculations to original HU values vs calculations to over-written (material-assigned) HU values. Conclusion: The investigated 4D dose calculation method agrees well with measurement.

60 APPLIED LIFE SCIENCES↗

Managing a Real-Time Embedded Linux Platform with Buildroot

Developers of real-time embedded software often need to build the operating system, kernel, tools and supporting applications from source to work with the differences in their hardware configuration. The first attempts to introduce Linux-based real-time embedded systems into the Fermilab accelerator controls system used this approach but it was found to be time-consuming, difficult to maintain and difficult to adapt to different hardware configurations. Buildroot is an open source build system with a menu-driven configuration tool (similar to the Linux kernel build system) that automates this process. A customized Buildroot [1] system has been developed for use in the Fermilab accelerator controls system that includes several hardware configuration profiles (including Intel, ARM and PowerPC) and packages for Fermilab support software. A bootable image file is produced containing the Linux kernel, shell and supporting software suite that varies from 3 to 20 megabytes large – ideal for network booting. The result is a platform that is easier to maintain and deploy in diverse hardware configurations

Diamond, J.↗

Supply Chain-Based Solution to Prevent Fuel Tax Evasion: Final Report

The primary source of funding for the United States transportation system is derived from motor fuel and other highway use taxes. Loss of revenue attributed to fuel-tax evasion has been assessed to be somewhere between 1 billion and 3 billion per year. Any solution that addresses this problem needs to include not only the tax-collection agencies and auditors, but also the carriers transporting oil products and the carriers customers. This report presents a system developed by the Oak Ridge National Laboratory (ORNL) for the Federal Highway Administration which has the potential to reduce or eliminate many fuel-tax evasion schemes. The solution balances the needs of tax-auditors and those of the fuel-hauling companies and their customers. The system has three main components. The on-board subsystem combined sensors, tracking and communication devices, and software (the on-board Evidential Reasoning System, or obERS) to detect, monitor, and geo-locate the transfer of fuel among different locations. The back office sub-system (boERS) used self-learning algorithms to determine the legitimacy of the fuel loading and offloading (important for tax auditors) and detect potential illicit operations such as fuel theft (important for carriers and their customers, and may justify the deployment costs). The third sub-system, the Fuel Distribution and Auditing System or FDAS, is a centralized database, which together with a user interface allows tax auditors to query the data submitted by the fuel-hauling companies and correlate different parameters to quickly identify any anomalies. Industry partners included Barger Transport of Weber City, Virginia (fleet); Air-Weigh, of Eugene, Oregon (and their wires and harnesses); Liquid Bulk Tank (LBT) of Omaha, Nebraska (three five-compartment trailers); and Innovative Software Engineering (ISE) of Coralville, Iowa(on-board telematics device and back-office system). ORNL conducted a pilot test with the three instrumented vehicles collecting real-world data during an eight-month period (October 2014 to June 2015). The vehicles logged a total of 375,000 miles and transported more than 7.5 million gallons of fuel. The drivers entered information on the telematics device about 77% of the time, with about 15% (10%) of the data showing at least one compartment with loading (offloading) valve activity, but not driver information entry for that compartment. Seven events (or about 0.7% of the trips) were identified by the boERS as having a missing fuel diversion number. The solution developed and tested in the pilot test had federal- and state-level tax auditors as its main audience. Nonetheless, the technology has to be adopted by the fuel carriers and therefore the solution had to address the needs of fuel-hauling companies and their customers (i.e., fuel theft and cocktailing). Sensors in the hatches allowed the obERS recognize when one of those were opened (always a suspicious activity, unless it happens at locations where maintenance is performed on the vehicle). The fuel-theft issue was addressed by a self-learning algorithm deployed on the boERS that continuously processed the data from the field to construct probability distributions of measures such as elapsed time of fuel loading and offloading by driver, vehicle, and compartment, valve actuation sequence, elapsed time between the first two valve actuations (by driver, compartment, and location), and other parameters. Probability thresholds, which can be set up by the carrier, determined how to classify the observed events. The boERS also kept track of valve sequencing at a given location and analyzed these actuations to help identify any suspicious activities. Technical and economic recommendations include: a) simplification of the driver data-entry task; b) reduction of system deployment cost; c) addition of capabilities that make the system more appealing to industry (i.e., identify and avoid fuel misdelivery); and d) incorporation of additional capabilities to the FDAS. Path to commercialization recommendations include: a) identify a company willing to further develop, test, and certify a hardened system with a price point the market will bear; b).transfer ERS to licensee; c) establish a fleet or fleets who want this technology for carrier benefits or trailer manufacturer who wants to offer it as optional technology; and d) identify a location for the deployment of FDAS.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Reinforcement Learning for Distribution Grid Optimization (PyCIGAR) v0.1

PyCIGAR is a python software package that merges off-the-shelf reinforcement learning libraries (RLLib and Ray) with electric power distribution system simulation tools (OpenDSS and a custom power flow solver built by LBL). PyCIGAR enables the training of neural networks to optimize the behavior of different components in the electric distribution grid, such as control systems in photovoltaic rooftop solar inverters and electric battery storage systems. The software package has been used to train neural networks to update settings in photovoltaic rooftop solar inverter control systems to mitigate cyber attacks on other solar photovoltaic rooftop devices.

Arnold, Daniel↗

PB-AM: An open-source, fully analytical linear poisson-boltzmann solver

Here, we present the open source distributed software package Poisson-Boltzmann Analytical Method (PB-AM), a fully analytical solution to the linearized PB equation, for molecules represented as non-overlapping spherical cavities. The PB-AM software package includes the generation of outputs files appropriate for visualization using visual molecular dynamics, a Brownian dynamics scheme that uses periodic boundary conditions to simulate dynamics, the ability to specify docking criteria, and offers two different kinetics schemes to evaluate biomolecular association rate constants. Given that PB-AM defines mutual polarization completely and accurately, it can be refactored as a many-body expansion to explore 2- and 3-body polarization. Additionally, the software has been integrated into the Adaptive Poisson-Boltzmann Solver (APBS) software package to make it more accessible to a larger group of scientists, educators, and students that are more familiar with the APBS framework.

97 MATHEMATICS AND COMPUTING↗

Non-invasive airway health measurement using synchrotron x-ray microscopy of high refractive index glass microbeads

Cystic fibrosis (CF) is caused by a gene defect that compromises the ability of the mucociliary transit (MCT) system to clear the airways of debris and pathogens. To directly characterise airway health and the effects of treatments we have developed a synchrotron X-ray microscopy method that non-invasively measures the local rate and patterns of MCT behaviour. Although the nasal airways of CF mice exhibit the CF pathophysiology, there is evidence that nasal MCT is not altered in CF mice1. The aim of this experiment was to determine if our non-invasive local airway health assessment method could identify differences in nasal MCT rate between normal and CF mice, information that is potentially lost in bulk MCT measurements. Experiments were performed on the BL20XU beamline at the SPring-8 Synchrotron in Japan. Mice were anaesthetized, a small quantity of micron-sized marker particles were delivered to the nose, and images of the nasal airways were acquired for 15 minutes. The nasal airways were treated with hypertonic saline or mannitol to increase surface hydration and MCT. Custom software was used to locate and track particles and calculate individual and bulk MCT rates. No statistically significant differences in MCT rate were found between normal and CF mouse nasal airways or between treatments. However, we hope that the improved sensitivity provided by this technique will accelerate the ability to identify useful CF lung disease-modifying interventions in small animal models, and enhance the development and efficacy of proposed new therapies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Enhancing the physical modeling capability of open-source MFIX-DEM software for handling particle size polydispersity: Implementation and validation

Multiphase flows are ubiquitous in many industrial processes. The inherent coupling of different phases poses many unique challenges in predicting and effectively controlling these processes. Hence, computational modeling and simulation offers a viable approach to overcome these challenges. Here in this study, we present recent development efforts for enhancing the physical modeling capabilities of an open-source computational modeling tool for real life industrial multiphase processes by enabling particle-size polydispersity and demonstrating with an associated validation study. The proposed implementation was performed in MFIX open-source framework due to its unique feature of tightly integrated computational fluid dynamics and discrete element method solvers for simulating coupled continuum fluid and granular flows. We have implemented the polydispersity feature in a minimally invasive way and provided means to allow easy specification of an arbitrary particle size distribution function, which also enables the user to easily handle an arbitrary number of solid phases, each possessing a distinct arbitrary particle-size distribution. To establish the credibility of improvements, we have carried out a preliminary verification and validation (V&V) study for the polydispersity feature by employing a hopper bin discharge problem, which is frequently encountered in industrial applications. Specifically, two types of micro-glass beads with distinct size distributions are used to fill the hopper in two possible packing arrangements, i.e., well-mixed and layered configurations, with varying mass (particle number) ratios. The experimentally obtained discharge dynamics (e.g., normalized discharge mass fraction for one of the phases versus the overall discharge mass fraction) for different systems is found to be in excellent agreement with the corresponding simulation results.

42 ENGINEERING↗