Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed method”

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Microtomographic Measurements of Total Air‐Water Interfacial Areas for Soils

Abstract Synchrotron X‐ray microtomography (XMT) was used to measure total air‐water interfacial areas (A aw ) as a function of water saturation (S w ) for several soils that comprise a range of physical and geochemical properties. Measurements were also conducted for glass beads and quartz sands for comparison. Apparent near‐linearA aw ‐S w relationships are observed for the three sands and the three sandy soils. In contrast, the measured interfacial areas for two soils that contain greater proportions of silt and clay are strongly nonlinear functions of water saturation. The greater degree of nonlinearity observed for these two soils is due to their much greater particle‐size distributions (i.e., uniformity coefficients) and their concomitant greater range in pore sizes. Interfacial areas determined with the thermodynamic method were used to benchmark the XMT measurements. XMT‐measured interfacial areas compare well to the thermodynamic‐determined values for the sands and sandy soils. In contrast, the XMT‐measured interfacial areas for the two soils with larger particle‐size distributions are not fully congruent with the thermodynamic‐determined values. Both of these soils have large fractions of pore space comprising nominal pore diameters smaller than the resolution of the XMT imaging. These results suggest that air‐water interfacial area may not always be fully characterized by standard XMT for soils with large particle‐size distributions.

Environmental Sciences & Ecology↗
DOE OSTI2024

Updimensioning strategy derived from synthetic equiaxed grain structures for approximating 3D grain size distributions from 2D visualizations with 1D parameters

We generated synthetic equiaxed grain structures using computer graphics software to explore the relationship between various grain size determination methods and true three-dimensional (3D) grain diameters. Mirroring grain measurement techniques, the synthetic 3D grain structures are imaged as 2D micrographs which are measured to yield 1D grain size parameters. Synthetic grain structures provide data at a mass scale and permit exploration of both polished and fractured surface micrographs, revealing one-to-one correspondence between exposed 2D grain cross-sections and individual 3D grains. Analysis of this correspondence yielded a procedure to approximate 3D equiaxed grain size and volume distributions based on the mode of the 2D fractograph grain size distribution. The 3D approximation procedure is shown to be less susceptible to different imaging conditions that affect small, undiscernible grains compared to the standard planimetric and linear intercept methods, which by design also tend to underestimate the 3D grain diameter. The procedure requires larger sample sizes to lower variance and a deeper analysis which could become more practical with machine learning (ML) models for grain boundary segmentation, which synthetic grain structures can help train. This work lays the foundation for analyzing other grain distributions such as columnar and composite grains in similar depth.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Distributed Conditions for Small-signal Stability of Power Grids and Local Control Design

Operating modern power grids with stability guarantees is markedly important. Typical methods for analyzing and certifying power grid stability are largely centralized relying on the ability of the system operator to gather network-wide information and accurately compute the system's eigenvalues. These methods are oftentimes not privacy-preserving and computationally burdensome. They are therefore, not well-suited to modern power grids where small-signal stability has to be evaluated timely, efficiently and in a privacy-preserving fashion. Herein, we introduce a distributed methodology for certifying small-signal stability of power grids and designing the local controllers. First, we analytically derive distributed conditions for network-wide stability that bus agents can inspect using local information. By leveraging these conditions, we then introduce a distributed control design algorithm (DCDA) that can guide the local control design so that stability of the interconnected system is guaranteed. The agents that adopt the proposed distributed algorithm are responsible for tuning their local controllers, producing their local control commands and ensuring that their local stability condition is met. The system operator is only responsible for verifying network-wide stability upon receiving affirmative responses from all agents and, announcing, that the overall system is stable. The proposed DCDA algorithm is numerically validated via simulations using the IEEE 39-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Designing an Intrusion Detection for an Adjustable Speed Drive System Controlling a Critical Process

