Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling methods”

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 559 records · Page 31

New signature method for identifying organic material in atmosphere and water samples

The levels of plastic contamination in the earth’s oceans have steadily increased over the past century, creating a threat to marine life and increasing climate change. Several methods have been developed over the past decade to try to determine the amount and accumulation rate of plastic impurities in the ecosystem. These methods have proven promising but suffer from long measurement times, bulky non-mobile systems, and sample degradation. We propose a smaller, faster system that employs an electron beam to examine seawater samples for polyethylene contamination levels. This process generates a full electron and photon energy spectrum analyses based on the subatomic particles generated from electron interactions. Our simulations have already been developed based on the known molecular composition of seawater and polyethylene, showing the possibility of discerning molecular plastics from microorganisms. With a system the size of a suitcase, this system would be easily deployable to remote locations around the world or installed as a monitoring system for manufacturing facilities. This concept could lead to a new methodology for the science of signatures and be applied to other biological hazard investigations, such as general virus detection.

43 PARTICLE ACCELERATORS↗

Facile fabrication of 2D material multilayers and vdW heterostructures with multimodal microscopy and AFM characterization

Reliable transfer processes that enable manipulation of two-dimensional (2D) materials, e.g., transition metal dichalcogenides (TMDCs) and MXenes, from one substrate to another has been a necessity for successful device fabrication. Here, with both mechanical exfoliation and chemical vapor deposition (CVD) widely used, a versatile, clean, deterministic, and yet simple transfer technique is highly needed. To address such need, we developed a transfer method that takes advantage of wettability contrast between interfaces without the use of sacrificial layers or chemical processes. More importantly, a setup was developed to carry out this transfer method with high sample selectivity and fine control of the position and orientation of transferred TMDC crystals, a feature required for fabrication of the devices based on vertical 2D heterostructures. Using both exfoliated and CVD grown materials and subsequent atomic force microscopy (AFM), photoluminescence (PL), confocal Raman and tip enhanced Raman spectroscopy (TERS) characterization, we ascertained the quality of interfaces resulting from the transfer process while preserving excellent 2D material integrity. PL and TERS maps revealed nanometer-scale heterogeneities in the interfaces of fabricated heterostructures, which should enable further perfection of the transfer technique. TERS/TEPL information were employed to identify areas suitable for nanodevice fabrication, making the reported transfer and characterization methods ideal for making high quality assembly of 2D heterostructure more accessible, which should facilitate exploration of vertical 2D heterostructures for applications in electronics, batteries, solar cells, and twistronics.

2D material transfer↗

SDSS-IV MaNGA: Identification and multiwavelength properties of Type-1 AGN in the DR15 sample

ABSTRACT We present a method to identify type-1 active galactic nuclei (AGNs) in the central 3 arcsec integrated spectra of galaxies in the MaNGA DR15 sample. It is based on flux ratios estimates in spectral bands flanking the expected h α broad component h αBC. The high signal-to-noise ratio obtained (mean S/N = 84) permits the identification of h αBC without prior subtraction of the host galaxy (HG) stellar component. A final sample of 47 type-1 AGN is reported out of 4700 galaxies at z < 0.15. The results were compared with those from other methods based on the SDSS DR7 and MaNGA data. Detection of type-1 AGN in those works compared to our method goes from 26 per cent to 81 per cent. Spectral indexes were used to classify the type-1 AGN spectra according to different levels of AGN-HG contribution, finding 9 AGN-dominated, 14 intermediate, and 24 HG-dominated objects. Complementary data in NIR-MIR allowed us to identify type I AGN-dominated objects as blue and HG-dominated as red in the WISE colours. From NVSS and FIRST radio continuum data, we identify 5 HERGs (high-excitation radio galaxies) and 4 LERGs (low-excitation radio galaxies), three showing evidence of radio-jets in the FIRST maps. Additional X-ray data from ROSAT allowed us to build [O iii] and h αBC versus X-ray, NIR-MIR, and radio continuum diagrams, showing that L(h αBC) and L([O iii]) provide good correlations. The range in h αBC luminosity is wide 38 < logL(h αBC) < 44, with log FWHM(h αBC) ∼ 3–4, covering a range of Eddington ratios of −5.15 < log Lbol/LEdd < 0.70. Finally, we also identify and report ten possible changing-look AGN candidates.

Cortes-Suárez, Edgar↗

