Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reproducibility”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Mapping TES Temperature Sensitivity and Current Sensitivity as a Function of Temperature, Current, and Magnetic Field with IV Curve and Complex Admittance Measurements

We have specialized astronomical applications for X-ray microcalorimeters with superconducting transition edge sensors (TESs) that require exceptionally good TES performance, but which operate in the small-signal regime. We have therefore begun a program to carefully characterize the entire transition surface of TESs with and without the usual zebra stripes to see if there are reproducible local “sweet spots” where the performance is much better than average. These measurements require precise knowledge of the circuit parameters. Here, we show how the Shapiro effect can be used to precisely calibrate the value of the shunt resistor. We are also investigating the effects of stress and external magnetic fields to better understand reproducibility problems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Host-specific and tissue-dependent orchestration of microbiome community structure in traditional rice paddy ecosystems

Abstract Background and aim Rice and duckweed are two monocotyledonous plants that naturally coexist in paddy fields. While the presence of duckweed in paddy fields significantly improves rice productivity, the interplay between soil microbes and the two plant hosts in this agroecosystem remains unexplored. Methods We compared the bacterial community structure between duckweed, rice and soil from multiple rice paddies. We also isolated bacteria from these communities and characterized their modes of bacterial colonization and plant growth-promotion using model plants. Results Our data indicate that host-specific and tissue-dependent factors reproducibly orchestrate the bacterial community structure associated with their plant hosts. This is corroborated by results from culture-dependent approaches in which the dominant genus Pantoea isolated from rice aerial tissues can strongly attach and colonize gnotobiotic duckweed in spite of the low representation of this bacterial genus in the natural duckweed microbiome. Our study identified a core of 254 bacterial taxa that are consistently found in all rice and duckweed tissue samples from rice paddy fields. Furthermore, characterizing auxin-producing bacteria isolates from both plant species identified potential plant growth-promoting bacteria that may improve growth for both duckweed and rice in paddy fields. Conclusions Results from this work provide evidence for the importance of the host tissue and species context in determining plant colonization by microbes in the paddy field system. The resources generated in this study could facilitate the agronomic deployment of microbes for more sustainable rice production.

Agriculture↗
DOE OSTI2019

Genome-wide association studies of bark texture in Populus trichocarpa

Plants have developed bark as a defense barrier to deal with environmental stresses, such as pathogen invasion, drought, and UV radiation. Bark is composed of dead differentiated cells and is formed when radial expansion pushes the cortex and epidermis outward creating secondary meristem. However, the genetic control of this complex phenotype is generally unknown. Here, we use association mapping to define the genomic regions associated with natural variation in bark texture (BT) in Populus trichocarpa. Clonally replicated provenance trials of P. trichocarpa were studied for BT collected across three sites, multiple years, and (2–3) biological replicates per site. Forty-one genomic intervals containing SNPs significantly associated with BT were detected that were highly reproducible across sites, years, and replicates. A list of candidate genes within these regions with related putative function was identified. A total of 98 genes were considered candidate genes due to significance and putative function with connection to the phenotype. Association mapping using low coverage sequencing allowed us to detect narrow genomic intervals (1–20 kb) with high reproducibility and shared candidate genes. For example, a membrane-associated apoptosis protein and a wall-associated receptor kinase (PR5K-like) protein, which both are involved in radial growth and tissue differentiation, lie within significant trait-associated region. Here, two copies of root hair defective3 genes, transmembrane protein in Populus differentiating xylem and phloem, were also significantly associated with BT and co-located with major quantitative trait locus (QTL) for BT, bark thickness, and diameter from our previous study suggesting an important role in radial growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

How variable are the volumetric measurements from gated perfusion SPECT when a one-day stress-rest protocol is used?