In this article, we address the cyber-security problem of industrial control systems (ICSs) when their sensor measurements may be compromised due to an attacker who has intercepted those measurements via a network. We introduce a general-purpose method “Dynamic Watermarking (DW)” to detect potential cyber-intrusions on speed sensor measurements within industrial control systems, which deploy an adjustable speed drive (ASD) to control a critical process. The DW method is injecting a random private low-amplitude signal with a zero mean Gaussian distribution, “watermark”, into one of the input phase voltages powering the ASD system. The watermark signal propagates through the system including pulse width modulation (PWM) power conversion stage and motor, then ultimately appears in the speed sensor measurements. By deploying two statistical DW tests with two proper thresholds, the system can detect potential cyber-intrusions or unobservable cyber-attacks such as replay attacks and false data injection attacks (FDIA). The DW method tested on a laboratory-scale ASD system experimentally to protect the system against cyber-intrusions. This system, powered by a commercial PWM drive operating at 208 V, 3-phase, and 3.7 kW, served as our experimental platform.

42 ENGINEERING↗
DOE OSTI2021

Review of Dew Point Evaporative Cooling Technology for Air Conditioning Applications

Indirect evaporative cooling has the potential to significantly improve the natural environment. It follows from a significant reduction in electricity consumption in the hot period, and hence lower operating costs for cooling systems. This paper presents the current state of knowledge and research directions on dew point indirect evaporative cooling. It was found that researchers focus on the development of dew point indirect evaporative coolers (DPIEC) by improving its design, geometry, water distribution, and new porous materials implementation. To evaluate the performance of new types of DPIEC, different methods are used by the scientists. Finally, optimized devices are studied in terms of their performance in different systems, like hybrid and desiccant systems, considering different climate conditions. Potential directions of development of evaporative technologies were indicated, such as increasing the coefficient of performance of solid desiccant evaporative cooling systems, developing novel geometry, and efficient water distribution, including development of porous materials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Cumulative short-lived photofission product yields for nuclear forensic application

Experiments have long been carried out with linear accelerator (linac) produced bremsstrahlung X-rays to study high-energy photon-induced fission, known as “photofission”. One benefit of photofission, when compared with neutron-induced fission, is the ability to investigate fundamental physics of fission at excitation energies lower than what is possible with neutrons [1]. Absorption of an incident neutron results in a compound nucleus with an excitation energy at least equal to the neutron binding energy; thus, the lowest energy region of the fission barrier cannot be studied. Fission product yield distributions are dependent on the incident particle energy and the target nucleus mass. They are generally represented by a double-humped curve, with low and high mass “peaks” and an intermediate mass “valley”. As incident particle energy increases, the excitation energy of the nucleus increases, thereby increasing the overall yield of isotopes in the valley with respect to the two maxima or peaks of the fission product curve [3]. Photofission is most probable for incident photons around 14 MeV due to the giant dipole resonance. Induced photofission around this energy results in a fission product yield distribution that resembles 14 MeV DT neutron fission more than thermal or fast (500 keV) neutron fission. In post-detonation nuclear forensics, the peak-to-valley ratio of the fission product yield curve could be used for characterizing device type by providing an indication of the neutron energy spectrum that had induced fission (e.g., 500 keV vs. 14 MeV). The goal of this work is to improve nuclear data for isotopes of interest in nuclear forensics and special nuclear material (SNM) detection, specifically the fission product yield data via photofission experiments for the isotopes of 238U and 232Th. Photofission has been proposed as a less costly and more flexible production method for the valley isotopes that are representative of DT neutron fission product yield distributions, which would benefit nuclear forensics exercises in testing for analysis techniques. Natural uranium and thorium targets were irradiated with bremsstrahlung photons (nominal endpoint energies of 8, 14, and 20 MeV), generated via impingement of high-energy electrons on a tungsten electron-to-photon radiator. Gamma-ray spectroscopy was used to observe and identify the decay of the resulting fission products. The photofission yields of 102 fission products were measured for the set of the 238U and 232Th targets at the three bremsstrahlung X-ray endpoint energies 8, 14, and 20 MeV. These fragments included masses ranging from A=84 to A=144, with half-lives as short as 1.07 seconds.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2021

