Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “projected gradient method”

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 325 records · Page 18

Gradient boosting machine for modeling the energy consumption of commercial buildings

Accurate savings estimations are important to promote energy efficiency projects and demonstrate their cost-effectiveness. The increasing presence of advanced metering infrastructure (AMI) in commercial buildings has resulted in a rising availability of high frequency interval data. These data can be used for a variety of energy efficiency applications such as demand response, fault detection and diagnosis, and heating, ventilation, and air conditioning (HVAC) optimization. This large amount of data has also opened the door to the use of advanced statistical learning models, which hold promise for providing accurate building baseline energy consumption predictions, and thus accurate saving estimations. The gradient boosting machine is a powerful machine learning algorithm that is gaining considerable traction in a wide range of data driven applications, such as ecology, computer vision, and biology. In the present work an energy consumption baseline modeling method based on a gradient boosting machine was proposed. To assess the performance of this method, a recently published testing procedure was used on a large dataset of 410 commercial buildings. The model training periods were varied and several prediction accuracy metrics were used to evaluate the model's performance. The results show that using the gradient boosting machine model improved the R-squared prediction accuracy and the CV(RMSE) in more than 80 percent of the cases, when compared to an industry best practice model that is based on piecewise linear regression, and to a random forest algorithm.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Demonstration of a Piston Plug feed System for Feeding Coal/Biomass Mixtures across a Pressure Gradient for Application to a Commercial CBTL System

Producing liquid transportation fuels and power via coal and biomass to liquids (CBTL) and integrated gasification combined cycle (IGCC) processes can significantly improve the nation's energy security. The Energy Independence and Security Act of 2007 mandates increasing renewable fuels nearly 10-fold to >2.3 million barrels per day by 2022. Coal is abundantly available and coal to liquids (CTL) plants can be deployed today, but they will not become sustainable without large scale CO 2 capture and storage. Co-processing of coal and biomass in CBTL processes in a 60 to 40 ratio is an attractive option that has the potential to produce 4 million barrels of transportation fuels per day by 2020 at the same level of CO 2 emission as petroleum. In this work, Southern Research Institute (Southern) has made an attempt to address one of the major barriers to the development of large scale CBTL processes - cost effective/reliable dry-feeding of coal-biomass mixtures into a high pressure vessel representative of commercial entrained-flow gasifiers. Present method for dry coal feeding involves the use of pressurized lock-hopper arrangements that are not only very expensive with large space requirements but also have not been proven for reliably feeding coal-biomass mixtures without the potential problems of segregation and bridging. The project involved the development of a pilot-scale 250 lb/h high pressure dry coal-biomass mixture feeder provided by TKEnergi and proven for feeding biomass at a scale up to 6 ton/day. The aim of this project is to demonstrate cost effective feeding of coal-biomass mixtures (50:50 to 70:30) made from a variety of coals (bituminous, lignite) and biomass (wood, corn stover, switch grass). The feeder uses a hydraulic piston-based approach to produce a series of plugs of the mixture that act as a seal against high back-pressure of the gasification vessel in to which the mixture is being fed. The plugs are then fed one by one via a plug breaker into the high pressure gasification vessel. A number of runs involving the feeding of coal and biomass mixtures containing 50 to 70 weight % coal into a high pressure gasification vessel simulator have shown that plugs of sufficient density can be formed to provide a seal against pressures up to 450 psig if homogeneity of the mixture can be maintained. However, the in-homogeneity of coal-biomass mixtures can occur during the mixing process because of density, particle size and moisture differences. Also, the much lower compressibility of coal as opposed to biomass can contribute to non-uniform plug formation which can result in weak plugs. Based on present information, the piston plug feeder offered marginal economic advantages over lock-hoppers. The results suggest a modification to the piston feeder that can potentially seal against pressure without the need for forming plugs. This modified design could result in lower power requirements and potentially better economics.

01 COAL, LIGNITE, AND PEAT↗

A multivariate variational objective analysis-assimilation method. Part 1: Development of the basic model

