Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “difficult”

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 775 records · Page 43

Simulation study of chemo-mechanical impacts of CO 2 injection in morrow b sandstone reservoir

Here, this study investigates the impacts of geomechanical and geochemical changes on carbon storage in a partially depleted oil reservoir, using results from four different coupled simulation models. Models were used to examine the relative importance of storage mechanisms, and how changing reservoir parameters might affect these mechanisms through time. The study uses data from a Morrowan sandstone reservoir in the Farnsworth Unit (FWU), Ochiltree County, Texas which is currently undergoing CO 2 enhanced oil recovery (EOR). Partially depleted oil reservoirs such as the FWU offer attractive carbon utilization and/or storage targets because of existing infrastructure and economic benefits from incremental oil recovery as well as tax credits. However, prediction of storage capacity or long-term fluid migration in these fields can be difficult because of the wide variation in formation fluids and operational histories that may have undergone. CO 2 injection can cause complex geomechanical and geochemical responses in a reservoir as a result of interplay between dynamic changes in pore pressure, reservoir temperature, fluid composition, and interactions between formation fluids, CO 2 , and reservoir rock. Thus, multiple coupled numerical simulation models must be developed and used to more precisely understand what CO 2 storage mechanisms are most significant, as well as the long-term fate of the stored CO 2 . Our study used results from hydrodynamic, coupled hydro-geomechanical, coupled hydro-geochemical, and coupled hydro-geomechanical-geochemical models to examine how changes in geomechanical and geochemical properties can impact the injectivity or storage capacity of CO 2 . Models simulated historical field operations and then forward-modeled a water-alternate-gas (WAG) operation for 20 years, followed by a 1000-year post-injection monitoring. The work demonstrates that in this specific reservoir, geomechanical impacts appear to be more significant than any geochemical processes, which showed minimal impact on reservoir properties. Furthermore, these simulation scenarios corroborated other studies showing more than 90% of CO 2 injected are being stored and will remain stored in the reservoir. The study was able to draw on a rich set of data from this currently operational project for the work. Because the Morrow is a widespread oil-bearing reservoir in the midcontinent of the United States and is fairly typical of many sandstone reservoirs elsewhere, the study has significance beyond the boundaries of our field area.

58 GEOSCIENCES↗

Spatial and temporal variation in forest transpiration across a forested boreal peatland complex

Transpiration is a globally important component of evapotranspiration. Careful upscaling of transpiration from point measurements is thus crucial for quantifying water and energy fluxes. In spatially heterogeneous landscapes common across the boreal biome, upscaled transpiration estimates are difficult to determine due to variation in local environmental conditions (e.g., basal area, soil moisture, permafrost). Here, we sought to determine stand-level attributes that influence transpiration scalars for a forested boreal peatland complex consisting of sparsely treed wetlands and densely treed permafrost plateaus as land cover types. The objectives were to quantify spatial and temporal variability in stand-level transpiration, and to identify sources of uncertainty when scaling point measurements to the stand-level. Using heat ratio method sap flow sensors, we determined sap velocity for black spruce and tamarack for 2-week periods during peak growing season in 2013, 2017 and 2018. We found greater basal area, drier soils, and the presence of permafrost increased daily sap velocity in individual trees, suggesting that local environmental conditions are important in dictating sap velocity. When sap velocity was scaled to stand-level transpiration using gridded 20 × 20 m resolution data across the ~10 ha Scotty Creek ForestGEO plot, we observed significant differences in daily plot transpiration among years (0.17–0.30 mm), and across land cover types. Daily transpiration was lowest in grid-cells with sparsely treed wetlands compared to grid-cells with well-drained and densely treed permafrost plateaus, where daily transpiration reached 0.80 mm, or 30% of the daily evapotranspiration. When transpiration scalars (i.e., sap velocity) were not specific to the different land cover types (i.e., permafrost plateaus and wetlands), scaled stand-level transpiration was overestimated by 42%. To quantify the relative contribution of tree transpiration to ecosystem evapotranspiration, we recommend that sampling designs stratify across local environmental conditions to accurately represent variation associated with land cover types, especially with different hydrological functioning as encountered in rapidly thawing boreal peatland complexes.

54 ENVIRONMENTAL SCIENCES↗

Soil water percolation and nutrient fluxes as a function of topographical, seasonal and soil texture variation in Central Amazonia, Brazil

