skip to main content

Title: Modeling root system growth around obstacles
Abstract State-of-the-Art models of Root System Architecture (RSA) do not allow simulating root growth around rigid obstacles. Yet, the presence of obstacles can be highly disruptive to the root system. We grew wheat seedlings in sealed petri dishes without obstacle and in custom 3D-printed rhizoboxes containing obstacles. Time-lapse photography was used to reconstruct the wheat root morphology network. We used the reconstructed wheat root network without obstacle to calibrate an RSA model implemented in the R-SWMS software. The root network with obstacles allowed calibrating the parameters of a new function that models the influence of rigid obstacles on wheat root growth. Experimental results show that the presence of a rigid obstacle does not affect the growth rate of the wheat root axes, but that it does influence the root trajectory after the main axis has passed the obstacle. The growth recovery time, i.e. the time for the main root axis to recover its geotropism-driven growth, is proportional to the time during which the main axis grows along the obstacle. Qualitative and quantitative comparisons between experimental and numerical results show that the proposed model successfully simulates wheat RSA growth around obstacles. Our results suggest that wheat roots follow patterns that could more » inspire the design of adaptive engineering flow networks. « less
; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Scientific Reports
Sponsoring Org:
National Science Foundation
More Like this
  1. Introduction The mechanical vulnerability of the atherosclerotic cap is a crucial risk factor in asymptomatic fibroatheromas. Our research group demonstrated using numerical modeling that microcalcifications (µCalcs) located in the fibrous cap can multiply the tissue background stress by a factor 2-7[1-3]. We showed how this effect depends on the size and the ratio of the gap between particles pairs (h) and their diameter (D) along the tensile axis. In this context, we studied the impact of micro-beads of varying diameters and concentration on the rupture of human fibroatheroma laboratory models. Methods We created silicone-based (DowsilEE-3200, Dow Corning) dumbbell-shaped models (80%-scaled ASTM D412-C) of arterial tissues. Samples were divided into three groups: (1) without μBeads (control, n=12), (2) with μBeads of varying diameter (D=30,50,100μm) at a constant concentration of 1% weight (n=36), (3) with μBeads of constant diameter (D=50μm) at different concentrations (3% and 5% weight) (n=24). Before testing, samples were scanned under Micro-CT, at a resolution of 4µm. Images were then reconstructed in NRecon (SkySCan, v.2014) and structural parameters obtained in CTan (SkyScan, v.2014). These data were used to calculate the number of beads and their respective h/D ratio in a custom-made MATLAB script. We tested the samples using amore »custom-made micro material testing system equipped with real-time control and acquisition software (LabVIEW, v. 2018, NI). The reaction force and displacement were measured by the system and images of the sample were recorded by a high-resolution camera. The true stress and strain profiles of each sample were obtained by means of Digital Image Correlation (DIC). Results Samples with and without μBeads exhibited a distinct hyperelastic behaviour typical of arterial tissues (Fig1). Comparison of the mean ultimate stress (UTS) between groups was performed by one-way ANOVA test followed by post-hoc pairwise comparison. Regardless of the group, the presence of μBeads determined a statistically significant reduction in UTS (Fig2). Increasing the μBeads concentration was also positively correlated with lower stresses at rupture as more clusters formed resulting in lower values of h/D (Table1). Discussions Our results clearly capture the influence of μBeads on the rupture threshold of a vascular tissue mimicking material. In fact, samples with μBeads exhibit levels of UTS that are around two times lower than the control group. This effect appears to be dependent on the μBeads proximity, as lower h/D correlates with higher UTS reductions. On the other hand, the effect of particle size is not apparent for the diameters considered in this study. The plausible explanation for the observed change in rupture threshold is the increase in stress concentration around spherical μBeads, which we have previously shown in analytical and numerical studies [1-3]. Our experimental observations support our previous studies suggesting that μCalcs located within the fibroatheroma cap may be responsible for significantly increasing the risk of cap rupture that precedes myocardial infarction and sudden death.« less
  2. Abstract
    Excessive phosphorus (P) applications to croplands can contribute to eutrophication of surface waters through surface runoff and subsurface (leaching) losses. We analyzed leaching losses of total dissolved P (TDP) from no-till corn, hybrid poplar (Populus nigra X P. maximowiczii), switchgrass (Panicum virgatum), miscanthus (Miscanthus giganteus), native grasses, and restored prairie, all planted in 2008 on former cropland in Michigan, USA. All crops except corn (13 kg P ha−1 year−1) were grown without P fertilization. Biomass was harvested at the end of each growing season except for poplar. Soil water at 1.2 m depth was sampled weekly to biweekly for TDP determination during March–November 2009–2016 using tension lysimeters. Soil test P (0–25 cm depth) was measured every autumn. Soil water TDP concentrations were usually below levels where eutrophication of surface waters is frequently observed (> 0.02 mg L−1) but often higher than in deep groundwater or nearby streams and lakes. Rates of P leaching, estimated from measured concentrations and modeled drainage, did not differ statistically among cropping systems across years; 7-year cropping system means ranged from 0.035 to 0.072 kg P ha−1 year−1 with large interannual variation. Leached P was positively related to STP, which decreased over the 7 years in all systems. These results indicate that both P-fertilized and unfertilized cropping systems mayMore>>
  3. This research presents an experimental program executed to understand the strength and stiffness properties of hollow built-up glass compression members that are intended for use in the modular construction of all glass, compression-dominant, shell-type structures. The proposed compression-dominant geometric form has been developed using the methods of form finding and three-dimensional graphical statics. This research takes the first steps towards a new construction methodology for glass structures where individual hollow glass units (HGU) are assembled using an interlocking system to form large, compression-dominant, shell-type structures, thereby exploiting the high compression strength of glass. In this study, an individual HGU has an elongated hexagonal prism shape and consists of two deck plates, two long side plates, and four short side plates, as is shown in Figure 1. Connections between glass plates are made using a two-sided transparent structural adhesive tape. The test matrix includes four HGUs, two each fabricated with 1 mm and 2 mm thick adhesive tape. All samples are dimensioned 64 cm on the long axis of symmetry, 51 cm on the short axis of symmetry, and are 10 cm in width. Glass plates are all 10 mm thick annealed float glass with geometric fabrication done using 5-axis abrasivemore »water jet cutting. HGU assembly is accomplished using 3D printed truing clips and results in a rigid three-dimensional glass frame. Testing was done with the HGU oriented such that load was introduced on the short side edges of the two deck plates, resulting in an asymmetric load-support condition. A soft interface material was used between the HGU and steel plates of the hydraulic actuator and support for the purpose of avoiding premature cracking from local stress concentrations on the glass edges at the load and support locations. Force was applied in displacement control at 0.25 mm/minute with a full array of displacement and strain sensors. Test results for load vs. center deck plate transverse deflection are shown in Figure 2. All samples failed explosively by flexural buckling with no premature cracking on the load and support edges of the deck plates. Strain and deformation data clearly show the presence of second-order behavior resulting from bending deformation perpendicular to the plane of the deck plates. In general, linear axial behavior transitions to nonlinear second-order behavior, with increasing rates in deflection and strain growth, ultimately ending in glass fracture on the tension surfaces of the buckled deck plates. The failure resulted in near-complete disintegration of the deck plates, but with no observable cracking in any side plates and a secure connection on all adhesive tape. Results of the experimental program clearly demonstrate the feasibility of using HGUs for modular construction of compression dominant all-glass shell-type structures. This method of construction can significantly reduce the self-weight of the structure, and it will inspire the use of sustainable materials in the construction of efficient structures.« less
  4. In this work, we explore how the emergence of collective motion in a system of particles is influenced by the structure of their domain. Using the Vicsek model to generate flocking, we simulate two-dimensional systems that are confined based on varying obstacle arrangements. The presence of obstacles alters the topological structure of the domain where collective motion occurs, which, in turn, alters the scaling behavior. We evaluate these trends by considering the scaling exponent and critical noise threshold for the Vicsek model, as well as the associated diffusion properties of the system. We show that obstacles tend to inhibit collective motion by forcing particles to traverse the system based on curved trajectories that reflect the domain topology. Our results highlight key challenges related to the development of a more comprehensive understanding of geometric structure's influence on collective behavior.
  5. mmWave communication in 60GHz band has been recognized as an emerging technology to support various bandwidth-hungry applications in indoor scenarios. To maintain ultra-high throughputs while addressing potential blockage problems for mmWave signals, maintaining line-of-sight (LoS) communications between client devices and access points (APs) is critical. To maximize LoS communications, one approach is to deploy multiple APs in the same room. In this paper, we investigate the optimal placement of multiple APs using both analytical methods and simulations. Considering the uncertainty of obstacles and clients, we focus on two typical indoor settings: random-obstacle-random-client (RORC) scenarios and fixed-obstacle-random-client (FORC) scenarios. In the first case, we analytically derive the optimal positions of APs by solving a thinnest covering problem. This analytical result is used to show that deploying up to 5 APs in a specific room brings substantial performance gains. For the FORC scenario, we propose the shadowing-elimination search (SES) algorithm based on an analytic model to efficiently determine the placement of APs. We show, through simulations, that with only a few APs, the network can achieve blockage-free operation in the presence of multiple obstacles and also demonstrate that the algorithm produces near-optimal deployments. Finally, we perform ns-3 simulations based on the IEEEmore »802.11ad protocol at mmWave frequency to validate our analytical results. The ns-3 results show that proposed multi-AP deployments produce significantly higher aggregate performance as compared to other common AP placements in indoor scenarios.« less