The variational method of undetermined multipliers is used to derive a multivariate model for objective analysis. The model is intended for the assimilation of 3-D fields of rawinsonde height, temperature and wind, and mean level temperature observed by satellite into a dynamically consistent data set. Relative measurement errors are taken into account. The dynamic equations are the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation. The model Euler-Lagrange equations are eleven linear and/or nonlinear partial differential and/or algebraic equations. A cyclical solution sequence is described. Other model features include a nonlinear terrain-following vertical coordinate that eliminates truncation error in the pressure gradient terms of the horizontal momentum equations and easily accommodates satellite observed mean layer temperatures in the middle and upper troposphere. A projection of the pressure gradient onto equivalent pressure surfaces removes most of the adverse impacts of the lower coordinate surface on the variational adjustment.

Achtemeier, Gary L.↗

Super Resolving Unrolled Neural Networks for Remote Sensing

In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.

47 OTHER INSTRUMENTATION↗

UV patternable thin film chemistry for shape and functionally versatile self-oscillating gels

We describe a new planar processing chemistry that allows the synthesis and patterning of dynamically self-actuating gels of diverse form. The chemical reaction covalently incorporates a methacrylate-modified ruthenium trisbipyridine monomer within a poly acrylamide (PAAm) gel, modified to provide a flexible chemistry for UV-curable, self-oscillating Belousov–Zhabotinsky (BZ) gels. The photoinitiated polymerization offers significant advantages over previously reported BZ gel chemistries – fast polymerization rates, broad flexibility in the structural and functional attributes that impact the BZ reaction driven mechanics, and a facility for net shape fabrication via projection mode photo lithographic exposure. The facility of the synthetic method enabled an opportunity to broadly explore the structure/property and structure/rate correlations characterizing the dynamics of these gels and develop means to both control and predict the detailed features of the self-oscillating chemical waves they embed. The PAAm systems initiate under the control of transport-dominated influences weighted heavily by the diffusion/depletion of reagents present in the ambient contacting fluid medium. These initiating modes of oscillation are generally unstable temporally, showing shape- and property-sensitive conversions to oscillating modes controlled by the attributes of the gel-phase environment. The present work establishes the nature of these sensitivities, developing insights related to such features in patterned gel objects or films as composition, size, thickness, shape, degree of crosslinking, and gradient forms of structure within the gel. We observed in this work that, across a broad design-rule and compositional space, the PAAm based BZ gel systems exhibit contractions in dimension for regions in which the Ru catalyst is present in the +3 oxidation state (PNIPAAm BZ gels expand in this state). The chemomechanics of the PAAm system are further differentiated from those of PNIPAAm in the generally (and relatively) fast oscillatory rates it can sustain. Furthermore, this work suggests the dynamics in the PAAm gels do not reflect a simple osmotic force driven by differences in the solvation of Ru(II/III) centers. Computational modeling supports a mechanism in which the charge state of Ru centers mediate a differential degree of intra/inter-chain segmental interactions – dynamics that mediate a temporally varying degree of reversibly formed association crosslinks within the gel.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-F-BRB-10: A Statistical Voxel Based Normal Organ Dose Prediction Model for Coplanar and Non-Coplanar Prostate Radiotherapy

Purpose: Knowledge learned from previous plans can be used to guide future treatment planning. Existing knowledge-based treatment planning methods study the correlation between organ geometry and dose volume histogram (DVH), which is a lossy representation of the complete dose distribution. A statistical voxel dose learning (SVDL) model was developed that includes the complete dose volume information. Its accuracy of predicting volumetric-modulated arc therapy (VMAT) and non-coplanar 4π radiotherapy was quantified. SVDL provided more isotropic dose gradients and may improve knowledge-based planning. Methods: 12 prostate SBRT patients originally treated using two full-arc VMAT techniques were re-planned with 4π using 20 intensity-modulated non-coplanar fields to a prescription dose of 40 Gy. The bladder and rectum voxels were binned based on their distances to the PTV. The dose distribution in each bin was resampled by convolving to a Gaussian kernel, resulting in 1000 data points in each bin that predicted the statistical dose information of a voxel with unknown dose in a new patient without triaging information that may be collectively important to a particular patient. We used this method to predict the DVHs, mean and max doses in a leave-one-out cross validation (LOOCV) test and compared its performance against lossy estimators including mean, median, mode, Poisson and Rayleigh of the voxelized dose distributions. Results: SVDL predicted the bladder and rectum doses more accurately than other estimators, giving mean percentile errors ranging from 13.35–19.46%, 4.81–19.47%, 22.49–28.69%, 23.35–30.5%, 21.05–53.93% for predicting mean, max dose, V20, V35, and V40 respectively, to OARs in both planning techniques. The prediction errors were generally lower for 4π than VMAT. Conclusion: By employing all dose volume information in the SVDL model, the OAR doses were more accurately predicted. 4π plans are better suited for knowledge-based planning than the VMAT plans that are strongly biased in its dose gradient orientation. This project is supported by Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