Canopy position influences the degree of light suppression of leaf respiration in abundant tree genera in the Amazon Forest

Leaf respiration in the dark (R dark ) and light (R day ) is poorly characterized in diverse tropical ecosystems, and little to no information exists on the degree of light suppression in common tree species within the Amazon basin, and their dependences upon plant functional traits and position within the canopy. We quantified R dark and apparent R day using the Kok method and measured key leaf traits in 26 tree individuals of different species distributed in three different canopy positions: canopy, lower canopy, and understory. To explore the relationships between the leaf traits we used the standardized major axis (SMA). We found that canopy trees had significantly higher rates of R dark and R day than trees in the understory. The difference between R dark and R day (the light suppression of respiration) was greatest in the understory (68 ± 9%, 95% CI) and lower canopy (49 ± 9%, 95% CI) when compared to the canopy (37 ± 10%, 95% CI). We also found that R day was significantly and strongly correlated with R dark (p < 0.001) for all the canopy positions. Also, leaf mass per area (LMA) and leaf Phosphorus concentration (P) had a significant relationship with R dark (p < 0.001; p = 0.003), respectively. In addition, a significant relationship was found for LMA in the canopy and lower canopy positions (p = 0.009; p = 0.048) while P was only significant in the canopy (p = 0.044). Finally, no significant relationship was found between R dark and nitrogen, sugars, and starch. Our results highlight the importance of including representation of the light suppression of leaf respiration in terrestrial biosphere models and also of accounting for vertical gradients within forest canopies and connections with functional traits.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Study of Quasielastic Barrier Distributions as a Step towards the Synthesis of Superheavy Elements with Hot Fusion Reactions

The excitation functions for quasielastic scattering of Ne 22 + Cm 248 , Mg 26 + Cm 248 , and Ca 48 + U 238 are measured using a gas-filled recoil ion separator. The quasielastic barrier distributions are extracted for these systems and are compared with coupled-channel calculations. The results indicate that the barrier distribution is affected dominantly by deformation of the actinide target nuclei, but also by vibrational or rotational excitations of the projectile nuclei, as well as neutron transfer processes before capture. From a comparison between the experimental barrier distributions and the evaporation residue cross sections for Sg ( Z = 106 ), Hs (108), Cn (112), and Lv (116), it is suggested that the hot fusion reactions take advantage of a compact collision, where the projectile approaches along the short axis of a prolately deformed nucleus. Overall, a new method is proposed to estimate the optimum incident energy to synthesize unknown superheavy nuclei using the barrier distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of cyber attacks on distributed compressive sensing based state estimation in power distribution grids

Modern power distribution grids suffer from multiple vulnerabilities due to the tight integration between the physical system and the cyber infrastructure. Sophisticated and malicious cyber attacks continue to adversely impact the grid operation leading to performance degradation, service interruption, and grid failure. State estimation plays an essential role in grid monitoring and advancing cyber-attack situational awareness. In this regard, this paper first proposes a distributed compressive sensing (CS) state estimation approach for an unobservable distribution grid. Further, the proposed distributed CS approach divides the distribution grid into sub-areas to perform local state estimation. Then an alternating direction method of multipliers (ADMM) based iterative information exchange among neighboring areas is employed to complete the estimation process. In this estimation process, the impact of loss of measurement data, false data injection (FDI), replay, and neighborhood cyber-attacks is analyzed. Extensive simulations are performed on the IEEE 37-bus and IEEE 123-bus standard networks to demonstrate the algorithm’s robustness to the aforementioned cyber-attacks. A quantitative analysis of computational complexity and simulation time of the distributed CS based approach is also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

