Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantitative structure”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Unveiling the purely young star formation history of the SMC’s northeastern shell from colour–magnitude diagram fitting

ABSTRACT We obtain a quantitative star formation history (SFH) of a shell-like structure (‘shell’) located in the northeastern part of the Small Magellanic Cloud (SMC). We use the Survey of the MAgellanic Stellar History to derive colour–magnitude diagrams (CMDs), reaching below the oldest main-sequence turnoff, from which we compute the SFHs with CMD-fitting techniques. We present, for the first time, a novel technique that uses red clump (RC) stars from the CMDs to assess and account for the SMC’s line-of-sight depth effect present during the SFH derivation. We find that accounting for this effect recovers a more accurate SFH. We quantify an $\sim$7 kpc line-of-sight depth present in the CMDs, in good agreement with depth estimates from RC stars in the northeastern SMC. By isolating the stellar content of the northeastern shell and incorporating the line-of-sight depth into our calculations, we obtain an unprecedentedly detailed SFH. We find that the northeastern shell is primarily composed of stars younger than $\sim$500 Myr, with significant star formation enhancements around $\sim$250 and $\sim$450 Myr. These young stars are the main contributors to the shell’s structure. We show synchronicity between the northeastern shell’s SFH with the Large Magellanic Cloud’s (LMC) northern arm, which we attribute to the interaction history of the SMC with the LMC and the Milky Way (MW) over the past $\sim$500 Myr. Our results highlight the complex interplay of ram pressure stripping and the influence of the MW’s circumgalactic medium in shaping the SMC’s northeastern shell.

Sakowska, Joanna D. (ORCID:0000000215941466)↗
DOE OSTI2021

Anomalously strong near-neighbor attraction in doped 1D cuprate chains

In the cuprates, one-dimensional (1D) chain compounds provide a distinctive opportunity to understand the microscopic physics, owing to the availability of reliable theories. However, progress has been limited by the challenge of controllably doping these materials. We report the synthesis and spectroscopic analysis of the 1D cuprate Ba 2-x Sr x CuO 3+δ over a wide range of hole doping. Our angle-resolved photoemission experiments reveal the doping evolution of the holon and spinon branches. We identify a prominent folding branch whose intensity fails to match predictions of the simple Hubbard model. An additional strong near-neighbor attraction, which may arise from coupling to phonons, quantitatively explains experiments for all accessible doping levels. Furthermore, considering structural and quantum chemistry similarities among cuprates, this attraction may play a similarly important role in high-temperature cuprate superconductors.

1D Chain↗
DOE OSTI2022

Multi-Defect Detection in Additively Manufactured Lattice Structures Using 3D Electrical Resistance Tomography

Cellular lattice structures possess high strength-to-weight ratios suitable for advanced lightweight engineering applications. However, their quality and mechanical performance can degrade because of defects introduced during manufacturing or in-service. Their complexity and small length scale features make defects difficult to detect using conventional nondestructive evaluation methods. Here we propose a current injection-based method, electrical resistance tomography (ERT), that can be used to detect damaged struts in conductive cellular lattice structures with their intrinsic electromechanical properties. The reconstructed conductivity distributions from ERT can reveal the severity and location of damaged struts without having to probe each strut. However, the low central sensitivity of ERT may result in image artifacts and inaccurate localization of damaged struts. To address this issue, this study introduces an absolute, high throughput, conductivity reconstruction algorithm for 3D ERT. The algorithm incorporates a strut-based normalized sensitivity map to compensate for lower interior sensitivity and suppresses reconstruction artifacts. Numerical simulations and experiments on fabricated representative cellular lattice structures were performed to verify the ability of ERT to quantitatively identify single and multiple damaged struts. The improved performance of this method compared with classical ERT was observed, based on greatly decreased imaging and reconstructed value errors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enabling Robust Exoplanet Atmospheric Retrievals with Gaussian Processes