Controlling distributed energy resources via deep reinforcement learning for load flexibility and energy efficiency

Behind-the-meter distributed energy resources (DERs), including building solar photovoltaic (PV) technology and electric battery storage, are increasingly being considered as solutions to support carbon reduction goals and increase grid reliability and resiliency. However, dynamic control of these resources in concert with traditional building loads, to effect efficiency and demand flexibility, is not yet commonplace in commercial control products. Traditional rule-based control algorithms do not offer integrated closed-loop control to optimize across systems, and most often, PV and battery systems are operated for energy arbitrage and demand charge management, and not for the provision of grid services. More advanced control approaches, such as MPC control have not been widely adopted in industry because they require significant expertise to develop and deploy. Recent advances in deep reinforcement learning (DRL) offer a promising option to optimize the operation of DER systems and building loads with reduced setup effort. However, there are limited studies that evaluate the efficacy of these methods to control multiple building subsystems simultaneously. Additionally, most of the research has been conducted in simulated environments as opposed to real buildings. This paper proposes a DRL approach that uses a deep deterministic policy gradient algorithm for integrated control of HVAC and electric battery storage systems in the presence of on-site PV generation. The DRL algorithm, trained on synthetic data, was deployed in a physical test building and evaluated against a baseline that uses the current best-in-class rule-based control strategies. Performance in delivering energy efficiency, load shift, and load shed was tested using price-based signals. The results showed that the DRL-based controller can produce cost savings of up to 39.6% as compared to the baseline controller, while maintaining similar thermal comfort in the building. The project team has also integrated the simulation components developed during this work as an OpenAIGym environment and made it publicly available so that prospective DRL researchers can leverage this environment to evaluate alternate DRL algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Measuring Diffusion of Liquids by Common-Path Interferometry

A method of observing the interdiffusion of a pair of miscible liquids is based on the use of a common-path interferometer (CPI) to measure the spatially varying gradient of the index refraction in the interfacial region in which the interdiffusion takes place. Assuming that the indices of refraction of the two liquids are different and that the gradient of the index of refraction of the liquid is proportional to the gradient in the relative concentrations of either liquid, the diffusivity of the pair of liquids can be calculated from the temporal variation of the spatial variation of the index of refraction. This method yields robust measurements and does not require precise knowledge of the indices of refraction of the pure liquids. Moreover, the CPI instrumentation is compact and is optomechanically robust by virtue of its common- path design. The two liquids are placed in a transparent rectangular parallelepiped test cell. Initially, the interface between the liquids is a horizontal plane, above which lies pure liquid 2 (the less-dense liquid) and below which lies pure liquid 1 (the denser liquid). The subsequent interdiffusion of the liquids gives rise to a gradient of concentration and a corresponding gradient of the index of refraction in a mixing layer. For the purpose of observing the interdiffusion, the test cell is placed in the test section of the CPI, in which a collimated, polarized beam of light from a low-power laser is projected horizontally through a region that contains the mixing layer.

Rashidnia, Nasser↗

QA/QC-ed Groundwater Level Time Series in PLM-1 and PLM-6 Monitoring Wells, East River, Colorado (2016-2022)