Abstract Understanding soil water dynamics and transport of nutrients is challenging in tropical rainforests due to the uniqueness of several properties related to soils, vegetation and seasonality that make relating patterns found in temperate environments to tropical sites difficult. We address the need for better edaphic characterization in tropical environments by investigating soil water percolation rates and chemistry across topographic, soil texture and seasonal gradients in a mature tropical rainforest in Central Amazonia, Brazil. We utilized a passive wick flux meter (e.g., drainage lysimeter) to directly measure real‐time percolation fluxes at 60‐cm depth, and to sample a suite of chemical species across plateau, slope and valley topographic positions. We found percolation flux volume and chemical exports generally increase with decreasing elevation and clay content, which was lowest in the valley. Daily percolation flux was observed to be 2.39 ± 0.44 in plateau, 3.01 ± 0.50 in slope and 6.16 ± 0.83 mm in valley. Most solutes were present in small amounts of <1 mg L −1 , such as PO₄ 3− , Fe 2+ /Fe 3+ and Mn 2+ ; however, NO 3 − concentrations were >20 mg L −1 , even exceeding 100 mg L −1 in the valley. Based on additional isotopic analysis, we speculate high NO 3 − concentrations are partially an artefact of root decomposition following installation of the flux meters. The empirical relationships we show among percolation volume and nutrient exports under varying topographies and soil textures can improve Earth System Model performance by better constraining ecohydrological relationships to nutrient fluxes, which can in‐turn better illuminate the important factors that govern their behaviour.

54 ENVIRONMENTAL SCIENCES↗

Hidden Features: How Subsurface and Landscape Heterogeneity Govern Hydrologic Connectivity and Stream Chemistry in a Montane Watershed

ABSTRACT Hydrologic connectivity is defined as the connection among stores of water within a watershed and controls the flux of water and solutes from the subsurface to the stream. Hydrologic connectivity is difficult to quantify because it is goverened by heterogeniety in subsurface storage and permeability and responds to seasonal changes in precipitation inputs and subsurface moisture conditions. How interannual climate variability impacts hydrologic connectivity, and thus stream flow generation and chemistry, remains unclear. Using a rare, four‐year synoptic stream chemistry dataset, we evaluated shifts in stream chemistry and stream flow source of Coal Creek, a montane, headwater tributary of the Upper Colorado River. We leveraged compositional principal component analysis and end‐member mixing to evaluate how seasonal and interannual variation in subsurface moisture conditions impacts stream chemistry. Overall, three main findings emerged from this work. First, three geochemically distinct end members were identified that constrained stream flow chemistry: reach inflows, and quick and slow flow groundwater contributions. Reach inflows were impacted by historic base and precious metal mine inputs. Bedrock fractures facilitated much of the transport of quick flow groundwater and higher‐storage subsurface features (e.g., alluvial fans) facilitated the transport of slow flow groundwater. Second, the contributions of different end members to the stream changed over the summer. In early summer, stream flow was composed of all three end members, while in late summer, it was composed predominantly of reach inflows and slow flow groundwater. Finally, we observed minimal differences in proportional composition in stream chemistry across all four years, indicating seasonal variability in subsurface moisture and spatial heterogeneity in landscape and geologic features had a greater influence than interannual climate fluctuation on hydrologic connectivity and stream water chemistry. These findings indicate that mechanisms controlling solute transport (e.g., hydrologic connectivity and flow path activation) may be resilient (i.e., able to rebound after perturbations) to predicted increases in climate variability. By establishing a framework for assessing compositional stream chemistry across variable hydrologic and subsurface moisture conditions, our study offers a method to evaluate watershed biogeochemical resilience to variations in hydrometeorological conditions.

