Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow path”

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 955 records · Page 53

Linear magnetic motor/generator

A linear magnetic motor/generator is disclosed which uses magnetic flux to provide mechanical motion or electrical energy. The linear magnetic motor/generator includes an axially movable actuator mechanism. A permament magnet mechanism defines a first magnetic flux path which passes through a first end portion of the actuator mechanism. Another permament magnet mechanism defines a second magnetic flux path which passes through a second end portion of the actuator mechanism. A drive coil defines a third magnetic flux path passing through a third central portion of the actuator mechanism. A drive coil selectively adds magnetic flux to and subtracts magnetic flux from magnetic flux flowing in the first and second magnetic flux path.

Studer, P. A.↗

Characteristic path analysis of confinement influence on steady two-dimensional detonation propagation

Steady detonation in multi-dimensional flow is controlled by the chemical energy release that occurs in a subsonic elliptic flow region known as the detonation driving zone (DDZ). It is the region encompassing the detonation shock and sonic flow locus (in the frame of the detonation shock). A detonation that is strongly confined by material surrounding the explosive has the shock and sonic locus separated at the material interface. Information about the material boundary is traditionally believed to influence the DDZ structure via the subsonic flow on the boundary ahead of the sonic locus. A detonation that is weakly confined has the detonation shock and sonic locus intersecting at the material boundary. The sonic nature of the flow at the intersection point on the boundary is believed to isolate the DDZ structure from the material properties of the confinement. We examine the paths of characteristics propagating information about the confinement through the supersonic hyperbolic flow region that exists beyond the sonic locus, and determine whether these paths may impinge on the sonic locus and consequently influence the DDZ structure. Our configuration consists of a solid wall boundary deflected through a specified angle on detonation shock arrival, so that the streamline turning angle of the wall at the explosive edge is unambiguously defined. By varying the wall deflection angle from small through large values, we can systematically capture the evolution of the DDZ structure and the characteristic flow regions that influence its structure for strongly to weakly confined detonations. In all strong and weak confinement cases examined, we find that a subset of characteristics from the supersonic flow regions always impinge on the sonic locus. Limiting characteristics are identified that define the boundary between characteristics that impinge on the sonic surface and those that propagate information downstream of the sonic surface. In combination with an oblique-shock polar analysis, we show that the effects on the DDZ of characteristic impingement can be significant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer

Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.

42 ENGINEERING↗

Gluing together modular flows with free fermions

We revisit the calculation of multi-interval modular Hamiltonians for free fermions using a Euclidean path integral approach. We show how the multi-interval modular flow is obtained by gluing together the single interval modular flows. Using this relation, we obtain an exact expression for the multi-interval modular Hamiltonian and entanglement entropy in agreement with existing results. An essential ingredient in our derivation is the introduction of the modular action. This determines the non-local field theory describing the free fermion reduced density matrix, and makes manifest it’s non-local conformal symmetry and U(1) Kacs-Moody symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Continuation in a parameter - Experience with viscous and free surface flows

The results of modifications in continuation methods applied to obtain solutions to the Navier-Stokes systems of equations for incompressible, two-dimensional, steady flows are reported. It is shown that parameter continuation permits prediction of accurate, initial estimates for iterative processing of nonlinear finite difference and finite element equations of motions. The new parameter steps are derived from values of the preceding parameter steps. The accuracy of the estimates is ensured with appropriate choices of the step size. The continuation predictor/iterative corrector is demonstrated to trace the branches of parameter space along which steady flow states are found, and techniques are available for tracing multiply branching paths. The techniques are applied to solving the Navier-Stokes equations for flow through a rotating square channel, the formation of a falling liquid curtain, and gyrostatic equilibria of rotating cylindrical drops.

Kheshgi, H. S.↗

SDN-NGenIA, a Software Defined Next Generation Integrated Architecture for HEP and Data Intensive Science