This data set contains QA/QC-ed (Quality Assurance and Quality Control) water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location at the East River Watershed, Colorado, USA. These wells are used to monitor subsurface water and carbon inventories and fluxes, and to determine the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry (see related references) can be used to estimate rates of solute export from the hillslope to the floodplain and river. QA/QC analysis of measured groundwater levels in monitoring wells PLM-1 and PLM-6 included identification and flagging of duplicated values of timestamps, gap filling of missing timestamps and water levels, removal of abnormal/bad and outliers of measured water levels. The QA/QC analysis also tested the application of different QA/QC methods and the development of regular (5-minute, 1-hour, and 1-day) time series datasets, which can serve as a benchmark for testing other QA/QC techniques, and will be applicable for ecohydrological modeling. The package includes a Readme file, one R code file used to perform QA/QC, a series of 8 data csv files (six QA/QC-ed regular time series datasets of varying intervals (5-min, 1-hr, 1-day) and two files with QA/QC flagging of original data), and three files for the reporting format adoption of this dataset (InstallationMethods, file level metadata (flmd), and data dictionary (dd) files).QA/QC-ed data herein were derived from the original/raw data publication available at Williams et al., 2020 (DOI: 10.15485/1818367). For more information about running R code file (10.15485_1866836_QAQC_PLM1_PLM6.R) to reproduce QA/QC output files, see README (QAQC_PLM_readme.docx). This dataset replaces the previously published raw data time series, and is the final groundwater data product for the PLM wells in the East River. Complete metadata information on the PLM1 and PLM6 wells are available in a related dataset on ESS-DIVE: Varadharajan C, et al (2022). https://doi.org/10.15485/1660962. These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales. 2022/09/09 Update: Converted data files using ESS-DIVE’s Hydrological Monitoring Reporting Format. With the adoption of this reporting format, the addition of three new files (v1_20220909_flmd.csv, V1_20220909_dd.csv, and InstallationMethods.csv) were added. The file-level metadata file (v1_20220909_flmd.csv) contains information specific to the files contained within the dataset. The data dictionary file (v1_20220909_dd.csv) contains definitions of column headers and other terms across the dataset. The installation methods file (InstallationMethods.csv) contains a description of methods associated with installation and deployment at PLM1 and PLM6 wells. Additionally, eight data files were re-formatted to follow the reporting format guidance (er_plm1_waterlevel_2016-2020.csv, er_plm1_waterlevel_1-hour_2016-2020.csv, er_plm1_waterlevel_daily_2016-2020.csv, QA_PLM1_Flagging.csv, er_plm6_waterlevel_2016-2020.csv, er_plm6_waterlevel_1-hour_2016-2020.csv, er_plm6_waterlevel_daily_2016-2020.csv, QA_PLM6_Flagging.csv). The major changes to the data files include the addition of header_rows above the data containing metadata about the particular well, units, and sensor description. 2023/01/18 Update: Dataset updated to include additional QA/QC-ed water level data up until 2022-10-12 for ER-PLM1 and 2022-10-13 for ER-PLM6. Reporting format specific files (v2_20230118_flmd.csv, v2_20230118_dd.csv, v2_20230118_InstallationMethods.csv) were updated to reflect the additional data. R code file (QAQC_PLM1_PLM6.R) was added to replace the previously uploaded HTML files to enable execution of the associated code. R code file (QAQC_PLM1_PLM6.R) and ReadMe file (QAQC_PLM_readme.docx) were revised to clarify where original data was retrieved from and to remove local file paths.

54 ENVIRONMENTAL SCIENCES↗

Final Technical Report for Department of Energy Award DE-SC0006625, “Predictability of the carbon-climate system on seasonal to decadal time scales.”

The project aims to investigate the feasibility of advancing our understanding of the carbon cycle, using a carbon-weather data assimilation system that updates the modeled carbon dioxide concentration and atmospheric circulation every six hours using CO 2 data (from the OCO 2 satellite) and weather data. At the core of the system is the DOE-NCAR-CAM5fv global circulation model coupled to the National Center for Atmospheric Research's Data Assimilation Testbed, running an ensemble of 30 models. This combination provides realistic vertical carbon dioxide gradients and conservation of dry air mass. A global four-dimensional distribution of atmospheric CO 2 concentration is produced. Our results show (1) that OCO 2 total precipitable water data are reliable and provide valuable uncertainty information for the OCO 2 data assimilation; and (2) that our approach is a promising method for monitoring national carbon dioxide emissions.

54 ENVIRONMENTAL SCIENCES↗

Commuting projections on graphs

