Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “target variability”

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 955 records · Page 53

Towards Ultra-high-resolution E3SM Land Modeling on Exascale Computers

Here we present an ultra-high-resolution E3SM land model (uELM) for high-fidelity land simulations targeting new Exascale computers. After considering modeling infrastructure compatibility and ELM software features, we designed a parallel model for the uELM development targeting hybrid architectures of new US Exascale computers. We also described a function unit test framework to expedite the piece-wise code porting (with compiler directives), verification, and global variable management. Furthermore, in this study, we report an early uELM model development using OpenACC within a function unit test framework on a pre-Exascale computer, demonstrate the performance of a uLEM submodel with a 3.0-time speedup, and summarize the code porting experience regarding global variable handling, deepcopy, memory reduction, and parallel loop reconstruction.

97 MATHEMATICS AND COMPUTING↗

An Improved Stochastic Unit Commitment Formulation to Accommodate Wind Uncertainty

The United States targets to supply 20% of its electricity generation using wind energy by 2030. The expansion of renewable resources, especially weather-based resources such as wind, creates more uncertainty and variability in the operation of the power grid. New methods and approaches in electricity market operations are needed to efficiently manage the continuing increase in variability and uncertainty caused by expanding intermittent wind. This paper proposes an improved stochastic programming approach for incorporating wind uncertainty into energy markets. The proposed formulation improves the two-stage stochastic unit commitment problem by introducing a dynamic decision making approach similar to a multi-stage formulation in the presence of wind power scenarios which are not well represented by a scenario tree. The numerical results present up to 1%-2% decrease in operational costs compared to the two-stage stochastic unit commitment formulation.

Electricity markets↗

Experimental evaluations of the accuracy of 3D and 4D planning in robotic tracking stereotactic body radiotherapy for lung cancers

Purpose: Due to the complexity of 4D target tracking radiotherapy, the accuracy of this treatment strategy should be experimentally validated against established standard 3D technique. This work compared the accuracy of 3D and 4D dose calculations in respiration tracking stereotactic body radiotherapy (SBRT). Methods: Using the 4D planning module of the CyberKnife treatment planning system, treatment plans for a moving target and a static off-target cord structure were created on different four-dimensional computed tomography (4D-CT) datasets of a thorax phantom moving in different ranges. The 4D planning system used B-splines deformable image registrations (DIR) to accumulate dose distributions calculated on different breathing geometries, each corresponding to a static 3D-CT image of the 4D-CT dataset, onto a reference image to compose a 4D dose distribution. For each motion, 4D optimization was performed to generate a 4D treatment plan of the moving target. For comparison with standard 3D planning, each 4D plan was copied to the reference end-exhale images and a standard 3D dose calculation was followed. Treatment plans of the off-target structure were first obtained by standard 3D optimization on the end-exhale images. Subsequently, they were applied to recalculate the 4D dose distributions using DIRs. All dose distributions that were initially obtained using the ray-tracing algorithm with equivalent path-length heterogeneity correction (3D{sub EPL} and 4D{sub EPL}) were recalculated by a Monte Carlo algorithm (3D{sub MC} and 4D{sub MC}) to further investigate the effects of dose calculation algorithms. The calculated 3D{sub EPL}, 3D{sub MC}, 4D{sub EPL}, and 4D{sub MC} dose distributions were compared to measurements by Gafchromic EBT2 films in the axial and coronal planes of the moving target object, and the coronal plane for the static off-target object based on the {gamma} metric at 5%/3mm criteria ({gamma}{sub 5%/3mm}). Treatment plans were considered acceptable if the percentage of pixels passing {gamma}{sub 5%/3mm} (P{sub {gamma}<1}) {>=} 90%. Results: The averaged P{sub {gamma}<1} values of the 3D{sub EPL}, 3D{sub MC}, 4D{sub EPL}, and 4D{sub MC} dose calculation methods for the moving target plans are 95%, 95%, 94%, and 95% for reproducible motion, and 95%, 96%, 94%, and 93% for nonreproducible motion during actual treatment delivery. The overall measured target dose distributions are in better agreement with the 3D{sub MC} dose distributions than the 4D{sub MC} dose distributions. Conversely, measured dose distributions agree much better with the 4D{sub EPL/MC} than the 3D{sub EPL/MC} dose distributions in the static off-target structure, resulting in higher P{sub {gamma}<1} values with 4D{sub EPL/MC} (91%) vs 3D{sub EPL} (24%) and 3D{sub MC} (25%). Systematic changes of target motion reduced the averaged P{sub {gamma}<1} to 47% and 53% for 4D{sub EPL} and 4D{sub MC} dose calculations, and 22% for 3D{sub EPL/MC} dose calculations in the off-target films. Conclusions: In robotic tracking SBRT, 4D treatment planning was found to yield better prediction of the dose distributions in the off-target structure, but not necessarily in the moving target, compared to standard 3D treatment planning, for reproducible and nonreproducible target motion. It is important to ensure on a patient-by-patient basis that the cumulative uncertainty associated with the 4D-CT artifacts, deformable image registration, and motion variability is significantly smaller than the cumulative uncertainty occurred in standard 3D planning in order to make 4D planning a justified option.

