SEARCH · Engineering Papers

Results for “flow conservation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Divergence study of a high-aspect ratio, forward-swept wing

A study has been conducted in the NASA-Langley Transonic Dynamics wind tunnel to determine the divergence characteristics of a high aspect ratio forward swept wing of rectangular planform, for sweep angles of zero, -15, -30, -45, and -60 deg. In addition to a rectangular wingtip, a tip whose geometry lay parallel with the freestream flow was tested at a sweep angle of -45 deg. The primary objective was to obtain data that could be used to verify the divergence prediction capabilities of an aeroelastic analysis code employing kernel function aerodynamics. The analytical predictions of convergence were found to be conservative for all forward sweep angles, and the effect of the two tip shapes on divergence dynamic pressure were accurately predicted.

Cole, S. R.↗

NASA NTRS1992

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗

NASA NTRS1992

Development of an improved model for runback water on aircraft surfaces

A computer simulation for 'running wet' and evaporative aircraft anti-icing systems is developed. The model is based on the analysis of the liquid water film which forms in the regions of direct impingement and, then, breaks up near the impingement limits into rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted using a stability analysis theory and the laws of mass energy conservation. The solid structure is modeled as a multiple layer wall. The anti-icing system modeled is of the thermal type utilizing hot air and/or electrical heating elements embedded within the wall layers. Experimental observations revealing some of the basic physics of the water flow on the surface are presented. Detailed qualitative documentation of the tests are given. Several numerical examples are considered, and the effect of some of the involved parameters on the system performance are investigated.

Al-Khalil, Kamel M.↗

NASA NTRS2021

Wilson Corners SWMU 001 2019-2020 Annual Long-Term Monitoring Report Kennedy Space Center, Florida

This report presents a summary of the long-term monitoring (LTM) activities that occurred in 2019 and 2020 at Wilson Corners, Solid Waste Management Unit 001, at the John F. Kennedy Space Center (KSC), Florida. The site is monitored under KSC’s Resource Conservation and Recovery Act Corrective Action Program. Adaptive site management is being utilized through ongoing assessment, design, and interim measures. This approach also meets the requirements of Chapter 62-780, Florida Administrative Code. The goal of LTM at this site is threefold: to determine groundwater flow characteristics, monitor the downgradient concentration trends, and monitor select locations internal to the groundwater plume. Every 5 years, upgradient and side-gradient monitoring wells are sampled to verify delineation. The last time this was performed was in 2015. Sampling of these wells in 2020 was replaced with the direct push technology investigations in October 2020 and April 2021. The activities presented in this report include three field events: (1) December 2019 -groundwater gauging of 38 monitoring wells and sampling of 36 monitoring wells; (2) May 2020 - groundwater sampling of 7 monitoring wells; and (3) December 2020 - groundwater gauging of 49 monitoring wells and sampling of 43 monitoring wells. During the December 2019 and May 2020 events, monitoring wells were sampled using passive diffusion bags and were analyzed for a “full list” of volatile organic compounds (VOCs), including Freon 113. During the December 2020 event, the low flow sampling method was used, and samples were analyzed for a “select list” of VOCs, including Freon 113. The following conclusions can be made based on the 2019 and 2020 LTM results: •In December 2019 and December 2020, groundwater flow for the site was generally to the west with the individual zones only varying by occasional northerly and southerly components. This is generally consistent with historical observations at the site. •The GCTL plume continues to extend both horizontally and vertically beyond our current monitoring well network. •The vertical extent of VOCs was historically generally delineated by monitoring wells screened greater than 48 feet bls. The results from the two vertical extent monitoring wells WILC-MW0078 and WILC-MW0130 that were sampled during the 2019 and 2020 LTM indicate that groundwater VOC concentrations in both wells were greater than the GCTL. •Freon 113 was not detected above GCTLs during the 2019 and 2020 LTM events. Based on groundwater sampling activities performed during this reporting period, including recent DPT groundwater sampling conducted for the implementation of an air sparge system, the following recommendations are provided: •Install 10 new monitoring wells, which will be sampled along with 49 existing monitoring wells, to assist with delineation of the low concentration plume. The wells are listed inTable 4-1, and the proposed locations are presented on Figures4 -1 through 4-5. (Figures 4-1, 4-2, 4-3, and 4-5 also show updated plumes and proposed performance monitoring well locations. This information will be discussed and presented under a separate cover.) •Modify the annual LTM sampling plan as presented inTable 4-1. •Sampling for future events at the site will be conducted using low flow pumping methods. The next annual LTM sampling event is currently planned for December 2021. Please note the 10 new monitoring wells will not be installed before this event but will be installed as part of the upcoming IM construction and will be sampled during baseline sampling, along with the 49 existing monitoring wells. Following this sampling event, it is recommended that the Annual LTM be combined with the performance monitoring under the IM implementation, and LTM be temporarily discontinued. Objectives of the 2021 LTM are to: (i) evaluate groundwater gradient and flow direction by collecting depth to water measurements from LTM wells; (ii) continue monitoring the peripheral VOC trends in the northern, southern, and western portions of the site by monitoring existing wells and installing and monitoring new wells as recommended; and (iii) monitor select internal plume wells.