This project aims to enable the LHC and other leading programs in high energy physics and other global science domains funded by the DOE to operate with a new level of efficiency and control, through the development of a Next Generation Integrated Architecture (NGenIA) based on intelligent software defined network (SDN)-driven network systems coupled to high throughput applications. While the initial focus is on the challenging LHC use case, the systems and products being developed are general, and apply to many fields of data intensive science ranging from astrophysical sky surveys to bioinformatics and earth observation, as well as other organizations facing the challenges of extracting knowledge from distributed multi-petabyte data stores. A central concept in this development program is a new paradigm ‘consistent network operations’ among widely distributed computing and storage facilities, where stable high throughput flows, at set rates, cross load-balanced network paths, up to flexible high water marks that are adjusted in real time to accommodate other network traffic. The large smooth flows are launched and managed by SDN services that act in concert with the experiments’ site-resident data distribution and management systems, to meet the expanding needs of the science programs. The technical goals include the construction of autonomous, intelligent site-resident services that interact dynamically with network-resident services, and with the science programs’ principal data distribution and management tools, to request or command network resources in support of high throughput petascale to exascale workflow Specific work items include: 1. Developing compact Data Transfer Nodes (DTNs) with auto-tuning functions that support data transfer rates in the 100 Gbps to the 1 Tbps range when used with high throughput data transfer applications 2. Deep site orchestration among virtualized clusters, storage subsystems and subnets to successfully coschedule CPU, storage and network resources 3. Science-program designed site architectures, operational modes, and policy and resource usage priorities, adjudicated across multiple network domains and virtual organizations, using the orchestration functions and methods 4. Seamlessly extending end-to-end operation across both extra-site and intra-site boundaries through the use of Open vSwitch (OVS), FireQoS + next generation Science DMZs 5. Novel methods of system integration that enable granular control of extreme scale long distance transfers through flow matching of scattered source-destination address pairs to multi-domain dynamic circuits 6. Funneling massive sets of streams to DTNs at the site edge hosting petascale buffer pools configured for flows of 100 Gbps and up, exploiting state of the art data transfers 7. Adaptive scheduling based on pervasive end-to-end monitoring, including DTN or compute-node resident agents providing comprehensive end-system profiling 8. Developing unsupervised and supervised machine learning and modeling methods to optimize the workflow involving terabyte to multi-petabyte datasets.

47 OTHER INSTRUMENTATION↗

SCC Aerosol Transport Model Summary Report

The purpose of the present work is to introduce a generic and reliable numerical model for prediction of aerosol transport, deposition, and plugging in leak paths similar to stress corrosion cracks while accounting for potential plugging formation. The model is dynamic (changing leak path geometry due to plugging) and it relies on the numerical solution of the aerosol transport equation in one dimension using finite differences. An extensive validation exercise (particle diameters: 0.01-10 μm and pressure difference up to 12 kPa) was conducted based on comparisons with experimental and theoretical data. The developed model is fairly general since it is based on a generic mechanistic description of the aerosol flow and particle deposition in the leak path, without requiring experimental fittings other than using the leak rate parameters commonly available in practical applications. Despite the approximate agreement of the model with the experimental data, the obtained predictions can be qualified as realistic. This is a major achievement of the model, which has been based on generic aerosol mechanics with no reference or adjustment to experimental data.

42 ENGINEERING↗

Cooled spool piston compressor

A hydraulically powered gas compressor receives low pressure gas and outputs a high pressure gas. The housing of the compressor defines a cylinder with a center chamber having a cross-sectional area less than the cross-sectional area of a left end chamber and a right end chamber, and a spool-type piston assembly is movable within the cylinder and includes a left end closure, a right end closure, and a center body that are in sealing engagement with the respective cylinder walls as the piston reciprocates. First and second annual compression chambers are provided between the piston enclosures and center housing portion of the compressor, thereby minimizing the spacing between the core gas and a cooled surface of the compressor. Restricted flow passageways are provided in the piston closure members and a path is provided in the central body of the piston assembly, such that hydraulic fluid flows through the piston assembly to cool the piston assembly during its operation. The compressor of the present invention may be easily adapted for a particular application, and is capable of generating high gas pressures while maintaining both the compressed gas and the compressor components within acceptable temperature limits.

