Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear potential flow”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Sequential Quadratic Programming Algorithm for Nonsmooth Problems with Upper- \({\boldsymbol{\mathcal{C}^2}}\) Objective

An optimization algorithm for nonsmooth nonconvex constrained optimization problems with upper- \({\boldsymbol{\mathcal{C}^2}}\) objective functions is proposed and analyzed. Upper- \({\boldsymbol{\mathcal{C}^2}}\) is a weakly concave property that exists in difference of convex (DC) functions and arises naturally in many applications, particularly certain classes of solutions to parametric optimization problems e.g., recourse of stochastic programming and projection onto closed sets. The algorithm can be viewed as an extension of sequential quadratic programming (SQP) to nonsmooth problems with upper- \({\boldsymbol{\mathcal{C}^2}}\) objectives or a simplified bundle method. It is globally convergent with bounded algorithm parameters that are updated with a trust-region criterion. The algorithm handles general smooth constraints through linearization and uses a line search to ensure progress. The potential inconsistencies from the linearization of the constraints are addressed through a penalty method. In conclusion, the capabilities of the algorithm are demonstrated by solving both simple upper- \({\boldsymbol{\mathcal{C}^2}}\) problems and a real-world optimal power flow problem used in current power grid industry practices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Overall chilled water system energy consumption modeling and optimization

The emergence of increasingly affordable variable-speed drive technology has changed the approach used to control chilled water systems equipped with these drives. The purpose of this research was to develop an integrated chilled water modeling technique that can determine the optimal system setpoints and estimate the energy saving potential of chiller system. The chiller system equipped with Variable Frequency Drives (VFDs) on cooling tower fans and condenser water pumps. To accomplish the objective, physical component models of the centrifugal chiller, cooling tower and condenser water pump were established with the goal of incorporating the system’s condenser water flow rate and cooling tower fan speeds as optimization variables. Furthermore, a cooling load prediction algorithm was developed using a multiple non-linear regression model to approximate the building’s cooling load subject to a range of environmental conditions. The inputs and outputs of the individual component models were linked to estimate how adjusting the cooling tower fan and condenser water pump speed would influence the system’s comprehensive performance. Here, the overall system model was then optimized using a generalized reduced gradient optimization algorithm to determine the potential energy savings through speed control with VFDs and to ascertain a control logic strategy for the building automation system to operate the heating and cooling system. A case-study was performed on a single chiller system at a museum and the model was calibrated according to logged data collected over four months. Results showed that for the system analyzed, the energy saving of optimizing the cooling tower fan system was found to be 12–15%, while the energy saving potential of optimizing the condenser water pump with the cooling tower fan was negligible. Additionally, comparing different cooling tower fan control strategies showed that a wet-bulb approach-based cooling tower control strategy was shown to have the highest correlation to the optimized fan speed with an R 2 of 0.924.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electrochemical mechanisms of leakage‐current in photovoltaic modules

Abstract This paper analyzes the mechanisms and pathways for leakage current flow observed in Si photovoltaic modules subjected to high temperature and humidity and a large voltage bias with respect to ground. The current inside of the frame is in the form of electron motion, but in the glass and polymer, it is in at least a large part attributable to the movement of ions. When the mode of current flow changes from electronic to ionic conduction, electrochemical reactions will take place at the interface. This can include reactions that produce volatile chemical species like H 2 , CO x , and O 2 , along with ionic species such as OH − and H 3 O + . Here, we see evidence of the importance of different charge carriers with different diffusion rates and the influence of electrochemical processes involved. The application of negative voltage to the cell circuit affects the resistivity of glass producing surfaces with poor conductivity but with some increases in the electrochemical potential producing complicated interactions that are important when the voltage is changed. In the polymer, there is the development of a space charge region and a chemical gradient providing two oppositional forces to current flow, which when released create a complicated discharge process. Here, we give a basic understanding of the charge/discharge of PV cells highlighting how the specific mechanisms are important in understanding some of the degradation processes in PV modules. We find that there is evidence of multiple significant charge carrier species with different diffusion time scales. The glass/polymer interface forms a depleted region of higher resistance after prolonged exposure to current. Charge also builds up at the polymer to cell antireflective coating interface and mostly flows to the gridlines to experience electrochemical reactions. These complexities result in non‐linear behavior where the apparent resistivities of the different layers change during charge/discharging processes, making the modeling of the current flow extremely difficult.