Atmospheric retrievals are essential tools for interpreting exoplanet transmission and eclipse spectra, enabling quantitative constraints on the chemical composition, aerosol properties, and thermal structure of planetary atmospheres. The James Webb Space Telescope (JWST) offers unprecedented spectral precision, resolution, and wavelength coverage, unlocking transformative insights into the formation, evolution, climate, and potential habitability of planetary systems. However, this opportunity is accompanied by challenges: modeling assumptions and unaccounted-for noise or signal sources can bias retrieval outcomes and their interpretation. To address these limitations, we introduce a Gaussian process (GP)-aided atmospheric retrieval framework that flexibly accounts for unmodeled features and correlated noise in exoplanet spectra. We validate this method on synthetic JWST observations, and show that GP-aided retrievals reduce bias in inferred abundances and better capture model–data mismatches than traditional approaches. We also introduce the concept of mean squared error to quantify the trade-off between bias and variance, arguing that this metric more accurately reflects retrieval performance than bias alone. We then reanalyze the NIRISS/SOSS JWST transmission spectrum of WASP-96 b, finding that GP-aided retrievals yield broader constraints on CO 2 and H 2 O, possibly alleviating tension between previous retrieval results and equilibrium predictions. Our GP framework provides precise and accurate constraints while highlighting regions where models fail to explain the data. As JWST matures and future facilities come online, a deeper understanding of the limitations of both data and models will be essential, and GP-enabled retrievals like the one presented here offer a principled path forward.