60 APPLIED LIFE SCIENCES↗

TH-A-9A-04: Incorporating Liver Functionality in Radiation Therapy Treatment Planning

Purpose: Liver SBRT patients have both variable pretreatment liver function (e.g., due to degree of cirrhosis and/or prior treatments) and sensitivity to radiation, leading to high variability in potential liver toxicity with similar doses. This work aims to explicitly incorporate liver perfusion into treatment planning to redistribute dose to preserve well-functioning areas without compromising target coverage. Methods: Voxel-based liver perfusion, a measure of functionality, was computed from dynamic contrast-enhanced MRI. Two optimization models with different cost functions subject to the same dose constraints (e.g., minimum target EUD and maximum critical structure EUDs) were compared. The cost functions minimized were EUD (standard model) and functionality-weighted EUD (functional model) to the liver. The resulting treatment plans delivering the same target EUD were compared with respect to their DVHs, their dose wash difference, the average dose delivered to voxels of a particular perfusion level, and change in number of high-/low-functioning voxels receiving a particular dose. Two-dimensional synthetic and three-dimensional clinical examples were studied. Results: The DVHs of all structures of plans from each model were comparable. In contrast, in plans obtained with the functional model, the average dose delivered to high-/low-functioning voxels was lower/higher than in plans obtained with its standard counterpart. The number of high-/low-functioning voxels receiving high/low dose was lower in the plans that considered perfusion in the cost function than in the plans that did not. Redistribution of dose can be observed in the dose wash differences. Conclusion: Liver perfusion can be used during treatment planning potentially to minimize the risk of toxicity during liver SBRT, resulting in better global liver function. The functional model redistributes dose in the standard model from higher to lower functioning voxels, while achieving the same target EUD and satisfying dose limits to critical structures. This project is funded by MCubed and grant R01-CA132834.

60 APPLIED LIFE SCIENCES↗

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗

Evaluation of the FIR Example using Xilinx Vivado High-Level Synthesis Compiler

Compared to central processing units (CPUs) and graphics processing units (GPUs), field programmable gate arrays (FPGAs) have major advantages in reconfigurability and performance achieved per watt. This development flow has been augmented with high-level synthesis (HLS) flow that can convert programs written in a high-level programming language to Hardware Description Language (HDL). Using high-level programming languages such as C, C++, and OpenCL for FPGA-based development could allow software developers, who have little FPGA knowledge, to take advantage of the FPGA-based application acceleration. This improves developer productivity and makes the FPGA-based acceleration accessible to hardware and software developers. Xilinx Vivado HLS compiler is a high-level synthesis tool that enables C, C++ and System C specification to be directly targeted into Xilinx FPGAs without the need to create RTL manually. The white paper [1] published recently by Xilinx uses a finite impulse response (FIR) example to demonstrate the variable-precision features in the Vivado HLS compiler and the resource and power benefits of converting floating point to fixed point for a design. To get a better understanding of variable-precision features in terms of resource usage and performance, this report presents the experimental results of evaluating the FIR example using Vivado HLS 2017.1 and a Kintex Ultrascale FPGA. In addition, we evaluated the half-precision floating-point data type against the double-precision and single-precision data type and present the detailed results.

97 MATHEMATICS AND COMPUTING↗

Computational Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. The computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

36 MATERIALS SCIENCE↗

Cloud types and the tropical Earth radiation budget, revised