14 SOLAR ENERGY↗
DOE OSTI2022

Irradiation of Sintered Mo Disks with Presence of Organic Impurities in He Gas Flow

Argonne National Laboratory (Argonne) is assisting NorthStar Medical Technologies in the development of a domestic supply of 99 Mo. Specifically, the present study focuses on the production of 99 Mo-feed solution used by the RadioGenix™ 99m Tc generator. During the target-irradiation phase of production, impurities can potentially be introduced into the feed, and can lead to disturbance of ligand- 99m Tc complexation chemistry and contamination of the final radiopharmaceutical that directly interacts with the patient. To address this issue, Argonne performed irradiations and chemical processing to identify whether the potential contamination of He flow with hydrocarbon oil during irradiation affects the radiochemical purity of the final K 2 MoO 4 (K 2 TcO 4 ) in 5M KOH solution. To mimic the conditions of real irradiation at NorthStar, Argonne used its electron linear accelerator and Van de Graaff facilities, heated the target and oil source to >800°C, and controlled oxygen in the presence of He during irradiation. Following irradiation, scanning electron microscopy (SEM) and carbon analysis (CA) were used to detect carbon contamination on the solid targets. The radiochemical purity of the dissolved targets was studied via thin-layer paper chromatography (TLC). As a result of these experiments, small regions of the surface of some irradiated disks were found to be high in carbon, but the total carbon content was still negligible in comparison to the reference sample, which had experienced no irradiation or contact with oil. The summarized results from the SEM, CA, and TLC tests lead to the conclusion that even an excess of oil and heating during irradiation do not affect the radiochemical purity of the final product.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

The Circular Economy Lifecycle Assessment and Visualization Framework: A Case Study of Wind Blade Circularity in Texas

Moving the current linear economy toward circularity is expected to have environmental, economic, and social impacts. Various modeling methods, including economic input-output modeling, life cycle assessment, agent-based modeling, and system dynamics, have been used to examine circular supply chains and analyze their impacts. This work describes the newly developed Circular Economy Lifecycle Assessment and Visualization (CELAVI) framework, which is designed to model how the impacts of supply chains might change as circularity increases. We first establish the framework with a discussion of modeling capabilities that are needed to capture circularity transitions; these capabilities are based on the fact that supply chains moving toward circularity are dynamic and therefore not at steady state, may encompass multiple industrial sectors or other interdependent supply chains and occupy a large spatial area. To demonstrate the capabilities of CELAVI, we present a case study on end-of-life wind turbine blades in the U.S. state of Texas. Our findings show that depending on exact process costs and transportation distances, mechanical recycling could lead to 69% or more of end-of-life turbine blade mass being kept in circulation rather than being landfilled, with only a 7.1% increase in global warming potential over the linear supply chain. We discuss next steps for framework development.

17 WIND ENERGY↗
DOE OSTI2026

Ultrasonic-Assisted Extrusion Processing for Enhancing Physical Properties of High-Density Polyethylene by Flow-Induced Crystallization

The evolution of crystallinity resulting from stress imposed on a melt, known as flow-induced crystallinity, can strongly influence the mechanical and physical properties of semicrystalline polymers. This study investigates shear-induced crystallization by applying an ultrasonic field to the melt flow as it passes through dies with various geometries. A custom-built sonication die is employed for controlling the dynamic temperature and shear environment, resulting in molecular alignment and potential for flow-induced crystallization. Application of both conventional and ultrasonic shear rates at the equilibrium melt temperature of high-density polyethylene (HDPE) was investigated to accelerate crystallinity and manipulate the crystal morphology across the film in pursuit of improved mechanical and gas barrier properties without the need for additives or other polymer layers. The relationships among ultrasonic-assisted extrusion processing, polymer structure, and performance were analyzed using wide- and small-angle X-ray scattering (WAXS and SAXS), tensile testing, and oxygen transmission rate (OTR) analysis. Multiple linear regression models were implemented to predict the correlation among HDPE structure, process, and properties. Structural analysis revealed that both conventional and ultrasonic shear rates had the most significant influence on lamellar spacing and redistribution of rigid and soft amorphous fractions within the crystalline domains, ultimately dictating the mechanical and physical properties of the films. The goal is to explore the potential of the ultrasonic-assisted high crystallinity monolayer that can replace some of the functionality of complex, heterogeneous multilayer packaging with a single-material film having enhanced oxygen barrier properties.

