Search NASA⌕ Search

SEARCH · Search NASA

Results for “Rigidity function”

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

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

At least 757 records · Page 42

Formations of Tethered Spacecraft as Stable Platforms for Far IR and Sub-mm Astronomy

In this paper we describe current research in tethered formations for interferometry, and a roadmap to demonstrating the required key technologies via on-ground and in-orbit testing. We propose an integrated kilometer-size tethered spacecraft formation flying concept which enables Far IR and Sub-mm astronomy observations from space. A rather general model is used to predict the dynamics, control, and estimation performance of formations of spacecraft connected by tethers in LEO and deep space. These models include the orbital and tethered formation dynamics, environmental models, and models of the formation estimator/controller/commander. Both centralized and decentralized control/sensing/estimation schemes are possible, and dynamic ranges of interest for sensing/control are described. Key component/subsystem technologies are described which need both ground-based and in-orbit demonstration prior to their utilization in precision space interferometry missions using tethered formations. Defining an orbiting formation as an ensemble of orbiting spacecraft performing a cooperative task, recent work has demonstrated the validity of the tethering the spacecraft to provide both the required formation rigidity and satisfy the formation reconfiguration needs such as interferometer baseline control. In our concept, several vehicles are connected and move along the tether, so that to reposition them the connecting tether links must vary in length. This feature enables variable and precise baseline control while the system spins around the boresight. The control architecture features an interferometer configuration composed of one central combiner spacecraft and two aligned collector spacecraft. The combiner spacecraft acts as the formation leader and is also where the centralized sensing and estimation functions reside. Some of the issues analyzed with the model are: dynamic modes of deformation of the distributed structure, architecture of the formation sensor, and sources of dynamical perturbation that need to be mitigated for precision operation in space. Examples from numerical simulation of an envisioned scenario in heliocentric orbit demonstrate the potential of the concept for space interferometry.

Quadrelli, Marco B.↗

Electron Microscopy Approaches to Unraveling the Structure of Amorphous Materials

Determining atomic structures in crystalline materials—where atoms are arranged in rigid, periodic lattices—has been highly successful using probes such as electrons, X-rays, and neutrons. In contrast, amorphous materials, despite their ubiquity and technological importance, remain far more challenging to characterize with comparable accuracy and precision. This review highlights existing, emerging, and potential (scanning) transmission electron microscopy ((S)TEM) techniques for probing short- and medium-range order in amorphous materials. Approaches ranging from high-resolution (S)TEM imaging and selected electron diffraction pattern to four-dimensional STEM (4D-STEM) based pair distribution function, fluctuation electron microscopy, tomography, ptychography, and spectroscopic methods are discussed, emphasizing their ability to provide complementary insights across multiple length scales—from sub-angstrom local environments to nanometer-scale correlations. Here, we further explore the promise of multimodal and correlative strategies, as well as the growing role of machine learning and physics-informed AI in enabling real-time, quantitative interpretation of complex structural signatures. Together, these advances point toward a future where electron microscopy not only reveals the hidden order in amorphous systems but also establishes robust structure–property relationships, paving the way for materials innovation in disordered matter.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Secondary cell wall biosynthesis

Secondary walls are synthesized in specialized cells, such as tracheary elements and fibers, and their remarkable strength and rigidity provide strong mechanical support to the cells and the plant body. The main components of secondary walls are cellulose, xylan, glucomannan and lignin. Biochemical, molecular and genetic studies have led to the discovery of most of the genes involved in the biosynthesis of secondary wall components. Cellulose is synthesized by cellulose synthase complexes in the plasma membrane and the recent success of in vitro synthesis of cellulose microfibrils by a single recombinant cellulose synthase isoform reconstituted into proteoliposomes opens new doors to further investigate the structure and functions of cellulose synthase complexes. Most genes involved in the glycosyl backbone synthesis, glycosyl substitutions and acetylation of xylan and glucomannan have been genetically characterized and the biochemical properties of some of their encoded enzymes have been investigated. The genes and their encoded enzymes participating in monolignol biosynthesis and modification have been extensively studied both genetically and biochemically. A full understanding of how secondary wall components are synthesized will ultimately allow us to produce plants with custom-designed secondary wall composition tailored to diverse applications.

59 BASIC BIOLOGICAL SCIENCES↗

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence E - Raw Data