Nimbus-7 cloud and Earth radiation budget data are compared in a study of the effects of clouds on the tropical radiation budget. The data consist of daily averages over fixed 500 sq km target areas, and the months of July 1979 and January 1980 were chosen to show the effect of seasonal changes. Six climate regions, consisting of 14 to 24 target areas each, were picked for intensive analysis because they exemplified the range in the tropical cloud/net radiation interactions. The normal analysis was to consider net radiation as the independent variable and examine how cloud cover, cloud type, albedo and emitted radiation varied with the net radiation. Two recurring themes keep repeating on a local, regional, and zonal basis: the net radiation is strongly influenced by the average cloud type and amount present, but most net radiation values could be produced by several combinations of cloud types and amount. The regions of highest net radiation (greater than 125 W/sq m) tend to have medium to heavy cloud cover. In these cases, thin medium altitude clouds predominate. Their cloud tops are normally too warm to be classified as cirrus by the Nimbus cloud algorithm. A common feature in the tropical oceans are large regions where the total regional cloud cover varies from 20 to 90 percent, but with little regional difference in the net radiation. The monsoon and rain areas are high net radiation regions.

Dhuria, Harbans L.↗

Personalized Risk Assesment for Adverse Drug Reactions and Treatment Failures

As NASA and other space agencies prepare for future long-term space missions beyond the LEO, the cumulative impact of risk factors encountered in space increases substantially rising concerns about astronauts health. Application of on-board medications to mitigate clinical symptoms associated with certain medical conditions and illnesses is the first line of response to ensure sustainable health and performance of crew. Unfortunately, very limited research has been conducted to determine efficacy of the earth-based pharmaceuticals in a microgravity environment. In some instances, orally administered medications taken during flight were reported to be less effective than expected. Evaluation of series of experiments involving astronauts from shuttle flights shows notable individual variability to several pharmaceuticals during flight. These data provide reasonable assumption of perturbation in CYP450 enzymes during spaceflight, which contribute to the hepatic metabolism of the majority of drugs and therefore may have significant effects on therapeutic efficacy and increase treatment-related toxicity. The genes encoding the CYP450 enzymes are highly variable in humans. Inheritable variations of CYP450 hepatic metabolizer enzymes and transport proteins play a crucial role in the inter-individual variability of drug efficiency and risks of adverse drug reactions. Additionally, there are some reports that document changes in the levels of production of drug-metabolizing enzymes in microgravity. Therefore, in order to provide a safe and effective pharmaceutical treatment in space, medications selection should be based not only on the specific efficacy of medications but also on the individual drug sensitivity and flight-induced changes in metabolism of astronauts chosen for a particular mission. To our knowledge, there was no pre-flight drug sensitivity testing on a genetic level for any of the previous manned NASA space missions. Therefore, technologies capable of predicting and managing medication efficacy, side effects, and toxicity of drugs based on individual genetic variability of crew members are increasingly needed. In this report, we present results of testing the market available Personalized Prescribing System (PPS), a comprehensive, non-invasive solution for safer, targeted medication management for every crew member. Statistical accuracy and simplicity of non-invasive sample analysis demonstrate the feasibility of drug sensitivity assessment and record-keeping tool for flight surgeons and astronauts in applying the recommended medications for situations arising in flight. The information on individual drug sensitivity will translate into personalized risk assessment for adverse drug reactions and treatment failures for each drug from the medication kit as well as predefined outcome. This will address the HHC’s raised “Concern of Clinically Relevant Unpredicted Effects of Medication” as recently updated.

medications↗

Gully Formation at the Haughton Impact Structure (Arctic Canada) Through the Melting of Snow and Ground Ice, with Implications for Gully Formation on Mars

