Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inverse Optimization”

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 1,423 records · Page 79

How consistent are top-down hydrocarbon emissions based on formaldehyde observations from GOME-2 and OMI?

The vertical columns of formaldehyde (HCHO) retrieved from two satellite instruments, the Global Ozone Monitoring Instrument-2 (GOME-2) on Metop-A and the Ozone Monitoring Instrument (OMI) on Aura, are used to constrain global emissions of HCHO precursors from open fires, vegetation and human activities in the year 2010. To this end, the emissions are varied and optimized using the ad-joint model technique in the IMAGESv2 global CTM (chem-ical transport model) on a monthly basis and at the model res-olution. Given the different local overpass times of GOME- 2 (09:30 LT) and OMI (13:30 LT), the simulated diurnal cy-cle of HCHO columns is investigated and evaluated against ground-based optical measurements at seven sites in Europe, China and Africa. The modeled diurnal cycle exhibits large variability, reflecting competition between photochemistry and emission variations, with noon or early afternoon max-ima at remote locations (oceans) and in regions dominated by anthropogenic emissions, late afternoon or evening max-ima over fire scenes, and midday minima in isoprene-rich re-gions. The agreement between simulated and ground-based columns is generally better in summer (with a clear after-noon maximum at mid-latitude sites) than in winter, and the annually averaged ratio of afternoon to morning columns is slightly higher in the model (1.126) than in the ground-based measurements (1.043).The anthropogenic VOC (volatile organic compound) sources are found to be weakly constrained by the inversions on the global scale, mainly owing to their generally minor contribution to the HCHO columns, except over strongly pol-luted regions, like China. The OMI-based inversion yields total flux estimates over China close to the bottom-up inven-tory (24.6 vs. 25.5 TgVOC yr -1 in the a priori) with, how-ever, pronounced increases in the northeast of China and re-ductions in the south. Lower fluxes are estimated based on GOME-2 HCHO columns (20.6 TgVOC yr -1 ), in particular over the northeast, likely reflecting mismatches between the observed and the modeled diurnal cycle in this region.

54 ENVIRONMENTAL SCIENCES↗

Optimization Of Ocean Color Algorithms: Application To Satellite And In Situ Data Merging

The objective of our program is to develop and validate a procedure for ocean color data merging which is one of the major goals of the SIMBIOS project (McClain et al., 1995). The need for a merging capability is dictated by the fact that since the launch of MODIS on the Terra platform and over the next decade, several global ocean color missions from various space agencies are or will be operational simultaneously. The apparent redundancy in simultaneous ocean color missions can actually be exploited to various benefits. The most obvious benefit is improved coverage (Gregg et al., 1998; Gregg & Woodward, 1998). The patchy and uneven daily coverage from any single sensor can be improved by using a combination of sensors. Beside improved coverage of the global ocean the merging of ocean color data should also result in new, improved, more diverse and better data products with lower uncertainties. Ultimately, ocean color data merging should result in the development of a unified, scientific quality, ocean color time series, from SeaWiFS to NPOESS and beyond. Various approaches can be used for ocean color data merging and several have been tested within the frame of the SIMBIOS program (see e.g. Kwiatkowska & Fargion, 2003, Franz et al., 2003). As part of the SIMBIOS Program, we have developed a merging method for ocean color data. Conversely to other methods our approach does not combine end-products like the subsurface chlorophyll concentration (chl) from different sensors to generate a unified product. Instead, our procedure uses the normalized waterleaving radiances (LwN( )) from single or multiple sensors and uses them in the inversion of a semianalytical ocean color model that allows the retrieval of several ocean color variables simultaneously. Beside ensuring simultaneity and consistency of the retrievals (all products are derived from a single algorithm), this model-based approach has various benefits over techniques that blend end-products (e.g. chlorophyll): 1) it works with single or multiple data sources regardless of their specific bands, 2) it exploits band redundancies and band differences, 3) it accounts for uncertainties in the LwN( ) data and, 4) it provides uncertainty estimates for the retrieved variables.

Maritorena, Stephane↗

Engineering CAM Photosynthetic Machinery into Bioenergy Crops for Biofuels Production in Marginal Environments

