Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large deviations”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Distributed Parareal in Time with Adaptive Coarse Solver for Large Scale Power System Simulations

This paper describes a distributed Parareal in time algorithm with the adaptive coarse solver. In the distributed algorithm, the propagation of a coarse solution is distributed among processors in addition to the fine solution propagation. In the adaptive coarse solver the complexity of the model is changed based on the size of generator rotor angle deviation. In addition, non-iterative approach to model a ZIP load is presented. The proposed method is tested on the large-scale 5617-machine 70285-bus Eastern Interconnection system.

Osipov, Denis↗
DOE OSTI2015

Technical Note: Study of the electron transport parameters used in PENELOPE for the Monte Carlo simulation of Linac targets

Purpose: The Monte Carlo simulation of electron transport in Linac targets using the condensed history technique is known to be problematic owing to a potential dependence of absorbed dose distributions on the electron step length. In the PENELOPE code, the step length is partially determined by the transport parameters C1 and C2. The authors have investigated the effect on the absorbed dose distribution of the values given to these parameters in the target. Methods: A monoenergetic 6.26 MeV electron pencil beam from a point source was simulated impinging normally on a cylindrical tungsten target. Electrons leaving the tungsten were discarded. Radial absorbed dose profiles were obtained at 1.5 cm of depth in a water phantom located at 100 cm for values of C1 and C2 in the target both equal to 0.1, 0.01, or 0.001. A detailed simulation case was also considered and taken as the reference. Additionally, lateral dose profiles were estimated and compared with experimental measurements for a 6 MV photon beam of a Varian Clinac 2100 for the cases of C1 and C2 both set to 0.1 or 0.001 in the target. Results: On the central axis, the dose obtained for the case C1 = C2 = 0.1 shows a deviation of (17.2% ± 1.2%) with respect to the detailed simulation. This difference decreases to (3.7% ± 1.2%) for the case C1 = C2 = 0.01. The case C1 = C2 = 0.001 produces a radial dose profile that is equivalent to that of the detailed simulation within the reached statistical uncertainty of 1%. The effect is also appreciable in the crossline dose profiles estimated for the realistic geometry of the Linac. In another simulation, it was shown that the error made by choosing inappropriate transport parameters can be masked by tuning the energy and focal spot size of the initial beam. Conclusions: The use of large path lengths for the condensed simulation of electrons in a Linac target with PENELOPE conducts to deviations of the dose in the patient or phantom. Based on the results obtained in this work, values of C1 and C2 larger than 0.001 should not be used in Linac targets without further investigation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets—Part I: Mesoscale simulations

Recently a mesoscale model was developed that includes thermomechanically consistent constitutive equations with dissipation due to porous compaction, inelastic distortional deformation rate, and slipping on persistent dry joint surfaces. In this Part I of a two-part paper, we describe enhancements of the mesoscale model which include: normal compliance with hysteresis, non-persistent joints with finite area, joint sets with angular deviation, and the weakening effect of fluid pressure in the joints. Examples study the influences of these effects on average strength. These computationally expensive simulations are used to inspire functional forms in a continuum model, described in Part II, which is needed for large scale simulations of real applications.

42 ENGINEERING↗
DOE OSTI2017

Statistical projections for multi-dimensional visual data exploration

When working with large, multidimensional and multivariate data, science users are frequently interested in understanding variation in data, as opposed to the actual data values. Our work focuses on exploring how a simple statistical metric, the Coefficient of Variation (or Cv), can be used in several different ways to facilitate understanding variation in large data. As a statistical measure, it offers a key advantage over more widely accepted measures like standard deviation, namely to its ability to capture local variation properties. As a multidimensional projection operator, Cv is an effective way of reducing data size while preserving the key variational signal. Visualizations produced from Cv that target conveying variation in data are highly informative, especially compared to those produced with more widely known methods. We demonstrate these ideas within the context of a two-part application case study focusing on understanding long-Term trends in the the changes in precipitation and winds in large-scale climate model ensemble output.

Nguyen, H↗
DOE OSTI2024

Project Dinos I: A joint lensing–dynamics constraint on the deviation from the power law in the mass profile of massive ellipticals