Root system architecture and environmental flux analysis in mature crops using 3D root mesocosms

Current methods of root sampling typically only obtain small or incomplete sections of root systems and do not capture their true complexity. To facilitate the visualization and analysis of full-sized plant root systems in 3-dimensions, we developed customized mesocosm growth containers. While highly scalable, the design presented here uses an internal volume of 45 ft 3 (1.27 m 3 ), suitable for large crop and bioenergy grass root systems to grow largely unconstrained. Furthermore, they allow for the excavation and preservation of 3-dimensional root system architecture (RSA), and facilitate the collection of time-resolved subterranean environmental data. Sensor arrays monitoring matric potential, temperature and CO 2 levels are buried in a grid formation at various depths to assess environmental fluxes at regular intervals. Methods of 3D data visualization of fluxes were developed to allow for comparison with root system architectural traits. Following harvest, the recovered root system can be digitally reconstructed in 3D through photogrammetry, which is an inexpensive method requiring only an appropriate studio space and a digital camera. We developed a pipeline to extract features from the 3D point clouds, or from derived skeletons that include point cloud voxel number as a proxy for biomass, total root system length, volume, depth, convex hull volume and solidity as a function of depth. Ground-truthing these features with biomass measurements from manually dissected root systems showed a high correlation. We evaluated switchgrass, maize, and sorghum root systems to highlight the capability for species wide comparisons. We focused on two switchgrass ecotypes, upland (VS16) and lowland (WBC3), in identical environments to demonstrate widely different root system architectures that may be indicative of core differences in their rhizoeconomic foraging strategies. Finally, we imposed a strong physiological water stress and manipulated the growth medium to demonstrate whole root system plasticity in response to environmental stimuli. Hence, these new “3D Root Mesocosms” and accompanying computational analysis provides a new paradigm for study of mature crop systems and the environmental fluxes that shape them.

59 BASIC BIOLOGICAL SCIENCES↗

GMRES with embedded ensemble propagation for the efficient solution of parametric linear systems in uncertainty quantification of computational models

In a previous work, embedded ensemble propagation was proposed to improve the efficiency of sampling-based uncertainty quantification methods of computational models on emerging computational architectures. It consists of simultaneously evaluating the model for a subset of samples together, instead of evaluating them individually. A first approach introduced to solve parametric linear systems with ensemble propagation is ensemble reduction. In Krylov methods for example, this reduction consists in coupling the samples together using an inner product that sums the sample contributions. Ensemble reduction has the advantages of being able to use optimized implementations of BLAS functions and having a stopping criterion which involves only one scalar. However, the reduction potentially decreases the rate of convergence due to the gathering of the spectra of the samples. In this paper, we investigate a second approach: ensemble propagation without ensemble reduction in the case of GMRES. This second approach solves each sample simultaneously but independently to improve the convergence compared to ensemble reduction. This raises two new issues which are solved in this paper: the fact that optimized implementations of BLAS functions cannot be used anymore and that ensemble divergence, whereby individual samples within an ensemble must follow different code execution paths, can occur. We tackle those issues by implementing a high-performing ensemble GEMV and by using masks. The proposed ensemble GEMV leads to a similar cost per GMRES iteration for both approaches, i.e. with and without reduction. For illustration, we study the performances of the new linear solver in the context of a mesh tying problem. Furthermore, this example demonstrates improved ensemble propagation speed-up without reduction.

BLAS↗

Proceedings for the Workshop on Applied Nuclear Data Activities 2024