Background: Using myocardial perfusion scintigraphy (MPS), an increase in left ventricular (LV) volumes or a decrease in ejection fraction (EF) from rest to stress may be clinically important. The variation in these measures between the low-dose stress acquisition and high-dose rest acquisition in a one-day stress-rest protocol has not been established. We assessed the reproducibility of gated volumetric indices between stress and rest and the normal variation in ungated TID ratio for a one-day stress-rest {sup 99m}Tc-tetrofosmin protocol. Methods: Two thousand and one hundred and fifty eight (2158) {sup 99m}Tc-tetrofosmin MPS patient studies were analyzed retrospectively. Studies were excluded for incomplete data, significant technical difficulties, or (for gated analysis but not for analysis of TID ratio) if the LV EF was > 75%. An analysis of gated data was undertaken to establish the reproducibility of ventricular volumes and EF between stress and rest scans. Ungated volume data were analyzed to determine the confidence limits of TID ratio according to ventricular volume. Results: Gated data were analyzed for 621 patients without inducible hypoperfusion. Mean EF at rest was slightly higher than after stress (62.4% ± 10.3% vs 61.2% ± 10.4%, P < 0.001), and the standard deviation of the difference was 5.2% (95% CI 4.9% to 5.5%). Ungated volumes were available for 992 non-ischaemic patients. The upper 95% CI for TID ratio was 1.23. This increased from 1.20 to 1.37 between the highest and lowest deciles of rest ungated volume. Conclusion: Using a one-day stress-rest {sup 99m}Tc-tetrofosmin protocol, a fall in LV EF between rest and stress of > 11.6% or a TID ratio of > 1.23 is likely to be clinically reliable. The upper limit of normal for TID ratio needs to be increased for patients with small LV chamber volumes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Capability of CAM5.1 in simulating maximum air temperature patterns over West Africa during boreal spring

This study classifies maximum air temperature patterns over West Africa into six groups and evaluates the capability of a global climate model (Community Atmospheric Model version 5.1; CAM) to simulate them. We analyzed 45-year (1961–2005) multi-ensemble (50 members) simulations from CAM and compared the results with those of the Climate Research Unit (CRU) and the twentieth Century Reanalysis data sets. Using Self Organizing Map algorithm to classify the spatial patterns of maximum air temperature during boreal spring, the study reveals the temperature patterns that CAM can simulate well and those the model struggles to reproduce. The results show that the best agreements between the composites of observation and CAM occur in the first temperature pattern group (which features positive temperatures anomalies over the Sahel) and Node 2 (which features near-normal temperature) pattern of the third group. CAM succeeded in reproducing some of the associated regional atmospheric dynamics and thermodynamic features in winds (horizontal and vertical), temperature fields, the cloud fractions, and the mean sea-level pressure. Although CAM struggles to capture the relationship between air temperature patterns and tele-connection indices during the boreal spring season over West Africa, it agrees with observations that temperature patterns over the sub-region cannot be associated with a single climate index. An ensemble member (SIM48) captures the inter-annual variation of the observed temperaure patterns with high sycronization (η > 44%), much better than that of ensembles mean (η < 30%). SIM48 also captures adequately four of the spatial patterns in comparison to three captured by the ensembles mean. This indicates that, for better seasonal forecasts and more reliable future climate projections, the practice whereby an ensemble mean is based on uniformly averaging the members rather than the performance of individual ensemble members needs to be reviewed. Finally, the results of the study may be used to improve the perfomance of CAM over West Africa, thereby strengthening the on-going efforts to include CAM as part of multi-model forecasting system over West Africa.

Lawal, Kamoru A. (ORCID:0000000281988844)↗
DOE OSTI2021

Propagation of a Stress Pulse in a Heterogeneous Elastic Bar

The propagation of a wave pulse due to low-speed impact on a one-dimensional, heterogeneous bar is studied here. Due to the dispersive character of the medium, the pulse attenuates as it propagates. This attenuation is studied over propagation distances that are much longer than the size of the microstructure. A homogenized peridynamic material model can be calibrated to reproduce the attenuation and spreading of the wave. The calibration consists of matching the dispersion curve for the heterogeneous material near the limit of long wavelengths. It is demonstrated that the peridynamic method reproduces the attenuation of wave pulses predicted by an exact microstructural model over large propagation distances.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase diagram of magnetic shape memory alloy Ni 50 Mn $50–x$ In $x$ , 0 < $x$ , 25 from first principles, via spin cluster expansion and phonon vibrational entropies