For a given (connected) graph, we consider vector spaces of (discrete) functions defined on its vertices and its edges. These two spaces are related by a discrete gradient operator, Grad and its adjoint, ₋Div, referred to as (negative) discrete divergence. We also consider a coarse graph obtained by aggregation of vertices of the original one. Then a coarse vertex space is identified with the subspace of piecewise constant functions over the aggregates. We consider the ℓ 2 -projection Q H onto the space of these piecewise constants. In the present paper, our main result is the construction of a projection π H from the original edge-space onto a properly constructed coarse edge-space associated with the edges of the coarse graph. The projections π H and Q H commute with the discrete divergence operator, i.e., we have div π H = Q H div. The respective pair of coarse edge-space and coarse vertexspace offer the potential to construct two-level, and by recursion, multilevel methods for the mixed formulation of the graph Laplacian which utilizes the discrete divergence operator. The performance of one two-level method with overlapping Schwarz smoothing and correction based on the constructed coarse spaces for solving such mixed graph Laplacian systems is illustrated on a number of graph examples.

97 MATHEMATICS AND COMPUTING↗

Studies using IMP, Voyager and Pioneer cosmic ray data to determine the size of the heliosphere

The purpose of this project was to use the cosmic ray data from the IMP, Voyager and Pioneer spacecraft in the heliosphere out to approximately 65 AU to estimate the size of the heliosphere. We used several techniques to develop a consistent picture of the size of the heliosphere. The first method used a response function approach which determined the intensity as a function of time by scaling the modulation effect as they move outward and eventually reach the boundary of the heliosphere. In this model the effects of transient cosmic ray disturbances is included. A second approach using the perturbation method in which drifts are considered as a perturbation to the standard diffusion-convection modulation models was not fully developed. In a third approach the location of the modulation boundary beyond the termination shock was estimated using observations of the intensity and radial gradients between Voyager 2 and Pioneer 10 along with new estimates of the interstellar intensity of more than 70 MeV galactic cosmic rays. Using this method we found that for 7 years, from 1983 to 1990, the modulation boundary remained constant at 83 +/- 5 AU. We infer from these studies that a modulation boundary can be estimated only by extrapolating the observed radial gradients when the solar magnetic field polarity is such that cosmic-ray particles are drifting in the heliosphere inward toward the Earth along the neutral sheet. The boundary distance is larger than the estimates of the location of the termination shock at 67 +/- 5 AU using the same method. Two other studies partially supported by this grant are attached. The first deals with the recovery period of the greater than 70 MeV cosmic rays in the outer heliosphere from 1992-1995. In the second paper we compare the rigidity dependence of the 11-year cosmic ray variation at the Earth in two cycles of opposite solar magnetic field polarity.

Lockwood, John A.↗

Microwave Imaging on Metal Objects

This final report for the project discusses the attempts to model, using different methods, microwave image reconstruction. Maximum Entropy Method was not successful. Attempts to use Singular Value Decomposition (SVD) got some good results after initial failure. SVD is based upon a theory of linear algebra, to the effect that any M X N Matrix A whose number of rows M is greater than or equal to its number of columns, N can be written as the product of an M X N column-orthogonal matrix U, an N X N diagonal Matrix, W, with m positive or zero elements (the singular values) and the transposition of an N X N orthogonal matrix V. In microwave imaging, the scattered fields can be expressed by the induced current distribution. The SVD method required more contiguous computer memory than was available. Work was also done on the Conjugate Gradient Method (CGM), which didn't work well when tried earlier. It was found that separation of the imaginary part and the real part during calculation may work. This work was considered incomplete as of the end of the grant period.

Tolliver, C. L.↗

EXCHANGE Campaign 1: A Community-Driven Baseline Characterization of Soils, Sediments, and Water Across Coastal Gradients

