Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic equivalent”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Application of Circulation Control Technology to Airframe Noise Reduction

This report is a summary of the work performed by Georgia Tech Research Institute (GTRI) under NASA Langley Grant NAG-1-2146, which was awarded as a part of NASA's Breakthrough Innovative Technologies (BIT) initiative. This was a three-year program, with a one-year no-cost extension. Each year's study has been an integrated effort consisting of computational fluid dynamics, experimental aerodynamics, and detailed noise and flow measurements. Year I effort examined the feasibility of reducing airframe noise by replacing the conventional wing systems with a Circulation Control Wing (CCW), where steady blowing was used through the trailing edge of the wing over a Coanda surface. It was shown that the wing lift increases with CCW blowing and indeed for the same lift, a CCW wing was shown to produce less noise. Year 2 effort dealt with a similar study on the role of pulsed blowing on airframe noise. The main objective of this portion of the study was to assess whether pulse blowing from the trailing edge of a CCW resulted in more, less, or the same amount of radiated noise to the farfield. Results show that a reduction in farfield noise of up to 5 dB is measured when pulse flow is compared with steady flow for an equivalent lift configuration. This reduction is in the spectral region associated with the trailing edge jet noise. This result is due to the unique advantage that pulsed flow has over steady flow. For a range of frequencies, more lift is experienced with the same mass flow as the steady case. Thus, for an equivalent lift and slot height, the pulsed system can operate at lower jet velocities, and hence lower jet noise. The computational analysis showed that for a given time-averaged mass flow rate, pulsed jets give a higher value of C(sub l) and a higher L/D than equivalent steady jets. This benefit is attributable to higher instantaneous jet velocities, and higher instantaneous C(sub mu) values for the pulsed jet. Pulsed jet benefits increase at higher frequencies. However, these advantages are somewhat offset by the unsteadiness in the loads, which will cause structural vibrations and fatigue. Additional studies must be done, perhaps with multiple jets on the upper and lower surfaces, to smooth out the fluctuations in lift while retaining the benefits. The rest of the effort was devoted to examining ways of reducing flap edge noise by blowing air through a Coanda nozzle over a rounded tip of the flap. In this case, we were successful in moving the tip vortex away from the tip, but the device producing the blowing was noisy and we were unable to examine the noise benefits, although we believe that the movement of the tip vortex far from the tip should provide noise benefits. It should be noted that in an effort to understand the fluid dynamics and the aeroacoustics of a jet blowing over a Coanda surface, we also carried out a very extensive study of the high aspect ratio slot jets. A first-ever set of far-field noise spectra were measured for jets exhausting from slots with aspect ratios in the range 100 to 3000. Parallel measurements of velocity profiles, length scales and convection velocities were measured to understand the noise generation of high aspect ratio jets. Attempts were also made to develop jet noise prediction schemes for such jets. Much of the work done under this effort has been described in five conference papers and two doctoral theses. The first year s work on the use of steady blowing was described in two AIAA papers presented at the 2001 AIAA Aerospace Sciences Meeting in Reno. Subsequent work was presented at the 9th AIMCEAS Aeroacoustics Conference and Exhibit held at Hilton Head May 12-13. Another paper is to be presented at the 2004 AIAA Aerospace Sciences Meeting in Reno in January 2004. All six papers are included with this report as Appendices. The bulk of the experimental work done in an effort to produce a pulsed flow that is free of upstream noise is also attached as an Appendix.

Ahuja, K. K.↗
DOE OSTI2020

Comparative modeling of the disregistry and Peierls stress for dissociated edge and screw dislocations in Al