The metamagnetic shape memory Heusler alloy Ni 50 Mn $50–x$ In $x$ exhibits a rich phase diagram featuring competing magnetic states, coupled magnetic–structural phase transitions, and strong compositional sensitivity. Existing first-principles approaches struggletocapturetheintertwinedchemical, magnetic, andvibrationaleffectsinthesealloys, necessitating a more integrated modeling framework. We develop a spin cluster expansion (spin-CE) framework augmented by a quasi-harmonic phonon model to capture both configurational (chemical and magnetic) and vibrational contributions to the free energy of Ni 50 Mn $50–x$ In $x$ over the full compositional range 0 ≤ x ≤25. The spin-CE includes both chemical clusters and composition-dependent Ising spin interactions, with parameters fit to a first-principles density functional theory (DFT) dataset. Using this approach, we predict the complete magnetostructural phase diagram and transformation temperatures of Ni 50 Mn $50–x$ In $x$ across the composition space. We find that vibrational entropy alone is insufficient to reproduce the martensitic transformation in the magnetic shape memory alloy regime, highlighting the essential role of magnetism. Incorporating both magnetic and vibrational contributions allows us to reproduce all experimentally known phases, including the disappearance of the stable martensite phase at a critical In concentration and the Curie temperature of the austenite phase. The method also captures the transition with increasing In in martensite from antiferromagnetic to ferromagnetic order and predicts re-entrant ferromagnetism, though the latter occurs at higher In content than reported experimentally. We discuss possible sources of this discrepancy and highlight the broader applicability of the method to other magnetostructurally complex materials, where it may offer mechanistic insight and predictive design capabilities.

Cluster expansion↗
DOE OSTI2024

Grain structure and texture selection regimes in metal powder bed fusion

Additive manufacturing (AM) offers opportunities to produce complex part geometries not possible with conventional processing and in some cases even improve part performance. However, adoption has been slowed by difficulties assessing microstructure variability and there is no straightforward approach to relate processing to grain structure characteristics. In this study, datasets from AdditiveFOAM heat transport simulations of laser powder bed fusion (LPBF) are used to drive ExaCA simulations of grain structure. The GPU utilization of ExaCA and an algorithmic update for modeling melt pool overlap region solidification enabled rapid and parallel simulation across a wider range of process conditions than previously explored with cellular automata-based solidification models. A texture selection angle $θ_s$ is defined based on melt pool overlap geometry, and the range of $θ_s$ over which a commonly observed texture transition occurs in characterized AM builds was well-reproduced by ExaCA simulations over a wide range of melt pool shape, hatch spacing, and layer height. ExaCA simulations with 90 degree rotation of the scan direction on every other layer reproduced a number of trends from the AM literature including grain refinement, the dominance of layers with larger melt pools on the final grain structure, and the weakening or strengthening of texture depending on odd and even layer melt pool overlap geometry. EBSD data from a benchmark AM part is used to validate the simulated mechanism of a layer rotation-induced texture strengthening effect. Importantly, these results expand the understanding of the mechanisms for texture selection in alloys with cubic crystal symmetry and offer an approach to easily evaluate processing conditions. With this new understanding, these modeling tools will enable anticipation of previously unexpected variations in grain structure and target specific microstructures and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Results of an interlaboratory study on the working curve in vat photopolymerization II: Towards a standardized method

