SEARCH · Engineering Papers

Results for “Trending analysis”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,729 records · Page 96

All sources NASA NTRS DOE OSTI

NASA NTRS1992

State-selected chemical reaction dynamics at the S matrix level - Final-state specificities of near-threshold processes at low and high energies

State-to-state reaction probabilities are found to be highly final-state specific at state-selected threshold energies for the reactions O + H2 yield OH + H and H + H2 yield H2 + H. The study includes initial rotational states with quantum numbers 0-15, and the specificity is especially dramatic for the more highly rotationally excited reactants. The analysis is based on accurate quantum mechanical reactive scattering calculations. Final-state specificity is shown in general to increase with the rotational quantum number of the reactant diatom, and the trends are confirmed for both zero and nonzero values of the total angular momentum.

Chatfield, David C.↗

NASA NTRS2004

Time Changes of the European Gravity Field from GRACE: A Comparison with Ground Measurements from Superconducting Gravimeters and with Hydrology Model Predictions

We investigate the time-variable gravity changes in Europe retrieved from the initial GRACE monthly solutions spanning a 18 month duration from April 2002 to October 2003. Gravity anomaly maps are retrieved in Central Europe from the monthly satellite solutions we compare the fields according to various truncation levels (typically between degree 10 and 20) of the initial fields (expressed in spherical harmonics to degree 120). For these different degrees, an empirical orthogonal function (EOF) decomposition of the time-variable gravity field leads us to its main spatial and temporal characteristics. We show that the dominant signal is found to be annual with an amplitude and a phase both in agreement with predictions in Europe modeled using snow and soil-moisture variations from recent hydrology models. We compare these GRACE gravity field changes to surface gravity observations from 6 superconducting gravimeters of the GGP (Global Geodynamics Project) European sub-network, with a special attention to loading corrections. Initial results suggest that all 3 data sets (GRACE, hydrology and GGP) are responding to annual changes in near-surface water in Europe of a few microGal (at length scales of approx.1000 km) that show a high value in winter and a summer minimum. We also point out that the GRACE gravity field evolution seems to indicate that there is a trend in gravity between summer 2002 and summer 2003 which can be related to the 2003 heatwave in Europe and its hydrological consequences (drought). Despite the limited time span of our analysis and the uncertainties in retrieving a regional solution from the network of gravimeters, the calibration and validation aspects of the GRACE data processing based on the annual hydrology cycle in Europe are in progress.

Hinderer, J.↗

NASA NTRS2000

Noise-Reduction Benefits Analyzed for Over-the-Wing-Mounted Advanced Turbofan Engines

As we look to the future, increasingly stringent civilian aviation noise regulations will require the design and manufacture of extremely quiet commercial aircraft. Also, the large fan diameters of modern engines with increasingly higher bypass ratios pose significant packaging and aircraft installation challenges. One design approach that addresses both of these challenges is to mount the engines above the wing. In addition to allowing the performance trend towards large diameters and high bypass ratio cycles to continue, this approach allows the wing to shield much of the engine noise from people on the ground. The Propulsion Systems Analysis Office at the NASA Glenn Research Center at Lewis Field conducted independent analytical research to estimate the noise reduction potential of mounting advanced turbofan engines above the wing. Certification noise predictions were made for a notional long-haul commercial quadjet transport. A large quad was chosen because, even under current regulations, such aircraft sometimes experience difficulty in complying with certification noise requirements with a substantial margin. Also, because of its long wing chords, a large airplane would receive the greatest advantage of any noise-shielding benefit.

Berton, Jeffrey J.↗

NASA NTRS2014

Modeling of Commercial Turbofan Engine With Ice Crystal Ingestion: Follow-On

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which is ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in flight. The computational tool was utilized to help guide a portion of the PSL testing, and was used to predict ice accretion could also occur at significantly lower altitudes. The predictions were qualitatively verified by subsequent testing of the engine in the PSL. In a previous study, analysis of select PSL test data points helped to calibrate the engine icing computational tool to assess the risk of ice accretion. This current study is a continuation of that data analysis effort. The study focused on tracking the variations in wet bulb temperature and ice particle melt ratio through the engine core flow path. The results from this study have identified trends, while also identifying gaps in understanding as to how the local wet bulb temperature and melt ratio affects the risk of ice accretion and subsequent engine behavior.