The mass distribution in massive elliptical galaxies encodes their evolutionary history, thus providing an avenue to constrain the baryonic astrophysics in their evolution. The power-law assumption for the radial mass profile in ellipticals has been sufficient to describe several observables to the noise level, including strong lensing and stellar dynamics. In this paper, we quantitatively constrained any deviation, or the lack thereof, from the power-law mass profile in massive ellipticals through joint lensing–dynamics analysis of a large statistical sample with 77 galaxy–galaxy lens systems. We performed an improved and uniform lens modelling of these systems from archival Hubble Space Telescope imaging using the automated lens modelling pipeline dolphin. We combined the lens model posteriors with the stellar dynamics to constrain the deviation from the power law after accounting for the line-of-sight lensing effects, a first for analyses on galaxy–galaxy lenses. We find that the Sloan Lens ACS Survey lens galaxies with a mean redshift of 0.2 are consistent with the power-law profile within 1.1σ (2.8σ) and the Strong Lensing Legacy Survey lens galaxies with a mean redshift of 0.6 are consistent within 0.8σ (2.1σ), for a spatially constant (Osipkov–Merritt) stellar anisotropy profile. We adopted the spatially constant anisotropy profile as our baseline choice based on previous dynamical observables of local ellipticals. However, spatially resolved stellar kinematics of lens galaxies are necessary to differentiate between the two anisotropy models. Future studies will use our lens models to constrain the mass distribution individually in the dark matter and baryonic components.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Assessing the costs of municipal solid waste treatment technologies in developing Asian countries

Highlights: • Municipal waste treatment costs are assessed for the first time in developing Asia. • A database with primary and secondary data from existing plants was prepared. • Capital and operational expenditures patterns were analyzed and discussed. • Incineration identified by large margin as the costliest approach for treating waste. • Context specific circumstances are key to explain deviations from the norm. - Abstract: The management of municipal solid waste (MSW) is one of the main costs incurred by local authorities in developing countries. According to some estimates, these costs can account for up to 50% of city government budgets. It is therefore of importance that policymakers, urban planners and practitioners have an adequate understanding of what these costs consist of, from collection to final waste disposal. This article focuses on a specific stage of the MSW value chain, the treatment of waste, and it aims to identify cost patterns associated with the implementation and operation of waste treatment approaches in developing Asian countries. An analysis of the capital (CAPEX) and operational expenditures (OPEX) of a number of facilities located in countries of the region was conducted based on a database gathering nearly 100 projects and which served as basis for assessing four technology categories: composting, anaerobic digestion (AD), thermal treatment, and the production of refuse-derived fuel (RDF). Among these, it was found that the least costly to invest, as a function of the capacity to process waste, are composting facilities, with an average CAPEX per ton of 21,493 USD{sub 2015}/ton. Conversely, at the upper end featured incineration plants, with an average CAPEX of 81,880 USD{sub 2015}/ton, with this treatment approach ranking by and large as the most capital intensive of the four categories assessed. OPEX figures of the plants, normalized and analyzed in the form of OPEX/ton, were also found to be higher for incineration than for biological treatment methods, although on this component differences amongst the technology groups were less pronounced than those observed for CAPEX. While the results indicated the existence of distinct cost implications for available treatment approaches in the developing Asian context, the analysis also underscored the importance of understanding the local context as a means to properly identify the cost structure of each specific plant. Moreover, even though CAPEX and OPEX figures are important elements to assess the costs of a waste treatment system, these should not be considered on a standalone basis for decision making purposes. In complement to this internal cost dimension, the broader impacts – to the economy, society and the environment – resulting from the adoption of a certain treatment approach should be properly understood and, ideally, measured and expressed in monetary terms.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Search for electroweak production of supersymmetric particles in final states with two τ -leptons in $\sqrt{s}$ = 13 TeV pp collisions with the ATLAS detector

Three searches for the direct production of τ-sleptons or charginos and neutralinos in final states with at least two hadronically decaying τ-leptons are presented. For chargino and neutralino production, decays via intermediate τ-sleptons or W and h bosons are considered. The analysis uses a dataset of pp collisions corresponding to an integrated luminosity of 139 fb -1 , recorded with the ATLAS detector at the Large Hadron Collider at a centre-of-mass energy of 13 TeV. No significant deviation from the expected Standard Model background is observed and supersymmetric particle mass limits at 95% confidence level are obtained in simplified models. For direct production of $\sim\atop{\chi}$ $^+_1$ $\sim\atop{\chi}$ $^-_1$, chargino masses are excluded up to 970 GeV, while $\sim\atop{\chi}$ $^±_1$ and $\sim\atop{\chi}$ $^0_2$ masses up to 1160 GeV (330 GeV) are excluded for $\sim\atop{\chi}$ $^±_1$ $\sim\atop{\chi}$ $^0_2$ / $\sim\atop{\chi}$ $^+_1$ $\sim\atop{\chi}$ $^-_1$ production with subsequent decays via τ-sleptons (W and h bosons). Masses of τ-sleptons up to 500 GeV are excluded for mass degenerate $\sim\atop{\tau}$ L,R scenarios and up to 425 GeV for $\sim\atop{\tau}$ L-only scenarios. Sensitivity to $\sim\atop{\tau}$ R -only scenarios from the ATLAS experiment is presented here for the first time, with $\sim\atop{\tau}$ R masses excluded up to 350 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Evaluation and development of flow condensation correlations using the data from low GWP refrigerants in an axial micro-fin aluminum tube