Johnson, Keira [College of Earth, Ocean, and Atmos↗

Misclassification of primary liver cancer in the Life Span Study of atomic bomb survivors

Primary liver cancer is difficult to diagnose accurately at death, due to metastases from nearby organs and to concomitant diseases, such as chronic hepatitis and cirrhosis. Trends in diagnostic accuracy could affect radiation risk estimates for incident liver cancer by altering background rates or by impacting risk modification by sex and age. We quantified the potential impact of death‐certificate inaccuracies on radiation risk estimates for liver cancer in the Life Span Study of atomic bomb survivors. True‐positive and false‐negative rates were obtained from a previous study that compared death‐certificate causes of death with those based on pathological review, from 1958 to 1987. We assumed various scenarios for misclassification rates after 1987. We obtained estimated true positives and estimated false negatives by stratified sampling from binomial distributions with probabilities given by the true‐positive and false‐negative rates, respectively. Poisson regression methods were applied to highly stratified person‐year tables of corrected case counts and accrued person years. During the study period (1958–2009), there were 1,885 cases of liver cancer, which included 383 death‐certificate‐only (DCO) cases; 1,283 cases with chronic liver disease as the underlying cause of death; and 150 DCO cases of pancreatic cancer among 105,444 study participants. Across the range of scenarios considered, radiation risk estimates based on corrected case counts were attenuated, on average, by 13–30%. Our results indicated that radiation risk estimates for liver cancer were potentially sensitive to death‐certificate inaccuracies. Additional data are needed to inform misclassification rates in recent years.

French, Benjamin↗

Unravelling Disorder in Aperiodic Crystals – Diffuse Scattering and Atomic Resolution Holography

The atomic–scale disorder of aperiodic crystals, and quasicrystals in particular, is inherently difficult to explore by experimental methods due to their complex atomic arrangements. Two advanced characterization techniques, a revived and an emerging one, offer direct experimental access even to such complex atomic structures: Diffuse Scattering and Atomic Resolution Holography. Finally, in this overview, we introduce their specific application to aperiodic crystals and discuss their merits and difficulties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation studies of polypeptoids using replica exchange with dynamical scaling and dihedral biasing

Polypeptoids differ from polypeptides in that the amide bond can more frequently adopt both cis and trans conformations. The transition between the two conformations requires overcoming a large energy barrier, making it difficult for conventional molecular simulations to adequately visit the cis and trans structures. A replica-exchange method is presented that allows for easy rotations of the amide bond and also an efficient linking to a high temperature replica. The method allows for just three replicas (one at the temperature and Hamiltonian of interest, a second high temperature replica with a biased dihedral potential, and a third connecting them) to overcome the amide bond sampling problem and also enhance sampling for other coordinates. Here our results indicate that for short peptoid oligomers, the conformations can range from all cis to all trans with an average cis/trans ratio that depends on side chain and potential model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Projector‐Based Quantum Embedding Study of Iron Complexes

Projection‐based embedding theory (PBET) is used to calculate and assess the challenging spin‐crossover energies for a selection of small Fe‐containing systems by embedding the metal center into the frozen potential of the ligands. MP2, CCSD, and CCSD(T) are embedded in potentials from the SCAN and r 2 SCAN functionals and compared with the canonical values for the constituent methods and previously reported reference values. Considering the PBET calculations as a correction for the underlying DFT, the embedding calculations are able to provided improvement for most cases. In some cases, the PBET methods are able to compensate for limitations in the wave function methods and produce results similar to more rigorous calculations from the literature. For the systems with spin‐crossover energies near zero, the current methodology fails to provide consistent improvement. In conclusion, the isolated recalculation of the electronic structure around the metal center when embedded into a DFT treatment of the ligand field shows promise as a pragmatic and lower cost treatment compared to the canonical treatment of the whole system of the difficult class of spin‐crossover complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hyper-Raman spectroscopy of CeO 2

Cerium dioxide (CeO 2 ) belongs to the fluorite structure group (Oh 5 point group) resulting in one Raman active T 2g mode and two IR active T 1u TO and LO modes. The T 2g Raman band in CeO 2 has been used to probe the symmetric Ce-O vibrational mode of different sized nanoparticulates. Far infrared studies of the Ce-O interactions are difficult to conduct due to low optical transmission, efficiency of the detector and water absorption bands interfering with the measurements. The selection rules of hyper-Raman spectroscopy (HRS) allow the study of IR bands in the far infrared spectral region, such as the motion of the Ce atoms. For the first time, HRS was used to investigate the T 1u TO and LO IR modes in CeO 2 . The position and full-width half-maximum (FWHM) of the TO and LO bands in the HRS spectrum were used to investigate the spectral response in relation to particle size (i.e., phonon confinement). The investigation of CeO2 with HRS provides a research platform for future work in the study of ThO 2 , UO 2 and PuO 2 IR modes.

74 ATOMIC AND MOLECULAR PHYSICS↗

Transitory sensitivity in automatic chemical kinetic mechanism analysis

Abstract Detailed chemical kinetic mechanisms are necessary for resolving many important chemical processes. As the chemistry of smaller molecules has become better grounded and quantum chemistry calculations have become cheaper, kineticists have become interested in constructing progressively larger kinetic mechanisms to model increasingly complex chemical processes. These large kinetic mechanisms prove incredibly difficult to refine and time‐consuming to interpret. Traditional sensitivity analysis on a large mechanism can range from inconvenient to practically impossible without special techniques to reduce the computational cost. We first present a new time‐local sensitivity analysis we term transitory sensitivity analysis. Transitory sensitivity analysis is demonstrated in an example to accurately identify traditionally sensitive reactions at an 18,000x speed up over traditional sensitivities. By fusing transitory sensitivity analysis with more traditional time‐local branching, pathway, and cluster analyses, we develop an algorithm for efficient automatic mechanism analysis. This automatic mechanism analysis at a time point is able to identify the reactions a target is most sensitive to using transitory sensitivity analysis and then propose hypotheses why the reaction might be sensitive using branching, pathway, and cluster analyses. We implement these algorithms within the reaction mechanism simulator (RMS) package, which enables us to report the automatic mechanism analysis results in highly readable text formats and in molecular flux diagrams.

Johnson, Matthew S.↗

Seasonal drivers of dissolved oxygen across a tidal creek–marsh interface revealed by machine learning

Abstract Dissolved oxygen (DO) is a key biogeochemical control in coastal systems, and its concentration and drivers vary markedly through time and space. This makes it difficult to accurately represent coastal DO and associated biogeochemical processes in models, limiting our ability to predict how these systems will respond to global change. We obtained high‐frequency (5‐min) in situ measurements of DO collected at three locations across the interface of a tidal creek and coastal marsh in the Pacific Northwest, USA. Random Forest machine learning models quantified the importance of three categories of environmental drivers (Aquatic, Climatic, and Terrestrial) of DO variability across the creek–marsh interface. We selected two 4‐month datasets representing Summer and Winter seasonal periods to test two hypotheses on the dominant drivers of DO at the coastal interface. We found that the Terrestrial driver—characterized by long periods of anaerobic conditions and episodic pulses in DO after floods—was most important during the Winter, whereas the Aquatic driver—characterized by variability over tidal, diel, and lunar cycles—was most important during the Summer. We explored how future climate change scenarios could alter the drivers of DO variability using a cumulative sums driver–response framework. Our results suggest that under climate change, Aquatic and Climatic drivers may increase in importance during the Summer, potentially linked to changing metabolic regimes and sea level, with Terrestrial driver importance potentially increasing during the Winter. Our approach highlights useful methods for understanding the spatiotemporal complexity of oxygen across coastal interfaces and quantifying the relative importance of distinct environmental drivers.

54 ENVIRONMENTAL SCIENCES↗

Self-diagnosis of model suitability for continuous measurements of stream-dissolved organic carbon derived from in situ UV–visible spectroscopy

Application of high-frequency monitoring of dissolved organic carbon (DOC) is difficult in instances where training datasets are challenging to develop (e.g., remote locations) and the relationship between optical features and DOC concentration changes due to environmental or landscape shifts (e.g., climate or land-use change). We developed and compared three partial least squares (PLS) models using in situ water level measurements, conductivity, and UV–Vis spectral attenuation to predict DOC. Two site-specific models were developed using data from a hillslope-dominated forest or a low-relief wetland-pond-dominated stream catchment. The third model, using data from both sites, exhibited the best performance (DOC range = 4–15.5 mg C L −1 , mean = 8.38 mg C L −1 , training RMSE = 0.34 mg C L −1 , internal validation RMSE = 0.50 mg C L −1 , external validation RMSE = 2.43 mg C L −1 ). We further demonstrate using PLS model statistics to monitor performance and elucidate when and how models should be updated. These statistics, Hotelling's T 2 and squared prediction errors, are useful consistency checks for the predictions made and detect underlying inconsistencies that, if undetected, can reduce the robustness of DOC prediction. For example, via the T 2 statistic, we identified the summer–autumn transition as a period when DOC composition differed from what was represented in the training dataset. We also determined that elevated SUVA 254 values contributed to the overall bias observed in predictions made during the subsequent year as part of the external validation. This enabled the application of a bias correction that reduced the RMSE from 2.43 to 0.89 mg C L −1 . The method presented here could be applied to future monitoring programs enabling model updates to monitor DOC fluxes accurately from optical datasets (e.g., attenuance or fluorescence) in the face of developing datasets in remote locations or environmental change. In conclusion, implementation of this approach may also identify possible regime shifts or landscape and hydrologic change associated with climate and other environmental changes relevant to terrestrial to aquatic fluxes.

UV-VIS spectroscopy↗

Noninvasive Near‐Field Spectroscopy of Single Subwavelength Complementary Resonators

Abstract Subwavelength metallic resonators provide a route to achieving strong light–matter coupling by means of tight confinement of resonant electromagnetic fields. Investigation of such resonators however often presents experimental difficulties, particularly at terahertz (THz) frequencies. A single subwavelength resonator interacts weakly with THz beams, making it difficult to probe it using far‐field methods, whereas arrays of resonators exhibit inter‐resonator coupling, which affect the resonator spectral signature and field confinement. Here, traditional far‐field THz spectroscopy is systematically compared with aperture‐type THz near‐field microscopy for investigating complementary THz resonators. While the far‐field method proves impractical for measuring single resonators, the near‐field technique gives high signal‐to‐noise spectral information, only achievable in the far‐field with resonator arrays. At the same time, the near‐field technique allows to analyze single resonators without significant interaction with the near‐field probe. Furthermore, the near‐field technique allows highly confined fields and surface waves to be mapped in space and time. This information gives invaluable insight into spectral response in resonator arrays. This near‐field microscopy and spectroscopy method enables investigations of strong light–matter coupling at THz frequencies in the single‐resonator regime.

Hale, Lucy L.↗

A Study on Recycled Polymers Recovered from Multilayer Plastic Packaging Films by Solvent–Targeted Recovery and Precipitation (STRAP)

Multilayer plastic film use increased in multiple packaging applications due to its versatility and overall increased performance over monolayer structures. However, the performance gains from multiple layers also make recycling difficult because they contain multiple polymers that can be immiscible and burdensome to traditional mechanical recycling operations. A possible solution is the solvent-targeted recovery and precipitation (STRAP) process, but the effect on the retrieved polymer is still unknown. Here, the STRAP process is applied to two different multilayer films and samples of recovered polymers are evaluated for physical, molecular, and thermal properties. Changes in the molecular weight are not significant, but differences in thermal properties are reported along with the coprecipitation of different polymers. Solvent retention in the polymer matrix from STRAP reduced the glass transition temperature of samples, but enhanced drying recovered it. Heavy metals, such as Cd, Cr, and Pb are not detected, indicating regulatory compliance for different applications.

36 MATERIALS SCIENCE↗

Denitrification and the challenge of scaling microsite knowledge to the globe

Here, our knowledge of microbial processes—who is responsible for what, the rates at which they occur, and the substrates consumed and products produced—is imperfect for many if not most taxa, but even less is known about how microsite processes scale to the ecosystem and thence the globe. In both natural and managed environments, scaling links fundamental knowledge to application and also allows for global assessments of the importance of microbial processes. But rarely is scaling straightforward: More often than not, process rates in situ are distributed in a highly skewed fashion, under the influence of multiple interacting controls, and thus often difficult to sample, quantify, and predict. To date, quantitative models of many important processes fail to capture daily, seasonal, and annual fluxes with the precision needed to effect meaningful management outcomes. Nitrogen cycle processes are a case in point, and denitrification is a prime example. Statistical models based on machine learning can improve predictability and identify the best environmental predictors but are—by themselves—insufficient for revealing process-level knowledge gaps or predicting outcomes under novel environmental conditions. Hybrid models that incorporate well-calibrated process models as predictors for machine learning algorithms can provide both improved understanding and more reliable forecasts under environmental conditions not yet experienced. Incorporating trait-based models into such efforts promises to improve predictions and understanding still further, but much more development is needed.

59 BASIC BIOLOGICAL SCIENCES↗

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗

Nonoverlapping block smoothers for the Stokes equations

Overlapping block smoothers efficiently damp the error contributions from highly oscillatory components within multigrid methods for the Stokes equations but they are computationally expensive. This paper is concentrated on the development and analysis of new block smoothers for the Stokes equations that are discretized on staggered grids. These smoothers are nonoverlapping and therefore desirable due to reduced computational costs. Traditional geometric multigrid methods are based on simple pointwise smoothers. However, using multigrid methods to efficiently solve more difficult problems such as the Stokes equations leads to computationally more expensive smoothers, for example, overlapping block smoothers. Nonoverlapping smoothers are less expensive, but have been considered less efficient in the literature. In this paper, we develop new nonoverlapping smoothers, the so-called triad-wise smoothers, and show their efficiency within multigrid methods to solve the Stokes equations. In addition, we compare overlapping and nonoverlapping smoothers by measuring their computational costs and analyzing their behavior by the use of local Fourier analysis.

97 MATHEMATICS AND COMPUTING↗

Toward efficient polynomial preconditioning for GMRES

Here, we present a polynomial preconditioner for solving large systems of linear equations. The polynomial is derived from the minimum residual polynomial (the GMRES polynomial) and is more straightforward to compute and implement than many previous polynomial preconditioners. Our current implementation of this polynomial using its roots is naturally more stable than previous methods of computing the same polynomial. We implement further stability control using added roots, and this allows for high degree polynomials. We discuss the effectiveness and challenges of root-adding and give an additional check for stability. In this article, we study the polynomial preconditioner applied to GMRES; however it could be used with any Krylov solver. This polynomial preconditioning algorithm can dramatically improve convergence for some problems, especially for difficult problems, and can reduce dot products by an even greater margin.

97 MATHEMATICS AND COMPUTING↗