Many elementary deformation processes in metals involve the motion of dislocations. The planes of glide and specific processes dislocations prefer depend heavily on their atomic core structures. Atomistic simulations are desirable for dislocation modeling but their application to even sub-micron scale problems is in general computationally costly. Accordingly, continuum-based approaches, such as the phase-field microelasticity, phase-field dislocation dynamics (PFDD), generalized Peierls–Nabarro (GPN) models, and the concurrent atomistic–continuum (CAC) method, have attracted increasing attention in the field of dislocation modeling because they well represent both short-range cores interactions and long-range stress fields of dislocations. To better understand their similarities and differences, it is useful to compare these methods in the context of benchmark simulations and predictions. In this paper, we apply the CAC method and different PFDD variants – one of them is equivalent to a GPN model – to simulate an extended (i.e., dissociated) dislocation in Al with initially pure edge or pure screw character in terms of the disregistry. CAC and discrete forms of PFDD are also employed to calculate the Peierls stress. By conducting comprehensive convergence studies, we quantify the dependence of these measures on time/grid resolution and simulation cell size. Several important but often overlooked differences between PFDD/GPN variants are clarified. In conclusion, our work sheds light on the advantages and limitations of each method, as well as the path towards enabling them to effectively model complex dislocation processes at larger length scales.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Identification of aerodynamic indicial functions using flight data

It is pointed out that the use of indicial function representation provides a model superior to the aerodynamic derivative model. Specific derivatives can be approximated from the indicial models. The model can also be used to compute equivalent stability and control parameters not usually available from flight data. It is shown that derivatives regarding the angle-of-attack and the side slip angle can be derived directly from the indicial functions without any identifiability problem. Attention is given to the pitch moment coefficient, linear indicial function representation, the identification problem for the pitch moment equation, the identifiability of linear systems, parametric representations of the indicial functions, an identification technique, angle-of-attack and pitch rate dynamics in the pitch plane, multivariate linear models, nonlinear aerodynamic indicial functions, measurement system accuracy, and poststall and spin-entry data from a scaled research vehicle.

Gupta, N. K.↗
NASA NTRS2012

Characterization of Deficiencies in the Frequency Domain Forced Response Analysis Technique for Supersonic Turbine Bladed Disks

Turbine blades in rocket and jet engine turbomachinery experience enormous harmonic loading conditions. These loads result from the integer number of upstream and downstream stator vanes as well as the other turbine stages. Assessing the blade structural integrity is a complex task requiring an initial characterization of whether resonance is possible and then performing a forced response analysis if that condition is met. The standard technique for forced response analysis in rocket engine turbines is to decompose a computational fluid dynamics (CFD).generated flow field into its harmonic components, and to then perform a frequency response analysis at the problematic natural frequencies using cyclically symmetric structural dynamic models. Recent CFD analysis and water-flow testing at NASA/MSFC, though, indicates that this technique may miss substantial harmonic and non ]harmonic excitation sources that become present in complex flows. This complex content can only be captured by a CFD flow field encompassing at least an entire revolution. A substantial development effort to create a series of software programs to enable application of the 360 degree forcing function in a frequency response analysis on cyclic symmetric models has been completed (to be described in a future paper), but the question still remains whether the frequency response analysis itself is capable of capturing the excitation content sufficiently. Two studies comparing frequency response analysis with transient response analysis, therefore, of bladed-disks undergoing this complex flow environment have been performed. The first is of a bladed disk with each blade modeled by simple beam elements and the disk modeled with plates (using the finite element code MSC/NASTRAN). The focus of this model is to be representative of response of realistic bladed disks, and so the dimensions are roughly equivalent to the new J2X rocket engine 1st stage fuel pump turbine. The simplicity of the model allows the CFD load to be able to be readily applied, along with analytical and experimental variations in both the temporal and spatial fourier components of the excitation. In addition, this model is a first step in identifying response differences between transient and frequency forced response analysis techniques. The second phase assesses this difference for a much more realistic solid model of a bladed-disk in order to evaluate the effect of the spatial variation in loading on blade dominated modes. Neither research on the accuracy of the frequency response method when used in this context or a comprehensive study of the effect of test-observed variation on blade forced response have been found in the literature, so this research is a new contribution to practical structural dynamic analysis of gas turbines. The primary excitation of the upstream nozzles interacts with the blades on fuel pump of the J2X causes the 5th Nodal diameter modes to be excited, as explained by Tyler and Sofrin1, so a modal analysis was first performed on the beam/plate model and the 5ND bladed-disk mode at 40167 hz was identified and chosen to be the one excited at resonance (see figure 1). The first forced response analysis with this model focuses on identifying differences between frequency and transient response analyses. A hypothesis going into the analysis was that perhaps the frequency response was enforcing a temporal periodicity that did not really exist, and so therefore it would overestimate the response. As high dynamic response was a considerable source of stress in the J2X, examining this concept could potentially be beneficial for the program.