To mitigate global warming, the world is transitioning to refrigerants with low global warming potential (GWP). Supporting this shift requires a model that can accurately predict the heat transfer and pressure drop of new refrigerants, crucial for designing efficient heat exchangers. Existing models, however, are largely based on currently deployed refrigerants and primarily developed for unexpanded micro-fin tubes with spiral angles of 6° to 30°. Their applicability to new refrigerants, especially in expanded micro-fin tubes, is uncertain. This study assesses the performance of four well-known condensation models for six emerging refrigerants—R-32, R-454B, R-454C, R-455A, R-1234yf, and R-1234ze(E)—against experimental data. Initially, the Han and Lee (2005) model shows the best prediction accuracy with a mean absolute deviation (MAD) of 22.1 %. To enhance the accuracy of heat transfer models for new refrigerants and geometries with large temperature glides, two approaches are proposed. Here, the first approach applies a simple correction factor, reducing the MAD of the Cavallini et al. (2009) model from 68.2 % to 15.4 %. The second approach uses the variable metric method for minimization, fitting new constants to the data. This optimization results in the Kedzierski and Goncalves (1997) model achieving the highest accuracy, with a MAD of 13.1 %. For pressure drop models, the Cavallini et al. (1997) model is the most accurate with a MAD of 6.4 %, followed by the Haraguchi et al. (1993) model with a MAD of 9.4 %. Due to its simplicity, the Haraguchi et al. (1993) model is a practical option for predicting frictional pressure drop.

Axial micro-fin tubes↗
DOE OSTI2016

Prospective Study Delivering Simultaneous Integrated High-dose Tumor Boost (≤70 Gy) With Image Guided Adaptive Radiation Therapy for Radical Treatment of Localized Muscle-Invasive Bladder Cancer

Purpose: Image guided adaptive radiation therapy offers individualized solutions to improve target coverage and reduce normal tissue irradiation, allowing the opportunity to increase the radiation tumor dose and spare normal bladder tissue. Methods and Materials: A library of 3 intensity modulated radiation therapy plans were created (small, medium, and large) from planning computed tomography (CT) scans performed at 30 and 60 minutes; treating the whole bladder to 52 Gy and the tumor to 70 Gy in 32 fractions. A “plan of the day” approach was used for treatment delivery. A post-treatment cone beam CT (CBCT) scan was acquired weekly to assess intrafraction filling and coverage. Results: A total of 18 patients completed treatment to 70 Gy. The plan and treatment for 1 patient was to 68 Gy. Also, 1 patient's plan was to 70 Gy but the patient was treated to a total dose of 65.6 Gy because dose-limiting toxicity occurred before dose escalation. A total of 734 CBCT scans were evaluated. Small, medium, and large plans were used in 36%, 48%, and 16% of cases, respectively. The mean ± standard deviation rate of intrafraction filling at the start of treatment (ie, week 1) was 4.0 ± 4.8 mL/min (range 0.1-19.4) and at end of radiation therapy (ie, week 5 or 6) was 1.1 ± 1.6 mL/min (range 0.01-7.5; P=.002). The mean D{sub 98} (dose received by 98% volume) of the tumor boost and bladder as assessed on the post-treatment CBCT scan was 97.07% ± 2.10% (range 89.0%-104%) and 99.97% ± 2.62% (range 96.4%-112.0%). At a median follow-up period of 19 months (range 4-33), no muscle-invasive recurrences had developed. Two patients experienced late toxicity (both grade 3 cystitis) at 5.3 months (now resolved) and 18 months after radiation therapy. Conclusions: Image guided adaptive radiation therapy using intensity modulated radiation therapy to deliver a simultaneous integrated tumor boost to 70 Gy is feasible, with acceptable toxicity, and will be evaluated in a randomized trial.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Test of lepton universality in beauty-quark decays