Morris, Brian G.↗

Forward rotor vortex effects on counter rotating propeller noise

Three configurations of a model counter rotating propeller manipulate the blade tip flow by: placing the CRP at angle of attack, installing shrouds, and turning the upstream blades to provide forward sweep. Flow visualization and flow measurements with thermal anemometry show no evidence of a tip vortex; however, a leading edge vortex was detected on aft swept blades. The modifications served to alter the strength and/or path of the leading edge vortex. The vortical flow is eliminated by forward sweep on the upstream propeller blades. Far field acoustic data from each test indicate only small influences on the level and directivity of the BPFs. The interaction tone at the sum of the two BPF's was significantly altered in a consistent manner. As the vortex system varied, the interaction tone was affected: far field noise levels in the forward quandrant increased and the characteristic noise minimum near the plane of rotation became less pronounced and in some cases were eliminated. If the forward propeller leading edge vortex system does not impact the rear propeller in the standard manner, a net increase in the primary interaction tone occurs for the model tested. If the leading edge vortex is removed, the interaction tone increases.

Laur, Michele↗

Validation of radiative transfer computation with Monte Carlo method for ultra-relativistic background flow

We developed a three-dimensional radiative transfer code for an ultra-relativistic background flow-field by using the Monte Carlo (MC) method in the context of gamma-ray burst (GRB) emission. For obtaining reliable simulation results in the coupled computation of MC radiation transport with relativistic hydrodynamics which can reproduce GRB emission, we validated radiative transfer computation in the ultra-relativistic regime and assessed the appropriate simulation conditions. The radiative transfer code was validated through two test calculations: (1) computing in different inertial frames and (2) computing in flow-fields with discontinuous and smeared shock fronts. The simulation results of the angular distribution and spectrum were compared among three different inertial frames and in good agreement with each other. If the time duration for updating the flow-field was sufficiently small to resolve a mean free path of a photon into ten steps, the results were thoroughly converged. The spectrum computed in the flow-field with a discontinuous shock front obeyed a power-law in frequency whose index was positive in the range from 1 to 10 MeV. The number of photons in the high-energy side decreased with the smeared shock front because the photons were less scattered immediately behind the shock wave due to the small electron number density. The large optical depth near the shock front was needed for obtaining high-energy photons through bulk Compton scattering. Even one-dimensional structure of the shock wave could affect the results of radiation transport computation. Although we examined the effect of the shock structure on the emitted spectrum with a large number of cells, it is hard to employ so many computational cells per dimension in multi-dimensional simulations. Therefore, a further investigation with a smaller number of cells is required for obtaining realistic high-energy photons with multi-dimensional computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Self-synchronizing Schlieren photography and interferometry for the visualization of unsteady transonic flows

The use of self synchronizing stroboscopic Schlieren and laser interferometer systems to obtain quantitative space time measurements of distinguished flow surfaces, steakline patterns, and the density field of two dimensional flows that exhibit a periodic content was investigated. A large field single path stroboscopic Schlieren system was designed, constructed and successfully applied to visualize four periodic flows: near wake behind an oscillating airfoil; edge tone sound generation; 2-D planar wall jet; and axisymmetric pulsed sonic jet. This visualization technique provides an effective means of studying quasi-periodic flows in real time. The image on the viewing screen is a spatial signal average of the coherent periodic motion rather than a single realization, the high speed motion of a quasi-periodic flow can be reconstructed by recording photographs of the flow at different fixed time delays in one cycle. The preliminary design and construction of a self synchronizing stroboscopic laser interferometer with a modified Mach-Zehnder optical system is also reported.