CAM provides an excellent opportunity to engineer enhanced photosynthetic performance and water-use efficiency (WUE) into bioenergy crops. Our vision was to understand the enzymatic and regulatory pathways required for CAM in order to engineer CAM photosynthetic machinery into Arabidopsis and Populus. Major advances were made in understanding the molecular “parts list” necessary for engineering CAM through the sequencing the genomes and transcriptomes of several model and agronomically important CAM species including Agave, Ananas, Clusia, Kalanchoe, Mesembryanthemum, and Opuntia. Comparison of diel co-expression network analyses in Agave, Ananas, and Kalanchoe with expression patterns in C 3 photosynthesis species revealed the rescheduling of gene expression associated with stomatal opening and closing and led to the identification of many candidate CAM-related genes for further functional testing. Both RNAi-mediated gene silencing and CRISPR/Cas genome editing for gene knockout were used to investigate over 70 CAM-related enzymes, regulatory kinases/phosphatases, and transporters and 25 transcription factors implicated in playing key regulatory roles in CAM. Knockdown of key C 4 cycle enzymes within the carboxylation and decarboxylation modules of CAM not only perturbed CAM itself, but also altered the expression of core circadian clock components demonstrating that loss of CAM function also perturbs the central circadian clock and alters the regulatory patterns of key guard cell signaling genes that are linked with the characteristic inverse pattern of stomatal opening and closing during CAM. An advanced system for joining D.AN fragments from a universal library was developed that automatically maintains open reading frames (ORFs) and does not require linkers, adaptors, regions of sequence homology, amplification (or mutation for domestication) of fragments to function systems for the construction of gene circuits. Alternatively, Gibson assembly was used to create gene circuits consisting of up to 15 unique gene expression cassettes expressing either the core C 4 cycle of CAM comprising carboxylation or decarboxylation modules or both modules together. The expression of each gene cassette was driven by a drought-stress inducible 5’ regulatory region or promoter (and cognate 3’ terminator region) from A. thaliana so that CAM gene expression would be enhanced under conditions of water-deficit stress. These gene circuits were introduced into the C 3 photosynthesis model species A. thaliana. Expression of the nocturnal carboxylation module resulted in increased biomass and seed yield arising from net CO 2 fixation at night and malate accumulation typical of a CAM plant along with daytime CO 2 fixation typical of a C 3 photosynthesis plant. To create a leaf anatomy suited for the optimal performance of CAM in a C 3 photosynthesis plant, increased tissue succulence was engineered in A. thaliana (and Nicotiana sylvestris), which also improved biomass and seed yield and drought attenuation and improved salinity tolerance. The basic CAM Biodesign principles described in this report can be extended to increase the productivity and WUE of food and bioenergy crops and are expected to have broad applicability and ensure sustainable biofuel feedstock production in the face of the global climate crisis.

59 BASIC BIOLOGICAL SCIENCES↗

Energy efficient integrated photonic systems based on inverse design

The energy footprint of modern information processing and communications systems is immense and utilizes a significant fraction of total global energy usage. Data centers alone consume over 70 billion kilowatt-hours per year. Much of this energy usage is intrinsic to the use of electronic wiring, making optical-based technologies a necessary and promising route to mitigating energy consumption in short- to medium-distance communication links. In this project, we developed and implemented a framework for the design of optical components, based on machine learning, which enables optical components relevant to optical information processing to be realized at their physical performance limits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Core-shell and egg-shell zeolite catalysts for enhanced hydrocarbon processing

Developing structure-performance relationships with the underlying goal of optimizing known zeolite catalysts involves the manipulation of their physicochemical properties. Here, we systematically assessed the impact of mesoscopic gradients in acid site concentration, which has generally received little attention in the design of zeolite catalysts for hydrocarbon upgrading. A series of core–shell MEL-type zeolites were synthesized with catalytically active ZSM-11 cores and passivated silicalite-2 shells of varying thickness. Our findings revealed that ZSM-11@silicalite-2 particles with ultrathin shells (<10 nm) have enhanced mass transport, characteristic of relatively smaller particles, compared to the corresponding ZSM-11 core. Additionally, catalytic testing using the methanol-to-hydrocarbon (MTH) reaction showed that core–shell zeolites exhibit longer lifetimes, higher total turnovers, and an unexpected promotion of the aromatic cycle in the hydrocarbon pool mechanism. Time-resolved acid titration of core and core–shell catalysts confirmed that the siliceous shell introduces a hydrophobic exterior that impacts molecular diffusion. In comparison, prepared MFI core-shells (ZSM-5@silicalite-1) showed similar enhancement in catalyst performance. Moreover, we prepared egg-shell configurations of each zeolite, silicalite-2@ZSM-11 and silicalite-1@ZSM-5, comprised of an inert core and catalytically active shell. This inverse design of the egg-shell created pseudo nanosheets with total turnovers that were markedly higher than their homogeneous counterparts. Collectively, this study demonstrated that mesoscopic gradients in acid concentration via the design of core–shell and egg-shell zeolites significantly improve catalyst performance over conventional analogues for hydrocarbon upgrading.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗

Energy breakdown in capacitive deionization

We explored the energy loss mechanisms in capacitive deionization (CDI). We hypothesize that resistive and parasitic losses are two main sources of energy losses. We measured contribution from each loss mechanism in water desalination with constant current (CC) charge/discharge cycling. Resistive energy loss is expected to dominate in high current charging cases, as it increases approximately linearly with current for fixed charge transfer (resistive power loss scales as square of current and charging time scales as inverse of current). On the other hand, parasitic loss is dominant in low current cases, as the electrodes spend more time at higher voltages. We built a CDI cell with five electrode pairs and standard flow between architecture. We performed a series of experiments with various cycling currents and cut-off voltages (voltage at which current is reversed) and studied these energy losses. To this end, we measured series resistance of the cell (contact resistances, resistance of wires, and resistance of solution in spacers) during charging and discharging from voltage response of a small amplitude AC current signal added to the underlying cycling current. We performed a separate set of experiments to quantify parasitic (or leakage) current of the cell versus cell voltage. We then used these data to estimate parasitic losses under the assumption that leakage current is primarily voltage (and not current) dependent. Our results confirmed that resistive and parasitic losses respectively dominate in the limit of high and low currents. We also measured salt adsorption and report energy-normalized adsorbed salt (ENAS, energy loss per ion removed) and average salt adsorption rate (ASAR). As a result, we show a clear tradeoff between ASAR and ENAS and show that balancing these losses leads to optimal energy efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uncertainty Quantification of GEOS-5 L-band Radiative Transfer Model Parameters Using Bayesian Inference and SMOS Observations

Uncertainties in L-band (1.4 GHz) radiative transfer modeling (RTM) affect the simulation of brightness temperatures (Tb) over land and the inversion of satellite-observed Tb into soil moisture retrievals. In particular, accurate estimates of the microwave soil roughness, vegetation opacity and scattering albedo for large-scale applications are difficult to obtain from field studies and often lack an uncertainty estimate. Here, a Markov Chain Monte Carlo (MCMC) simulation method is used to determine satellite-scale estimates of RTM parameters and their posterior uncertainty by minimizing the misfit between long-term averages and standard deviations of simulated and observed Tb at a range of incidence angles, at horizontal and vertical polarization, and for morning and evening overpasses. Tb simulations are generated with the Goddard Earth Observing System (GEOS-5) and confronted with Tb observations from the Soil Moisture Ocean Salinity (SMOS) mission. The MCMC algorithm suggests that the relative uncertainty of the RTM parameter estimates is typically less than 25 of the maximum a posteriori density (MAP) parameter value. Furthermore, the actual root-mean-square-differences in long-term Tb averages and standard deviations are found consistent with the respective estimated total simulation and observation error standard deviations of m3.1K and s2.4K. It is also shown that the MAP parameter values estimated through MCMC simulation are in close agreement with those obtained with Particle Swarm Optimization (PSO).

MCMC↗

Demonstration of Decentralized Physics-Driven Learning