The standard model of particle physics currently provides our best description of fundamental particles and their interactions. The theory predicts that the different charged leptons, the electron, muon and tau, have identical electroweak interaction strengths. Previous measurements have shown that a wide range of particle decays are consistent with this principle of lepton universality. This article presents evidence for the breaking of lepton universality in beauty-quark decays, with a significance of 3.1 standard deviations, based on proton–proton collision data collected with the LHCb detector at CERN’s Large Hadron Collider. The measurements are of processes in which a beauty meson transforms into a strange meson with the emission of either an electron and a positron, or a muon and an antimuon. If confirmed by future measurements, this violation of lepton universality would imply physics beyond the standard model, such as a new fundamental interaction between quarks and leptons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Interference of electrical double layers: Confinement effects on structure, dynamics, and screening of ionic liquids

Ionic liquids are widely used as electrolytes in electronic devices in which they are subject to nanoconfinement within nanopores or nanofilms. Because the intrinsic width of an electrical double layer is on the order of several nanometers, nanoconfinement is expected to fundamentally alter the double layer properties. Furthermore, in confined systems, a large portion of the ions are interfacial, e.g., at the electrode interface, leading to significant deviations of electrostatic screening and ion dynamics as compared to bulk properties. In this article, we systematically investigate the interference between electrical double layers for nanoconfined ionic liquids and the resulting influence on the structure, dynamics, and screening behavior. We perform molecular dynamics simulations for the ionic liquids [BMIm + ][BF$^-_4$] and [BMIm + ][PF$^-_6$] confined between two flat electrodes at systematic separation distances between 1.5 nm and 4.5 nm for both conducting and insulating boundary conditions. We find that while ion dynamics is expectedly slower than in the bulk (by ~2 orders of magnitude), there is an unexpected non-linear trend with the confinement length that leads to a local maximum in dynamic rates at ~3.5–4.5 nm confinement. We show that this nonlinear trend is due to the ion correlation that arises from the interference between opposite double layers. We further evaluate confinement effects on the ion structure and capacitance and investigate the influence of electronic polarization of the ionic liquid on the resulting properties. This systematic evaluation of the connection between electrostatic screening and structure and dynamics of ionic liquids in confined systems is important for the fundamental understanding of electrochemical supercapacitors.

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

Statistics of anisotropies in inflation with spectator vector fields

We study the statistics of the primordial power spectrum in models where massless gauge vectors are coupled to the inflaton, paying special attention to observational implications of having fundamental or effective horizons embedded in a bath of infrared fluctuations. As quantum infrared modes cross the horizon, they classicalize and build a background vector field. We find that the vector experiences a statistical precession phenomenon. Implications for primordial correlators and the interpretation thereof are considered. Firstly, we show how in general two, not only one, additional observables, a quadrupole amplitude and an intrinsic shape parameter, are necessary to fully describe the correction to the curvature power spectrum, and develop a unique parametrization for them. Secondly, we show that the observed anisotropic amplitude and the associated preferred direction depend on the volume of the patch being probed. We calculate non-zero priors for the expected deviations between detections based on microwave background data (which probes the entire Hubble patch) and large scale structure (which only probes a fraction of it)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Cosmology based on f(R) gravity with O(1) eV sterile neutrino

We address the cosmological role of an additional O(1) eV sterile neutrino in modified gravity models. We confront the present cosmological data with predictions of the FLRW cosmological model based on a variant of f(R) modified gravity proposed by one of the authors previously. This viable cosmological model which deviation from general relativity with a cosmological constant Λ decreases as R{sup −2n} for large, but not too large values of the Ricci scalar R (while no Λ is introduced by hand at small R) provides an alternative explanation of present dark energy and the accelerated expansion of the Universe (the case n=2 is considered in the paper). Various up-to-date cosmological data sets exploited include measurements of the cosmic microwave background (CMB) anisotropy, the CMB lensing potential, the baryon acoustic oscillations (BAO), the cluster mass function and the Hubble constant. We find that the CMB+BAO constraints strongly restrict the sum of neutrino masses from above. This excludes values of the model parameter λ∼ 1 for which distinctive cosmological features of the model are mostly pronounced as compared to the ΛCDM model, since then free streaming damping of perturbations due to neutrino rest masses is not sufficient to compensate their extra growth occurring in f(R) modified gravity. Thus, in the gravity sector we obtain λ>8.2 (2σ) with the account of systematic uncertainties in galaxy cluster mass function measurements and λ>9.4 (2σ) without them. At the same time in the latter case we find for the sterile neutrino mass 0.47 eV < m{sub ν, sterile} < 1 eV (2σ) assuming that the sterile neutrinos are thermalized and the active neutrinos are massless, not significantly larger than in the standard ΛCDM with the same data set: 0.45 eV < m{sub ν, sterile} < 0.92 eV (2σ). However, a possible discovery of a sterile neutrino with the mass m{sub ν, sterile} ≈ 1.5 eV motivated by various anomalies in neutrino oscillation experiments would favor cosmology based on f(R) gravity rather than the ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