The Workshop for Applied Nuclear Data Activities (WANDA) is designed to increase communication among nuclear data (ND) users in multidisciplinary federal programs, ND producers, ND funders, and other ND experts. It also presents an opportunity to cross-pollinate ideas as well as introduce ND gaps identified by federal programs to ND experts and ND capabilities to the various federal ND users. WANDA 2024 included five technical sessions, three of which focused on Fusion Energy Sciences (FES)—FES Fusion Neutronics, FES Tritium Production, and FES Material Damage—and two stand-alone sessions—Isotopes and Targetry for Nuclear Data and Uncertainty Quantification. The FES sessions successfully brought new voices to the WANDA discussions, expanding the application space in which nuclear data are critical. FES programs need accurate nuclear data with realistic uncertainty quantification to properly estimate, for example, shielding, activation, tritium production, helium production, structural material integrity, and superconducting magnet operation. This includes a variety of projectile (neutrons, photons, charged particles) and target atoms. One of the action items common to all the FES sessions was a need to perform sensitivity studies to identify the prioritization of nuclear data needs. The Isotopes and Targetry session highlighted the many capabilities available to produce high-quality targets for nuclear data measurements, including 3D printing with spherical powders, combustion synthesis coupled with spin coating & electrospraying, inkjet printing, and isotopic doping. These new methods open doors for more accurate measurement, but it was also stressed that sample characterization following any method of fabrication is of the highest importance to accurately interpret nuclear data measurement results that used that sample. The Uncertainty Quantification (UQ) session was broken into two categories: nuclear data uncertainty quantification and the use of that uncertainty quantification. Thematic to the UQ session was the loss of information when going from nuclear data measurement, to evaluation, to evaluated file, and finally to neutron transport calculations. Current evaluated ND libraries typically only contain covariances, which assume that the probability distributions are Gaussian. Beyond being a simplified assumption for many evaluations, this can lead to negative values on many observables when attempting to sample the covariance. The covariance format, however, is very efficient in that a simple set of linear equations can transform uncertainty from parameters or cross sections to the application of interest. Focused collaboration is needed between nuclear data evaluators and nuclear data users to ensure that needs are being met.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bay Area Regional Energy: Network Integrated Commercial Retrofits (BRICR) Project. Final Report

The BRICR project applied large-scale building energy modeling concepts with the aim of reducing the cost of energy efficiency targeting, design, and project development, and measurement of energy savings for energy efficiency programs implemented by local governments that serve small and medium commercial buildings (SMB). The project leveraged the services and resources of existing local government energy programs serving disadvantaged and hard-to-reach SMB customers. In contrast to programs run by utilities, local government programs generally do not have direct access to energy billing records for an entire class of customers in a geographic area, which prior research demonstrated useful for large-scale building energy model baseline development and calibration. , However, local governments are rich in public records that offer important clues about physical attributes and uses that, along with behavior, determine energy use. Relying only on public records, BRICR demonstrated development of credible baseline energy models for 3,792 office, retail, and hotel buildings. Publicly disclosed annual energy use data from a local energy benchmarking program and anonymized data from the Building Performance Database, the nation’s largest dataset about energy-related characteristics of buildings, were utilized to validate and calibrate energy models via an innovative method comparing distributions of energy intensity by fuel type for portfolios of buildings of similar size, vintage, and use. Portfolio calibration does not provide certainty that an energy model fits an individual building; the method is useful when billing data is not accessible – a common situation for researchers, energy service providers and ESCOs, local governments, and any party other than a utility. A software component was developed, the BRICR gem, which automates simulation when relevant data is added or edited by the user to a file saved in the standardized BuildingSync XML schema for energy audit data. The component was demonstrated as a simplified means to generate a mass of energy models corresponding to public records containing basic attributes such as building scale, location, use, year built, and aspect ratio in combination with building energy code prototype data corresponding to use and vintage. The component was also demonstrated as a simplified means to automate energy simulation when attributes are revised; the intention was to enable iterative improvement of the baseline model and energy savings estimates for common energy conservation measures as users revise relevant attributes based on their observations. In the context of institutional change and uncertainty for the participating local government energy programs, 13 whole building retrofits were completed. Impacts were measured by applying the CalTRACK2.0 methods to standardize measurement of normalized metered energy consumption. The GRIDMeter methods of stratified sampling and individual load shape analysis were applied to adjust for impacts of the effect of COVID-19 on retrofitted buildings in the context of all local buildings of similar size and use. Excluding impacts of the pandemic, retrofitted buildings demonstrated between 1.6% and 25.1% reduction in energy use. The project contributed use cases and feedback that helped inform evolution of the software tools and data formats that were combined for the first time in the BRICR project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Attribute-Aware RBFs: Interactive Visualization of Time Series Particle Volumes Using RT Core Range Queries