The formation of gullies on Mars has been the topic of active debate and scientific study since their first discovery by Malin and Edgett in 2000. Several mechanisms have been proposed to account for gully formation on Mars, from “dry” mass movement processes, release of water or brine from subsurface aquifers, and the melting of near-surface ground ice or snowpacks. In their global documentation of martian gullies, report that gullies are confined to ~27–83°S and ~28–72°N latitudes and span all longitudes. Gullies on Mars have been documented on impact crater walls and central uplifts, isolated massifs, and on canyon walls, with crater walls being the most common situation. In order to better understand gully formation on Mars, we have been conducting field studies in the Canadian High Arctic over the past several summers, most recently in summer 2018 and 2019 under the auspices of the Canadian Space Agency-funded Icy Mars Analogue Program. It is notable that the majority of previous studies in the Arctic and Antarctica, including our recent work on Devon Island, have focused on gullies formed on slopes generated by regular endogenic geological processes and in “regular” bedrock. How-ever, as noted above, meteorite impact craters are the most dominant setting for gullies on Mars. Impact craters provide an environment with diverse lithologies – including impact-generated and impact-modified rocks – and slope angle, and thus greatly variable hill slope processes could occur within a localized area. Here, we investigate the formation of gullies within the Haughton impact structure and compare them to gullies formed in unimpacted target rock in the nearby Thomas Lee Inlet

Osinski, G. R.↗

Reflection and Timing Study of the Transient Black Hole X-Ray Binary MAXI J1803-298 with NuSTAR

The transient black hole X-ray binary MAXI J1803−298 was discovered on 2021 May 1, as it went into outburst from a quiescent state. As the source rose in flux it showed periodic absorption dips and fit the timing and spectral characteristics of a hard-state accreting black hole. We report on the results of a Target-of-Opportunity observation with NuSTAR obtained near the peak outburst flux beginning on 2021 May 13, after the source had transitioned into an intermediate state. MAXI J1803−298 is variable across the observation, which we investigate by extracting spectral and timing products separately for different levels of flux throughout the observation. Our timing analysis reveals two distinct potential quasiperiodic oscillations (QPOs) which are not harmonically related at 5.4 ± 0.2 Hz and 9.4 ± 0.3 Hz, present only during periods of lower flux. With clear relativistic reflection signatures detected in the source spectrum, we applied several different reflection models to the spectra of MAXI J1803−298. Here we report our results, utilizing high-density reflection models to constrain the disk geometry, and assess changes in the spectrum dependent on the source flux. With a standard broken power-law emissivity, we find a near-maximal spin for the black hole, and we are able to constrain the inclination of the accretion disk at 75° ± 2°, which is expected for a source that has shown periodic absorption dips. We also significantly detect a narrow absorption feature at 6.91 ± 0.06 keV with an equivalent width between 4 and 9 eV, which we interpret as the signature of a disk wind.

Benjamin M. Coughenour↗

Combustion mode switching with a turbocharged/supercharged engine

A method for switching between low- and high-dilution combustion modes in an internal combustion engine having an intake passage with an exhaust-driven turbocharger, a crankshaft-driven positive displacement supercharger downstream of the turbocharger and having variable boost controllable with a supercharger bypass valve, and a throttle valve downstream of the supercharger. The current combustion mode and mass air flow are determined. A switch to the target combustion mode is commanded when an operating condition falls within a range of predetermined operating conditions. A target mass air flow to achieve a target air-fuel ratio corresponding to the current operating condition and the target combustion mode is determined. The degree of opening of the supercharger bypass valve and the throttle valve are controlled to achieve the target mass air flow. The amount of residual exhaust gas is manipulated.

33 ADVANCED PROPULSION SYSTEMS↗

Experimental verification of a method to create a variable energy neutron beam from a monoenergetic, isotropic source using neutron elastic scatter and time of flight

An experiment was performed to determine the neutron energy of near-monoergetic deuterium–deu-terium (D–D) neutrons that elastically scatter in a hydrogenous target. The experiment used two liquid scintillators to perform time of flight (TOF) measurements to determine neutron energy, with the start detector also serving as the scatter target. The stop detector was placed 1.0 m away and at scatter angles of π/6, π/4, and π/3 rad, and 1.5 m at a scatter angle of π/4 rad. When discrete 1 ns increments were implemented, the TOF peaks had estimated errors between 21.2 and 3.6% relative to their expected locations. Full widths at half-maximum (FWHM) ranged between 9.6 and 20.9 ns, or approximately 0.56–0.66 MeV. Monte Carlo simulations were also conducted that approximated the experimental setup and had both D–D and deuterium–tritium (DT) neutrons. The simulated results had errors between -17.2 and 0.0% relative to their expected TOF peaks when 1 ns increments were applied. The largest D–D and D–T FWHMs were 26.7 and 13.7 ns, or approximately 0.85 and 4.98 MeV, respectively. These values, however, can be reduced through manipulation of the dimensions of the system components. Lastly, the results encourage further study of the neutron elastic scatter TOF system with particular interest in application to active neutron interrogation to search for conventional explosives.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Seasonal Impact of Regional Outdoor Biomass Burning on Air Pollution in Three Indian Cities: Delhi, Bengaluru, and Pune