OBJECTS APPEAR SMALLER AS THEY RECEDE: HOW PROPER MOTIONS CAN DIRECTLY REVEAL THE COSMIC EXPANSION, PROVIDE GEOMETRIC DISTANCES, AND MEASURE THE HUBBLE CONSTANT

Objects and structures gravitationally decoupled from the Hubble expansion will appear to shrink in angular size as the universe expands. Observations of extragalactic proper motions can thus directly reveal the cosmic expansion. Relatively static structures such as galaxies or galaxy clusters can potentially be used to measure the Hubble constant, and test masses in large scale structures can measure the overdensity. Since recession velocities and angular separations can be precisely measured, apparent proper motions can also provide geometric distance measurements to static structures. The apparent fractional angular compression of static objects is 15 μas yr{sup –1} in the local universe; this motion is modulated by the overdensity in dynamic expansion-decoupled structures. We use the Titov et al. quasar proper motion catalog to examine the pairwise proper motion of a sparse network of test masses. Small-separation pairs (<200 Mpc comoving) are too few to measure the expected effect, yielding an inconclusive 8.3 ± 14.9 μas yr{sup –1}. Large-separation pairs (200-1500 Mpc) show no net convergence or divergence for z < 1, –2.7 ± 3.7 μas yr{sup –1}, consistent with pure Hubble expansion and significantly inconsistent with static structures, as expected. For all pairs a 'null test' gives –0.36 ± 0.62 μas yr{sup –1}, consistent with Hubble expansion and excludes a static locus at ∼5-10σ significance for z ≅ 0.5-2.0. The observed large-separation pairs provide a reference frame for small-separation pairs that will significantly deviate from the Hubble flow. The current limitation is the number of small-separation objects with precise astrometry, but Gaia will address this and will likely detect the cosmic recession.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Information-theoretic astrophysical uncertainties in the effective theory of dark matter direct detection

The impact of astrophysical uncertainties in direct detection searches can vary significantly across particle dark matter models and detector targets, due to the different velocity and momentum dependencies of the scattering cross section. We address these uncertainties for all operators of the nonrelativistic effective field theory of dark-matter/nucleon interactions, making use of the Kullback-Leibler (KL) information divergence to measure the deviation of the true dark matter velocity distribution from the Maxwell-Boltzmann form. This approach quantifies how astrophysical uncertainties affect each operator in the effective theory, without assuming any specific functional form for the velocity distribution. While for some operators the uncertainties are smaller than 1 order of magnitude for entropically motivated deviations from the Maxwell-Boltzmann form, for other operators, these uncertainties can be as large as three orders of magnitude near threshold. Furthermore, we identify the dependence of the scattering rate for various operators of the effective theory with different velocity-weighted moments of the velocity distribution, functionally analogous to the mean, variance, or skewness. This provides new analytic insight into which features of the velocity distribution are most relevant to detect a given particle dark matter model. Our technique is general and could be applied to a broader class of physics problems where a physical observable depends on the statistical moments of an uncertain theoretical distribution.

Herrera, Gonzalo [MIT, MKI; Harvard U.; Virginia T↗
DOE OSTI2019

Downside risk in reservoir management

Downside risk, which refers to deviations below a threshold, is often important in water management decisions, especially in areas with large and skewed variations in precipitation patterns. In this paper, we present a model for a reservoir manager who is downside risk averse and who performs a dynamic allocation of irrigation water, taking into account the negative effects of droughts on farm profits and different environmental constraints. We analyse the water stock, flows and agricultural profits for alternative environmental restrictions and thresholds for irrigation levels and find that stricter environmental constraints increase total water supply and carryover stock, while higher penalty thresholds lead to their overall decrease. Furthermore, increasing penalty thresholds leads to a higher emphasis on avoiding shortages, at the expense of lower average profits.

Roseta‐Palma, Catarina↗