An Adaptive-Importance-Sampling-Enhanced Bayesian Approach for Topology Estimation in an Unbalanced Power Distribution System

The reliable operation of a power distribution system relies on a good prior knowledge of its topology and its system state. Although crucial, due to the lack of direct monitoring devices on the switch statuses, the topology information is often unavailable or outdated for the distribution system operators for real-time applications. Apart from the limited observability of the power distribution system, other challenges are the nonlinearity of the model, the complicated, unbalanced structure of the distribution system, and the scale of the system. To overcome the above challenges, we, in this paper, propose a Bayesian-inference framework that allows us to simultaneously estimate the topology and the state of a three-phase, unbalanced power distribution system. Specifically, by using the very limited number of measurements available that are associated with the forecast load data, we efficiently recover the full Bayesian posterior distributions of the system topology under both normal and outage operation conditions. This is performed through an adaptive importance sampling procedure that greatly alleviates the computational burden of the traditional Monte-Carlo (MC)-sampling-based approach while maintaining a good estimation accuracy. The simulations conducted on the IEEE 123-bus test system and an unbalanced 1282-bus system reveal the excellent performances of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

CoSHA: Code for Stellar Properties Heuristic Assignment—for the MaStar Stellar Library

We introduce CoSHA: a Code for Stellar properties Heuristic Assignment. In order to estimate the stellar properties, CoSHA implements a Gradient Tree Boosting algorithm to label each star across the parameter space (T eff , $\mathrm{log}g$, [Fe/H], and [α/Fe]). We use CoSHA to estimate the stellar atmospheric parameters of 22,000 unique stars in the MaNGA Stellar Library (MaStar). To quantify the reliability of our approach, we run internal tests, using both the Göttingen Stellar Library (a theoretical library) and the first data release of MaStar, and external tests, by comparing the resulting distributions in the parameter space with the APOGEE estimates of the same properties. In summary, our parameter estimates span the ranges T eff = [2900, 12,000] K, $\mathrm{log}g=[-0.5,5.6]$, [Fe/H] = [-3.74, 0.81], and αM = [-0.22, 1.17]. We report internal (external) uncertainties of the properties of ${\sigma }_{{T}_{\mathrm{eff}}}\sim 43(240)$ K, ${\sigma }_{\mathrm{log}g}\sim 0.2(0.4)$, σ [Fe/H] ~ 0.16(0.24), and σ [α/Fe] ~ 0.09(0.08). These uncertainties are comparable to those of other methods with similar objectives. Despite the fact that CoSHA is not aware of the spatial distributions of these physical properties in the Milky Way, we are able to recover the main trends known in the literature. The catalog of physical properties for MaStar can be accessed online.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Hybrid Gradient Method to Designing Bayesian Experiments for Implicit Models

Bayesian experimental design (BED) aims at designing an experiment to maximize the information gathering from the collected data. The optimal design is usually achieved by maximizing the mutual information (MI) between the data and the model parameters. When the analytical expression of the MI is unavailable, e.g.,having implicit models with intractable data distributions, a neural network-based lower bound of the MI was recently proposed and a gradient ascent method was used to maximize the lower bound [1]. However, the approach in [1] requires a pathwise sampling path to compute the gradient of the MI lower bound with respect to the design variables, and such a pathwise sampling path is usually inaccessible for implicit models. In this work, we propose a hybrid gradient approach that leverages recent advances in variational MI estimator and evolution strategies (ES)combined with black-box stochastic gradient ascent (SGA) to maximize the MI lower bound. This allows the design process to be achieved through a unified scalable procedure for implicit models without sampling path gradients. Several experiments demonstrate that our approach significantly improves the scalability of BED for implicit models in high-dimensional design space.

Zhang, Jiaxin↗
DOE OSTI2024

Proton and neutron density distributions at supranormal density in low- and medium-energy heavy-ion collisions. II. Central Pb + Pb collisions