Jennifer Lynn Joyal↗

NASA NTRS2012

A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing

Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.

Bellan, Josette↗

NASA NTRS2002

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗

DOE OSTI2014

On the interaction between turbulence and a planar rarefaction

The modeling of turbulence, whether it be numerical or analytical, is a difficult challenge. Turbulence is amenable to analysis with linear theory if it is subject to rapid distortions, i.e., motions occurring on a timescale that is short compared to the timescale for nonlinear interactions. Such an approach (referred to as rapid distortion theory) could prove useful for understanding aspects of astrophysical turbulence, which is often subject to rapid distortions, such as supernova explosions or the free-fall associated with gravitational instability. As a proof of principle, a particularly simple problem is considered here: the evolution of vorticity due to a planar rarefaction in an ideal gas. Analytical solutions are obtained for incompressive modes having a wave vector perpendicular to the distortion; as in the case of gradient-driven instabilities, these are the modes that couple most strongly to the mean flow. Vorticity can either grow or decay in the wake of a rarefaction front, and there are two competing effects that determine which outcome occurs: entropy fluctuations couple to the mean pressure gradient to produce vorticity via baroclinic effects, whereas vorticity is damped due to the conservation of angular momentum as the fluid expands. Whether vorticity grows or decays depends upon the ratio of entropic to vortical fluctuations at the location of the front; growth occurs if this ratio is of order unity or larger. In the limit of purely entropic fluctuations in the ambient fluid, a strong rarefaction generates vorticity with a turbulent Mach number on the order of the rms of the ambient entropy fluctuations. The analytical results are shown to compare well with results from two- and three-dimensional numerical simulations. Analytical solutions are also derived in the linear regime of Reynolds-averaged turbulence models. This highlights an inconsistency in standard turbulence models that prevents them from accurately capturing the physics of rarefaction-turbulence interaction. In addition to providing physical insight, the solutions derived here can be used to verify algorithms of both the Reynolds-averaged and direct numerical simulation variety. Finally, dimensional analysis of the equations indicates that rapid distortion of turbulence can give rise to two distinct regimes in the turbulent spectrum: a distortion range at large scales where linear distortion effects dominate, and an inertial range at small scales where nonlinear effects dominate.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Insulated Concrete Form Walls Integrated With Mechanical Systems in a Cold Climate Test House

Transitioning from standard light frame to a thermal mass wall system in a high performance home will require a higher level of design integration with the mechanical systems. The much higher mass in the ICF wall influences heat transfer through the wall and affects how the heating and cooling system responds to changing outdoor conditions. This is even more important for efficient, low-load homes with efficient heat pump systems in colder climates where the heating and cooling peak loads are significantly different from standard construction. This report analyzes a range of design features and component performance estimates in an effort to select practical, cost-effective solutions for high performance homes in a cold climate. Of primary interest is the influence of the ICF walls on developing an effective air sealing strategy and selecting an appropriate heating and cooling equipment type and capacity. The domestic water heating system is analyzed for costs and savings to investigate options for higher efficiency electric water heating. A method to ensure mechanical ventilation air flows is examined. The final solution package includes high-R mass walls, very low infiltration rates, multi-stage heat pump heating, solar thermal domestic hot water system, and energy recovery ventilation. This solution package can be used for homes to exceed 2012 International Energy Conservation Code requirements throughout all climate zones and achieves the DOE Challenge Home certification.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1986

Transonic flow solutions using a composite velocity procedure for potential, Euler and RNS equations

Solutions for transonic viscous and inviscid flows using a composite velocity procedure are presented. The velocity components of the compressible flow equations are written in terms of a multiplicative composite consisting of a viscous or rotational velocity and an inviscid, irrotational, potential-like function. This provides for an efficient solution procedure that is locally representative of both asymptotic inviscid and boundary layer theories. A modified conservative form of the axial momentum equation that is required to obtain rotational solutions in the inviscid region is presented and a combined conservation/nonconservation form is applied for evaluation of the reduced Navier-Stokes (RNS), Euler and potential equations. A variety of results is presented and the effects of the approximations on entropy production, shock capturing, and viscous interaction are discussed.

