Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Input Sensitivity”

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

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

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗
DOE OSTI2017

Mineral properties, microbes, transport, and plant-input profiles control vertical distribution and age of soil carbon stocks

The role of organo-mineral interactions, microbial dynamics, and vertical plant input profiles are hypothesized to be important in controlling soil organic matter (SOM) stocks and dynamics. In this paper, to test this hypothesis, we enhanced and applied a model (Biotic and Abiotic Model of SOM – BAMS1) that represents microbial dynamics and organo-mineral interactions integrated with a multiphase reactive transport solver for variably saturated porous media. The model represents aqueous chemistry, aqueous advection and diffusion, gaseous diffusion, sorption processes, bacterial and fungal activity, and archetypal monomer- and polymer-carbon substrate groups, including dead cell wall material. This model structure is fundamentally different, and produces different SOM dynamics, than the pseudo-first order relationships that are used in most site- and global-scale terrestrial SOM models. We simulated two grasslands, a seven-site chronosequence (~3.9–240 Ky) in Northern California, and a Russian Chernozem site where soils were sampled 100 years apart. We calibrated the model's vertically-resolved soil bulk specific surface area (SBSSA) using observed bulk SOM content, and then tested the model against observed Δ 14 C profiles. The modeled microbial processes, organo-mineral interactions, and vertical aqueous transport produced realistic vertically-resolved predictions of bulk SOM content, Δ 14 C values of SOM, lignin content, and fungi-to-aerobic bacteria biomass ratios. Using sensitivity analyses, we found that vertical carbon input profiles were important controls over the Δ 14 C depth distribution. Shallower carbon input profiles lead to older carbon at depth. In addition, the SBSSA was the dominant control over the magnitude and vertical distribution of SOM stocks. Finally, the findings of this study demonstrate the value of explicitly incorporating microbial activity, sorption, and vertical transport into land models to predict SOM dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Solar parameters in long-baseline accelerator neutrino oscillations

Long-baseline (LBL) accelerator neutrino oscillation experiments, such as NOvA and T2K in the current generation, and DUNE-LBL and HK-LBL in the coming years, will measure the remaining unknown oscillation parameters with excellent precision. These analyses assume external input on the so-called “solar parameters,” θ 12 and $Δm^{2}_{21}$, from solar experiments such as SNO, SK, and Borexino, as well as reactor experiments like KamLAND. Here we investigate their role in long-baseline experiments. We show that, without external input on $Δm^{2}_{21}$ and θ 12 , the sensitivity to detecting and quantifying CP violation is significantly, but not entirely, reduced. Thus long-baseline accelerator experiments can actually determine $Δm^{2}_{21}$ and θ 12 , and thus all six oscillation parameters, without input from any other oscillation experiment. In particular, $Δm^{2}_{21}$ can be determined; thus DUNE-LBL and HK-LBL can measure both the solar and atmospheric mass splittings in their long-baseline analyses alone. While their sensitivities are not competitive with existing constraints, they are very orthogonal probes of solar parameters and provide a key consistency check of a less probed sector of the three-flavor oscillation picture. Furthermore, we also show that the true values of $Δm^{2}_{21}$ and θ 12 play an important role in the sensitivity of other oscillation parameters such as the CP violating phase δ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Sensitivity Analysis of MFiX-PIC Parameters Using Nodeworks, PSUADE, and DAKOTA

The study presented in this report was aimed to demonstrate UQ analysis performed not only with Nodeworks, but also two other well-established UQ software tools from the U.S. DOE’s National Laboratories (PSUADE from Lawrence Livermore National Laboratory and DAKOTA from Sandia National Laboratory). It is important to emphasize that the motivation of this study was not to determine the best UQ software, but to verify if the global sensitivity analyses from the end-to-end workflow in Nodeworks are consistent with the results of other two UQ software. The components of Nodeworks from Python’s ecosystem have been tested as standalone libraries. However, an assessment study for the complete workflow targeting a specific UQ analysis has not been performed for Nodeworks. Hence, this study is expected to serve as an equivalent of solution verification for Nodeworks using other established UQ tools as reference solution. For this purpose, three distinct flow configurations (i.e., settling bed, bubbling fluidized, and circulating fluidized bed) have been used as representative multiphase flow problems of interest. The results of the systematic simulation campaigns performed in an earlier study using the particle-in-cell (PIC) approach in the Multiphase Flow with Interphase eXchanges (MFIX) suite of solvers (i.e., MFiX-PIC) was utilized. The same set of tabulated results was provided as input to the different UQ software for global sensitivity analysis. Results for the three cases indicate that based on the Sobol’ Sensitivity Indices method the order of importance ranking determined by Nodeworks for the Sobol’ Total Sensitivity Indices is consistent with PSUADE and DAKOTA in each case for the five model parameters considered. The input files for Nodeworks for the three cases are also shared through NETL’s Gitlab repository for the reader interested in reproducibility and further analysis (See Section 1.2).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