Jorgenson, Philip C. E.↗

NASA NTRS2014

Modeling of Commercial Turbofan Engine with Ice Crystal Ingestion; Follow-On

Jorgenson, Philip C. E.↗

NASA NTRS2013

Virtual Construction of Space Habitats: Connecting Building Information Models (BIM) and SysML

Current trends in design, construction and management of complex projects make use of Building Information Models (BIM) connecting different types of data to geometrical models. This information model allow different types of analysis beyond pure graphical representations. Space habitats, regardless their size, are also complex systems that require the synchronization of many types of information and disciplines beyond mass, volume, power or other basic volumetric parameters. For this, the state-of-the-art model based systems engineering languages and processes - for instance SysML - represent a solid way to tackle this problem from a programmatic point of view. Nevertheless integrating this with a powerful geometrical architectural design tool with BIM capabilities could represent a change in the workflow and paradigm of space habitats design applicable to other aerospace complex systems. This paper shows some general findings and overall conclusions based on the ongoing research to create a design protocol and method that practically connects a systems engineering approach with a BIM architectural and engineering design as a complete Model Based Engineering approach. Therefore, one hypothetical example is created and followed during the design process. In order to make it possible this research also tackles the application of IFC categories and parameters in the aerospace field starting with the application upon the space habitats design as way to understand the information flow between disciplines and tools. By building virtual space habitats we can potentially improve in the near future the way more complex designs are developed from very little detail from concept to manufacturing.

space architecture↗

NASA NTRS2018

Influence of Load Profile on Biomechanics of the Squat and Deadlift

Long duration space travel will expose astronauts to extended periods of reduced gravity. Since gravity is not present to assist loading, astronauts will use resistive and aerobic exercise regimes for the duration of the space flight to counteract the effect reduced gravity has on the body. Astronauts will exercise on a flywheel based device on the second Orion Exploration Mission (EM2). The effect that the flywheel load profile has on biomechanics is unknown when compared to free weights or a simulated free weight profile. The purpose of this evaluation is to compare the differences in lower body kinematics and kinetics between the flywheel and free weight profile. Test subjects were instrumented with reflective markers for motion capture data collection while exercising on the Device for Aerobic and Resistive Training (DART) prototype developed by TDA Research, Inc. of Wheat Ridge, CO. Exercises performed were either the squat while wearing a harness or deadlift while grasping a T-bar, both of which interfaced with the DART through a rope cable. The DART, a motorized device controlled via software, is capable of producing different load profiles. Profiles include simulated free weights with adjustable eccentric overload capability and flywheel. Test variables included the applied device load and the load profile, both set through the DART software interface. Motion capture data was collected with a 12 camera system (Smart-DX, BTS Bioengineering, Brooklyn, NY). Bilateral ground reaction force data were collected with force plates (P6000, BTS Bioengineering). DART cable force was recorded through an internal load cell. Data was collected from a total of four test subjects. The three repetition maximum was determined for each test subject for both squat and deadlift. High, medium and low loads were determined based on this maximum. The test subject performed five repetitions of each exercise at each load and each test trial was repeated twice. Cadence was controlled during exercising. Biomechanical data were used to drive the models utilizing the OpenSim software platform. Subject specific models were scaled to match the anthropometrics of the test subjects and used to estimate the peak joint angle, joint range of motion, the peak joint moment and the joint moment impulse of the lower extremity joints and the lumbar joint.These are preliminary results because the data analysis is ongoing. There was a lot of inter-subject variability, however, in general, there was a greater peak lumbar flexion angle for the flywheel squat exercise. There was a trend towards a greater range of motion at the hip for the squat exercise with free weight. There was also a greater hip extension, lumbar extension and ankle plantar flexion impulse moment during the squat exercise with free weight. During the deadlift exercise, the peak knee flexion angle and peak knee extension moment were greater when using the free weight profile. For deadlift, the hip extension, lumbar extension and ankle plantar flexion impulse moment tended to be greater with the free weight profile. Overall, the kinematic and kinetic outcomes for the flywheel profile were either statistically the same as free weight profile in many cases, significantly reduced in a few cases, and rarely higher.

Bioengineering↗

NASA NTRSDate not supplied

Improved Characterization of PSC Processes Derived from a Third-Generation CALIOP and MLS Detection and Composition Classification Algorithm