Gordnier, R. E.↗

NASA NTRS1986

Higher-order flux difference splitting schemes for the Euler equations using upstream interpolations

A class of explicit two time-level, 2p + 1 space-point, (2p 1)-th order, upwind-biased flux difference splitting schemes are proposed for the numerical advection based on Lagrange's interpolation, and the method is an accord with the physical domain of dependence. A normalized Jacobian coefficient matrix is introduced to convert the schemes to hyperbolic systems of conservation laws, and approaches to make the higher-order schemes total variation stable are discussed. Accuracy and stability of the present schemes are examined, and implicit total variation diminishing schemes are developed for steady-state calculations.Application to gasdynamic problems for both steady and unsteady flows covering a wide range of Mach numbers is considered, and results for a blast wave passing a cylinder, and head-on collision of two blast waves over a circular arc, are presented. The flow patterns were found to be symmetric, and good resolution of flow structures was obtained.

Yang, J. Y.↗

NASA NTRS1989

Transonic flow solutions using a composite velocity procedure for potential, Euler and RNS equations

Gordnier, R. E.↗

NASA NTRS1992

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Goldberg, Louis F.↗

NASA NTRS2004

Part 1 of a Computational Study of a Drop-Laden Mixing Layer

This first of three reports on a computational study of a drop-laden temporal mixing layer presents the results of direct numerical simulations (DNS) of well-resolved flow fields and the derivation of the large-eddy simulation (LES) equations that would govern the larger scales of a turbulent flow field. The mixing layer consisted of two counterflowing gas streams, one of which was initially laden with evaporating liquid drops. The gas phase was composed of two perfect gas species, the carrier gas and the vapor emanating from the drops, and was computed in an Eulerian reference frame, whereas each drop was tracked individually in a Lagrangian manner. The flow perturbations that were initially imposed on the layer caused mixing and eventual transition to turbulence. The DNS database obtained included transitional states for layers with various liquid mass loadings. For the DNS, the gas-phase equations were the compressible Navier-Stokes equations for conservation of momentum and additional conservation equations for total energy and species mass. These equations included source terms representing the effect of the drops on the mass, momentum, and energy of the gas phase. From the DNS equations, the expression for the irreversible entropy production (dissipation) was derived and used to determine the dissipation due to the source terms. The LES equations were derived by spatially filtering the DNS set and the magnitudes of the terms were computed at transitional states, leading to a hierarchy of terms to guide simplification of the LES equations. It was concluded that effort should be devoted to the accurate modeling of both the subgridscale fluxes and the filtered source terms, which were the dominant unclosed terms appearing in the LES equations.

Okong'o, Nora A.↗

NASA NTRS2012

High Order Finite Difference Methods with Subcell Resolution for 2D Detonation Waves

In simulating hyperbolic conservation laws in conjunction with an inhomogeneous stiff source term, if the solution is discontinuous, spurious numerical results may be produced due to different time scales of the transport part and the source term. This numerical issue often arises in combustion and high speed chemical reacting flows.

Wang, W.↗

DOE OSTI2014

Building America Case Study: Lancaster County Career and Technology Center Green Home 3, Mt Joy, Pennsylvania

Transitioning from standard light frame to a thermal mass wall system in a high performance home will require a higher level of design integration with the mechanical systems. The much higher mass in the ICF wall influences heat transfer through the wall and affects how the heating and cooling system responds to changing outdoor conditions. This is even more important for efficient, low-load homes with efficient heat pump systems in colder climates where the heating and cooling peak loads are significantly different from standard construction.This report analyzes a range of design features and component performance estimates in an effort to select practical, cost-effective solutions for high performance homes in a cold climate. Of primary interest is the influence of the ICF walls on developing an effective air sealing strategy and selecting an appropriate heating and cooling equipment type and capacity. The domestic water heating system is analyzed for costs and savings to investigate options for higher efficiency electric water heating. A method to ensure mechanical ventilation air flows is examined. The final solution package includes high-R mass walls, very low infiltration rates, multi-stage heat pump heating, solar thermal domestic hot water system, and energy recovery ventilation. This solution package can be used for homes to exceed 2012 International Energy Conservation Code requirements throughout all climate zones and achieves the DOE Challenge Home certification.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS2002