In typical artificial neural networks, neurons adjust according to global calculations of a central processor, but in the brain, neurons and synapses self-adjust based on local information. Contrastive learning algorithms have recently been proposed to train physical systems, such as fluidic, mechanical, or electrical networks, to perform machine-learning tasks from local evolution rules. However, to date, such systems have only been implemented in silico due to the engineering challenge of creating elements that autonomously evolve based on their own response to two sets of global boundary conditions. Here, we introduce and implement a physics-driven contrastive learning scheme for a network of variable resistors, using circuitry to locally compare the response of two identical networks subjected to the two different sets of boundary conditions. Using this method, our system effectively trains itself, optimizing its resistance values without the use of a central processor or external information storage. Once the system is trained for a specified allostery, regression, or classification task, the task is subsequently performed rapidly and automatically by the physical imperative to minimize power dissipation in response to the given voltage inputs. Here, we demonstrate that, unlike typical computers, such learning systems are robust to extreme damage (and thus manufacturing defects) due to their decentralized learning. Our twin-network approach is therefore readily scalable to extremely large or nonlinear networks, where its distributed nature will be an enormous advantage; a laboratory network of only 500 edges will already outpace its in silico counterpart.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A cloud computing approach to massively-parallel distributed stochastic optimization with application to geomechanical reservoir characterization

Many search strategies used to optimize parameters during hydrologic model calibration are inherently limited in that simulations must be performed sequentially in order for the algorithm to converge. As a result, the overall time required to reach a solution increases linearly with the number of steps necessary. In contrast, sampling strategies that attempt to approximate the posterior distribution of the parameters do not require the model simulations to be run in sequence and therefore invite the use of a massively parallel approach across heterogeneous computing environments. We apply a cloud computing methodology to this problem where a head node running on Amazon Web Services continuously produces candidate models that are distributed to arbitrary compute environments to run the simulation, in our case Clemson’s Palmetto cluster computer. Thus the time required to obtain a sufficient number of samples to produce an adequate approximation of the posterior distribution of the parameters is limited by the number of available compute sites. Monte Carlo sampling of the posterior space would be duplicative and inefficient, thus we combine Latin hypercube methods with genetic algorithms to ensure that no large regions of the parameter space are left unexplored, and that high probability regions are investigated thoroughly. We independently assess the fit to different types of data sets without the need to construct a single objective function. Sets of Pareto optimal models are particularly valuable when the quality and information content of data is poorly known or conceptual model errors exist in our simulation. We illustrate our approach for a problem where in situ strain data are used to calibrate the parameters of an oil reservoir in the Avant oil field in Oklahoma. We use well injection and shut-in to stress a deep (500m) confined formation, and observe the geomechanical response using borehole strainmeters installed in the shallow subsurface (30m). These data were then used to calibrate the poroelastic properties of the reservoir and caprock along with the geometry of an elliptic channel associated with the formation. Using our approach we are able to identify issues with model error as well as substantial tradeoffs that occur within the parameters.

inversion, stochastic, geomechanics, poroelastic, ↗

Characterization of Climate Change and Variability with GPS

We compared zonal mean specific humidity derived from the 21 June-4 July 1995 Global Positioning System (GPS)/MET occultation observations with that derived from the European Center for Medium-Range Weather Forecasts (ECMWF) global analyses. The GPS/MET results indicate a drier troposphere, especially near the subtropical tradewind inversion. A small, moist bias in the GPS/MET upper northern-hemisphere troposphere compared to ECMWF may be due to a small radiosonde temperature bias. A diagram shows the difference (g/kg) between the GPS/MET zonal mean specific humidity and that for June-August derived from 1963-1973 radiosondes. Although the observing period is short, GPS and ECMWF results both indicate a significantly wetter boundary layer at most latitudes consistent with decadal trends observed in radiosonde data. GPS/MET results exhibit higher tropical convective available potential energy (CAPE), suggesting a more vigorous tropical Hadley circulation. Drier, free troposphere air in the descending branches of the Hadley circulation is due in part to a moist radiosonde bias but may also reflect some negative moisture feedback. Using 1992-1997 ground GPS observations and recent advancements in GPS technology, we removed an apparent altimetric drift (-1.2 +/- 0.4 mm/yr) due to columnar water vapor from the Topography (Ocean) Experiment (TOPEX) microwave radiometer, which brought the TOPEX mean sea level change estimates into better agreement with historical tide gauge records, suggesting global mean sea level is rising at a rate of 1.5-2.0 mm/yr. We can also discern a statistically significant increase of 0.2 +/- 0.1 kg/square m/yr in mean columnar water vapor over the ocean from 1992-1997. Optimal fingerprinting can be used for the detection and attribution of tropospheric warming due to an anthropogenic greenhouse. Optimal fingerprinting distinguishes between different types of signals according to their spatial and temporal patterns, while minimizing the influence of natural climate variability. S. Leroy concludes that the signal-to-noise ratio of global warming detection increases by unity approximately every 10 years if a single oceanic region is chosen. Less time for detection is likely when many global regions are considered simultaneously. GPS occultation constellations allow the possibility of detecting small changes in upper air temperature with inconsequential calibration errors, making occultation an ideal data type for global warming detection studies. Our initial study of a 22-GHz satellite-satellite occultation system predicts upper troposphere moisture sensitivities of 3-5 ppmv and 1-2 percent in the middle and lower troposphere. Additional information contained in original.