Smoothed-particle hydrodynamics (SPH) is a mesh-free method used to simulate volumetric media in fluids, astrophysics, and solid mechanics. Visualizing these simulations is problematic because these datasets often contain millions, if not billions of particles carrying physical attributes and moving over time. Radial basis functions (RBFs) are used to model particles, and overlapping particles are interpolated to reconstruct a high-quality volumetric field; however, this interpolation process is expensive and makes interactive visualization difficult. Existing RBF interpolation schemes do not account for color-mapped attributes and are instead constrained to visualizing just the density field. To address these challenges, we exploit ray tracing cores in modern GPU architectures to accelerate scalar field reconstruction. We use a novel RBF interpolation scheme to integrate per-particle colors and densities, and leverage GPU-parallel tree construction and refitting to quickly update the tree as the simulation animates over time or when the user manipulates particle radii. We also propose a Hilbert reordering scheme to cluster particles together at the leaves of the tree to reduce tree memory consumption. Finally, we reduce the noise of volumetric shadows by adopting a spatially temporal blue noise sampling scheme. Our method can provide a more detailed and interactive view of these large, volumetric, time-series particle datasets than traditional methods, leading to new insights into these physics simulations.

Particle Volumes↗

Continuous polymer films deposited on top of porous substrates using plasma-enhanced atomic layer deposition and molecular layer deposition

Molecular layer deposition (MLD) is able to produce ultrathin polymer films with control over thickness, cross-linking, and chemical composition. With these capabilities, MLD should be useful in the fabrication of novel polymer membranes on porous supports. However, confining a continuous MLD film to the surface of porous substrates is difficult because of MLD film growth in the pores. The deposition in the pores lowers the conductance of the porous support. This paper presents a method to deposit continuous polymer films on top of porous substrates. In this method, Al2O3 plasma-enhanced atomic layer deposition (PE-ALD) using trimethylaluminum and oxygen plasma as the reactants was first used to cap the pores of the substrate. Subsequently, a polyamide MLD film was deposited on the Al2O3 PE-ALD capping layer using m-phenylenediamine and trimesoyl chloride as the reactants. The Al2O3 pore caps were then removed from the porous substrate by etching from the backside using a timed exposure to a dilute sodium hydroxide solution. This method was demonstrated using anodic aluminum oxide (AAO) and polyethersulfone (PES) porous substrates. Al2O3 PE-ALD film growth was limited to the top of the porous substrate, resulting in rapid surface recombination or high sticking coefficients for the reactive plasma species within the pores. Gas permeance measurements confirmed the pore capping of the AAO substrates. The reopening of the pores by dissolving the Al2O3 pore caps with a sodium hydroxide solution was monitored using gas permeance versus etch time. The removal of the Al2O3 pore caps from the PES substrates could also dissolve the Al2O3 layer underneath the MLD film. The loss of this Al2O3 layer led to the detachment of the MLD film from the PES substrate. However, the MLD film could be anchored to the PES support at fractures located in the Al2O3 film prior to the MLD. The Al2O3 film fracture allowed the MLD film to anchor firmly to the PES substrate by MLD in the pores of the PES porous substrate. The distance between the anchor points was a function of fracture density. This distance could be controlled by applying a tensile stress to the Al2O3 PE-ALD film to fracture the film through sample bending. This method produced firmly anchored polymer MLD films on top of the PES porous substrates.

Welch, Brian C. (ORCID:000000029747045X)↗

eDNAjoint: An R package for interpreting paired or semi‐paired environmental DNA and traditional survey data in a Bayesian framework

Abstract Environmental DNA (eDNA) sampling is increasingly used in surveys of species distribution as a potentially sensitive and efficient monitoring method. Yet access to modelling tools designed specifically for interpreting this new data type lags behind its ubiquity. While occupancy modelling software has dominated the analytical landscape for eDNA data analysis of single species, this type of model may not always be the most appropriate. The rate of eDNA detection often corresponds to species density, rather than just occupancy, and researchers often have access to observations from non‐genetic sampling methods at the same sites. To provide users access to a modelling framework designed to maximize the use of all available data, we developed an R package, eDNAjoint . The package provides an easy‐to‐use interface for fitting a ‘joint’ model that integrates data from paired or semi‐paired eDNA and traditional surveys in a Bayesian framework. The model can be used to estimate parameters like the probability of a false positive eDNA detection and mean catch rate at a site, and the package allows access to multiple model variations and Bayesian prior customization. Additional functionality can be used for model selection, summarising posteriors and comparing the relative sensitivities of the two survey methods. We demonstrate the use of eDNAjoint by fitting a variation of the model with site‐level covariates that scale the sensitivity of eDNA sampling relative to traditional sampling. The example workflow uses binary eDNA and seine count data for the endangered tidewater goby ( Eucyclogobius newberryi ) from a study by Schmelzle and Kinziger (2016). This use case includes a prior sensitivity analysis and an evaluation of the relationship between detection rates and environmental variables. eDNAjoint has the potential to greatly increase the range of users who will be able to rigorously analyse eDNA and traditional survey data in a Bayesian framework, understand if and how eDNA can improve monitoring practices, and gain confidence in the interpretability of eDNA data.