Brown, Andrew M.↗
NASA NTRS1981

Transient dynamic analysis of high-speed lightly loaded cylindrical roller bearings. 2: Computer program and results

The governing differential equations of motion for a high speed cylindrical roller bearing are programmed for numerical solution and plotted output. This computer program has the capability of performing a two dimensional or three dimensional simulation. Two numerical solutions of the governing differential equations were obtained to simulate the motion of a roller bearing, one for the two dimensional equations of motion and one for the three dimensional equations of motion. Computer generated plots were obtained and present such data as roller/cage interaction forces, roller/race traction forces, roller/race relative slip velocities and cage angular speed over a nondimensional time equivalent to 1.2 revolutions of the inner race. Roller axial displacement, roller skew angle, and skew moment are also plotted for the three dimensional solution. The trajectory of the cage center is plotted for both the two dimensional and three dimensional solutions.

Conry, T. F.↗
NASA NTRS1992

Using adaptive grid in modeling rocket nozzle flow

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which cannot be solved analytically. However, this system of equations called the Navier-Stokes equations can be solved numerically. The accuracy and the convergence of the solution of the system of equations will depend largely on how precisely the sharp gradients in the domain of interest can be resolved. With the advances in computer technology, more sophisticated algorithms are available to improve the accuracy and convergence of the solutions. An adaptive grid generation is one of the schemes which can be incorporated into the algorithm to enhance the capability of numerical modeling. It is equivalent to putting intelligence into the algorithm to optimize the use of computer memory. With this scheme, the finite difference domain of the flow field called the grid does neither have to be very fine nor strategically placed at the location of sharp gradients. The grid is self adapting as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzles by taking the refinement part of grid generation out of the hands of computational fluid dynamics (CFD) specialists and place it into the computer algorithm itself.

Chow, Alan S.↗
NASA NTRS1993

Hybrid state-space self-tuning control of uncertain linear systems

The paper presents a hybrid state-space self-tuner using a new dual-rate sampling scheme for digital adaptive control of continuous-time uncertain linear systems. A state-space-based recursive least-squares algorithm, together with a variable forgetting factor, is used for direct estimations of both the equivalent discrete-time uncertain linear system parameters and the associated discrete-time state of a continuous-time uncertain linear system from the sampled input and output data. An analogue optimal regional pole-placement design method is used for designing an optimal observer-based analogue controller. A suboptimal observer-based digital controller is then designed from the designed analogue controller using digital redesign technique. To enhance the robustness of parameter identification and state estimation algorithms, a dynamic bound for a class of uncertain bilinear parameters and a fast-rate digital controller are developed at each fast-sampling period. Also, to accommodate computation loads and computation delay for developing the advanced hybrid self-tuner, the designed analogue controller and observer gains are both updated at each slow-sampling period. This control technique has been successfully applied to benchmark control problems.

Shieh, L. S.↗
NASA NTRS2004

Induced Power of the Helicopter Rotor

A simplified rotor model was used to explore fundamental behavior of lifting rotor induced power at moderate and high advance ratios. Several rotor inflow theories, including dynamic inflow theory and prescribed-wake vortex theory, together with idealized notional airfoil stall models were employed. A number of unusual results were encountered at high advance ratios including trim control reversal and multiple trim solutions. Significant increases in rotor induced power (torque) above the ideal minimum were observed for moderately high advance ratio. Very high induced power was observed near and above unity advance ratio. The results were sensitive to the stall characteristics of the airfoil models used. An equivalent wing analysis was developed to determine induced power from Prandtl lifting line theory and help interpret the rotor induced power behavior in terms of the spanwise airload distribution. The equivalent wing approach was successful in capturing the principal variations of induced power for different configurations and operating conditions. The effects blade root cutout were found to have a significant effect on rotor trim and induced power at high advance ratios.

Ormiston, Robert A.↗
NASA NTRS2018