Aeroacoustic Flow Phenomena Accurately Captured by New Computational Fluid Dynamics Method

One of the challenges in the computational fluid dynamics area is the accurate calculation of aeroacoustic phenomena, especially in the presence of shock waves. One such phenomenon is "transonic resonance," where an unsteady shock wave at the throat of a convergent-divergent nozzle results in the emission of acoustic tones. The space-time Conservation-Element and Solution-Element (CE/SE) method developed at the NASA Glenn Research Center can faithfully capture the shock waves, their unsteady motion, and the generated acoustic tones. The CE/SE method is a revolutionary new approach to the numerical modeling of physical phenomena where features with steep gradients (e.g., shock waves, phase transition, etc.) must coexist with those having weaker variations. The CE/SE method does not require the complex interpolation procedures (that allow for the possibility of a shock between grid cells) used by many other methods to transfer information between grid cells. These interpolation procedures can add too much numerical dissipation to the solution process. Thus, while shocks are resolved, weaker waves, such as acoustic waves, are washed out.

Blech, Richard A.↗

DOE OSTI2020

Performance evaluation of an air conditioning system based on quasi isothermal dehumidifcation

This work presents a first study of novel air conditioning system based on quasi isothermal dehumidification combined with dew-point evaporative cooling through the Maisotsenko cycle (M-Cycle). The dehumidification process is realized in desiccant coated heat exchanger with heat rejection system. The key idea behind the solution is to maximize the cooling effectiveness by creating favourable conditions for dehumidification and evaporative cooling process. This is achieved by rejecting the heat of sorption in the desiccant system and providing cold and dry air for M-Cycle air cooler, which allows the system to be powered by using low-grade heat energy at 50 °C. The analysis is performed on the basis of experimentally validated numerical models of each component combined in a unified model of the system. The goal of the study is to determine the performance of the solution under different operational conditions. In addition two types of M-Cycle air coolers (i.e. cross-flow and regenerative unit) are compared to establish which unit is more suitable as the source of cooling energy for the system. It was established that proposed system is able to provide comfortable room parameters (supply air temperature equal 17 °C and humidity ratio equal 0.010 kg/kg-da) for outdoor air temperatures ranging from 28 to 36 °C and outdoor air humidity ranging from 0.010 to 0.018 kg/kg-da. This means that the proposed solution is able to provide comfortable conditions in temperate and humid climates. Results also indicate that the system consisting of the regenerative M-Cycle HMX allows better cooling load management, offers higher thermal Coefficient of Performance (COP th ) and uses less water compared to the system with the cross-flow HMX. However, the system with cross-flow heat and mass exchanger (HMX) has a relatively greater electrical Coefficient of Performance (COP el ) in most cases.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2016

High-frequency monitoring of water fluxes and nutrient loads to assess the effects of controlled drainage on water storage and nutrient transport

High nitrogen (N) and phosphorus (P) fluxes from upstream agriculture threaten aquatic ecosystems in surface waters and estuaries, especially in areas characterized by high agricultural N and P inputs and densely drained catchments like the Netherlands. Controlled drainage has been recognized as an effective option to optimize soil moisture conditions for agriculture and to reduce unnecessary losses of fresh water and nutrients. This is achieved by introducing control structures with adjustable overflow levels into subsurface tube drain systems. A small-scale (1 ha) field experiment was designed to investigate the hydrological and chemical changes after introducing controlled drainage. Precipitation rates and the response of water tables and drain fluxes were measured in the periods before the introduction of controlled drainage (2007–2008) and after (2009–2011). For the N and P concentration measurements, auto-analyzers for continuous records were combined with passive samplers for time-averaged concentrations at individual drain outlets. The experimental setup enabled the quantification of changes in the water and solute balance after introducing controlled drainage. The results showed that introducing controlled drainage reduced the drain discharge and increased the groundwater storage in the field. To achieve this, the overflow levels have to be elevated in early spring, before the drain discharge stops due to dryer conditions and falling groundwater levels. The groundwater storage in the field would have been larger if the water levels in the adjacent ditch were controlled as well by an adjustable weir. The N concentrations and loads increased, which was largely related to elevated concentrations in one of the three monitored tube drains. The P loads via the tube drains reduced due to the reduction in discharge after introducing controlled drainage. Furthermore, this may be counteracted by the higher groundwater levels and the larger contribution of N- and P-rich shallow groundwater and overland flow to the surface water.

54 ENVIRONMENTAL SCIENCES↗