Air pollution in many of India's cities exceeds national and international standards, and effective pollution control strategies require knowledge of the sources that contribute to air pollution and their spatiotemporal variability. In this study, we examine the influence of a single pollution source, outdoor biomass burning, on particulate matter (PM) concentrations, surface visibility, and aerosol optical depth (AOD) from 2007 to 2013 in three of the most populous Indian cities. We define the upwind regions, or "airsheds," for the cities by using atmospheric back trajectories from the HYSPLIT model. Using satellite fire radiative power (FRP) observations as a measure of fire activity, we target pre-monsoon and post-monsoon fires upwind of the Delhi National Capital Region and pre-monsoon fires surrounding Bengaluru and Pune. We find varying contributions of outdoor fires to different air quality metrics. For the post-monsoon burning season, we find that a subset of local meteorological variables (air temperature, humidity, sea level pressure, wind speed and direction) and FRP as the only pollution source explained 39% of variance in Delhi station PM(sub 10) anomalies, 77% in visibility, and 30% in satellite AOD; additionally, per unit increase in FRP within the daily airshed (1000 MW), PM(sub 10) increases by 16.34 micrograms per cubic meter, visibility decreases by 0.097 km, and satellite AOD increases by 0.07. In contrast, for the pre-monsoon burning season, we find less significant contributions from FRP to air quality in all three cities. Further, we attribute 99% of FRP from post-monsoon outdoor fires within Delhi's average airshed to agricultural burning. Our work suggests that although outdoor fires are not the dominant air pollution source in India throughout the year, post-monsoon fires contribute substantially to regional air pollution and high levels of population exposure around Delhi. During 3-day blocks of extreme PM(sub 2.5) in the 2013 post-monsoon burning season, which coincided with statistically significant high fire activity, concentrations in Delhi averaged 304 micrograms per cubic meter, or more than 1000% above the 24-h PM(sub 2.5) guideline (25 micrograms per cubic meter) of the World Health Organization. These results suggest that providing viable alternatives to agricultural residue burning could help improve post-monsoon air quality for a growing population of 63 million (39% in urban areas) within Delhi's airshed.

Air quality; Outdoor fires; Crop residue burning; ↗

Intraspecific variability in root traits and edaphic conditions influence soil microbiomes across 12 switchgrass cultivars

Microbial communities help plants access nutrients and tolerate stress. Some microbiomes are specific to plant genotypes and, therefore, may contribute to intraspecific differences in plant growth and be a promising target for plant breeding. Switchgrass (Panicum virgatum L.) is a potential bioenergy crop with broad variation in yields and environmental responses; recent studies suggest that associations with distinct microbiomes may contribute to variation in cultivar yields. We used a common garden experiment to investigate variation in 12 mature switchgrass cultivar soil microbiomes and, further, to examine how root traits and soil conditions influence microbiome structure. Here, we found that average root diameter varied up to 33% among cultivars and that they associated with distinct soil microbiomes. Cultivar had a larger effect on the soil bacterial than fungal community, but both were strongly influenced by soil properties. Root traits had a weaker effect on microbiome structure, but root length contributed to variation in the fungal community. Unlike the soil communities, the root bacterial communities did not group by cultivar, based on a subset of samples. Microbial biomass carbon and nitrogen and the abundance of several dominant bacterial phyla varied between ecotypes, but overall the differences in soil microbiomes were greater among cultivars than between ecotypes. Our findings show that there is not one soil microbiome that applies to all switchgrass cultivars, or even to each ecotype. These subtle but significant differences in root traits, microbial biomass, and the abundance of certain soil bacteria could explain differences in cultivar yields and environmental responses.

59 BASIC BIOLOGICAL SCIENCES↗

NEON AOP Imaging Spectroscopy Survey of Upper East River Colorado Watersheds: Raw-Space Radiance and Observational Variable Dataset