Kadlec, R.↗

Rarefied-continuum gas dynamics transition for SUMS project

This program is to develop an analytic method for reducing SUMS data for the determination of the undisturbed atmosphere conditions ahead of the shuttle along its descending trajectory. It is divided into an internal flow problem, an external flow problem and their matching conditions. Since the existing method of Direct Simulation Monte Carlo (DSMC) failed completely for the internal flow problem, the emphasis is on the internal flow of a highly non-equilibrium, rarefied air through a short tube of a diameter much less than the gaseous mean free path. A two fluid model analysis of this internal flow problem has been developed and studied with typical results illustrated. A computer program for such an analysis and a technical paper published in Lecture Notes in Physics No. 323 (1989) are included as Appendices 3 and 4. A proposal for in situ determination of the surface accommodation coefficients sigma sub t and sigma e is included in Appendix 5 because of their importance in quantitative data reduction. A two fluid formulation for the external flow problem is included as Appendix 6 and a review article for AIAA on Hypersonic propulsion, much dependent on ambient atmospheric density, is also included as Appendix 7.

Cheng, Sin-I↗

Shifts in hydropower operation to balance wind and solar will modify effects on aquatic biota

To avoid negative consequences to freshwater biota from climate change, society must complete the transition from fossil to renewable electricity sources. However, temporal patterns in hydropower generation (and flow releases that affect aquatic biota) may change with increased wind and solar penetration. We used power cost modeling to characterize current and future within-day and seasonal patterns in hydropower generation across the Eastern Interconnection in a wet and a dry year. Compared to the baseline, future hydropower generation across the grid decreased during the day and increased before dawn and after dusk. At a project level, such a pattern would suggest 'double peaking' operation (up- and down-ramping before dawn and after dusk, with lower releases midday). Variation in generation was higher in wet years than dry years, foreshadowing possible flow constraints on hydropower flexibility. At the grid scale, projected ramping rates were higher in all seasons. A review of the ecological literature suggests that these changes would shift the timing of invertebrate drift and elevate the risk of nest scouring during up-ramping and the risk of stranding or dewatering during down ramping. Thermal conditions may be moderated by increased ramping. Strategies for adapting to future shifts in the renewable portfolio range from re-regulation in reservoir cascades to providing flow refuge (structures and vegetation) below individual projects. Coordinated basin-scale operation can distribute peaking operation to maintain grid support while restricting local ramping at critical ecological times. In addition, research to design hybrid renewable systems that add battery storage is needed to understand how we can mitigate future risks to aquatic communities while promoting the use of renewable energy. This study, which is among the first to examine ecological side-effects of the shift to renewable energy in freshwater ecosystems, lays out a path toward understanding and navigating changes to flow regimes under the energy transition.

54 ENVIRONMENTAL SCIENCES↗

Mercury Speciation via Direct Mercury Analyzer

