Results for “dissipative systems” — Page 72

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2025

Microstructure-sensitive mechanical behavior of an additively manufactured psuedoelastic shape memory alloy

The additive manufacturing of shape memory alloys into complex geometries enables fabrication of advanced functional systems across a variety of fields and domains. This work presents results focused on the mechanical behavior of additively manufactured shape memory pseudoelastic NiTi. The deformation induced solid state phase transformation from austenite to martensite allows this system to accommodate large recoverable strains. This deformation behavior is fundamentally driven by crystal-scale transformation physics. Laser powder bed fusion processing reveals that the resulting microstructure, both grain morphology and crystallographic texture, is strongly dependent on the manufacturing processing history. Exhaustive mechanical testing demonstrates that these microstructural factors strongly impact both tensile and cyclic stress–strain behavior. Cyclic dissipative behavior, however, is similar across all tested microstructures following an initial transient period. Remarkably, analysis of spatial strain fields during tensile loading reveals two distinctly different localization “modes”. The first is initiation of localized deformation bands which continuously propagate through the tensile bar during loading. In the second mode localization is observed but lacks propagation; instead additional localization cites nucleate during subsequent loading. The latter phenomena is suspected to be driven by grain-scale deformation physics as the localized band morphologies coincide with grain morphologies. These phenomena strongly impact the resulting aggregate stress–strain behavior. Hence, manufacturers and designers of psuedoelastic functional components must at the very least consider the potential variability in properties when considering additive manufacturing processing. More ideally the process–structure–property relations can be used to further tailor and optimize final functional performance.

Additive manufacturing↗

DOE OSTI2025

Geometric transport signatures of strained multi-Weyl semimetals

The minimal coupling of strain to Dirac and Weyl semimetals, and its modeling as a pseudogauge field has been extensively studied, resulting in several proposed topological transport signatures. In this work, we study the effects of strain on higher winding number Weyl semimetals and show that strain is not a pseudogauge field for any winding number larger than one. Here we focus on the double-Weyl semimetal as an illustrative example to show that the application of strain splits the higher winding number Weyl nodes and produces an anisotropic Fermi surface. Specifically, the Fermi surface of the double-Weyl semimetal acquires nematic order. By extending chiral kinetic theory for such nematic fields, we determine the effective gauge fields acting on the system and show how strain induces anisotropy and affects the geometry of the semiclassical phase space of the double-Weyl semimetal. Further, the strain-induced deformation of the Weyl nodes results in transport signatures related to the covariant coupling of the strain tensor to the geometric tensor associated with the Weyl nodes giving rise to strain-dependent dissipative corrections to the longitudinal as well as the Hall conductance. Thus, by extension, we show that in multi-Weyl semimetals, strain produces geometric signatures rather than topological signatures. Further, we highlight that the most general way to view strain is as a symmetry-breaking field rather than a pseudogauge field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1974

The stability of motion of satellites with cavities partially filled with liquid

The purpose of the present paper is to determine rigorously the effect of viscous dissipation on the stability and time-dependent motion of a spinning rigid body with a spherical cavity partially filled with liquid. The cavity is such that its center does not lie on the body axis of rotation. The work includes the problem formulation, consisting of the boundary-value problem for the liquid and moment equations for the entire system. Because of large Reynold's numbers involved, viscosity effects are negligible everywhere except for a thin boundary layer near the wetted surface. Using a boundary-layer analysis, the liquid problem is reduced to the solution of the inviscid equations subject to modified boundary conditions.

Nayfeh, A. H.↗

NASA NTRS1977

Solar array maximum power tracking with closed-loop control of a 30-centimeter ion thruster

A new solar array/ion thruster system control concept has been developed and demonstrated. An ion thruster beam load is used to automatically and continuously operate an unregulated solar array at its maximum power point independent of variations in solar array voltage and current. Preliminary tests were run which verified that this method of control can be implemented with a few, physically small, signal level components dissipating less than two watts.

Gruber, R. P.↗

NASA NTRS1986

Magnetic energy storage and conversion in the solar atmosphere