Sequence E: Yaw Releases (P) This test sequence used a downwind, rigid turbine with a 3.4° cone angle. The wind speeds ranged from 7 m/s to 17 m/s. Initial yaw angles of ±90° were achieved. The blade tip pitch was 3°. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. The turbine was positioned at each specified yaw angle. A 30- second data set was collected with the yaw brake engaged. These points used a letter for the last digit beginning with A and proceeding through the alphabet as repeat points were needed. These points were plotted and compared to ascertain the functionality of the instrumentation. The files ending in numbers beginning with 0 and increasing with each repetition represent yaw release points. Once the fixed-position yaw test was complete, the yaw drive was engaged to hold the turbine. The yaw brake was released, and the yaw drive was disabled allowing the turbine to yaw freely. The yaw drive was disabled about 5 seconds into the campaign, and the turbine was allowed to yaw freely for the rest of the 30-second duration. When the yaw drive is disabled, no torque is applied to the motor, and the inertia and yawing force of the nacelle and rotor overcomes the friction and inertia of the yaw gear and motor. The yaw releases were repeated five times at each condition. Four additional campaigns were collected to determine the turbine’s natural yaw error by releasing the brake at 0° yaw error. These file names use the E designation, followed by two digits for wind speed, followed by XXXX, followed by the repetition digit.

17 WIND ENERGY↗

Principal Effects of Axial Load on Moment-Distribution Analysis of Rigid Structures

This thesis presents the method of moment distribution modified to include the effect of axial load upon the bending moments. This modification makes it possible to analyze accurately complex structures, such as rigid fuselage trusses, that heretofore had to be analyzed by approximate formulas and empirical rules. The method is simple enough to be practicable even for complex structures, and it gives a means of analysis for continuous beams that is simpler than the extended three-moment equation now in common use. When the effect of axial load is included, it is found that the basic principles of moment distribution remain unchanged, the only difference being that the factors used, instead of being constants for a given member, become functions of the axial load. Formulas have been developed for these factors, and curves plotted so that their applications requires no more work than moment distribution without axial load. Simple problems have been included to illustrate the use of the curves.

James, Benjamin Wylie↗

CO 2 Adsorption on PIMs Studied with 13 C NMR Spectroscopy

Polymers of intrinsic microporosity (PIMs) can be functionalized with nitrogen-bearing groups, such as a nitrile group for PIM-1 and a methylated tetrazole group for MTZ-PIM, to physisorb CO 2 . Static 13 C NMR spectroscopy was used to explore CO 2 dynamics by analysis of one-dimensional spectra and T 2 relaxation times at temperatures ranging from 6 to 295 K. Each polymer had two different loadings of CO 2 to confirm the dominating relaxation mechanism at lower temperatures and deduce the identity of the weaker secondary adsorption site on the rigid polymer backbone. Both polymers, with a high loading of CO 2 , experience a narrowing of the 13 C NMR line width as the samples are warmed to 80 K, which is the onset of the weakly adsorbed CO 2 hopping on and off the polymer backbone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding Pore Filling Processes and Adsorption/Desorption Hysteresis in Nanoporous Metal–Organic Frameworks: Insights from Grand Canonical Monte Carlo Simulations and Free Energy Calculations

Grand canonical Monte Carlo (GCMC) simulations were used to investigate pore filling and hysteresis in nanoporous metal-organic frameworks (MOFs). Adsorption and desorption isotherms were calculated for argon at 87 K in 1866 MOFs from the CoRE MOF database and for short n-alkanes in selected MOFs, keeping the adsorbent structure rigid. Analysis of the molecular configurations showed two different mechanisms and origins of hysteresis: one involving a transition of the adsorbate arrangement in the pores similar to a gas-to-liquid transition associated with a large change in the loading and one more similar to a liquid-to-solid transition associated with a relatively small change in the loading. Our GCMC simulations in MOFs with diverse pore topologies indicate exceptions to an empirical relationship for the minimum diameter of a cylindical pore required for hysteresis as a function of the adsorbate diameter and reduced temperature. The simulations reveal some structures where isotherms exhibit two steps in the adsorption branch and only one step in the desorption branch. Hysteresis loops with a different number of adsorption and desorption steps are not common. Here, to better understand why hysteresis is observed in the GCMC simulations, the concept of the transition probability for observing a step in the adsorption isotherm at a given pressure in a GCMC simulation is introduced. We used two different methods to calculate the transition probabilities and find that these yield comparable results. Furthermore, the transition probability provides a measure for the length of GCMC simulations to yield reliable results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure of the lutein-binding domain of human StARD3 at 1.74 Å resolution and model of a complex with lutein