FREYA Developments 2015-2016

We have been studying the sensitivity of neutron observables to the FREYA input distributions. Patrick Talou generated 1000 different yield functions, Y(A,Z,TKE), for 252 Cf(sf) based on fits to available data on the yields as a function of fragment mass, Y(A), and the total kinetic energy distribution as a function of heavy fragment mass, TKE(AH), and its width, σTKE(AH), to obtain Y(A,TKE). These results are augmented by Wahl systematics to include the fragment charge. We have so far found that the neutron observables do not depend strongly on the yield functions. We also found that the yield functions, constrained by data, are not sufficient to constrain the average neutron multiplicity within the evaluated uncertainty. A paper on this is in progress. We began working on a paper on photofission with the University of Michigan while on a Consortium for Verification Technology fellowship. This is an important paper because they have measured the neutron multiplicity distribution, the prompt fission neutron spectrum, and the energy correlation of two neutrons as a function of angle of separation for the first time in photofission of 235 U. We have made FREYA calculations to compare with the data, but the paper is still in progress since the student left and other researchers have had to take over the analysis. We expect to submit the paper in the next several months, after the analysis is complete.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

NASA rotor systems research aircraft: Fixed-wing configuration flight-test results

The fixed-wing, airplane configuration flight-test results of the Rotor System Research Aircraft (RSRA), NASA 740, at Ames/Dryden Flight Research Center are documented. Fourteen taxi and flight tests were performed from December 1983 to October 1984. This was the first time the RSRA was flown with the main rotor removed; the tail rotor was installed. These tests confirmed that the RSRA is operable as a fixed-wing aircraft. Data were obtained for various takeoff and landing distances, control sensitivity, trim and dynamics stability characteristics, performance rotor-hub drag, and acoustics signature. Stability data were obtained with the rotor hub both installed and removed. The speed envelope was developed to 261 knots true airspeed (KTAS), 226 knots calibrated airspeed (KCAS) at 10,000 ft density altitude. The airplane was configured at 5 deg. wing incidence with 5 deg. wing flaps as a normal configuration. Level-flight data were acquired at 167 KCAS for wing incidence from 0 to 10 deg. Step inputs and doublet inputs of various magnitudes were utilized to acquire dynamic stability and control sensitivity data. Sine-wave inputs of constantly increasing frequency were used to generate parameter identification data. The maximum load factor attained was 2.34 g at 206 KCAS.

Erickson, R. E.↗
DOE OSTI2015

A revised model for microbially induced calcite precipitation: Improvements and new insights based on recent experiments

Abstract The model for microbially induced calcite precipitation (MICP) published by Ebigbo et al. (2012) has been improved based on new insights obtained from experiments and model calibration. The challenge in constructing a predictive model for permeability reduction in the underground with MICP is the quantification of the complex interaction between flow, transport, biofilm growth, and reaction kinetics. New data from Lauchnor et al. (2015) on whole‐cell ureolysis kinetics from batch experiments were incorporated into the model, which has allowed for a more precise quantification of the relevant parameters as well as a simplification of the reaction kinetics in the equations of the model. Further, the model has been calibrated objectively by inverse modeling using quasi‐1D column experiments and a radial flow experiment. From the postprocessing of the inverse modeling, a comprehensive sensitivity analysis has been performed with focus on the model input parameters that were fitted in the course of the model calibration. It reveals that calcite precipitation and concentrations of and are particularly sensitive to parameters associated with the ureolysis rate and the attachment behavior of biomass. Based on the determined sensitivities and the ranges of values for the estimated parameters in the inversion, it is possible to identify focal areas where further research can have a high impact toward improving the understanding and engineering of MICP.