The new 3-year CloudSat and CALIPSO Science Team project described in this poster will use a unique combination of data from the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) instrument on CALIPSO and the Microwave Limb Sounder (MLS) on Aura, in conjunction with supporting meteorological information and detailed modeling studies, to advance our understanding of polar stratospheric cloud (PSC) processes and their role in ozone depletion. We will develop a third-generation (Gen3) PSC detection and composition algorithm that incorporates a new, more robust two-dimensional, multi-channel CALIOP feature detection scheme (2D-McDA). We will also devise and implement an improved two-dimensional PSC composition classification scheme that utilizes multiple parameters (e.g., CALIOP 532-nm parallel and perpendicular scattering ratios, CALIOP 1064-nm total scattering ratio, MLS HNO3 and H2O, ambient temperature, and temperature histories) in a Bayesian approach to determine the most likely PSC composition and help constrain solid PSC particle number density and size/shape. The combined CALIOP/MLS analyses will allow us to study in detail the full life cycle of PSCs and their resulting impact on gas-phase HNO3 and H2O, which should lead to improved parameterizations of PSC microphysics in global CCMs where detailed particle information is not available. The Gen3 CALIOP PSC algorithm will be a natural stepping-stone toward the analysis of data collected during future spaceborne lidar missions, such as NASA’s Atmosphere Observation System (AtmOS) mission currently scheduled for launch late in this decade. We will also investigate possible trends in PSC occurrence and composition over the entire CALIOP data record and through further comparisons with the Stratospheric Aerosol Measurement (SAM) II solar occultation PSC record from 1979-1989. Finally, we will validate the mountain-wave parameterization and PSC schemes used in the UM-UKCA (Unified Model coupled to the United Kingdom Chemistry and Aerosol module) chemistry-climate model through detailed comparisons with earlier CALIOP PSC data products and those developed under this proposal.

CALIPSO↗

NASA NTRSDate not supplied

Improved Characterization of PSC Processes Derived from a Third-Generation CALIOP and MLS Detection and Composition Classification Algorithm

The new 3-year CloudSat and CALIPSO Science Team project described in this poster will use a unique combination of data from the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) instrument on CALIPSO and the Microwave Limb Sounder (MLS) on Aura, in conjunction with supporting meteorological information and detailed modeling studies, to advance our understanding of polar stratospheric cloud (PSC) processes and their role in ozone depletion. We will develop a third-generation (Gen3) PSC detection and composition algorithm that incorporates a new, more robust two-dimensional, multi-channel CALIOP feature detection scheme (2D-McDA). We will also devise and implement an improved two-dimensional PSC composition classification scheme that utilizes multiple parameters (e.g., CALIOP 532-nm parallel and perpendicular scattering ratios, CALIOP 1064-nm total scattering ratio, MLS HNO3 and H2O, ambient temperature, and temperature histories) in a Bayesian approach to determine the most likely PSC composition and help constrain solid PSC particle number density and size/shape. The combined CALIOP/MLS analyses will allow us to study in detail the full life cycle of PSCs and their resulting impact on gas-phase HNO3 and H2O, which should lead to improved parameterizations of PSC microphysics in global chemistry-climate models (CCMs) where detailed particle information is not available. The Gen3 CALIOP PSC algorithm will be a natural stepping-stone toward the analysis of data collected during future spaceborne lidar missions, such as NASA’s Atmosphere Observation System (AOS) mission currently scheduled for launch late in this decade. We will also investigate possible trends in PSC occurrence and composition over the entire CALIOP data record and through further comparisons with the Stratospheric Aerosol Measurement (SAM) II solar occultation PSC record from 1979-1989. Finally, we will validate the mountain-wave parameterization and PSC schemes used in the UM-UKCA (Unified Model coupled to the United Kingdom Chemistry and Aerosol module) CCM through detailed comparisons with earlier CALIOP PSC data products and those developed under this proposal.

CALIPSO↗

NASA NTRS2009

The Duration of Chemical Weathering of Gusev Crater's Wishstone-Watchtower Sequence