According to the approach employed in this investigation, particularly important simple configurations of magnetic field and plasma are identified, and it is attempted to achieve an understanding of the large-scale dynamic processes and transformations which these systems can undergo. Fundamental concepts are discussed, taking into account aspects of magnetic energy generation, ideal MHD theory, non-MHD properties, the concept of 'anomalous' resistivity, and global electrodynamic coupling. Questions of magnetically controlled energy conversion are examined, giving attention to magnetic modifications of plasma transport, the transition region structure and flows, channeling and acceleration of plasma, channeling and dissipation of MHD waves, and anomalous dissipation of field-aligned currents. A description of the characteristics of magnetohydrodynamic energy conversion is also provided, and outstanding questions are discussed.

Spicer, D. S.↗

NASA NTRS1986

Daytime lidar measurements of tidal winds in the mesospheric sodium layer at Urbana, Illinois

For more than 15 years lidar systems have been used to study the chemistry and dynamics of the mesospheric sodium layer. Because the layer is an excellent tracer of atmospheric wave motions, sodium lidar has proven to be particularly useful for studying the influence of gravity waves and tides on mesospheric dynamics. These waves, which originate in the troposphere and stratosphere, propagate through the mesosphere and dissipate their energy near the mesopause making important contributions to the momentum and turbulence budget in this region of the atmosphere. Recently, the sodium lidar was modified for daytime operation so that wave phenomena and chemical effects could be monitored throughout the complete diurnal cycle. The results of continuous 24 hour lidar observations of the sodium layer structure are presented alond with measurement of the semidiurnal tidal winds.

Kwon, K. H.↗

NASA NTRS1995

A turbulence module for the NPARC code

A turbulence module is developed for the 2D version of the NPARC code which is currently restricted to planar or axisymmetric flows without swirling. Four turbulence models have been built into the module: Baldwin-Lomax, Chien, Shih-Lumley and CMOTT models. The first is a mixing-length eddy-viscosity model which is mainly used for initialization of computational fields and the last three are the low Reynolds number two-equation models. Unlike chien's model, both the Shih-Lumley and CMOTT models do not involve the dimensionless wall distance y(sup +), an advantage for separated flow calculations. Contrary to the NPARC and most other compressible codes, the non-delta form of transport equations is used which leads to a simpler linearization and is more effective than using the delta form in ensuring the positiveness of the turbulent kinetic energy and its dissipation rate. To reduce numerical diffusion while maintaining necessary stability, a second-order accurate and bounded scheme is used for the convective terms of the turbulent transport equations. This scheme is implemented in a deferred correction manner so that the main coefficients of the resulting difference equations are always positive, thus making the numerical solutions process unconditionally stable. The system of equations is solved via a decoupled method and by the alternating direction TDMA of Thomas. The module can be easily linked to the NPARC code for turbulent flow calculations.

Zhu, J.↗

NASA NTRS2009

Enceladus: An Active Cryovolcanic Satellite

Enceladus is one of the most remarkable satellites in the solar system, as revealed by Cassini's detection of active plumes erupting from warm fractures near its south pole. This discovery makes Enceladus the only icy satellite known to exhibit ongoing internally driven geological activity. The activity is presumably powered by tidal heating maintained by Enceladus 2:1 mean-motion resonance with Dione, but many questions remain. For instance, it appears difficult or impossible to maintain the currently observed radiated power (probably at least 6 GW) in steady state. It is also not clear how Enceladus first entered its current self-maintaining warm and dissipative state initial heating from non-tidal sources is probably required. There are also many unanswered questions about Enceladus interior. The silicate fraction inferred from its density of 1.68 g per cubic centimeter is probably differentiated into a core, though we have only indirect evidence for differentiation. Above the core there is probably a global or regional water layer, inferred from several models of tidal heating, and an ice shell thick enough to support the ~1 kilometer amplitude topography seen on Enceladus. It is possible that dissipation is largely localized beneath the south polar region. Enceladus surface geology, ranging from moderately cratered terrain to the virtually crater-free active south polar region, is highly diverse, tectonically complex, and remarkably symmetrical about the rotation axis and the direction to Saturn. South polar activity is concentrated along the four tiger stripe fractures, which radiate heat at temperatures up to at least 167 K and are the source of multiple plumes ejecting ~200 kilograms per second of H2O vapor along with significant N2 (or C2H4), CO2, CH4, NH3, and higher-mass hydrocarbons. The escaping gas maintains Saturn's neutral gas torus, and the plumes also eject a large number of micron-sized H2O ice grains that populate Saturn's E-ring. The mechanism that powers the plumes is not well understood, and whether liquid water is involved is a subject of active debate (but likely nonetheless). Enceladus provides a promising potential habitat for life in the outer solar system, and the active plumes allow the unique opportunity for direct sampling of that zone. Enceladus is thus a prime target for Cassini's continued exploration of the Saturn system, and will be a tempting target for future missions.