The EXploration of Coastal Hydrobiogeochemistry Across a Network of Gradients and Experiments (EXCHANGE) program is a consortium of scientists working together to improve our understanding of how the two-way exchange of water between estuaries or large lake lacustuaries and the terrestrial landscape influence the state and function of ecosystems across the coastal interface. EXCHANGE Campaign 1 (EC1) focuses on the spatial variation in biogeochemical structure and function at the coastal terrestrial-aquatic interface (TAI). In the Fall of 2021, the EXCHANGE Consortium gathered samples from 52 TAIs. Samples collected from EC1 were analyzed for bulk geochemical parameters, bulk physicochemical parameters, organic matter characteristics, and redox-sensitive elements.Please download ec1_README.pdf for a complete list of available data in each .zip folder, package version history, and detailed information about the project. This README will serve as the central place for EC1 Data Package updates. Experimental setup and v1 methods are documented in Myers-Pigg and Pennington et al., 2023 (https://doi.org/10.1038/s41597-023-02548-7).EC1 Data Package Structure:ec1_README.pdfec1_methods.pdfec1_metadata_v3.zip...ec1_dd.csv...ec1_flmd.csv...ec1_sample_catalog.csv...ec1_metadata_kitlevel.csv...ec1_metadata_collectionlevel.csv...ec1_data_collectionlevel.csv...ec1_igsn_metadata.csvec1_soil_v3.zipec1_sediment_v3.zipec1_water_v3.zipec1_processingscripts_v3.zipThis data package is on v3 and was originally published May 2023 (v1). Subsequent updates will be published here with new version numbers. Please see the Change History section in ec1_README.pdf for detailed changes.---Acknowledging EXCHANGE: General Support and Data Product UseWe ask that users of EXCHANGE data add the following acknowledgement when publishing data in scholarly articles and data repositories:"This research is based on work supported by COMPASS-FME, a multi-institutional project supported by the U.S. Department of Energy, Office of Science, Biological and Environmental Research as part of the Environmental System Science Program."

54 ENVIRONMENTAL SCIENCES↗

AC and DC Fault Management for Megawatt Electrified Aircraft Electrical Powertrains Task 3: Lifetime and Reliability of Electrical Insulators

This research project was a collaborative investigation between researchers at the RTX Technology Research Center (RTRC) and the University of Texas at Austin and made a significant contribution to enabling electric aircraft. The transport of electric power between the points of generation and use requires power cables. These cables must be smaller, lighter and provide a more predictable life than power cables used in stationary applications. Consequently, this investigation provided heretofore unavailable information supporting the safety and reliability of smaller lighter power cables for electrified aircraft. In addition, the research identified key additional engineering data needed to support quantitative reliability assessments. Important advances included: • Demonstrated that at least one manufacturer can make a novel, smaller, lighter power cable that is free from serious defects. • Developed and published an appropriate analytical construct to describe the life of this novel cable. This is a necessary step for use in aviation where the understanding of remaining life is critical. • Demonstrated thermal-mechanical aging that suggested 1000+ flights before the thermal-mechanical processes produced defects large enough that the defect growth was accelerated electrically. • Showed that electrical aging took place at two rates. The first possibly lasting weeks to months and the second possibly days to weeks. If robust, this provides a good diagnostic for cable replacement. • Demonstrated that the traditional electrical testing of cable materials using manufactured voids can be misleading due to the size of the voids. Emerging laser drilling technology permitted demonstration that the physics of failure in realistically small voids is different from that in the unrealistically large voids used in earlier research, which is very important for high-quality, high-performance, small aircraft cables. Although this project represents a significant contribution to the specifics of cable aging in the aircraft environment, important additional research remains to be completed, including: • Non-uniform thermal cycling by applying the heat from the center conductor to maximize thermal stress next to the core area where the electric gradient is the strongest. This builds on the uniform thermal cycling that has been completed. • The augmentation of the thermal-mechanical failure rate by electrical processes. Better understanding of these time constants strongly affects the ability to predict life. • Termination design: Terminations provide not only electrical reflection potential, but a location for a series arc fault and an area where ozone can diffuse into the center conductor and negatively affect cable insulation. • The abrasion and ozone resistance of the cable jacket. • Pressure cycling as an accelerant of thermal, mechanical, and/or electrical aging. • Possible methods for online PD detection and offline PD localization

model↗

Offshore application of landslide susceptibility mapping using gradient-boosted decision trees: a Gulf of Mexico case study

Abstract Among natural hazards occurring offshore, submarine landslides pose a significant risk to offshore infrastructure installations attached to the seafloor. With the offshore being important for current and future energy production, there is a need to anticipate where future landslide events are likely to occur to support planning and development projects. Using the northern Gulf of Mexico (GoM) as a case study, this paper performs Landslide Susceptibility Mapping (LSM) using a gradient-boosted decision tree (GBDT) model to characterize the spatial patterns of submarine landslide probability over the United States Exclusive Economic Zone (EEZ) where water depths are greater than 120 m. With known spatial extents of historic submarine landslides and a Geographic Information System (GIS) database of known topographical, geomorphological, geological, and geochemical factors, the resulting model was capable of accurately forecasting potential locations of sediment instability. Results of a permutation modelling approach indicated that LSM accuracy is sensitive to the number of unique training locations with model accuracy becoming more stable as the number of training regions was increased. The influence that each input feature had on predicting landslide susceptibility was evaluated using the SHapely Additive exPlanations (SHAP) feature attribution method. Areas of high and very high susceptibility were associated with steep terrain including salt basins and escarpments. This case study serves as an initial assessment of the machine learning (ML) capabilities for producing accurate submarine landslide susceptibility maps given the current state of available natural hazard-related datasets and conveys both successes and limitations.

Dyer, Alec S. (ORCID:0000000219813904)↗

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗

Utilization of the St. Peter Sandstone in the Illinois Basin for CO2 Sequestration

This project is part of a larger project co-funded by the United States Department of Energy (US DOE) under cooperative agreement DE-FE0002068 from 12/08/2009 through 9/31/2014. The study is to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon Sandstone as potential targets for carbon dioxide (CO2) sequestration in the Illinois and Michigan Basins. This report evaluates the potential injectivity of the Ordovician St. Peter Sandstone. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data acquired through funding in this project as well as existing data from two additional, separately funded projects: the US DOE funded Illinois Basin – Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium (MGSC) in Macon County, Illinois, and the Illinois Industrial Carbon Capture and Sequestration (ICCS) Project funded through the American Recovery and Reinvestment Act (ARRA), which received a phase two award from DOE. This study addresses the question of whether or not the St. Peter Sandstone may serve as a suitable target for CO2 sequestration at locations within the Illinois Basin where it lies at greater depths (below the underground source of drinking water (USDW)) than at the IBDP site. The work performed included numerous improvements to the existing St. Peter reservoir model created in 2010. Model size and spatial resolution were increased resulting in a 3 fold increase in the number of model cells. Seismic data was utilized to inform spatial porosity distribution and an extensive core database was used to develop porosity-permeability relationships. The analysis involved a Base Model representative of the St. Peter at “in-situ” conditions, followed by the creation of two hypothetical models at in-situ + 1,000 feet (ft.) (300 m) and in-situ + 2,000 ft. (600 m) depths through systematic depthdependent adjustment of the Base Model parameters. Properties for the depth shifted models were based on porosity versus depth relationship extracted from the core database followed by application of the porosity-permeability relationship. Each of the three resulting models were used as input to dynamic simulations with the single well injection target of 3.2 million tons per annum (MTPA) for 30 years using an appropriate fracture gradient based bottom hole pressure limit for each injection level. Modeling results are presented in terms of well bottomhole pressure (BHP), injection rate profiles, and three-dimensional (3D) saturation and differential pressure volumes at selected simulation times. Results suggest that the target CO2 injection rate of 3.2 MTPA may be achieved in the St. Peter Sandstone at in-situ conditions and at the in-situ +1,000 ft. (300 m) depth using a single injector well. In the latter case the target injection rate is achieved after a ramp up period which is caused by multi-phase flow effects and thus subject to increased modeling uncertainty. Results confirm that the target rate may not be achieved at the in-situ +2,000 ft. (600 m) level even with multiple wells. These new modeling results for the in-situ case are more optimistic than previous modeling results. This difference is attributed to the difference in methods and data used to develop model permeability distributions. Recommendations for further work include restriction of modeling activity to the in-situ +1,000 ft. (300 m) and shallower depth interval, sensitivity and uncertainty analysis, and refinement of porosity and permeability estimates through depth and area selective querying of the available core database. It is also suggested that further modeling efforts include scope for evaluating project performance in terms of metrics directly related to the Environmental Protection Agency (EPA) Class VI permit requirements for the area of review (AoR) definition and post injection site closure monitoring.

Will, Robert↗