Mineralogical abundance of primary minerals versus secondary minerals, chemical mixing relationships, and elemental ratios have been used to assess the degree of aqueous alteration at Gusev Crater and Meridiani Planum. However, limited work has used Ti-normalized mass-balance analysis chemical data to quantify gains and losses of elements from altered materials as well as estimate the duration of aqueous alteration on Mars. The Ti-normalized mass-balance approach accounts for volumetric changes associated with geochemical alteration. If volumetric changes are not considered, observed geochemical trends based on un-normalized data have the potential to be misleading. Assessing gains and losses from altered materials can indicate the geochemistry of fluids involved in the alteration. Furthermore, elemental losses can be combined with dissolution rates to estimate the duration of chemical weathering. Knowledge of the duration of aqueous alteration will provide insight into the climate history of Mars as well as indicate the potential for microbial habitability. The Wishstone-Watchtower materials in Gusev Crater are suitable candidates for Ti-normalized mass-balance analysis because mixing relationships of these two materials indicate that Watchtower materials may be derived from Wishstone-like materials. The objectives of this work are to (1) employ Ti-normalized mass-balance to assess gains and losses from the Wishstone-Watchtower sequence and (2) to combine losses with laboratory dissolution rates to estimate alteration times of the Watchtower material.

Sutter, Brad↗

NASA NTRS2000

Comparison of On-Orbit Performance of Rate Sensing Gyroscopes

This work presents results from the study of a large volume of spacecraft flight data pertaining to gyroscope performance. We have examined long and short term trends of gyroscope biase of first, an exact solution for the time-dependence of the attitude part of the state error covariance, averaged over the three spacecraft axes, is given. This solution is more complete than the usual cubic polynomial for the variance in that it includes coupling among the three gyroscope axes and includes an important term arising from the initial correlations. Second, several continuous 24-hour spans of gyroscope data are examined to verify the short-term statistical model. This analysis demonstrates that in-flight data can be used to determine the strength of the white noise driving the random walk of the gyroscope bias. This may be useful for postlaunch improvement to the noise model and for diagnosing the health of the gyroscope. Third, the long-term trends in gyroscope biases show a nearly linear systematic variation over time scales of years. This has been found on three different missions. While the random walk model is adequate as a basis for onboard Kalman filters or for state estimation using relatively short time spans, these trends indicate that some applications could benefit by accounting for the secular changes in the biases. One example is a new gyroscope calibration method that is under development that allows for multi-epoch bias solutions.

Sedlak, Joseph↗

DOE OSTI2021

Reconstructing porous media using generative flow networks

One area of intense scientific interest for the study of sandstones, carbonates, and shale at the pore scale is the use of limited image and petrophysical data to generate multiple realizations of a rock’s pore structure. Such images aid efforts to quantify uncertainty in petrophysical properties, including porosity–permeability transforms. We develop and evaluate a deep learning-based method to synthesize porous media volumes using a so-called generative flow model trained on x-ray microscope images of rock texture and pore structure. These models are optimized on a log-likelihood objective and they synthesize large and realistic three-dimensional images. Further, we demonstrate the rapid generation of sandstone image volumes that display realism as gauged by quantitative comparison of topological features using Minkowski functionals of porosity, specific surface area, and the Euler–Poincaré characteristic (i.e., pore connectivity). We also evaluate the single-phase permeability using Navier–Stokes and lattice Boltzmann methods and show that transport properties of the generated samples match measured trends.

58 GEOSCIENCES↗

NASA NTRS1993

Active tectonics in southern Xinjiang, China: Analysis of terrace riser and normal fault scarp degradation along the Hotan-Qira fault system

The northern piedmont of the western Kunlun mountains (Xinjiang, China) is marked at its easternmost extremity, south of the Hotan-Qira oases, by a set of normal faults trending N50E for nearly 70 km. Conspicuous on Landsat and SPOT images, these faults follow the southeastern border of a deep flexural basin and may be related to the subsidence of the Tarim platform loaded by the western Kunlun northward overthrust. The Hotan-Qira normal fault system vertically offsets the piedmont slope by 70 m. Highest fault scarps reach 20 m and often display evidence for recent reactivations about 2 m high. Successive stream entrenchments in uplifted footwallls have formed inset terraces. We have leveled topographic profiles across fault scarps and transverse abandoned terrace risers. The state of degradation of each terrace edge has been characterized by a degradation coefficient tau, derived by comparison with analytical erosion models. Edges of highest abandoned terraces yield a degradation coefficient of 33 +/- 4 sq.m. Profiles of cumulative fault scarps have been analyzed in a similar way using synthetic profiles generated with a simple incremental fault scarp model.