The working curve measurement in photopolymer additive manufacturing is a ubiquitous measure of the cure depth of a printing resin as a function of radiant exposure of light. The fit parameters from this measurement (the depth of light penetration D p and the critical exposure E c ) are used to evaluate and report a resin’s printability, optimize processing parameters, and inform print and resin quality control. Despite its widespread use, the working curve lacks a standard measurement method. Here, following up on our paper “Results of an Interlaboratory Study on the Working Curve in Vat Photopolymerization” from last year, an interlaboratory study on the working curve was performed using calibrated, reproducible, bandpass filtered light sources. With these light sources, the variability between labs in measured working curves was dramatically reduced from the initial interlaboratory study. Aggregate data from this experiment produced reliable D p and E c measurements at 385 nm of 39.2 ± 3.7 µm and 12.3 ± 3.0 mJ cm −2 , respectively. At 405 nm the values of D p and E c are 69.3 ± 3.8 µm and 17.9 ± 2.3 mJ cm −2 , respectively. The results are agnostic to the thickness measurement tool utilized by participants, ensuring broad applicability across laboratories. We also tested the generalizability of the proposed method of using a filtered light source by filtering a commercial 405 nm light source and obtaining a working curve in agreement with the aggregate data from the interlaboratory study. This interlaboratory study provides a basis for a documentary standard for the working curve, so that the entire photopolymer additive manufacturing industry can share reproducible and interoperable working curve data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Noise-aware optimization in nominally identical manufacturing and measuring systems for high-throughput parallel workflows

Device-to-device variability in experimental noise critically impacts reproducibility, especially in automated, high-throughput systems like additive manufacturing farms. While manageable in small labs, such variability can escalate into serious risks at larger scales, such as architectural 3D printing, where noise may cause structural or economic failures. This contribution presents a noise-aware decision-making algorithm that quantifies and models device-specific noise profiles to manage variability adaptively. It uses distributional analysis and pairwise divergence metrics with clustering to choose between single-device and robust multi-device Bayesian optimization strategies. Unlike conventional methods that assume homogeneous devices or enforce generic robustness, the proposed framework explicitly determines whether shared optimization across devices is appropriate based on the degree of inter-device noise heterogeneity. This enables improved performance, reproducibility, and efficiency. An experimental case study involving three nominally identical 3D printers (same brand, model, and close serial numbers) demonstrates reduced redundancy, lower resource usage, and improved reliability, along with improved convergence stability and solution quality through the selection of the appropriate optimization strategy based on the degree of inter-device noise heterogeneity. Overall, this framework establishes a general approach for precision- and resource-aware optimization in scalable, automated experimental platforms, demonstrated here on a representative multi-device 3D printing case study.

Schenk, Christina↗
DOE OSTI2021

Quantifying carbon budget, crop yields and their responses to environmental variability using the ecosys model for U.S. Midwestern agroecosystems

As one of the major agricultural production areas in the world, the United States (U.S.) Midwest plays a vital role in the global food supply and agricultural ecosystem services. Although significant efforts have been made in modeling the carbon cycle dynamics over this area, large uncertainty still exists in the previous simulations in terms of reproducing individual components of the carbon cycle and their responses to environmental variability. Here we evaluated the performance of an advanced agroecosystem model, ecosys, in simulating carbon budgets over the U.S. Midwest, considering both the magnitude of carbon flux/yield and its response to environmental (climate and soil) variability. We conducted model simulations and evaluations at 7 cropland eddy-covariance sites as well as over 293 counties of Illinois, Indiana, and Iowa in the U.S. Midwest. The site-level simulations showed that ecosys captured both the magnitude and seasonal patterns of carbon fluxes (i.e., net ecosystem carbon exchange (NEE), ecosystem gross primary production (GPP), and ecosystem respiration (Reco)), leaf area index (LAI), and dynamic plant carbon allocation processes, with R2 equal to 0.92, 0.87, 0.87, and 0.78 for GPP, NEE, Reco, and LAI, respectively across all the sites compared with the observations. For regional scale simulations, ecosys reproduced the spatial distribution and interannual variability of corn and soybean yields with the constraints of observed yields and a new remotely sensed GPP product, with R2 of multi-year averaged simulated and observed yield equal 0.83 and 0.80 for corn and soybean, respectively. The simulated responses of carbon cycle dynamics to environmental variability were consistent with that from the empirical observations at both site and regional scales. Furthermore, our results demonstrated the applicability of ecosys in simulating the carbon cycle and soil carbon dynamics of the U.S. Midwestern agroecosystems under different climate and soil conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Hybrid super homogenization and discontinuity factor method for continuous finite element diffusion