A crystal structure of the lutein-binding domain of human StARD3 (StAR-related lipid-transfer protein 3; also known as MLN64) has been refined to 1.74 Å resolution. A previous structure of the same protein determined to 2.2 Å resolution highlighted homology with StARD1 and shared cholesterol-binding character. StARD3 has since been recognized as a carotenoid-binding protein in the primate retina, where its biochemical function of binding lutein with specificity appears to be well suited to recruit this photoprotective molecule. The current and previous structures correspond closely to each other (r.m.s.d. of 0.25 Å), especially in terms of the helix-grip fold constructed around a solvent-filled cavity. Regions of interest were defined with alternate conformations in the current higher-resolution structure, including Arg351 found within the cavity and Ω1, a loop of four residues found just outside the cavity entrance. Models of the complex with lutein generated by rigid-body docking indicate that one of the ionone rings must protrude outside the cavity, and this insight has implications for molecular interactions with transport proteins and enzymes that act on lutein. Interestingly, models with the ϵ-ionone ring characteristic of lutein pointing towards the bottom of the cavity were associated with fewer steric clashes, suggesting that steric complementarity and ligand asymmetry may play a role in discriminating lutein from the other ocular carotenoids zeaxanthin and meso -zeaxanthin, which only have β-ionone rings.

Biochemistry & Molecular Biology↗

Contaminated Solar Array Handrail Samples Retrieved From Mir Analyzed

In January 1998 during the shuttle STS 89 mission, an eight-section Russian solar array panel was retrieved after more than 10 years of exposure to the orbital space environment on Mir. The array was deployed June 16, 1987, and removed on November 3, 1997. It had been actively used as a source of electrical power for 8 years. This operational array had been located on the Mir core module, located directly above the Kvant-2 module. Its retrieval provided a unique opportunity to study the effects of the low-Earth-orbit environment on a functional solar array. The intact solar array underwent scientific inspections and preliminary tests by a joint team of U.S. and Russian investigators to evaluate the effects of long-term space exposure. Upon initial examination, significant contamination was observed over most components of the array. One panel, panel 8, was provided to the U.S. scientists for further evaluation. As part of the U.S. investigations, two solar array handrail samples from panel 8 were evaluated for contamination at the NASA Glenn Research Center at Lewis Field. One is a section of a rigid handrail, and the other is a section of woven fabric tape that was overwrapped around a flexible handhold. Both the flexible handhold woven fabric and the rigid handrail were significantly darkened after 10 years of space exposure. They were evaluated with optical microscopy, field emission scanning electron microscopy (FESEM), and energy-dispersive spectroscopy. Solar absorptance and room-temperature emittance values also were obtained. The returned contaminated solar array segment is very similar in design to the solar arrays being supplied by the Russians for the International Space Station. Therefore, it was desirable to determine what the contaminants on various surfaces are and what the sources of the contamination were.

deGroh, Kim K.↗

Circumferential pressure probe

A probe for measuring circumferential pressure inside a body cavity is disclosed. In the preferred embodiment, a urodynamic pressure measurement probe for evaluating human urinary sphincter function is disclosed. Along the length of the probe are disposed a multiplicity of deformable wall sensors which typically comprise support tube sections with flexible side wall areas. These are arranged along the length of the probe in two areas, one just proximal to the tip for the sensing of fluid pressure inside the bladder, and five in the sensing section which is positioned within the urethra at the point at which the urinary sphincter constricts to control the flow of urine. The remainder of the length of the probe comprises multiple rigid support tube sections interspersed with flexible support tube sections in the form of bellows to provide flexibility.

Holmes, Harlan K.↗

Robonaut 2 - IVA Experiments On-Board ISS and Development Towards EVA Capability