Avouac, Jean-Philippe↗

DOE OSTI2020

Probing the Atomic-Scale Structure of Amorphous Aluminum Oxide Grown by Atomic Layer Deposition

Atomic layer deposition (ALD) is a well-established technique for depositing nanoscale coatings with pristine control of film thickness and composition. The trimethylaluminum (TMA) and water (H 2 O) ALD chemistry is inarguably the most widely used and yet to date, we have little information about the atomic-scale structure of the amorphous aluminum oxide (AlO x ) formed by this chemistry. This lack of understanding hinders our ability to establish process-structure-property relationships and ultimately limits technological advancements employing AlOx made via ALD. In this work, we employ synchrotron high-energy X-ray diffraction (HE-XRD) coupled with pair distribution function (PDF) analysis to characterize the atomic structure of amorphous AlO x ALD coatings. We combine ex situ and in operando HE-XRD measurements on ALD AlO x and fit these experimental data using stochastic structural modeling to reveal variations in the Al-O bond length, Al and O coordination environment, and extent of Al vacancies as a function of growth conditions. In particular, the local atomic structure of ALD AlO x is found to change with the substrate and number of ALD cycles. The observed trends are consistent with the formation of bulk Al 2 O 3 surrounded by an O-rich surface layer. We deconvolute these data to reveal atomic-scale structural information for both the bulk and surface phases. Overall, this work demonstrates the usefulness of HE-XRD and PDF analysis in improving our understanding of the structure of amorphous ALD thin films and provides a pathway to evaluate how process changes impact the structure and properties of ALD films.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Root anatomy helps to reconcile observed root trait syndromes in tropical tree species

Premise Studying the organization of functional traits in plant leaves and stems has revealed notable patterns linking function and form; however, evidence of similarly robust organization in root tissues remains controversial. We posit that anatomical traits in roots can provide insight on the overall organization of the root system. We hypothesized that size variation in the tissue outside the stele is related in a nonlinear fashion with functional traits associated with direct resource uptake, including a negative relationship with root architectural traits, and that similar relationships detected in tropical areas also hold true in other biomes. Methods We addressed our hypotheses using empirical data from 24 tropical tree species in French Guiana, including anatomical measurements in first order roots and functional trait description for the entire fine root system. In addition, we compiled a global meta‐analysis of root traits for 500+ forest species across tropical, subtropical, and temperate forests. Results Our results supported the expected nonlinear relationships between cortical size and morphological traits and a negative linear trend with architectural traits. We confirmed a global negative relationship among specific root length (SRL), diameter, and tissue density, suggesting similar anatomical constraints in root systems across woody plants. However, the importance of factors varies across biomes, possibly related to the unequal phylogenetic representation across latitudes. Conclusions Our findings imply that the rhizocentric hypothesis can be a valuable approach to understand fine root trait syndromes and the evolution of absorptive roots in vascular plants.

Valverde‐Barrantes, Oscar J.↗

DOE OSTI2025

Implications of Changes in Water Stress and Precipitation Extremes for Cocoa Production in Côte D'ivoire and Ghana