Spencer, J. R.↗

NASA NTRS2017

STURM: Resuspension Mesocosms with Realistic Bottom Shear Stress and Water Column Turbulence for Benthic-Pelagic Coupling Studies: Design and Applications

Shear TUrbulence Resuspension Mesocosm (STURM) tanks, with high instantaneous bottom shear stress and realistic water column mixing in a single system, allow realistic benthic-pelagic coupling studies with sediment resuspension. The 1 m3 tanks can be programmed to produce tidal or episodic sediment resuspension for extended time periods (e.g. 4 weeks), over muddy sediments with a variety of benthic organisms. A resuspension paddle produces uniform bottom shear stress across the sediment surface while gently mixing a 1 m deep overlying water column. The STURM tanks can be programmed to different magnitudes, frequencies, and durations of bottom shear stress (and thus resuspension) with proportional water column turbulence levels over a wide range of mixing settings for benthic-pelagic coupling experiments. Over ten STURM calibration settings, RMS turbulent velocity ranged from 0.26 to 4.52 cm s−1, energy dissipation rate from 0.0016 to 2.65 cm2 s−3, the average bottom shear stress from 0.0035 to 0.19 Pa, and the instantaneous maximum bottom shear stress from 0.07 to 1.7 Pa. We have performed four 4-week benthic-pelagic coupling ecosystem experiments with tidal resuspension and stepwise erosion experiments (both with and without infaunal bivalves), carried out experiments on oyster biodeposit resuspension, mimicked storms overlain on tidal resuspension, and studied the effects of varying frequency and duration of resuspension on sedimentary contaminant release. The large size of the tanks allows water quality and particle measurements using standard oceanographic instrumentation. The realistic scale and complexity of the contained ecosystems has revealed indirect feedbacks and responses that are not observable in smaller, less complex experimental systems.

Porter, Elka T.↗

NASA NTRS2000

Computation of Feedback Aeroacoustic System by the CE/SE Method

It is well known that due to vortex shedding in high speed flow over cutouts, cavities, and gaps, intense noise may be generated. Strong tonal oscillations occur in a feedback cycle in which the vortices shed from the upstream edge of the cavity convect downstream and impinge on the cavity lip, generating acoustic waves that propagate upstream to excite new vortices. Numerical simulation of such a complicated process requires a scheme that can: (1) resolve acoustic waves with low dispersion and numerical dissipation, (2) handle nonlinear and discontinuous waves (e.g. shocks), and (3) have an effective (near field) nonreflecting boundary condition (NRBC). The new space time conservation element and solution element method, or CE/SE for short, is a numerical method that meets the above requirements.

Loh, Ching Y.↗

DOE OSTI2025

Cryogenically enhanced quasi-optical resonator for megawatt pulsed millimeter-wave sources

We report the design, cryogenic optimization, and performance modeling of a compact quasi-optical ring resonator intended to compress microwaves pulses at 170 and 250 GHz to the megawatt level. By combining ultra-low-loss CVD diamond and gold-doped silicon wafers with high-RRR copper mirrors, the calculated unloaded quality factor exceeds 4.3 × 10 5 at 20 K and yields simulated gains up to $\mathscr{G}$ = 4.1 × 10 3 . Coupling the resonator with a laser-driven semiconductor switch described by an extended Vogel model shows that 1 MW, nanosecond pulses can be generated from only 445 W of microwave drive power while dissipating 272 W into the cryostat. A practical cooling architecture using two Gifford–McMahon stages (20 and 80 K) is proposed, demonstrating that high-repetition-rate (10–20 kHz) operation is feasible with commercially available cryocoolers. The results outline a clear path toward cost-effective, table-top sources for extreme-ultraviolet lithography, dynamic nuclear polarization, and fusion systems.