Kursinski, R.↗

NASA Tech Briefs, November 2010

Topics covered include: Portable Handheld Optical Window Inspection Device; Salience Assignment for Multiple-Instance Data and Its Application to Crop Yield Prediction; Speech Acquisition and Automatic Speech Recognition for Integrated Spacesuit Audio Systems ; Predicting Long-Range Traversability from Short-Range Stereo-Derived Geometry; Browser-Based Application for Telemetry Monitoring of Robotic Assets; Miniature Low-Noise G-Band I-Q Receiver; Methods of Using a Magnetic Field Response Sensor Within Closed, Electrically Conductive Containers; Differential Resonant Ring YIG Tuned Oscillator; Microfabricated Segmented-Involute-Foil Regenerator for Stirling Engines; Reducing Seal Adhesion in Low Impact Docking Systems; Optimal Flow Control Design; Corrosion-Resistant Container for Molten-Material Processing; Reusable Hot-Wire Cable Cutter; Deployment of a Curved Truss; High-Volume Airborne Fluids Handling Technologies to Fight Wildfires; Modeling of Alkane Oxidation Using Constituents and Species; Fabrication of Lanthanum Telluride 14-1-11 Zintl High-Temperature Thermoelectric Couple; A Computer Model for Analyzing Volatile Removal Assembly; Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature; Optical Sidebands Multiplier; Single Spatial-Mode Room-Temperature-Operated 3.0 to 3.4 micrometer Diode Lasers; Self-Nulling Beam Combiner Using No External Phase Inverter; Portable Dew Point Mass Spectrometry System for Real-Time Gas and Moisture Analysis; Maximum Likelihood Time-of-Arrival Estimation of Optical Pulses via Photon-Counting Photodetectors; Handheld White Light Interferometer for Measuring Defect Depth in Windows; Decomposition Algorithm for Global Reachability on a Time-Varying Graph; Autonomous GN and C for Spacecraft Exploration of Comets and Asteroids; Efficient Web Services Policy Combination; Using CTX Image Features to Predict HiRISE-Equivalent Rock Density; Isolation of the Paenibacillus phoenicis, a Spore-Forming Bacterium; Monolithically Integrated, Mechanically Resilient Carbon-Based Probes for Scanning Probe Microscopy; Cell Radiation Experiment System; Process to Produce Iron Nanoparticle Lunar Dust Simulant Composite; Inversion Method for Early Detection of ARES-1 Case Breach Failure; Use of ILTV Control Laws for LaNCETS Flight Research;and Evaluating Descent and Ascent Trajectories Near Non-Spherical Bodies.

Source record↗

Electron-Spin Filters Based on the Rashba Effect

Semiconductor electron-spin filters of a proposed type would be based on the Rashba effect, which is described briefly below. Electron-spin filters more precisely, sources of spin-polarized electron currents have been sought for research on, and development of, the emerging technological discipline of spintronics (spin-based electronics). There have been a number of successful demonstrations of injection of spin-polarized electrons from diluted magnetic semiconductors and from ferromagnetic metals into nonmagnetic semiconductors. In contrast, a device according to the proposal would be made from nonmagnetic semiconductor materials and would function without an applied magnetic field. The Rashba effect, named after one of its discoverers, is an energy splitting, of what would otherwise be degenerate quantum states, caused by a spin-orbit interaction in conjunction with a structural-inversion asymmetry in the presence of interfacial electric fields in a semiconductor heterostructure. The magnitude of the energy split is proportional to the electron wave number. The present proposal evolved from recent theoretical studies that suggested the possibility of devices in which electron energy states would be split by the Rashba effect and spin-polarized currents would be extracted by resonant quantum-mechanical tunneling. Accordingly, a device according to the proposal would be denoted an asymmetric resonant interband tunneling diode [a-RITD]. An a-RITD could be implemented in a variety of forms, the form favored in the proposal being a double-barrier heterostructure containing an asymmetric quantum well. It is envisioned that a-RITDs would be designed and fabricated in the InAs/GaSb/AlSb material system for several reasons: Heterostructures in this material system are strong candidates for pronounced Rashba spin splitting because InAs and GaSb exhibit large spin-orbit interactions and because both InAs and GaSb would be available for the construction of highly asymmetric quantum wells. This mate-rial system affords a variety of energy-band alignments that can be exploited to obtain resonant tunneling and other desired effects. The no-common-atom InAs/GaSb and InAs/AlSb interfaces would present opportunities for engineering interface potentials for optimizing Rashba spin splitting.