Hommel, Johannes↗
DOE OSTI2022

A general framework for quantifying uncertainty at scale

Abstract In many fields of science, comprehensive and realistic computational models are available nowadays. Often, the respective numerical calculations call for the use of powerful supercomputers, and therefore only a limited number of cases can be investigated explicitly. This prevents straightforward approaches to important tasks like uncertainty quantification and sensitivity analysis. This challenge can be overcome via our recently developed sensitivity-driven dimension-adaptive sparse grid interpolation strategy. The method exploits, via adaptivity, the structure of the underlying model (such as lower intrinsic dimensionality and anisotropic coupling of the uncertain inputs) to enable efficient and accurate uncertainty quantification and sensitivity analysis at scale. Here, we demonstrate the efficiency of this adaptive approach in the context of fusion research, in a realistic, computationally expensive scenario of turbulent transport in a magnetic confinement tokamak device with eight uncertain parameters, reducing the effort by at least two orders of magnitude. In addition, we show that this refinement method intrinsically provides an accurate surrogate model that is nine orders of magnitude cheaper than the high-fidelity model.

Farcaş, Ionuţ-Gabriel (ORCID:0000000238841580)↗
DOE OSTI2022

Estimating Fast Charging Infrastructure Requirements to Fully Electrify Ride-Hailing Fleets Across the United States

Ride-hailing electrification is made challenging by a variety of factors, including limited driver access to overnight charging and a nascent network of public direct-current fast charging (DCFC). This article projects the DCFC infrastructure needed to electrify ride-hailing across 384 cities in the United States. Research leveraging a novel infrastructure projection tool, Electric Vehicle Infrastructure Tool for On-Demand Mobility Services model (EVI-OnDemand), estimates that 25,400 DCFC ports capable of charging vehicles at a mean power of 150 kW are required to electrify 1.45 million ride-hailing vehicles in the baseline scenario. Results show that ride-hailing vehicles require 17.5 DCFC ports per 1000 vehicles, three times that of personal electric vehicles (EVs). Infrastructure demand is also shown to vary across cities; 63% more DCFC ports are needed to support 1000 battery EVs in Chicago versus Miami. The analysis also includes sensitivity studies that explore variation in charging demand across input assumptions, finding that charging power, vehicle class, overnight charging access rates, and port utilization influence the network size needed. Projections vary from 18,600 to 41,300 DCFC ports across sensitivity scenarios, indicating that ride-hailing electrification will require significant expansion of access to fast charging even in the most optimistic scenarios, and investments into supporting assets can reduce the scale of the public charging expansion needed.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1975

Numerical procedure for analyzing Langmuir probe data

A numerical procedure is proposed for analyzing Langmuir probe data in the presence of a two-group plasma containing both primary and Maxwellian electrons. The procedure is known as a least-squares differential-correction technique for determining the unknown coefficients of the governing equation. It is shown that for a given set of input data the results of the analysis are unique and independent of the initial estimate of electron temperature, that convergence is fastest when electron temperature is overestimated, that the results are sensitive to the region of the curve used as input and also to the voltage increment, and that plasma properties determined by the proposed numerical procedure are either consistent with those determined graphically or closer to expected values. With a suitable data acquisition system, the Fortran IV program worked out for this procedure could be used to provide real-time plasma diagnostic information for an operating ion thruster.

Beattie, J. R.↗
NASA NTRS2008

Uncertainty and Sensitivity Analyses of a Two-Parameter Impedance Prediction Model