Rotman, Yoav [Arizona State Univ., Tempe, AZ (Unit↗
DOE OSTI2020

Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data

Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-Resolution Comonomer Sequencing of Blocky Brominated Syndiotactic Polystyrene Copolymers Using 13C NMR Spectroscopy and Computer Simulations

This work demonstrates the first high-resolution comonomer sequencing of Blocky brominated syndiotactic polystyrene (sPS-co-sPS-Br) copolymers based on pentad assign- ments of the quaternary carbon region of the nuclear magnetic resonance spectrum. Copolymers containing p-bromostyrene (Br-Sty) units were prepared in matched sets using postpolymerization bromination methods carried out in the heterogeneous gel state (Blocky) and homogeneous solution state (Random). Quantitative information from the quaternary carbon spectra, heteronuclear multiple bond correlation spectroscopy, electronic structure calculations, and simulated statistically random copolymers was correlated to confirm the carbon resonance assignments for all 20 possible pentad comonomer sequences. Using the experimental pentad sequence prevalences, a computer code was developed to simulate chains with microstructures typical of each sample as a means to visually represent the copolymer blockiness with quantitative precision. Based on the microstructure and distribution of run lengths in these chains, the simulations revealed that the Blocky copolymers contain a high degree of blockiness. By comparing the run lengths in the simulated chains to the average number of styrene units in a crystalline segment of sPS (found by small-angle X-ray scattering), copolymer crystallizability was predicted. For the simulated Blocky B-21% (21 mol % Br-Sty) chain, the probability of randomly selecting a styrene unit in a crystallizable block was 25.8%, while that in the simulated Random R-18% was zero, in excellent agreement with the experimental crystallization behavior measured by differential scanning calorimetry. Additionally, these predictions confirmed that the simulated chains accurately represent the ensemble of chains in their respective copolymer samples. Furthermore, each simulated Blocky chain contained one or more long sPS blocks that paralleled the measured 38-40 styrene units spanning a crystalline segment within the sPS/CCl4 gel. This finding affirmed that the long sPS segments originated from the precise lamellar structure within the heterogeneous gel morphology (i.e., block length is correlated with lamellar thickness). Overall, the ability to tailor the copolymer microstructure through control of the semicrystalline gel morphology opens the door to synthesizing ordered copolymers by postpolymerization functionalization processes with unprecedented levels of compositional control.

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

Intermolecular correlations of liquid and glassy CS 2 studied by synchrotron radiation x-ray diffraction

How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. Here, we have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS 2 , which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16–25.7 Å -1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS 2 were also recorded in a wide temperature range of 160–300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS 2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .

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

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

State-dependent motion of a genetically encoded fluorescent biosensor

Genetically encoded biosensors can measure biochemical properties such as small-molecule concentrations with single-cell resolution, even in vivo. Despite their utility, these sensors are “black boxes”: Very little is known about the structures of their low- and high-fluorescence states or what features are required to transition between them. We used LiLac, a lactate biosensor with a quantitative fluorescence-lifetime readout, as a model system to address these questions. X-ray crystal structures and engineered high-affinity metal bridges demonstrate that LiLac exhibits a large interdomain twist motion that pulls the fluorescent protein away from a “sealed,” high-lifetime state in the absence of lactate to a “cracked,” low-lifetime state in its presence. Understanding the structures and dynamics of LiLac will help to think about and engineer other fluorescent biosensors.

Rosen, Paul C. (ORCID:000000017414454X)↗
DOE OSTI2020

Machine learning bridges local static structure with multiple properties in metallic glasses

A long-standing challenge in the metallic glass (MG) community has been how to quantitatively gauge the influence of the intricate local packing environment on the response (such as the propensity for atomic rearrangement) of the atomic configuration to external stimuli. Here we establish this structure–property relation by representing the complex amorphous structure using a single, flexibility-orientated structural quantity. This structural flexibility (SF) couples to a bona fide structural representation, the pair distribution function (PDF) of individual atoms, through a weighting function that reflects what matters in the static atomic configuration to dynamic responses. Machine learning is used, employing microscopic flexibility volume as the supervisory signal, to establish via direct regression an optimized weighting vector, which is proven robust for all quenching rates, deformation conditions, and different compositions in a given (e.g., Cu x Zr 100-x ) alloy system. Additionally, the SF is evaluated solely from the particle positions (PDF), for any structure variation, from the atomic scale up to sample average. Strong correlations are demonstrated between SF and a broad range of properties, including vibrational, diffusional, as well as elastic and plastic relaxation responses.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Using mixed methods to construct and analyze a participatory agent-based model of a complex Zimbabwean agro-pastoral system

Complex social-ecological systems can be difficult to study and manage. Simulation models can facilitate exploration of system behavior under novel conditions, and participatory modeling can involve stakeholders in developing appropriate management processes. Participatory modeling already typically involves qualitative structural validation of models with stakeholders, but with increased data and more sophisticated models, quantitative behavioral validation may be possible as well. In this study, we created a novel agent-basedmodel applied to a specific context: Zimbabwean non-governmental organization the Muonde Trust has been collecting data on their agro-pastoral system for the last 35 years and had concerns about land-use planning and the effectiveness of management interventions in the face of climate change. We collaboratively created an agent-based model of their system using their data archive, qualitatively calibrating it to the observed behavior of the real system without tuning any parameters to match specific quantitative outputs. We then behaviorally validated the model using quantitative community-based data and conducted a sensitivity analysis to determine the relative impact of underlying parameter assumptions, Indigenous management interventions, and different rainfall variation scenarios. We found that our process resulted in a model which was successfully structurally validated and sufficiently realistic to be useful for Muonde researchers as a discussion tool. The model was inconsistently behaviorally validated, however, with some model variables matching field data better than others. We observed increased model system instability due to increasing variability in underlying drivers (rainfall), and also due to management interventions that broke feedbacks between the components of the system. Interventions that smoothed year-to-year variation rather than exaggerating it tended to improve sustainability. The Muonde trust has used the model to successfully advocate to local leaders for changes in land-use planning policy that will increase the sustainability of their system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Bottom-Up Simulation, Reconstruction, and Quantification of Macromolecule Sequences from Experimental Polymerizations

Motivated by the canonical sequence–structure–function paradigm, tools to characterize chemical patterning in natural biomacromolecules, from proteins to nucleic acids, have grown exponentially in recent years. However, analogous strategies for synthetic macromolecules remain in nascent stages, complicated by sequence polydispersity and analytical limitations. To address this, we have developed a comprehensive and open-source Python package, PRISM (polymer rate insights and sequence modeling), an end-to-end workflow that provides a path from experimental kinetics measurements to quantitative and qualitative metrics for describing chemical patterning in stochastic polymers. First, a numerical integration strategy was constructed to simulate and fit experimental data from reversible addition–fragmentation chain transfer (RAFT) polymerization kinetics, enabling the facile estimation of relevant reactivity ratios. These ratios were then used in a mechanism-specific stochastic kinetic simulation strategy to simulate sequence ensembles corresponding to model systems spanning experimental copolymers, classes of statistical polymers (e.g., alternating, block, and gradient), and multiblock copolymers. Lastly, inspired by sequence homology metrics from bioinformatics, we introduce visualization strategies and quantitative metrics to facilitate comparisons of different sequence ensembles. As the sequence–structure–function paradigm becomes increasingly central in de novo design of synthetic macromolecules, this toolkit provides a first step toward accurate and representative sequence description and featurization.

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

Scalable method for micro-CT analysis enables large scale quantitative characterization of brain lesions and implants

Anatomic evaluation is an important aspect of many studies in neuroscience; however, it often lacks information about the three-dimensional structure of the brain. Micro-CT imaging provides an excellent, nondestructive, method for the evaluation of brain structure, but current applications to neurophysiological or lesion studies require removal of the skull as well as hazardous chemicals, dehydration, or embedding, limiting their scalability and utility. Here we present a protocol using eosin in combination with bone decalcification to enhance contrast in the tissue and then employ monochromatic and propagation phase-contrast micro-CT imaging to enable the imaging of brain structure with the preservation of the surrounding skull. Instead of relying on descriptive, time-consuming, or subjective methods, we develop simple quantitative analyses to map the locations of recording electrodes and to characterize the presence and extent of hippocampal brain lesions.

3-D reconstruction↗
DOE OSTI2024

Towards an understanding of particle-scale flaws and microstructure evolution in cold spray via accumulation of single particle impacts

Cold spray coatings are the sum of countless individual bonding events between single particles impacting on top of one another at high velocities. Thus, the collective behavior of microparticles must be considered to elucidate the origins of coating flaws at the scale of the particles and larger, or the dynamic evolution of the overall coating microstructure. Laser-induced particle impact testing (LIPIT) has been extensively used to study single-particle impacts, and in this work is adapted to study the accumulation of numerous particles with knowledge of each individual particle's impact parameters (particle size, velocity). The method reproducibly deposits stacks of gold particles (>20 particles) with different characteristic spectra of impact velocity. The observation of impact-induced erosion lets us define a critical velocity for material-build-up that is higher than that for single-particle bonding. The quantitative single-particle data are analyzed in a correlative manner to the structure and flaws in the resulting stacks, providing some first statistical connections between, e.g., strain and recrystallization, or aberrant particle characteristics and defects. The results highlight opportunities for the study of many-particle phenomena in microparticle impact—from interaction of particles in cold spray to multi-step erosion processes—with a quantitative view of the behavior of single particles.

Erosion↗
DOE OSTI2023

Helical Organic and Inorganic Polymers

Despite being a staple of synthetic plastics and biomolecules, helical polymers are scarcely studied with Gaussian-basis-set ab initio electron-correlated methods on an equal footing with molecules. This article introduces an ab initio second-order many-body Green’s-function [MBGF(2)] method with nondiagonal, frequency-dependent Dyson self-energy for infinite helical polymers using screw-axis-symmetry-adapted Gaussian-spherical-harmonics basis functions. Together with the Gaussian-basis-set density-functional theory for energies, analytical atomic forces, translational-period force, and helical-angle force, it can compute correlated energy, quasiparticle energy bands, structures, and vibrational frequencies of an infinite helical polymer, which smoothly converge at the corresponding oligomer results. These methods can handle incommensurable structures, which have an infinite translational period and are hard to characterize by any other method, just as efficiently as commensurable structures. We apply them to polyethylene (2/1 helix), polyacetylene (Peierls’ system), and polytetrafluoroethylene (13/6 helix) to establish the quantitative accuracy of MBGF(2)/cc-pVDZ in simulating their (angle-resolved) ultraviolet photoelectron spectra, and of B3LYP/cc-pVDZ or 6-31G** in reproducing their structures, infrared and Raman band positions, phonon dispersions, and (coherent and incoherent) inelastic neutron scattering spectra. We then predict the same properties for infinitely catenated chains of nitrogen or oxygen and discuss their possible metastable existence under ambient conditions. They include planar zigzag polyazene (N2)x (Peierls’ system), 11/3-helical isotactic polyazane (NH)x, 9/4-helical isotactic polyfluoroazane (NF)x, and 7/2-helical polyoxane (O)x as potential high-energy-density materials.

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

Identifying Topological Defects in Lamellar Phases through Contour Analysis of Complex Wave Fields

Lamellar phases frequently contain structural imperfections that significantly affect their behaviors and properties. Our previous research successfully reconstructed real-space configurations of defective lamellar phases from diffuse scattering patterns, indicating the presence of phase vortices as a potential method for identifying topological defects disrupting the smectic ordering. Here, this report presents a mathematical framework using regularized wave fields to represent defective lamellar structures in real space. Phase singularities, resulting from the interference of random waves and indicating lamellar order disruption, are identified through a contour integral. These wave fields, derived from coherent scattering in reciprocal space, were validated via computational benchmarks analyzing small-angle neutron scattering data from AOT surfactant solutions, facilitating further statistical analysis of the defects. Our study highlights the potential to extract meaningful information about topological defects in lyotropic phases by inversely analyzing experimentally measured two-point static correlations. Our method allows for detailed structural analysis of various lyotropic phases, both particulate and nonparticulate, in their quiescent states and facilitates quantitative investigation of defects’ role in phase transitions. By integrating small-angle scattering, deep learning, and vortex tangle analysis, our comprehensive approach shows promise in addressing complex challenges in the structural analysis of soft matter systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Orientation-dependent structural properties during growth and growth mechanism of CoO films

The orientation-dependent local structural properties of CoO films on sapphire substrates during growth were investigated through linearly-polarized extended X-ray absorption fine structure (EXAFS) measurements. Specifically, CoO(1 1 1) and (100) crystals with a rock-salt structure (Fm3m) were epitaxially grown on alpha-Al2O3(0001) and (1012) substrates, respectively, at 700 degrees C using a radio-frequency sputtering system. The local structural properties of CoO films in the in-plane and out-of-plane orientations were quantitatively determined using linearly-polarized EXAFS at the Co K-edge during growth. The EXAFS analysis revealed that during the initial stages of growth, the local structural properties exhibit significant differences compared to thick films, with short atomic distances and large (small) Debye-Waller factors observed in the out-of-plane (in-plane) orientations. The local structural strain mostly diminished when approximately 20 CoO layers accumulated on the substrate. Density functional theory (DFT) calculations further supported these findings by confirming that cobalt atoms initially form stable bonds with the sapphire surface, leading to the simultaneous growth of oxygen and cobalt layers in a coordinated manner through layer-by-layer growth.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A VOI Web Application for Distinct Geothermal Domains: Statistical Evaluation of Different Data Types within the Great Basin

The Great Basin region contains different domains that have different structural and hydrothermal flow patterns. Depending on the characteristics of these patterns, certain data types may be more successful at detecting hidden geothermal resources. In this paper, we quantitatively evaluate if certain data types are more successful in certain domains. Given different aquifer, strain and structural conditions, we explore which data types statistically reveal positively labeled geothermal sites. We utilize value of information (VOI) metrics to help quantify the reliability of data types to discriminate against "positive" and "negative" labeled geothermal sites. We also evaluate how kernel density estimation can help generalize the statistics that inform VOI, which is necessary given the limited data in geothermal exploration. Except for the Carbonate Aquifer, the highest ranking of the Vimperfect is the Local Structural Setting. Next, the slip and dilation tendency is first for Carbonate Aquifer and second for Central Nevada Seismic Belt and Western Great Basin. For the Carbonate Aquifer, heat flow is has the lowest Vimperfect value compared to the other three domains, which is consistent with the understanding of how heat flow measurements are masked by regional groundwater flow.

Bayesian analysis↗