Mercury has no known biological role, but is present in every living thing and widespread in the environment. As mercury enters an aqueous system, it is subject to methylation or demethylation. Methylmercury selectively partitions into biota and biomagnifies through the food chain. Nearly all methylmercury exposures in the U.S. occur through eating fish and shellfish that contain higher levels of methylmercury. Methyl mercury and other organomercury compounds are more toxic than inorganic mercury. Thus, monitoring and understanding mercury speciation are key to effective environmental management. It is particularly important to be able to separate and measure bioaccumulative-toxic species of mercury such as methyl mercury. We are developing a protocol for utilizing the Direct Mercury Analyzer (DMA-80) for mercury speciation of water samples and to monitor methylmercury (meHg). Objectives: - Advance the understanding of the performance of the direct mercury analysis (DMA) under a range of DOE relevant water chemistry conditions to quantify total mercury: - Demonstrate the applicability of DMA to analysis of typical DOE high level liquid 'tank' waste samples using nonradioactive simulants; Test innovative sorbents to extend the range of the direct mercury analyzer to environmental concentration levels. - Perform scoping tests to demonstrate proof of principle for the chemical reduction amalgamation modular microcolumns for mercury speciation: Examine various microcolumn configurations and document speciation performance for typical DOE high level liquid 'tank' waste samples using nonradioactive simulants; Perform tests to examine extending the speciation paradigm to environmental concentration levels. General Method: all samples analyzed using a direct mercury analyzer (DMA 80) which automates the following steps: 1) pyrolysis to release the mercury, 2) catalysis to convert all mercury to elemental, 3) amalgamation to trap the mercury, 4) desorption to release the mercury as a pulse and 5) spectroscopy to detect the mercury.. Instrument analyzes all mercury species in solids and liquids. The DMA 'eliminates' sample preparation, reduces labor, decreases waste, and results in more practical and cost effective methods. To test the performance of the DMA for DOE relevant conditions, we analyzed samples representative of at high (mg/L) tank waste concentrations in water and in diluted tank waste simulant and analyzed samples representative of environmental (ng/L levels) in water. High concentration samples were analyzed directly in the instrument. Low concentration samples were analyzed by concentrating the mercury onto sorbents such as silica thiol or a mixture of silica thiol and copper. The concentration was performed using microcolumns and placing the columns into the DMA (see below) or by mixing the sorbent into the bulk sample and then filtering ad analyzing the filter in the DMA. Speciation Methods: The speciation tests were based on the selective reduction of inorganic mercury ions - Hg{sup (+} {sup or} {sup 2+)} - to elemental mercury - Hg(0). The elemental mercury is then trapped by amalgamation into the metal. These reactions do not impact methyl mercury allowing for separation. All reactions were performed using powdered elemental copper. To demonstrate proof-of-principle for the selective collection of inorganic mercury, samples were exposed to different amounts of powdered copper (and a control with no copper). This was done using microcolumns and a vacuum box. The columns were placed directly on the DMA for analysis. All materials and conditions were tested using solutions of both inorganic mercury and methyl mercury and all trials were in triplicate. Results and Observations: Validated DMA for high concentration Hg samples such as DOE tank waste: - Liquid inHg - samples were stable on autosampler carousel and results were accurate; - Liquid meHg - samples required a stabilizer. Si-Thiol (0.02 g) in the sample boat generated reliable meHg results. Stabilizer can be added for routine operations. Successful 'proof of principal' in extending DMA range to near environmental levels - using both microcolumns and filtration; - Filtration - Adding Si-Thiol (0.5 g) to sample extracted both inHg and meHg and preserved sample. Handling filters was labor intensive. - Microcolumns - best performance using a combination of copper and Si-Thiol - columns clogged and flowed slowly in some cases; - Practical limit may be in the range of 1 to 5 ng/L due to blanks; - Path Forward - Test additional packing materials and methods to improve flow and consistency, and approaches to minimize blanks. Successful 'proof of principal' for Hg speciation using chemical reduction and amalgamation; - Copper microcolumns (0.5 g Cu) selectively trapped inHg for high concentration mercury samples; - Copper microcolumns (0.5 g Cu selectively trapped inHg for large volume samples at environmental concentrations; - Some sorption of me Hg so that separations not quantitative; - Path Forward - test additional packing materials and alternate eluents to better trap inHg and minimize meHg retention in column.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inertially enhanced mass transport using 3D-printed porous flow-through electrodes with periodic lattice structures

Significance The efficient utilization of electrical energy is an increasingly important challenge, especially as renewable energy sources become cheaper and demand increases. Electrochemical reactors utilizing flow-through electrodes (FTEs) provide an attractive path toward the efficient utilization of electrical energy. Their commercial viability and ultimate adoption hinge on attaining high current densities to drive cost competitiveness. There are limited opportunities for engineering FTE materials, as these are often random, disordered media. Alternatively, three-dimensional (3D)–printed FTEs provide the opportunity to quickly explore the impact of engineered electrode architectures on device performance. We demonstrate that 3D-printed FTEs have the potential to exceed the performance of conventional materials by using the expanded design freedom to engineer the internal flow.

30 DIRECT ENERGY CONVERSION↗

Modifying Ship Air-Wake Vortices for Aircraft Operations

Columnar-vortex generators (CVG) have been proposed as means to increase the safety of takeoffs and landings of aircraft on aircraft or helicopter carriers and other ships at sea. According to the proposal, CVGs would be installed at critical edge locations on ships to modify the vortices in the air wakes of the ships. The desired effects of modifications are to smooth airflows over takeoff and landing deck areas and divert vortices from takeoff and landing flight paths. With respect to aircraft operations, the wake flows of primary interest are those associated with the bow and side edges of aircraft-carrier decks and with superstructures of ships in general (see Figure 1). The bow and deck-edge vortices can adversely affect airplane and helicopter operations on carriers, while the superstructure wakes can primarily affect operations of helicopters. The concept of the CVG is not new; what is new is the proposed addition of CVGs to ship structures to effect favorable modifications of air wakes. Figure 2 depicts a basic CVG, vertical and horizontal CVGs installed on a simple superstructure, and horizontal CVGs installed on the bow and deck edges. The vertical CVGs would be closed at the deck but open at the top. Each horizontal CVG would be open at both ends. The dimensions of the CVGs installed on the aft edges of the superstructure would be chosen so that the portion of the flow modified by the vertical CVGs would interact synergistically with the portion of the flow modified by the horizontal CVG to move the air wake away from the takeoff-and-landing zone behind the superstructure. The deck-edge CVGs would be mounted flush with, and would extend slightly ahead of the bow of, the flight deck. The overall length of each tube would exceed that of the flight deck. Each deck-edge CVG would capture that portion of the airflow that generates a deck-edge vortex and would generate a columnar vortex of opposite sense to that of the unmodified vortex. The vortex generated by the CVG could be dispersed at its base, thereby removing unwanted turbulence in the path of an approaching airplane. The deck-edge CVGs would promote smooth flow over the entire flight deck. In the case of a Nimitz-class aircraft carrier like that of Figure 1, there would be a CVG on each of the outer edges of the two left portions of the flight deck and a single CVG on the right side of the flight deck. The forward-most CVG on the left side would take the generated vortex underneath the angled flight deck. A CVG could also be installed on the bow of the flight deck to smooth the flow of air onto the flight deck. In the case of wind incident on the deck from an azimuth other than straight ahead, the vortex generated by the bow CVG could, perhaps, be used to feed the CVG(s) of the leeward side edge of the flight deck.

Lamar, John E.↗

Tri-Sectional Approximation of the Shortest Path to Long-Term Voltage Stability Boundary with Distributed Energy Resources

Ensuring long-term voltage stability is critical for reliable operations of power grids. High share of distributed energy resources (DERs) can create complicated system operation modes that may invalidate the traditional long-term voltage stability analysis based on typical operation modes. To address this challenge, this paper investigates how to compute the shortest path to the voltage stability boundary in the DER aggregated load space with large dispersion. Instead of working in the Euclidean space, we establish the analysis and computations on the algebraic power flow manifold to better capture the curvature change of the shortest path along the direction of losing stability. A modified optimal control framework is presented for obtaining the ground-truth of the smooth shortest path on the manifold. To efficiently and accurately solve for the shortest path, we further leverage the geometric features of the power flow manifold and propose a tri-sectional approximation model that is scalable for large-scale systems. Several numerical examples, up to the 1354-bus system, with different DER penetration levels and high dimensional renewable power injection variations are evaluated. The simulation results demonstrate that the tri-sectional approximation achieves high accuracy and efficiency to approximate the shortest path to the voltage stability boundary.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