This paper presents a novel homogenization equivalence technique aiming to simultaneously leveragethe simplicity of the Super Homogenization (SPH) method to reproduce reference reaction rates andthe ability to preserve reference leakage rates at desired surfaces through Discontinuity Factors (DF).The need for this new class of methods arises from the inability of the current state-of-the-art SPH technology to properly reproduce the reactor multiplication factor for problems with significant leakage. This work shows that this defect lies in the use of normalization factors in the SPH algorithm: while they are necessary with purely reflecting problems to ensure uniqueness of the solution, they introduce homogenization inconsistencies if at least one vacuum boundary condition is present. Two solutions to this problem are presented in this work: (i) to simply remove the normalization factor or (ii) to introduce additional degrees of freedom in the form of DFs in such a way that the normalization factors can still constrain the problem. While the former clearly offers unrivaled simplicity and can lead to very satisfactory results for leakage-dominated cores or if spatial restriction of the SPH regions is applied – as demonstrated for the High Temperature Test Reactor, the latter is more robust and does not require to scale the fission terms with the reference eigenvalue.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Advances in the calculation of optical properties in superlattices; novel insights derived from the theory of finite periodic systems

Highlights: • Truly quantum calculation of optical properties in superlattices. • An approach that goes beyond the abilities of the current approaches in the literature. • Accurate calculations based on bona fide superlattice eigenvalues and eigenfunctions. • Accurate calculation of optical-transition matrix elements and new selection rules. • Optical spectra of AlGaAs and InGaN superlattices agree with experimental results. • Optical spectra features of blue diodes faithfully reproduced and fully explained. • Unexplained lines in high-resolution spectra assigned to surface and edge states. - Abstract: Using the genuine superlattice eigenvalues and eigenfunctions, and their parity symmetries, new transition selection rules are derived, and a novel approach for the explicit calculation of the optical spectra of periodic structures presented. The optical response and infrared spectra published in the last 40 years for (In, Ga)N and (Al, Ga)As superlattices, with different number of unit cells n (varying from ∼10 to ∼400) and barrier and valley widths in the range of 2.5 nm to 40 nm, and different properties have been revisited. Among the various examples discussed here, we study the narrow peaks clustered in groups reported in high resolution spectra of blue-emitting devices, that could not be explained until now. These properties are faithfully reproduced and fully understood. We show that these effects are related to and depend on the surface energy levels and their tuneable detachment. We shown also that the optical transitions experimentally observed, but forbidden by the current optical physics, are really allowed transitions. Since the number of energy eigenvalues E{sub μ,ν}{sup r,v} and eigenfunctions Ψ{sub μ{sup ′},ν{sup ′r,v}}(z) is now much greater, the number of allowed matrix-elements can be extremely large (of the order of n{sup 2} n{sub c}n{sub v}∕2, with n{sub c} and n{sub v} the number of subbands in the conduction and valence bands). We show that this number can be reduced by using the eigenfunctions’ symmetries, which are behind the parity symmetry selection rules and the leading order selection rules. These rules reduce the number of matrix-elements evaluations from ∼n{sup 2}n{sub c}n{sub v}∕2 to ∼nn{sub c}n{sub v}∕2, i.e., depending on the superlattice, from about 1000 to 100. We discuss also a third rule, that collects the contributions of the surface and edge states, reducing the number of transitions even more to N{sub s}≤n{sub c}n{sub v}. With these rules, the main peaks are conserved and their number practically matches that of the real spectrum. Excellent agreements are found with experimental results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Numerical analysis of experimental studies of methane hydrate dissociation induced by depressurization in a sandy porous medium