Generation of Performance Model for the Aeolian Wind Tunnel (AWT) Rotor at Reduced Pressure

The NASA Jet Propulsion Laboratory (JPL) designed the Mars Helicopter (MH) in collaboration with AeroVironment Inc., NASA Ames Research Center, and NASA Langley Research Center to explore the possibility of a vertical takeoff and landing (VTOL) Unmanned Aerial Vehicle (UAV) for flight on Mars. A 40-inch-diameter Aeolian Wind Tunnel (AWT) rotor, roughly approximating the proposed MH design by JPL, was tested in forward flight at Mars atmospheric pressure at the NASA Ames Planetary Aeolian Laboratory (PAL) in support of MH research efforts. This report describes the generation of the rotor model used to correlate with that experimental effort as reported by Ament and Koning. The 40-inch-diameter rotor was 3D-scanned and transformed into an airfoil deck. The scanned rotor airfoil sections are analyzed using C81 Generator (C81Gen) to generate the sectional aerodynamic coefficients for comprehensive analyses. A mid-fidelity computational fluid dynamics (CFD) simulation using Rotorcraft CFD (RotCFD) is pursued to efficiently estimate rotor hover and forward flight performance. Simulations at two pressures, 7 mbar (approximate Martian atmospheric pressure) and 1018 mbar (1 atmosphere), are performed to gain an understanding of the performance differences and Reynolds number effects observed. Experimental 1-atmosphere thrust for single- and dual-rotor isolated hover cases correlate well with the modeled rotor. Performance results at reduced pressure (7 mbar) show a drastic decrease in lift for equivalent RPMs tested at 1 atmosphere. Although this is primarily due to pressure reduction, Reynolds number effects also contribute to this decrease, as airfoil lift and drag coefficients are affected when compared with 1-atmosphere results. Further, simulated rotor power coefficient shows drastic increases at reduced pressures, attributed to laminar boundary layer separation, as described in Koning et al. for the MH rotor analysis. PAL experimental Martian Surface Wind Tunnel (MARSWIT) results are presented in the paper by Ament and Koning. The very low Reynolds number range is currently not well understood and presents various challenges for both experimentation and simulation.

Koning, Witold J. F.↗
DOE OSTI2022

The influence of intake flow and coolant temperature on gasoline spray morphology during early-injection DISI engine operation

Multi-hole gasoline injectors operating at conditions spanning throttled early-intake stroke operation produce spray plumes that either remained separated or merge and collapse due to flash boiling. Flash boiling occurs due to the sudden expansion of gas bubbles in the liquid fuel at high fuel temperature and low ambient pressure. This study records high-speed images of spray-morphology changes due to in-cylinder flow, thereby revealing operating conditions that do and do not affect the self-induced morphology observed in quiescent vessels. Specifically, in a central-injection, four-valve, high-tumble engine, where the thermodynamic state and in-cylinder cross flow are dynamic. Motivated by cold start and hot restart operation, the fuel pressure, coolant temperature, in-cylinder air pressure, and engine rpm were systematically varied over relevant operating conditions, which bracketed the range from non- to flash-boiling sprays. The results reveal the operating conditions at which the in-cylinder cross flow disrupts the spray morphology as well as the extent of the disruption. At 650 rpm, the spray morphology was similar to that observed in quiescent vessels at nominally equivalent fuel temperature and in-cylinder pressure, indicating that the spray’s self-induced entrainment flow dominated the in-cylinder flow. However, for fuel temperature and ambient pressure near the transition between non- and flash-boiling, the intake cross flow at higher engine speed (1950 rpm) significantly disrupted the spray morphology. The high cross-flow velocity appears to induce plume merging and collapse, whereas none was evident at low rpm (650 rpm). This study led to the postulate that the spray merging and collapse are governed by the rate of atomization near the nozzle exit, presumed to be controlled by either or both aerodynamic atomization and flash-boiling intensity. It would then follow that spray modeling in CFD requires atomization models that blend the effects of both physical processes.

42 ENGINEERING↗
NASA NTRS2001

Relativistic Iron Emission and Disk Reflection in Galactic Microquasar XTE J1748-288