Keller, Abigail G. [Department of Environment Scie↗

Fast HARDI Uncertainty Quantification and Visualization with Spherical Sampling

In this paper, we study uncertainty quantification and visualization of orientation distribution functions (ODF), which corresponds to the diffusion profile of high angular resolution diffusion imaging (HARDI) data. The shape inclusion probability (SIP) function is the state‐of‐the‐art method for capturing the uncertainty of ODF ensembles. The current method of computing the SIP function with a volumetric basis exhibits high computational and memory costs, which can be a bottleneck to integrating uncertainty into HARDI visualization techniques and tools. We propose a novel spherical sampling framework for faster computation of the SIP function with lower memory usage and increased accuracy. In particular, we propose direct extraction of SIP isosurfaces, which represent confidence intervals indicating spatial uncertainty of HARDI glyphs, by performing spherical sampling of ODFs. Our spherical sampling approach requires much less sampling than the state‐of‐the‐art volume sampling method, thus providing significantly enhanced performance, scalability, and the ability to perform implicit ray tracing. Our experiments demonstrate that the SIP isosurfaces extracted with our spherical sampling approach can achieve up to 8164× speedup, 37282× memory reduction, and 50.2% less SIP isosurface error compared to the classical volume sampling approach. We demonstrate the efficacy of our methods through experiments on synthetic and human‐brain HARDI datasets.

97 MATHEMATICS AND COMPUTING↗

Strong parallel evidence of selection during switchgrass sward establishment in hybrid and lowland ecotypes

Switchgrass sward establishment results in up to 90% seedling mortality. The degree of selection during sward establishment has not been reported using modern genetic methods. Pooled leaf samples were sequenced from replicated swards of 46 half-sib families from two breeding groups (lowland and hybrid) before and through 3 years of stand establishment. Pooled allele frequencies were then assessed using fixation indices (Fst) and an independent data set was used to predict the polygenic impact of establishment selection on two traits (heading date and winter survivorship). Last, the DNA pools were assigned survival rankings to predict the sward survival genomically estimated breeding values within the training data set. Strong and parallel selection occured in both breeding groups. Five genomic regions exceeded the significant threshold of 99.9% in >10 families, indicating consistent selection across families and breeding groups. Polygenic trait predictions determined that establishment selection was partially associated with winter survivorship but resulted in variable heading date alterations. The genomewide variation is consistent with selection for a small number of related parental lines. This study observed strong selection for a small number of hybrid and coastal ecotype individuals which are promising germplasm sources for improved sward survival. This confirms prior reports of sward selection during grassland establishment and highlights the strength of pooled DNA sequencing for survival traits.

54 ENVIRONMENTAL SCIENCES↗

Continuous assessment of metabolic activity of mitochondria using resonance Raman microspectroscopy

Dysfunctional mitochondrial activity can lead to a variety of different diseases. As such, there exists a need to quantify changes in mitochondria function as it relates to these specific diseased states. Here, we present the use of resonance Raman (RR) spectroscopy as a tool to determine changes in isolated mitochondrial activity. RR spectroscopy, using 532 nm as the excitation source, specifically provides information on the reduction and oxidation (RedOx) state of cytochrome c, which is determined by the activity of protein complexes in the electron transport chain. In this model, injection of the substrate succinate into the mitochondrial sample is used to drive the electron transport chain, which causes a subsequent change in cytochrome c RedOx state. This change in RedOx state is tracked by RR spectroscopy. We report this tool gives real-time information on the rise and fall of the amount of reduced cytochrome c within the mitochondrial sample, providing a method for rapid assessment of mitochondrial metabolism that has broad applications in both basic science and medical research.

59 BASIC BIOLOGICAL SCIENCES↗

A stochastic biomass blending problem in decentralized supply chains

Blending biomass materials of different physical or chemical properties provides an opportunity to adjust the quality of the feedstock to meet the specifications of the conversion platform. We propose a model which identifies the right mix of biomass to optimize the performance of the thermochemical conversion process at the mini-mum cost. This is a chance-constraint programming (CCP) model which takes into account the stochastic nature of biomass quality. The proposed CCP model ensures that process requirements, which are impacted by physical and chemical properties of biomass, are met most of the time. We consider two problem settings, a centralized and a decentralized supply chain. We propose a mixed-integer linear program to model the blending problem in the centralized setting and a bilevel program to model the blending problem in the decentralized setting. We use the sample average approximation method to approximate the chance constraints, and propose solution algorithms to solve this approximation. We develop a case study for South Carolina using data provided by the Billion Ton Study. Based on our results, the blends identified consist mainly of pine and softwood residues. The blends identified and the suppliers selected by both models are different. The cost of the centralized supply chain is 2%–6% lower. The implications of these results are twofold. First, these results could lead to improved collaborations in the supply chain. Second, these results provide an estimate of the approximation error from assuming centralized decision making in the supply chain.

09 BIOMASS FUELS↗

Establishing the accuracy of position–specific carbon isotope analysis of propane by GC–pyrolysis–GC–IRMS

Rationale: Position-specific (PS) δ 13 C values of propane have proven their ability to provide valuable information on the evolution history of natural gases. Two major approaches to measure PS δ 13 C values of propane are isotopic 13 C nuclear magnetic resonance (NMR) and gas chromatography-pyrolysis-gas chromatography-isotope ratio mass spectrometry (GC-Py-GC-IRMS). Measurement accuracy of the isotopic 13 C NMR has been verified, but the requirements of large sample size and long experimental time limit its applications. GC-Py-GC-IRMS is a more versatile method with a small sample size, but its accuracy has not been demonstrated. Methods: We measured the PS δ 13 C values of propane from nine natural gases using both 13 C NMR and GC-Py-GC-IRMS, then evaluated the accuracy of the GC-Py-GC-IRMS method. Results: The results show that large carbon isotope fractionations occurred for both terminal and central carbons within propane during pyrolysis. The isotope fractionations during the pyrolysis are reproducible at optimum conditions, but vary between the two GC-Py-GC-IRMS systems tested, affected by experimental conditions (e.g., pyrolysis temperature, flow rate, and reactor conditions). Conclusions: It is necessary to evaluate and calibrate each GC-Py-GC-IRMS system using propane gases with accurately determined PS δ 13 C values. Finally, this study also highlights a need for PS isotope standards for propane and other molecules (e.g., butane and acetic acid).

03 NATURAL GAS↗

Visualization and efficient generation of constrained high-dimensional theoretical parameter spaces

We describe a set of novel methods for efficiently sampling high-dimensional parameter spaces of physical theories defined at high energies, but constrained by experimental measurements made at lower energies. Often, theoretical models such as supersymmetry are defined by many parameters, $\mathcal{O}$(10 – 100), expressed at high energies, while relevant experimental constraints are often defined at much lower energies, preventing them from directly ruling out portions of the space. Instead, the low-energy constraints define a complex, potentially non-contiguous subspace of the theory parameters. Naive scanning of the theory space for points which satisfy the low-energy constraints is hopelessly inefficient due to the high dimensionality, and the inverse problem is considered intractable. As a result, many theoretical spaces remain under-explored. We introduce a class of modified generative autoencoders, which attack this problem by mapping the high-dimensional parameter space to a structured low-dimensional latent space, allowing for easy visualization and efficient generation of theory points which satisfy experimental constraints. An extension without dimensional compression, which focuses on limiting potential information loss, is also introduced.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimization models for integrated biorefinery operations

Variations of physical and chemical characteristics of biomass lead to an uneven flow of biomass in a biorefinery, which reduces equipment utilization and increases operational costs. Uncertainty of biomass supply and high processing costs increase the risk of investing in the US’s cellulosic biofuel industry. We propose a stochastic programming model to streamline processes within a biorefinery. A chance constraint models system’s reliability requirement that the reactor is operating at a high utilization rate given uncertain biomass moisture content, particle size distribution, and equipment failure. The model identifies operating conditions of equipment and inventory level to maintain a continuous flow of biomass to the reactor. Furthermore, the sample average approximation method approximates the chance constraint and a bisection search-based heuristic solves this approximation. A case study is developed using real-life data collected at Idaho National Laboratory’s biomass processing facility. An extensive computational analysis indicates that sequencing of biomass bales based on moisture level, increasing storage capacity, and managing particle size distribution, increases utilization of the reactor and reduces operational costs.

09 BIOMASS FUELS↗