Methane Hydrates (MHs) are a promising energy source abundantly available in nature. Understanding the complex processes of MH formation and dissociation is critical for the development of safe and efficient technologies for energy recovery. Many laboratory and numerical studies have investigated these processes using synthesized MH-bearing sediments. A near-universal issue encountered in these studies is the spatial heterogeneous hydrate distribution in the testing apparatus. In the absence of direct observations (e.g. using X-ray computed tomography) coupled with real time production data, the common assumption made in almost all numerical studies is a homogeneous distribution of the various phases. In an earlier study (Yin et al., 2018) that involved the numerical description of a set of experiments on MH-formation in sandy medium using the excess water method, we showed that spatially heterogeneous phase distribution is inevitable and significant. In the present study, we use as a starting point the results and observations at the end of the MH formation and seek to numerically reproduce the laboratory experiments of depressurization-induced dissociation of the spatially-heterogeneous MH distribution. This numerical study faithfully reproduces the geometry of the laboratory apparatus, the initial and boundary conditions of the system, and the parameters of the dissociation stimulus, capturing accurately all stages of the experimental process. Using inverse modelling (history-matching) that minimized deviations between the experimental observations and numerical predictions, we determined the values of all the important flow, thermal, and kinetic parameters that control the system behaviour, which yielded simulation results that were in excellent agreement with the measurements of key monitored variables, i.e. pressure, temperature, cumulative production of gas and water over time. We determined that at the onset of depressurization (when the pressure drop – the driving force of dissociation – is at its maximum), the rate of MH dissociation approaches that of an equilibrium reaction and is limited by the heat transfer from the system surroundings. Furthermore, as the effect of depressurization declines over time, the dissociation reaction becomes kinetically limited despite significant heat inflows from the boundaries, which lead to localized temperature increases in the reactor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Isospin symmetry holds beyond the proton drip line: First observation of excited states in 55Cu

Isospin symmetry, combined with shell-model description, has long been successful in reproducing nuclear structure and predicting mirror energy differences. However, these descriptions have been primarily limited to mirror nuclei with well-bound ground states. As one approaches the proton drip line, coupling to the continuum becomes important. Despite this, little is known about medium-mass nuclei at the drip line. Here we report the measurement of the excitation energies of several unbound states in a nucleus beyond the proton drip line, 55Cu, performed via γ-ray spectroscopy at the Radioactive Isotope Beam Factory (RIKEN). Several low-lying yrast states are observed to decay predominantly by γ-ray emission. In contrast, the state–predicted to be strongly populated in the one-proton knockout reaction–is notably absent, suggesting that it decays primarily via prompt proton emission. We show that even for such an unbound system, shell-model calculations based on a harmonic oscillator potential reproduce the observed mirror energy differences remarkably well in all the observed states. Based on this result, we propose that a general pattern might emerge where, in medium-mass systems at the proton drip line, the stronger Coulomb barrier together with the centrifugal barrier effectively confines the wave functions, so that the shell model is able to describe low-lying states, enabling accurate predictions of mirror energy differences between isobaric analogue states.

Pigliapoco, S [Dipartimento di Fisica e Astronomia↗
DOE OSTI2018

MYPT1 is targeted by miR-145 inhibiting viability, migration and invasion in 2D and 3D HeLa cultures

Highlights: • miR-143/145 affect viability, migration and invasion on monolayer and spheroid cultures. • miR-145 produces a stronger inhibitory effect and increases pMLC levels through MYPT1 targeting. • MYPT1 silencing by siRNAs reproduces miR-145 effects. The miR-143/145 cluster is down-regulated in cervical tumor cells suggesting a role in tumorigenesis including cytoskeleton remodeling, a key event for tumor progression. The aim of the present work was to determine the role of miR-143/145 in the modulation of the myosin regulator phospho-myosin light chain (pMLC). HeLa monolayer and tridimensional cultures were transfected with miR-143 or miR-145 mimics inhibiting cell viability, proliferation, migration and invasion, mainly through miR-145. MiR-145 transfection increased pMLC levels by targeting the MYPT1 subunit of the regulatory myosin phosphatase. MYPT1 knockdown by siRNAs reproduced miR-145 effects suggesting miR-145 as a tumor suppressor through MYPT1 targeting, leading to a subsequent increase of pMLC levels with implications for cervical cell viability, migration and invasion.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