This paper presents comparisons of predicted impedance uncertainty limits derived from Monte-Carlo-type simulations with a Two-Parameter (TP) impedance prediction model and measured impedance uncertainty limits based on multiple tests acquired in NASA Langley test rigs. These predicted and measured impedance uncertainty limits are used to evaluate the effects of simultaneous randomization of each input parameter for the impedance prediction and measurement processes. A sensitivity analysis is then used to further evaluate the TP prediction model by varying its input parameters on an individual basis. The variation imposed on the input parameters is based on measurements conducted with multiple tests in the NASA Langley normal incidence and grazing incidence impedance tubes; thus, the input parameters are assigned uncertainties commensurate with those of the measured data. These same measured data are used with the NASA Langley impedance measurement (eduction) processes to determine the corresponding measured impedance uncertainty limits, such that the predicted and measured impedance uncertainty limits (95% confidence intervals) can be compared. The measured reactance 95% confidence intervals encompass the corresponding predicted reactance confidence intervals over the frequency range of interest. The same is true for the confidence intervals of the measured and predicted resistance at near-resonance frequencies, but the predicted resistance confidence intervals are lower than the measured resistance confidence intervals (no overlap) at frequencies away from resonance. A sensitivity analysis indicates the discharge coefficient uncertainty is the major contributor to uncertainty in the predicted impedances for the perforate-over-honeycomb liner used in this study. This insight regarding the relative importance of each input parameter will be used to guide the design of experiments with test rigs currently being brought on-line at NASA Langley.

Jones, M. G.↗
DOE OSTI2015

Modeling stream temperature in the Anthropocene: An earth system modeling approach

A new large-scale stream temperature model has been developed within the Community Earth System Model (CESM) framework. The model is coupled with the Model for Scale Adaptive River Transport (MOSART) that represents river routing and a water management model (WM) that represents the effects of reservoir operations and water withdrawals on flow regulation. The coupled models allow the impacts of reservoir operations and withdrawals on stream temperature to be explicitly represented in a physically based and consistent way. The models have been applied to the Contiguous United States driven by observed meteorological forcing. It is shown that the model is capable of reproducing stream temperature spatiotemporal variation satisfactorily by comparison against the observed streamflow from over 320 USGS stations. Including water management in the models improves the agreement between the simulated and observed streamflow at a large number of stream gauge stations. Both climate and water management are found to have important influence on the spatiotemporal patterns of stream temperature. More interestingly, it is quantitatively estimated that reservoir operation could cool down stream temperature in the summer low-flow season (August – October) by as much as 1~2oC over many places, as water management generally mitigates low flow, which has important implications to aquatic ecosystems. In conclusion, sensitivity of the simulated stream temperature to input data and reservoir operation rules used in the WM model motivates future directions to address some limitations in the current modeling framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Smoothing by cubic spline modified applied to solve inverse thermal problem

This paper presents an alternative to cubic spline regularization and its weighted form applied in solving inverse thermal problems. The inverse heat transfer problems are classified as ill-posed, that is, the solution may become unstable, mainly because they are sensitive to random errors deriving from the input data, necessitating a regularization method to soften these effects. The smoothing technique proposed by cubic spline regularization ensures that the global data tend to be more stable, with fewer data oscillations and dependent on a single arbitrary parameter input. It also shows that the weighted cubic spline is able to enhance filter action. The methods have been implemented in order for the search engine to optimize the choice of parameters and weight and, thus, the smoothing gains more flexibility and accuracy. The simulated and experimental tests confirm that the techniques are effective in reducing the amplified noise by inverse thermal problem presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Technical and economic feasibility of solar flat-plate collector thermal energy systems for small and medium manufacturers

Manufacturing industries require process heating for a wide variety of operations, such as drying, cooking, washing, and sterilizing. Studies of industrial energy consumption have shown that more than 50% of industrial heating demand is in the low (<60 °C), medium (60-150 °C), and medium-high (<250 °C) temperature ranges. A significant fraction of this low temperature thermal energy could be generated by commercially available solar flat-plate collector thermal energy systems, which would reduce the use of fossil fuels, energy cost volatility, and carbon emissions. However, manufacturing industries have not been able to take advantage of solar flat-plate collector systems for two main reasons: First, information is not readily available about when and what flat-plate collector system sizes can be economically applied and, second, manufacturers often do not have appropriate personnel or analysis tools to make decisions about installation of solar flat-plate collector systems. The goal of this research is to provide design decision support to small- and medium-sized manufacturers in conducting feasibility analysis for the application of solar flat-plate collector systems in their operations. The $\bar{φ}, f$-chart method and modified $\bar{φ}, f$-chart design method are used to predict the amount of energy collected by solar flat-plate collectors. A simplified model is proposed to assess the cost of solar flat-plate collector energy systems. To demonstrate the presented feasibility analysis methodology, a decision support software application is developed. What-if scenarios are presented for several geographic locations of varying climates. A sensitivity analysis of savings-to-investment ratio (SIR) was performed by changing selected input parameters by ±30% from the base case. It was found that SIR is sensitive to natural gas price, followed by collector area dependent cost, discount rate, cost of the auxiliary gas heater, and cost of the preheating system.