crystallinity↗
DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗
DOE OSTI2022

A generic scenario analysis of end-of-life plastic management: Chemical additives

Plastic growing demand and the increment in global plastics production have raised the number of spent plastics, out of which over 90% are either landfilled or incinerated. Both methods for handling spent plastics are susceptible to releasing toxic substances, damaging air, water, soil, organisms, and public health. Improvements to the existing infrastructure for plastics management are needed to limit chemical additive release and exposure resulting from the end-of-life (EoL) stage. This article analyzes the current plastic waste management infrastructure and identifies chemical additive releases through a material flow analysis. Additionally, we performed a facility-level generic scenario analysis of the current U.S. EoL stage of plastic additives to track and estimate their potential migration, releases, and occupational exposure. Potential scenarios were analyzed through sensitivity analysis to examine the merit of increasing recycling rates, using chemical recycling, and implementing additive extraction post-recycling. Our analyses identified that the current state of plastic EoL management possesses high mass flow intensity toward incineration and landfilling. Although maximizing the plastic recycling rate is a reasonably straightforward goal for enhancing material circularity, the conventional mechanical recycling method requires improvement because major chemical additive release and contamination routes act as obstacles to achieving high-quality plastics for future reuse and should be mitigated through chemical recycling and additive extraction. In conclusion, the potential hazards and risks identified in this research create an opportunity to design a safer closed-loop plastic recycling infrastructure to handle additives strategically and support sustainable materials management efforts to transform the US plastic economy from linear to circular.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Physics-constrained deep learning of nonlinear normal modes of spatiotemporal fluid flow dynamics

In this study, we present a physics-constrained deep learning method to discover and visualize from data the invariant nonlinear normal modes (NNMs) which contain the spatiotemporal dynamics of the fluid flow potentially containing strong nonlinearity. Specifically, we develop a NNM-physics-constrained convolutional autoencoder (NNM-CNN-AE) integrated with a multi-temporal-step dynamics prediction block to learn the nonlinear modal transformation, the NNMs containing the spatiotemporal dynamics of the flow, and reduced-order reconstruction and long-time future-state prediction of the flow fields, simultaneously. In test cases, we apply the developed method to analyze different flow regimes past a cylinder, including laminar flows with low Reynolds number in transient and steady states (RD = 100) and high Reynolds number flow (RD = 1000), respectively. The results indicate that the identified NNMs are able to reveal the nonlinear spatiotemporal dynamics of these flows, and the NNMs-based reduced-order modeling consistently achieves better accuracy with orders of magnitudes smaller errors in construction and prediction of the nonlinear velocity and vorticity fields, compared to the linear proper orthogonal decomposition (POD) method and the Koopman-constrained-CNN-AE using the same number or dimension of modes. We perform an analysis of the modal energy distribution of NNMs and find that compared to POD modes, the few fundamental NNMs capture a very high level of total energy of the flow, which is advantageous for reduced-order modeling and representation of the complex flows. Finally, we discuss the potentials and limitations of the presented method.

Mechanics↗
DOE OSTI2023

Design and Development of a Particle Flow Control Mechanism for Particle-Based CSP System Applications

The submission is a two-page extended abstract. Within particle-based Concentrated Solar Power Systems (CSP), particles are the heat transfer medium, and therefore, the heat transfer effectiveness of these systems is reliant on their ability to effectively control particle mass flow rate. The vast majority of approaches for particle flow control in high temperature CSP is based on a slide-gate approach where a gate is linearly actuated across a tear-drop shaped opening. The slide-gate mechanism has a sensitivity of 0.354 g/s-mm, which is defined as the mass flow rate over the position of the gate. In the present work, we have designed and tested a new particle flow control mechanism with the potential for faster response times and different sensitivities than a slide gate approach enabling more control of particle flow.

14 SOLAR ENERGY↗
DOE OSTI2025

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗
DOE OSTI2023

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S=αN/f, whereα,N, andfare the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases withS, asymptotically matching Lentz and Chapman forS< 1. When the solution conserves interior potential vorticity, the momentum transport is proportional toS 2 forS< 1 and is in better agreement with numerical simulations.