Robonaut 2 (R2) has completed its fixed base activities on-board the ISS and is scheduled to receive its climbing legs in early 2014. In its continuing line of firsts, the R2 torso finished up its on-orbit activities on its stanchion with the manipulation of space blanket materials and performed multiple tasks under teleoperation control by IVA astronauts. The successful completion of these two IVA experiments is a key step in Robonaut's progression towards an EVA capability. Integration with the legs and climbing inside the ISS will provide another important part of the experience that R2 will need prior to performing tasks on the outside of ISS. In support of these on-orbit activities, R2 has been traversing across handrails in simulated zero-g environments and working with EVA tools and equipment on the ground to determine manipulation strategies for an EVA Robonaut. R2 made significant advances in robotic manipulation of deformable materials in space while working with its softgoods task panel. This panel features quarter turn latches that secure a space blanket to the task panel structure. The space blanket covers two cloth cubes that are attached with Velcro to the structure. R2 was able to open and close the latches, pull back the blanket, and remove the cube underneath. R2 simulated cleaning up an EVA worksite as well, by replacing the cube and reattaching the blanket. In order to interact with the softgoods panel, R2 has both autonomously and with a human in the loop identified and localized these deformable objects. Using stereo color cameras, R2 identified characteristic elements on the softgoods panel then extracted the location and orientation of the object in its field of view using stereo disparity and kinematic transforms. R2 used both vision processing and supervisory control to successfully accomplish this important task. Teleoperation is a key capability for Robonaut's effectiveness as an EVA system. To build proficiency, crewmembers have attempted increasingly difficult tasks using R2 inside the Station. After donning motion capture equipment and a virtual reality visor, Expedition 34/35 flight engineer Tom Marshburn began operations with simple hand movements. Having gained confidence, Marshburn guided R2's arms in a leader-follower exercise with crewmate Chris Cassidy. He was also able to use the hand to grab a tumbling roll of tape, a task only demonstrable in microgravity. Later efforts saw Cassidy handle softgoods through shared control with ground operators, mimicking an activity previously achieved using only autonomy. Robotic climbing through the ISS on handrails requires both precision motion and compliant grasps in order to both position grippers on handrails/seat track and prevent large internal forces. R2 climbs using actively controlled compliance and torque limiting to meet both the precision and softness requirements. During a step, the attached leg is controlled to be strong and stiff in order to maintain precision trajectory tracking. The swing leg is controlled to be stiff but weak to minimize unintentional impact forces while maintaining precision. During a simulated dual limb grasp (as shown in Figure 1), the R2 controller maintains one limb rigid and one limb soft to prevent large internal forces from building up. R2's grippers also use a form of force control to limit grip force while not fully closed on either a handrail or seat track thus limiting unintentional forces on cables/objects that may be present in R2's translational path. The on-board torso R2 safety system relies on a single end-effector velocity limit to prevent potential impact forces from exceeding Station maximum load requirements. R2's mobile configuration required modifications to the velocity limiting safety function due to its large, dynamic inertia. R2's legs maneuver the robot's mass creating configuration dependent, joint-relative inertias. A single all-encompassing velocity limit to cover worst case inertia is prohibitively low. The upgraded R2 control and safety systems solve this problem using momentum limiting, momentum control, and kinetic energy minimization. Momentum and kinetic energy take the robot mass into account relieving low velocity restrictions on low inertia end-effectors while ensuring that the overall mass of R2 is limited from hazardous velocities. The momentum of R2's five safety nodes (each of the four end-effectors and the body) is monitored and compared to a single momentum limit. If any of the five nodes exceeds the safety limit, the motor power is removed and the robot comes to a stop. Momentum control/limiting also provides a simple, reliable method to integrate hand held tools into the safety system by providing the tool mass to the control system thus automatically reducing the allowable velocity of the end-effector with the tool. Work on the ground continues to build the skill set for an EVA Robonaut. Recent experiments (Figure 2) demonstrate how a teleoperator can use R2 to manipulate a tether hook, an important safety precaution on spacewalks. Another task displayed Robonaut's ability to pull back a protective jacket over a hose and search for damage, as well as inspect a quick-disconnect fitting for debris. Demonstrations such as these are indicative of EVA work done on ISS, specifically seen during a series of spacewalks over 2012 and 2013 where astronauts searched for an ammonia leak in one of the external cooling loops. Through experiments both on ISS and on the ground, R2 is evolving and providing the information needed to plan out the upgrades that will make an EVA Robonaut an effective tool. With the addition of legs, R2 will start climbing inside the space station and supply invaluable information on how the climbing strategies and task stabilization techniques must be refined. Ground R2 systems will continue to work with additional EVA tools and equipment in preparation for onboard IVA testing and future EVA applications.

Diftler, Myron↗

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

New Variational Formulations of Hybrid Stress Elements