We report evidence for an Fe K-alpha fluorescence line feature and disk reflection in the very high, high-, and low-state X-ray spectra of the Galactic microquasar XTE J1748 - 288 during its 1998 June outburst. Spectral analyses are made on data gathered throughout the outburst by the Rossi X-Ray Timing Explorer Proportional Counter Array. Gaussian line, relativistic disk emission line, and ionized disk reflection models are fitted to the data. In the very high state the line profile appears strongly redshifted, consistent with disk emission from the innermost stable orbits around a maximally rotating Kerr black hole. In the high state the line profile is less redshifted and increasingly prominent. The low-state line profile is very strong (approx. 0.5 keV equivalent width) and centered at 6.7 +/- 0.10 keV; disk line emission model fits indicate that the inner edge of the disk fluctuates between approx. 20R(sub g) and - approx. 100R(sub g) in this state. The disk reflection fraction is traced through the outburst; reflection from an ionized disk is preferred in the very high and high states, and reflection from a relatively neutral disk is preferred in the low state. We discuss the implications of our findings for the binary system dynamics and accretion flow geometry in XTE J1748 - 288.

Miller, J. M.↗
NASA NTRS2004

Test of the Equivalence Principle in an Einstein Elevator

The scientific goal of the experiment is to test the equality of gravitational and inertial mass (i.e., to test the Principle of Equivalence) by measuring the independence of the rate of fall of bodies from their compositions. The measurement is accomplished by measuring the relative displacement (or equivalently acceleration) of two falling bodies of different materials which are the proof masses of a differential accelerometer spinning about an horizontal axis to modulate a possible violation signal. A non-zero differential acceleration appearing at the signal frequency will indicate a violation of the Equivalence Principle. The goal of the experiment is to measure the Eotvos ratio og/g (differential acceleration/common acceleration) with a targeted accuracy that is about two orders of magnitude better than the state of the art (presently at several parts in 10(exp 13). The analyses carried out during this first grant year have focused on: (1) evaluation of possible shapes for the proof masses to meet the requirements on the higher-order mass moment disturbances generated by the falling capsule; (2) dynamics of the instrument package and differential acceleration measurement in the presence of errors and imperfections; (3) computation of the inertia characteristic of the instrument package that enable a separation of the signal from the dynamics-related noise; (4) a revised thermal analysis of the instrument package in light of the new conceptual design of the cryostat; (5) the development of a dynamic and control model of the capsule attached to the gondola and balloon to define the requirements for the leveling mechanism (6) a conceptual design of the leveling mechanism that keeps the capsule aligned before release from the balloon; and (7) a new conceptual design of the customized cryostat and a preliminary valuation of its cost. The project also involves an international cooperation with the Institute of Space Physics (IFSI) in Rome, Italy. The group at IFSI is in charge of prototyping the differential accelerometer and carrying out precursor laboratory measurements. During this grant year, our partners analyzed and then designed a new prototype of differential accelerometer that has several characteristics in common with the flight accelerometer at this point of the instrument development. The highlights of these activities are documented in a section of this report.

Shapiro, Irwin I.↗
NASA NTRS2008

Interactive Inverse Design Optimization of Fuselage Shape for Low-Boom Supersonic Concepts

This paper introduces a tool called BOSS (Boom Optimization using Smoothest Shape modifications). BOSS utilizes interactive inverse design optimization to develop a fuselage shape that yields a low-boom aircraft configuration. A fundamental reason for developing BOSS is the need to generate feasible low-boom conceptual designs that are appropriate for further refinement using computational fluid dynamics (CFD) based preliminary design methods. BOSS was not developed to provide a numerical solution to the inverse design problem. Instead, BOSS was intended to help designers find the right configuration among an infinite number of possible configurations that are equally good using any numerical figure of merit. BOSS uses the smoothest shape modification strategy for modifying the fuselage radius distribution at 100 or more longitudinal locations to find a smooth fuselage shape that reduces the discrepancies between the design and target equivalent area distributions over any specified range of effective distance. For any given supersonic concept (with wing, fuselage, nacelles, tails, and/or canards), a designer can examine the differences between the design and target equivalent areas, decide which part of the design equivalent area curve needs to be modified, choose a desirable rate for the reduction of the discrepancies over the specified range, and select a parameter for smoothness control of the fuselage shape. BOSS will then generate a fuselage shape based on the designer's inputs in a matter of seconds. Using BOSS, within a few hours, a designer can either generate a realistic fuselage shape that yields a supersonic configuration with a low-boom ground signature or quickly eliminate any configuration that cannot achieve low-boom characteristics with fuselage shaping alone. A conceptual design case study is documented to demonstrate how BOSS can be used to develop a low-boom supersonic concept from a low-drag supersonic concept. The paper also contains a study on how perturbations in the equivalent area distribution affect the ground signature shape and how new target area distributions for low-boom signatures can be constructed using superposition of equivalent area distributions derived from the Seebass-George-Darden (SGD) theory.

Li, Wu↗
DOE OSTI2022

Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states

Here, we develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems, including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input–output and the hierarchical equations of motion formalisms into a double hierarchy of density matrix equations. We show analytically that under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an N-photon Fock state input induces no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon with or without the influence of phonons. For short pulses (or, equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter Δ and an exciton–light coupling constant Γ. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyze the absorption probability in the opposite long pulse (narrow bandwidth) regime. We then derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective vs independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy equations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2014

The Broad-Band X-Ray Spectrum of IC 4329A from a Joint NuSTAR/Suzaku Observation

We have obtained a deep, simultaneous observation of the bright, nearby Seyfert galaxy IC 4329A with Suzaku and NuSTAR (Nuclear Spectroscopic Telescope Array). Through a detailed spectral analysis, we are able to robustly separate the continuum, absorption, and distant reflection components in the spectrum. The absorbing column is found to be modest (approximately 6 times 10 (sup 21) per square centimeter), and does not introduce any significant curvature in the Fe K band. We are able to place a strong constraint on the presence of a broadened Fe K alpha line (E (sub rest) equals 6.46 (from plus 0.08 to minus 0.07) kiloelectronvolts with sigma equals 0.33 (from plus 0.08 to minus 0.07) kiloelectronvolts and EW equals 34 (from plus 8 to minus 7) electronvolts), though we are not able to constrain any of the parameters of a relativistic reflection model. These results highlight the range in broad Fe K line strengths observed in nearby, bright, active galactic nuclei (roughly an order of magnitude), and imply a corresponding range in the physical properties of the inner accretion disk in these sources. We have also updated our previously reported measurement of the high-energy cutoff of the hard X-ray emission using both observatories rather than just NuSTAR alone: E (sub cut) equals 186 plus or minus 14 kiloelectronvolts. This high-energy cutoff acts as a proxy for the temperature of the coronal electron plasma, enabling us to further separate this parameter from the plasmas optical depth and to update our results for these parameters as well. We derive kT 50 (from plus 6 to minus 3) kiloelectronvolts with tau equals 2.34 (from plus +0.16 to minus 0.11) using a spherical geometry, kT equals 61 plus or minus1 kiloelectronvolt with tau equals 0.68 plus or minus 0.02 for a slab geometry, with both having an equivalent goodness-of-fit.

Dynamics - stars↗
DOE OSTI2020

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↗
DOE OSTI2020

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

Spatially resolved reflectivities of Saturn - 3000-6000 A

Absolute reflectivity reductions are presented for spectrophotometric observations of seven spatially-resolved regions on the disk of Saturn in the 3000-6000 A region. An asymmetry is found between northern and southern hemisphere data, in the equivalent width of the CH4, 5430-A band and in the continuum reflectivity. In addition, the spectra are found not to include the 4700 A reflectivity peak found by Bergstralh et al (1981), and to exhibit general agreement with the Podolack and Danielson (1977) model. The blue and visual spectral region continuum reflectivity of a planetary atmosphere provides constraints for determining aerosol particle single-scattering albedo and vertical distribution, which are in turn important to an understanding of planetary energy balance, atmospheric dynamics, and planetary molecular absorption spectrum interpretation.

Cochran, W. D.↗