14 SOLAR ENERGY↗
DOE OSTI2026

Performance and transport in the ARC tokamak

The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓⁡𝑢⁢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓⁡𝑢⁢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓⁡𝑢⁢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Seasonal Evolution of the Subglacial Hydrologic System Modified by Supraglacial Lake Drainage in Western Greenland

The impact of summer surface melt on Greenland Ice Sheet dynamics is modulated by the state of the subglacial hydrologic system. Studies of ice motion indicate that efficiency of the subglacial system increases over the melt season, decreasing the sensitivity of ice motion to surface melt inputs. However, the behavior of the subglacial hydrologic system is complex and some characteristics are still poorly constrained. Here we investigate the coevolution of subglacial hydrology and ice motion in the Pâkitsoq region of western Greenland during the 2011 melt season. We analyze measurements from 11 Global Positioning System stations, from which we derive ice velocity, longitudinal strain rates, and basal uplift, alongside observations of surface ablation and supraglacial lake drainages. We observe ice acceleration after the onset of local surface melting, followed by gradual ice deceleration, consistent with increasing subglacial efficiency. In the study area, supraglacial lake drainages cooccur with a change in regional strain rate patterns and ice deceleration, suggesting that lake drainages contribute to rapid subglacial reorganization. At lower ice surface elevations (below ~900 m above sea level), ice motion is correlated with both total basal uplift and its rate of change, while at higher elevations (~900–1,100 m above sea level), ice motion correlated only with the basal uplift rate. This pattern suggests that continued cavity growth or subglacial sediment dynamics may be important in the apparent increase in subglacial drainage efficiency at higher elevations in the ablation zone. Furthermore, our results further suggest that transient subglacial behavior is important in the seasonal evolution of ice motion.

58 GEOSCIENCES↗
DOE OSTI2019

Can Convolution and Deconvolution Be Used as Tools for Modeling Multicomponent, Mixing-Limited Reaction Networks?

Transfer functions (i.e., convolutions) have been powerful tools for both forward modeling and reconstructing past states of transport systems in a wide range of subsurface flow problems. The advantage of the transfer function is that it is a simple alternative to complicated, distributed parameter models of flow and transport, but the majority of applications of transfer functions in hydrology have been limited to relatively simple cases, like passive tracers or first-order decay. The central question evaluated in this note is whether or not multicomponent mixing-limited reactive transport can be represented within a transfer function framework. Our examples consider forward-in-time predictions and backward-in-time reconstructions of a carbonate system that represents the intrusion of seawater into a freshwater aquifer. The main result is that accurate forward-in-time and backward-in-time models are developed by posing the problem in terms of conservative components. As with all convolution-based methods, the results are sensitive to errors and/or noise in the input functions, but we show that smoothed approximations of the requisite functions provide good representations of transport. Given the vast unknowns in any subsurface transport problem, such generalized, reactive transfer function models may have yet unexplored advantages when the trade-offs between overall computational cost, accuracy, and uncertainty are explored in more detail.

42 ENGINEERING↗
DOE OSTI2017

Results of the GABLS3 diurnal-cycle benchmark for wind energy applications

We present results of the GABLS3 model intercomparison benchmark revisited for wind energy applications. The case consists of a diurnal cycle, measured at the 200-m tall Cabauw tower in the Netherlands, including a nocturnal low-level jet. The benchmark includes a sensitivity analysis of WRF simulations using two input meteorological databases and five planetary boundary-layer schemes. A reference set of mesoscale tendencies is used to drive microscale simulations using RANS k-ϵ and LES turbulence models. The validation is based on rotor-based quantities of interest. Cycle-integrated mean absolute errors are used to quantify model performance. The results of the benchmark are used to discuss input uncertainties from mesoscale modelling, different meso-micro coupling strategies (online vs offline) and consistency between RANS and LES codes when dealing with boundary-layer mean flow quantities. Altogether, all the microscale simulations produce a consistent coupling with mesoscale forcings.

17 WIND ENERGY↗