Oceanography↗
DOE OSTI2023

Variation of Aliphatic Diisocyanates in Bio-Based TPUs

The utilization of bio-based materials for polymer production poses a challenge to both the industrial and academic sectors due to the availability and production costs of the necessary raw materials. Diisocyanates necessary for polyurethane synthesis have posed a particular challenge, given the paucity of natural diamine precursors for traditional phosgenation. We recently developed a phosgene-free flow chemistry methodology that allows for the safe, efficient, and scalable preparation of diisocyanates from naturally produced diacids. This chemistry broadens the potential for the development of renewable diisocyanates for a broad variety of applications. Here, we expand upon this work by demonstrating the scaled production of a panel of bio-based linear, aliphatic diisocyanates, with chain lengths ranging from four to eight carbons in good yields of 57–81% and high purity (~97%). Furthermore, these are applied to the synthesis of thermoplastic polyurethanes (TPUs) containing up to a 100% renewable carbon content. TPUs formulated using shorter carbon-chained diisocyanates displayed a higher tensile strength compared to those formulated using longer chains. Ready access to diisocyanates of varying chain length affords the ability to tailor TPU properties through careful selection of the diisocyanate, polyol, and chain extender, offering hard and soft TPUs for multiple end-use applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modified equilibrium distributions for Cooper-Frye particlization

We introduce a positive-definite single-particle distribution that is suitable for describing the transition from a macroscopic hydrodynamic to a microscopic kinetic description during the late stages of heavy-ion collisions in the presence of moderately large viscous corrections. The modified equilibrium distribution function can be constructed with hydrodynamic input from either relativistic viscous fluid dynamics or anisotropic fluid dynamics. We test the modified equilibrium distribution's hydrodynamic output for a stationary hadron resonance gas subject to either shear stress, bulk pressure, or baryon diffusion current at a given freeze-out temperature and baryon chemical potential. While it does not reproduce all components of the net baryon current and energy-momentum tensor exactly, it significantly improves upon the customary linearized approximations for the nonequilibrium correction δf n which typically lead to unphysical negative distribution functions at large particle momenta. A comparison of particle spectra and p T –differential elliptic flow coefficients from the Cooper–Frye formula computed with the modified equilibrium distribution and with linearized δf n corrections is presented, for two different (2 + 1)–dimensional hypersurfaces corresponding to central and noncentral Pb + Pb collisions at the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

The value of in-reservoir energy storage for flexible dispatch of geothermal power

Geothermal systems making use of advanced drilling and well stimulation techniques have the potential to provide tens to hundreds of gigawatts of clean electricity generation in the United States by 2050. With near-zero variable costs, geothermal plants have traditionally been envisioned as providing “baseload” power, generating at their maximum rated output at all times. However, as variable renewable energy sources (VREs) see greater deployment in energy markets, baseload power is becoming increasingly less competitive relative to flexible, dispatchable generation and energy storage. Herein we conduct an analysis of the potential for future geothermal plants to provide both of these services, taking advantage of the natural properties of confined, engineered geothermal reservoirs to store energy in the form of accumulated, pressurized geofluid and provide flexible load-following generation. We develop a linear optimization model based on multi-physics reservoir simulations that captures the transient pressure and flow behaviors within a confined, engineered geothermal reservoir. We then optimize the investment decisions and hourly operations of a power plant exploiting such a reservoir against a set of historical and modeled future electricity price series. Further, we find that operational flexibility and in-reservoir energy storage can significantly enhance the value of geothermal plants in markets with high VRE penetration, with energy value improvements of up to 60% relative to conventional baseload plants operating under identical conditions. Across a range of realistic subsurface and operational conditions, our modeling demonstrates that confined, engineered geothermal reservoirs can provide large and effectively free energy storage capacity, with round-trip storage efficiencies comparable to those of leading grid-scale energy storage technologies. Optimized operational strategies indicate that flexible geothermal plants can provide both short- and long-duration energy storage, prioritizing output during periods of high electricity prices. Sensitivity analysis assesses the variation in outcomes across a range of subsurface conditions and cost scenarios.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