In the variational formulations of finite elements by the Hu-Washizu and Hellinger-Reissner principles the stress equilibrium condition is maintained by the inclusion of internal displacements which function as the Lagrange multipliers for the constraints. These versions permit the use of natural coordinates and the relaxation of the equilibrium conditions and render considerable improvements in the assumed stress hybrid elements. These include the derivation of invariant hybrid elements which possess the ideal qualities such as minimum sensitivity to geometric distortions, minimum number of independent stress parameters, rank sufficient, and ability to represent constant strain states and bending moments. Another application is the formulation of semiLoof thin shell elements which can yield excellent results for many severe test cases because the rigid body nodes, the momentless membrane strains, and the inextensional bending modes are all represented.

Pian, T. H. H.↗

Mathematical modelling of flow distribution in the human cardiovascular system

The paper presents a detailed model of the entire human cardiovascular system which aims to study the changes in flow distribution caused by external stimuli, changes in internal parameters, or other factors. The arterial-venous network is represented by 325 interconnected elastic segments. The mathematical description of each segment is based on equations of hydrodynamics and those of stress/strain relationships in elastic materials. Appropriate input functions provide for the pumping of blood by the heart through the system. The analysis employs the finite-element technique which can accommodate any prescribed boundary conditions. Values of model parameters are from available data on physical and rheological properties of blood and blood vessels. As a representative example, simulation results on changes in flow distribution with changes in the elastic properties of blood vessels are discussed. They indicate that the errors in the calculated overall flow rates are not significant even in the extreme case of arteries and veins behaving as rigid tubes.

NASA Center JSC↗

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases: Preprint

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Structures of honeybee-infecting Lake Sinai virus reveal domain functions and capsid assembly with dynamic motions

Understanding the structural diversity of honeybee-infecting viruses is critical to maintain pollinator health and manage the spread of diseases in ecology and agriculture. We determine cryo-EM structures of $T$ = 4 and $T$ = 3 capsids of virus-like particles (VLPs) of Lake Sinai virus (LSV) 2 and delta-N48 LSV1, belonging to tetraviruses, at resolutions of 2.3–2.6 Å in various pH environments. Structural analysis shows that the LSV2 capsid protein (CP) structural features, particularly the protruding domain and C-arm, differ from those of other tetraviruses. The anchor loop on the central β-barrel domain interacts with the neighboring subunit to stabilize homo-trimeric capsomeres during assembly. Delta-N48 LSV1 CP interacts with ssRNA via the rigid helix α1’, α1’–α1 loop, β-barrel domain, and C-arm. Cryo-EM reconstructions, combined with X-ray crystallographic and small-angle scattering analyses, indicate that pH affects capsid conformations by regulating reversible dynamic particle motions and sizes of LSV2 VLPs. C-arms exist in all LSV2 and delta-N48 LSV1 VLPs across varied pH conditions, indicating that autoproteolysis cleavage is not required for LSV maturation. The observed linear domino-scaffold structures of various lengths, made up of trapezoid-shape capsomeres, provide a basis for icosahedral $T$ = 4 and $T$ = 3 architecture assemblies. These findings advance understanding of honeybee-infecting viruses that can cause Colony Collapse Disorder.

59 BASIC BIOLOGICAL SCIENCES↗

Few-cycle atomic layer deposition to nanoengineer polybenzimidazole for H 2 /CO 2 separation

Atomic layer deposition (ALD) creates uniform sub-nanometer films on a variety of surfaces and nanopore walls and has been used to modify polymers to improve surface affinity towards specific molecules, solvent resistance, and barrier properties to gases and vapors. Here, for the first time, we demonstrate that few-cycle ALD can be used to engineer functional polymers at a sub-nanometer scale to improve both molecular size-sieving ability and counterintuitively, gas permeability. Particularly, 1-cycle ALD treatment of polybenzimidazole (PBI) by sequential exposure to trimethylaluminum (TMA) and water vapor remarkably increases H 2 permeability by 120% - 270% and H 2 /CO 2 selectivity by 30% at 35–200 °C. The ALD not only deposits an AlOx layer on the surface but also enables the TMA to infiltrate and react with the bulk PBI to form an AlO x network, disrupting polymer chain packing and increasing chain rigidity. The membrane exhibits excellent stability when challenged with simulated syngas, overcoming the permeability/selectivity tradeoff for H 2 /CO 2 separation. In conclusion, this study showcases a facile and scalable way of engineering polymeric membranes at a sub-nanometer level to improve molecular separation performance.

36 MATERIALS SCIENCE↗

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↗