This paper represents a continuation of our investigation into the limitations on total, proton, and neutron particle number densities, as well as the asymmetry of proton and neutron density distributions achievable in central heavy-ion collisions. We explore these aspects at low and medium energies within the framework of the Boltzmann-Uhlenbeck-Uehling (pBUU) transport and time-dependent Hartree-Fock (TDHF) models. Previous studies, focusing on symmetric and asymmetric collisions of Ca and Sn nuclei [1], and initial results on Pb-nuclei collisions [2], emphasized the role of the Coulomb interaction in these events. Our findings indicated that: (i) the highest total densities predicted at 𝐸 beam = 800 MeV/nucleon were on the order of ≈ 2.5⁢𝜌 0 (𝜌 0 = 0.16 fm −3 ), (ii) the proton-neutron asymmetry for maximal densities, 𝛿=(𝜌$^{max}_{n}$− 𝜌$^{max}_{p}$)/(𝜌$^{max}_{n}$+𝜌$^{max}_{p}$), did not generally exceed the asymmetry in the initial state of the collision at all beam energies and tended to decrease during the reaction, and (iii) a significant portion of this asymmetry had its microscopic origin in Coulomb forces, masking the pure nuclear contribution. These new findings, particularly relevant in the astrophysical context, are further examined in this work, focusing on the heaviest target-projectile combination 212,208 Pb accessible in an experiment. We introduce the SkT3 Skyrme force model, not previously used for the Pb system, and compare it to the SV-bas and SV-sym34 models to explore the symmetry-energy dependence of the results. Contour plots of nucleonic densities are presented, contrasting the time evolution of the density distributions in low (TDHF) and high (pBUU) models. We also present the evolution of normalized maximal proton, neutron, and total nucleon number density with increasing beam energy in the full pBUU model and the Vlasov approximation, aiming to explore the impact of correlations in the reaction. The time evolution of the proton and neutron density distributions in the plane transverse to the beam direction is illustrated at both low and high beam energy. In conclusion, our detailed examination of the Pb system in this work provides further essential evidence that the aforementioned findings (i)–(iii) are only weakly dependent on system size and a symmetry-energy model, and thus, they are of more general importance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Standardizing Performance Metrics for Building-Level Electrical Distribution Systems

Building-level electrical distribution systems comprise a myriad of current-carrying equipment, conversion devices, and protection devices that deliver power from the utility or local distributed energy resources to end-use building loads. Electric power has traditionally been generated, transmitted, and distributed in alternating current (AC). However, the last decade has seen a significant increase in the integration of native direct current (DC) equipment that has elevated the importance of DC distribution systems. Numerous studies have comparatively examined the performance of various electrical distribution systems in buildings but have failed to achieve uniform conclusions, primarily because of a lack of consistent and analogous performance evaluation methods. This paper aims to fill this gap by providing a standard set of metrics and measurement boundaries to consistently evaluate the performance of AC, DC, or hybrid AC/DC electrical distribution systems. The efficacy of the proposed approach is evaluated on a representative medium-sized commercial office building model with AC distribution and an equivalent hybrid AC/DC and DC distribution model, wherein the AC distribution model is concluded to be the most efficient. The simulation results show variation in computed metrics with different selected boundaries that verify the effectiveness of the proposed approach in ensuring consistent computation of the performance of building-level electrical distribution systems. This paper provides an initial set of guidelines for building energy system stakeholders to adopt appropriate solutions, thus leading to more efficient energy systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Signal-preserving CMB component separation with machine learning