Panisset, Constant [Ecole Polytechnique Federale L↗

NASA NTRS1985

A parameterization technique for nonlinear spectral models

An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure of hydrodynamical systems, regardless of their complexity, is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. As a prototype of the generic problem containing this typical structure, the parameterization technique is applied to various three component subsets of a five component nonlinear spectral model of forced, dissipative quasi-geostrophic flow in a channel. The results obtained lead to specification of the necessary data coverage requirements for applying the technique in general.

Gelaro, R.↗

NASA NTRS1986

A parameterization technique for nonlinear spectral models

An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure of hydrodynamical systems, regardless of their complexity, is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. As a prototype of the generic problem containing this typical structure, the parameterization technique is applied to various three component subsets of a five component nonlinear spectral model of forced, dissipative quasi-geostrophic flow in a channel. The results obtained lead to specification of the necessary data coverage requirements for applying the technique in general.

Gelaro, R.↗

NASA NTRS1989

Vibration suppression for precision segmented reflector backup structure

An ongoing effort at the JPL on vibration suppression for the Precision Segmented Reflector backup structure is centered on the vibration damping augmentation through a system consisting of active and passive damping members. An active member is a structural member with built-in piezoelectric actuator and sensors. An efficient method was developed for optimal placement of active and passive damping members in the truss-type backup structure. A simple synergistic model between the active and passive damping was proposed based on a weighted energy dissipation criterion. A baseline passive member design with constrained viscoelastic material treatment was used as the source of passive damping.

Chen, G.-S.↗

NASA NTRS1999

Particle Segregation in Collisional Shearing Flows

The size segregation of flowing or shaken grains is a commonly observed phenomenon in industrial processes and in nature. In systems that do not involve much agitation of the grains, several mechanisms that involve gravity have been identified as leading to such segregation. In highly agitated flows, there is a mechanism independent of gravity that is available to drive separation of different grains. This is associated with spatial gradients in the energy of their velocity fluctuations. Because collisions between grains inevitably dissipate energy, collisional granular shear flows are usually of limited extent in the direction transverse to the flow. One consequence of this is that shear flows are strongly influenced by their boundaries. Because grains, on average, slip relative to boundaries, a bumpy or frictional boundary can convert slip energy into fluctuation energy. However, because each collision between a grain and the boundary dissipates fluctuation energy, there is a competition between production and dissipation. In principle, it is possible to design the geometry of the boundary - for example, the size and spacing of the bumps - so that the boundary either produces or dissipates fluctuation energy. This permits the control of the component of the spatial gradient of the fluctuation energy that is normal to the boundary. The gradients in fluctuation energy established by such boundaries may be exploited to drive the separation by size or other properties in a binary mixture of spherical grains. Microgravity makes the visual observations possible by permitting us to employ moderate rates of shear. On earth, the effects of gravity can be minimized by shearing so rapidly that the particle pressure overwhelms gravity. However, in this event, separation takes place too rapidly for visual observation, buoyancy and/or condensation associated with the centripetal acceleration must be accounted for, and the particles can be severely damaged. Because, in the absence of gravity, the only available time scale is proportional to the speed of the moving boundary, this speed can be made arbitrarily slow to permit observations and to avoid particle damage, without altering the phenomenon under study. The primary goal of this research is to carry out a physical experiment in which particle segregation is induced and maintained in a collisional flow of a binary mixture of two different types of spheres.

Jenkins, J. T.↗

NASA NTRS1998

Model Refinement Using Eigensystem Assignment

This paper presented a novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, respectively. Once, the free elements of the structural matrices have been defined, the problem of model refinement reduces to obtaining position, velocity, and acceleration gain matrices, which reassign a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to those obtained from system identification test data. A sequential procedure is used to assign one self-conjugate pair of closed-loop eigenvalues at each step using symmetric output feedback gain matrices, and the closed-loop eigenvectors are partially assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure can also impose that gain matrices be dissipative in order to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley (CSI Evolutionary Model) is presented to demonstrate the feasibility of the proposed approach.

Maghami, Pieman G.↗