ABSTRACT Climate change induces high variability in drought patterns and extreme precipitation indices in rainfed cocoa farming, impacting cocoa production. This study evaluated water stress, meteorological and agricultural drought conditions, and critical extreme precipitation indices in the world's two largest cocoa‐producing nations from 1981 to 2022. The results revealed a significant reduction in total annual precipitation (PRCPTOT), in the last three decades, with the greatest decline in the 1991–2000 and 2011–2022 periods. Ghana experienced the most significant reduction up to 15% (200 mm/year) in the last decade, attributed to a substantial decrease in wet days number (RR1) up to 25 days per year, a reduction in maximum consecutive wet days (CWD) up to 6 days per year, and an increase in maximum consecutive dry days (CDD) up to 15 days per year. Moreover, there was a notable decline in the Simple Daily Intensity Index (SDII), with reductions of up to 4 mm/day in certain areas, contributing to increased drought frequency, severity, and duration. In the most recent decade (2011–2022), particularly during the extremely dry years of 2013 and 2015, cocoa‐growing regions in Ghana (GHA) and eastern Côte d'Ivoire (CIV) experienced prolonged agricultural drought expressed by soil moisture deficit, typically extending from May to September. Additionally, large portions of central and eastern Ghana, as well as northeastern Côte d'Ivoire experienced sustained water stress, with over three consecutive months of total monthly precipitation falling below 100 mm, negatively impacting cocoa productivity. The decrease in the yield in the range of 2.5% to 37% was noted in the dry years and the following years, varying according to the country depending on the severity of the drought. Sensitivity analysis highlights cocoa yield's responsiveness to drought and water stress, particularly in specific years when water stress occurred, such as 1984,1985, 1989, 1995, 1999, 2000, and 2008. Considering the observed trends in precipitation patterns and their impact on cocoa production, it is crucial to acknowledge the inherent uncertainty of future precipitation patterns due to climate change. To address this challenge effectively, our study underscores the importance of identifying and closely monitoring regions currently facing water stress, as determined by precipitation and drought indicators. Over the analysed period (1981–2022), we have noted shifts in the distribution of water‐stressed areas, highlighting the dynamic nature of this issue. Consequently, we advocate for a targeted approach to implement cocoa supplementary irrigation in consistently water‐stressed regions.

Obahoundje, Salomon [International Water Managemen↗

DOE OSTI2020

Interlaboratory Reproducibility of Contour Method Data Analysis and Residual Stress Calculation

While the contour method for residual stress assessment has developed rapidly, no published study documents its interlaboratory reproducibility. Objective: Here we report an initial reproducibility experiment focused on contour method data analysis and residual stress calculation. The experiment uses surface topography data from a physical process simulation of elastic-plastic beam bending. The simulation provides surface topography, for input to the contour method data analysis, as well as a known residual stress field with 130 MPa peak magnitude. Additionally, to increase realism, noise and specific artifacts are added to the topography data. A group of participants received the topography data (without the known residual stress), independently analyzed the data, and submitted results as a two-dimensional residual stress field. Analysis of submissions provides a group average residual stress field and the spatial distribution of reproducibility standard deviation. The group average residual stress agrees with the known stress in magnitude and spatial trend. The reproducibility standard deviation ranges from 2 to 54 MPa over the measurement plane, with an average of 5.4 MPa. Reproducibility standard deviation is smaller in the cross-section interior (≤ 5 MPa), modest near local extrema in the stress field (5 to 10 MPa), and larger near the cross-section boundaries (10 to 30 MPa). Overall, the largest values of reproducibility standard deviation (up to 54 MPa) occur in limited areas where artifacts had been added to the topography data; while some participants identified and removed these artifacts, some did not, leading to systematic differences that elevated the standard deviation.

42 ENGINEERING↗

DOE OSTI2021

Emerging contact force heterogeneity in ordered soft granular media

Under external perturbations, inter-particle forces in disordered granular media are well known to form a heterogeneous distribution with filamentary patterns. Better understanding these forces and the distribution is important for predicting the collective behavior of granular media, the media second only to water as the most manipulated material in global industry. However, studies in this regard so far have been largely confined to granular media exhibiting only geometric heterogeneity, leaving the dimension of mechanical heterogeneity a rather uncharted area. In this work, through a FEM contact mechanics model, we show that a heterogeneous inter-particle force distribution can also emerge from the dimension of mechanical heterogeneity alone. Specifically, we numerically study inter-particle forces in packing of mechanically heterogeneous disks arranged over either a square or a hexagonal lattice and under quasi-static isotropic compression. Our results show that, at the system scale, a hexagonal packing exhibit a more heterogeneous inter-particle force distribution than a square packing does; At the particle scale, for both packing lattices, preliminary analysis shows the consistent coexistence of outliers (i.e., softer disks sustaining larger forces while stiffer disks sustaining smaller forces) in comparison to their homogeneous counterparts, which implies the existence of nonlocal effect. Further analysis on the portion of outliers and on spatial contact force correlations suggest that the hexagonal packing shows more pronounced nonlocal effect over the square packing under small mechanical heterogeneity. However, such trend is reversed when assemblies becomes more mechanically heterogeneous. Lastly, we confirm that, in the absence of particle reorganization events, contact friction merely plays the role of packing stabilization while its variation has little effect on inter-particle forces and their distribution.

36 MATERIALS SCIENCE↗