Analysis of microwave sky signals, such as the cosmic microwave background, often requires component separation using multifrequency methods, whereby different signals are isolated according to their different frequency behaviors. Many so-called blind methods, such as the internal linear combination (ILC), make minimal assumptions about the spatial distribution of the signal or contaminants, and only assume knowledge of the frequency dependence of the signal. The ILC produces a minimum-variance linear combination of the measured frequency maps. In the case of Gaussian, statistically isotropic fields, this is the optimal linear combination, as the variance is the only statistic of interest. However, in many cases the signal we wish to isolate, or the foregrounds we wish to remove, are non-Gaussian and/or statistically anisotropic (in particular for the case of Galactic foregrounds). In such cases, it is possible that machine learning (ML) techniques can be used to exploit the non-Gaussian features of the foregrounds and thereby improve component separation. However, many ML techniques require the use of complex, difficult-to-interpret operations on the data. We propose a hybrid method whereby we train an ML model using only combinations of the data that , and combine the resulting ML-predicted foreground estimate with the ILC solution to reduce the error from the ILC. We demonstrate our methods on simulations of extragalactic temperature and Galactic polarization foregrounds and show that our ML model can exploit non-Gaussian features, such as point sources and spatially varying spectral indices, to produce lower-variance maps than ILC—e.g., reducing the variance of the B-mode residual by factors of up to 5—while preserving the signal of interest in an unbiased manner. Moreover, we often find improved performance even when applying our ML technique to foreground models on which it was not trained. Published by the American Physical Society 2025

McCarthy, Fiona (ORCID:0000000253893565)↗
DOE OSTI2024

Impact of Floating PV (FPV) System on Water Evaporation: Case Study of an Open Loop Reservoir in North Carolina

Fresh water scarcity is becoming one of the important concerns worldwide. This situation of water scarcity is further exacerbated by the increase in water evaporation due to climate change. This manuscript will address this concern by considering the impact of floating photovoltaics (FPVs) systems on water evaporation. The installation of PV modules on the water bodies can help reduce the water evaporation rate from reservoirs by reducing the evaporating area of the water surface. This work evaluates the effect of FPV on water evaporation rates in the TuckerTown reservoir located in North Carolina state of United States. Meteorological data, such as mean temperature and dew-point, are used to calculate the rate of the water evaporation by Linacre method. The simulation results for the TuckerTown case study showed that the aggregated avoided evaporated volume for one year is 9 million gallons.

14 SOLAR ENERGY↗
DOE OSTI2020

Unsaturated and Saturated Hydraulic Parameters of Soils at the East River Watershed, Colorado Calculated using the Pedotransfer Function Method.

The purpose of the data package is to present the results of calculations of the following unsaturated-saturated hydraulic parameters of soils of the East River Watershed, Colorado: fully saturated water content (Qs), residual saturation (Qr), van genuchten parameters alpha and n, field capacity (FC), and saturated hydraulic conductivity (Ks). Calculations were conducted for the following soil texture types: sandy_loam, loam, silty_clay_loam, clay_loam, silt loam; the particle size distribution included the percentage of sand, silt, clay, and the bulk density (BD) of soils. Calculations were conducted for the following areas- Lower_Subalpine (6 samples),- Lower_Montane_lower_hillslope (7 samples),- Floodplain_at_base_of_lower_montane_hillslope, (11 samples),- PLM-0 (10 samples),- PLM-1 (10 samples), and - PLM-3 (13 samples)Calculations were conducted using the Pedotransfer Function Method developed by Tóth et al. (2014) with the application of the code ”euptf” in Rstudio. The data are being summarized in the Excel file ”Unsaturated_Saturated_Hydraulic_Parameters_PTF_calculated_EASt_EIver.xls" including two worksheets:“full_data_EAST_RIVER” which includes (a) the input data used for calculations, the soil texture triangle, and the plots of the water retention curves and unsaturated hydraulic conductivity for several types of soil texture, and (b) the summary table with the averaged parameters for each area. The the texture triangle was plotted using the Rstudio "soiltexture." Microsoft Excel software is required to open the Excel file. The metadata and a list of references to the Rstudio libraries used for calculations and data visualization are giving in the Word Readme file.

54 ENVIRONMENTAL SCIENCES↗