Ting, David Z.-Y.↗

Numerical and experimental studies on effects of moisture content on combustion characteristics of simulated municipal solid wastes in a fixed bed

Highlights: • The effects of moisture content on the burning process of MSW are investigated. • A two-dimensional mathematical model was built to simulate the combustion process. • Temperature distributions, process rates, gas species were measured and simulated. • The The conversion ratio of C/CO and N/NO in MSW are inverse to moisture content. - Abstract: In order to reveal the features of the combustion process in the porous bed of a waste incinerator, a two-dimensional unsteady state model and experimental study were employed to investigate the combustion process in a fixed bed of municipal solid waste (MSW) on the combustion process in a fixed bed reactor. Conservation equations of the waste bed were implemented to describe the incineration process. The gas phase turbulence was modeled using the k–ε turbulent model and the particle phase was modeled using the kinetic theory of granular flow. The rate of moisture evaporation, devolatilization rate, and char burnout was calculated according to the waste property characters. The simulation results were then compared with experimental data for different moisture content of MSW, which shows that the incineration process of waste in the fixed bed is reasonably simulated. The simulation results of solid temperature, gas species and process rate in the bed are accordant with experimental data. Due to the high moisture content of fuel, moisture evaporation consumes a vast amount of heat, and the evaporation takes up most of the combustion time (about 2/3 of the whole combustion process). The whole bed combustion process reduces greatly as MSW moisture content increases. The experimental and simulation results provide direction for design and optimization of the fixed bed of MSW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

Effects of Vegetation and of Heat and Vapor Fluxes from Soil on Snowpack Evolution and Radiobrightness

The radiobrightness of a snowpack is strongly linked to the snow moisture content profile, to the point that the only operational inversion algorithms require dry snow. Forward dynamic models do not include the effects of freezing and thawing of the soil beneath the snowpack and the effect of vegetation within the snow or above the snow. To get a more realistic description of the evolution of the snowpack, we reported an addition to the Snow-Soil-Vegetation-Atmosphere- Transfer (SSVAT) model, wherein we coupled soil processes of the Land Surface Process (LSP) model with the snow model SNTHERM. In the near future we will be adding a radiobrightness prediction based on the modeled moisture, temperature and snow grain size profiles. The initial investigations with this SSVAT for a late winter and early spring snow pack indicate that soil processes warm the snowpack and the soil. Vapor diffusion needs to be considered whenever the ground is thawed. In the early spring, heat flow from the ground into a snow and a strong temperature gradient across the snow lead to thermal convection. The buried vegetation can be ignored for a late winter snow pack. The warmer surface snow temperature will affect radiobrightness since it is most sensitive to snow surface characteristics. Comparison to data shows that SSVAT provides a more realistic representation of the temperature and moisture profiles in the snowpack and its underlying soil than SNTHERM. The radiobrightness module will be optimized for the prediction of brightness when the snow is moist. The liquid water content of snow causes considerable absorption compared to dry snow, and so longer wavelengths are likely to be most revealing as to the state of a moist snowpack. For volumetric moisture contents below about 7% (the pendular regime), the water forms rings around the contact points between snow grains. Electrostatic modeling of these pendular rings shows that the absorption of these rings is significantly higher than a sphere of the same volume. The first implementation of the radiobrightness module will therefore be a simple radiative transfer model without scattering.

soil↗