Lawrence Berkeley National Laboratory (LBNL) contracted the National Ecological Observatory Network Airborne Observation Platform (NEON AOP) to observe watersheds of interest near Crested Butte, CO with remotely sensed data including imaging spectroscopy. The flight box design encompassed the watersheds, and flight design targeted a minimum of 30% overlap between adjacent flight lines. The total area surveyed was 334 km2 across 72 lines. Flight lines were oriented north-south for consistency in solar illumination geometry. Seven days were required to cover the flight area. Some areas were flown additional times to improve data quality or temporal coincidence between airborne acquisitions and field sampling. In total, 102 lines were collected and can be downloaded here: https://portal.nersc.gov/wfsfa/doi-10-15485-1617204/. The instrument used for the acquisition was an Airborne Visible-Infrared Imaging Spectrometer - Next Generation (AVIRIS-NG), which records radiance within 426 contiguous spectral bands in the wavelength range from 380 to 2510 nm. The instrument was flown at a target aboveground altitude of 1000 m for a processed data spatial resolution of 1 m. These data are non-orthorectified, radiance data that require further processing prior to use as a surface reflectance data set.

2018 NEON and 2025 CHESS Campaigns↗

Survey on stochastic distribution systems: A full probability density function control theory with potential applications

Complex systems seen either in general engineering practice or economics are subjected to ever increased uncertainties that are mostly represented as random variables or parameters, and the characteristics of random variables are represented by their probability density functions (PDFs). Controlling their PDFs means to shape their stochastic distributions and in general it would provide a full treatment for system analysis and operational control and optimization. This leads to the development of stochastic distribution control (SDC) systems theory in the past decades, where the original aim of the controller design is to realize a shape control of the distributions of certain random variables in their PDFs sense for some engineering processes. Indeed, once the PDFs of these random variables or parameters are used to describe their distribution characters, the control task is to obtain control signals so that the output PDFs of stochastic systems are made to follow their target PDFs. The subject of SDC was initially originated for non-Gaussian stochastic control systems design but has found a wide spectrum of applications in general systems in terms of data-driven modeling, analysis, signal processing (filtering), data mining via multivariable statistics, decision-making (optimization) for systems subjected to uncertainties and even in economics. In this context, SDC constitutes an effective primer tool for complex system analysis, control and operational optimizations. In this review paper, a detailed survey of the developments on the research of SDC systems will be made together with their wide spectrum applications and future perspectives.

42 ENGINEERING↗

Wettability variation and its impact on CO 2 storage capacity at the Wyoming CarbonSAFE storage hub: An experimental approach

Meeting global and national net zero carbon emission targets will require geologic carbon disposal. The U.S. Department of Energy (DOE) has accordingly funded significant research in this area, including the Wyoming CarbonSAFE project at Dry Fork Station (DFS) in Campbell County, Wyoming. This work studied wettability on micro- and macro-scales, CO 2 storage potential, and the correlation between the two to support the Wyoming CarbonSAFE project’s subsurface assessment. During the study, a target formation’s wettability was found to affect how much CO 2 can be stored in a given formation. Here, in this study, representative rock samples were selected from the target storage formations— Lakota, Hulett, and Minnelusa—based on the heterogeneity of the lithology, permeability, and porosity of the respective formations. The rock samples are all fine-grained sandstone with variable cementation and bedding structure, including different scales of laminated bedding. The porosity and permeability vary within the range of 9.0–14.3% and 0.1–28.9 mD, respectively. These rock samples were prepared for the micro-scale wettability (contact angle measurement), macro-scale wettability (wettability index derived from unsteady-state flow characterization for the core plugs), and CO 2 storage evaluation. The macro-scale experiments suggested that wettability appeared to dominate the CO 2 storage potential performance during the drainage process, where less water-wet behavior promoted higher CO 2 storage potential. The micro-scale wettability tests showed that the rock samples at the studied reservoir conditions behaved water-wet and became more water-wet as pressure increased. This kind of wettability change discourages further CO 2 storage potential yet benefits the CO 2 residual trapping as the CO 2 injection proceeds for the studied area. The results allow the recommendation of the best reservoir candidate for storage based on wettability that affects CO 2 storage. The work presented in this study provides valuable insights into wettability’s effect on the CO 2 storage capacity and wettability’s importance when identifying the optimal CO 2 storage formation to meet the project’s goals.

58 GEOSCIENCES↗