Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Electric Potential and Electric Field Imaging with Dynamic Applications: 2017 Research Award Innovation

The technology and methods for remote quantitative imaging of electrostatic potentials and electrostatic fields in and around objects and in free space is presented. Electric field imaging (EFI) technology may be applied to characterize intrinsic or existing electric potentials and electric fields, or an externally generated electrostatic field may be used for illuminating volumes to be inspected with EFI. The baseline sensor technology (e-Sensor) and its construction, optional electric field generation (quasi-static generator), and current e- Sensor enhancements (ephemeral e-Sensor) are discussed. Critical design elements of current linear and real-time two-dimensional (2D) measurement systems are highlighted, and the development of a three dimensional (3D) EFI system is presented. Demonstrations for structural, electronic, human, and memory applications are shown. Recent work demonstrates that phonons may be used to create and annihilate electric dipoles within structures. Phonon induced dipoles are ephemeral and their polarization, strength, and location may be quantitatively characterized by EFI providing a new subsurface Phonon-EFI imaging technology. Initial results from real-time imaging of combustion and ion flow, and their measurement complications, will be discussed. These new EFI capabilities are demonstrated to characterize electric charge distribution creating a new field of study embracing areas of interest including electrostatic discharge (ESD) mitigation, crime scene forensics, design and materials selection for advanced sensors, combustion science, on-orbit space potential, container inspection, remote characterization of electronic circuits and level of activation, dielectric morphology of structures, tether integrity, organic molecular memory, atmospheric science, and medical diagnostic and treatment efficacy applications such as cardiac polarization wave propagation and electromyography imaging.

Generazio, Ed↗
DOE OSTI2025

Transportability of exogenous microbial community correlates with interwell connectivity in deep aquifers

Subsurface resource engineering operations often utilize continuous injection of externally-sourced water into geological reservoirs for formation pressure maintenance, resource recovery or energy/waste storage. Such injected water generally contains naturally occurring microbes. Little is known, however, about how the injectate microbes transport through geological media as a community, how such transportability is affected by injector-producer connectivity, and whether such knowledge can be utilized for flowpath characterization. In this study, we analyzed daily-to-weekly timeseries microbial community data from the injected- and produced-fluids of a ten-month flow test at a deep, well-characterized engineered aquifer. We found that the injectate microbial community was distinct from the indigenous community at the amplicon sequence variant (ASV) level, and that the transportability of injectate community towards a given producer, quantified by an “nASV-Overlap” metric we propose, had strong and significant positive correlation with known injector-producer connectivities at our site. This suggests that the better the connectivity, the higher the probability for more injectate species to flow through the interwell region and arrive at a producer. Because interwell connectivity is an important yet usually unknown parameter in subsurface resource engineering, such correlation in turn points to nASV-Overlap as a useful indicator of interwell connectivity for aquifer characterization and long-term monitoring. Based on our findings, an nASV-Overlap-based microbial tracing approach was developed for characterizing and monitoring the relative connectivities across multiple producers with a given injector. A side-by-side comparison between the new nASV-Overlap approach and traditional artificial tracer methods is presented, and their respective strengths and limitations are discussed.

Deep biosphere↗
DOE OSTI2020

Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions

Understanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo–Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High–resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first–order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. Here, the results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.

58 GEOSCIENCES↗
DOE OSTI2022

Geological activity shapes the microbiome in deep-subsurface aquifers by advection

Subsurface environments host diverse microorganisms in fluid-filled fractures; however, little is known about how geological and hydrological processes shape the subterranean biosphere. Here, we sampled three flowing boreholes weekly for 10 mo in a 1478-m-deep fractured rock aquifer to study the role of fracture activity (defined as seismically or aseismically induced fracture aperture change) and advection on fluid-associated microbial community composition. We found that despite a largely stable deep-subsurface fluid microbiome, drastic community-level shifts occurred after events signifying physical changes in the permeable fracture network. The community-level shifts include the emergence of microbial families from undetected to over 50% relative abundance, as well as the replacement of the community in one borehole by the earlier community from a different borehole. Null-model analysis indicates that the observed spatial and temporal community turnover was primarily driven by stochastic processes (as opposed to deterministic processes). We, therefore, conclude that the observed community-level shifts resulted from the physical transport of distinct microbial communities from other fracture(s) that outpaced environmental selection. Given that geological activity is a major cause of fracture activity and that geological activity is ubiquitous across space and time on Earth, our findings suggest that advection induced by geological activity is a general mechanism shaping the microbial biogeography and diversity in deep-subsurface habitats across the globe.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

"Hidden" hydrothermal technical potential & technoeconomics: Revealing permeability & fluids with more data

Historical hydrothermal estimates have largely relied on temperature or heat flow estimates ignoring the need for natural flowing fluids. More accurate hydrothermal estimates require some indication of permeability and fluids that naturally exist in the subsurface. This paper describes a novel approach that includes proxies of permeability and fluids in hydrothermal estimates by leveraging the relatively data-rich Great Basin. Specifically, nameplate capacities (megawatts) of operating geothermal plants, negative (0 megawatt) locations and 48 geophysical and geologic features are used to used in eXtreme Gradient Boosting (XGBoost) regression to make hydrothermal capacity predictions. Additionally, this work inputs the XGBoost-based hydrothermal predictions into the Renewable Energy Potential (reV) model to quantify technical capacity, its uncertainty and techno-economics. Compared to historical hydrothermal estimates, these predictions adhere to the 37 operating geothermal plants and negative locations. We present a method for subsampling the negative sites to bring the labels into balance that uses the geologic domain knowledge to proportionally represent negatives. Overall, the distributions of the hydrothermal technical capacity and the site levelized cost of energy are respectively much tighter, lower and more accurate than the previous estimates for the Great Basin, as they include geological and geophysical surrogates for permeability and fluids. Percentile (50th and 90th, median and high estimate, respectively) models provide bookends for these metrics.

13 HYDRO ENERGY↗
NASA NTRS1993

Sunspot dynamics

This report describes recent results of our theoretical and observational work on dynamical phenomena in sunspots. The overall goal of this research has been a better understanding of the various oscillatory, transient, and steady motions in a sunspot and their relation to the basic structure of the sunspot. The principal topics of the research reported here are the following: (1) sunspot seismology, i.e., the study of the interaction of solar p-modes with a sunspot as a probe of the subsurface structure of a sunspot; (2) local sources of acoustic waves in the solar photosphere; and (3) siphon flows in isolated magnetic flux tubes and their relation to the photospheric Evershed flow and to intense magnetic elements outside of sunspots.

Thomas, John H.↗
NASA NTRS1992

Thermographic imaging of cracks in thin metal sheets

The presence of cracks significantly decreases the structural integrity of thin metal sheets used in aerospace applications. Thermographic detection of surface temperature variations due to these cracks is possible after external heating. An approximate line source of heat is used to produce an inplane flow of heat in the sheet. A crack in the sheet perturbs the inplane flow of heat and can be seen in an image of the surface temperature of the sheet. An effective technique for locating these perturbations is presented which reduces the surface temperature image to an image of variations in the inplane heat flow. This technique is shown to greatly increase the detectability of the cracks. This thermographic method has advantages over other techniques in that it is able to remotely inspect a large area in a short period of time. The effectiveness of this technique depends on the shape, position and orientation of the heat source with respect to the cracks as well as the extent to which the crack perturbs the surface heat flow. The relationship between these parameters and the variation in the heat flow is determined both by experimental and computational techniques. Experimental data is presented for through-the-thickness, subsurface and surface EDM notches. Data for through-the-thickness fatigue cracks are also presented.

Cramer, K. E.↗
NASA NTRS2001

The Evolution of the Martian Hydrosphere: Implications for the Fate of a Primordial Ocean and the Current State of the Northern Plains

In this paper we consider the hydraulic and thermal conditions that gave rise to the elevated source regions of the Late Hesperian outflow channels and explore their implications for the evolution of the Martian hydrosphere. We find that if the outflow channel floodwaters were derived from a subpermafrost aquifer, then it implies that, throughout the planet's first billion years of evolution, as much as one third of its surface was covered by standing bodies of water and ice. Following the development of the global dichotomy, the bulk of this water would have existed as an ice-covered ocean in the northern plains. We demonstrate that the progressive crustal assimilation of this early surface reservoir of H2O (punctuated by possible episodes of less extensive flooding) was a natural consequence of the planet's subsequent climatic and geothermal evolution-potentially cycling the equivalent of a km-deep global ocean of water through the atmosphere and subsurface every approx. 10(exp 9) years. In response to the long-term decline in planetary heat flow, the progressive coldtrapping of H2O into the growing cryosphere is expected to have significantly depleted the original inventory of groundwater-a development that could well explain the apparent decline in outflow channel activity observed during the Amazonian. Although primarily a theoretical analysis, our findings appear remarkably consistent with the geomorphic and topographic evidence that Mars once possessed a primordial ocean and that a substantial relic of that body continues to survive as massive ice deposits within the northern plains. Confirmation of the presence of such deposits, combined with the potential detection of a global-scale groundwater system, would provide persuasive support for the validity of this analysis.

Clifford, Stephen M.↗
NASA NTRS2012

Development of a Deep-Penetrating, Compact Geothermal Heat Flow System for Robotic Lunar Geophysical Missions

Geothermal heat flow measurements are a high priority for the future lunar geophysical network missions recommended by the latest Decadal Survey of the National Academy. Geothermal heat flow is obtained as a product of two separate measurements of geothermal gradient and thermal conductivity of the regolith/soil interval penetrated by the instrument. The Apollo 15 and 17 astronauts deployed their heat flow probes down to 1.4-m and 2.3-m depths, respectively, using a rotary-percussive drill. However, recent studies show that the heat flow instrument for a lunar mission should be capable of excavating a ~3-m deep hole to avoid the effect of potential long-term changes of the surface thermal environment. For a future robotic geophysical mission, a system that utilizes a rotary/percussive drill would far exceed the limited payload and power capacities of the lander/rover. Therefore, we are currently developing a more compact heat flow system that is capable of 3-m penetration. Because the grains of lunar regolith are cohesive and densely packed, the previously proposed lightweight, internal hammering systems (the so-called moles ) are not likely to achieve the desired deep penetration. The excavation system for our new heat flow instrumentation utilizes a stem which winds out of a pneumatically driven reel and pushes its conical tip into the regolith. Simultaneously, gas jets, emitted from the cone tip, loosen and blow away the soil. Lab tests have demonstrated that this proboscis system has much greater excavation capability than a mole-based heat flow system, while it weighs about the same. Thermal sensors are attached along the stem and at the tip of the penetrating cone. Thermal conductivity is measured at the cone tip with a short (1- to 1.5-cm long) needle sensor containing a resistance temperature detector (RTD) and a heater wire. When it is inserted into the soil, the heater is activated. Thermal conductivity of the soil is obtained from the rate of temperature increase during the heating. By stopping during the excavation, it is possible to measure thermal conductivities at different depths. The gas jets are turned off when the penetrating cone reaches the target depth. Then, the stem pushes the needle sensor into the undisturbed soil at the bottom of the hole and carries out a thermal conductivity measurement. When the measurement is complete, the system resumes excavation. RTDs, placed along the stem at short (approx 30 cm) intervals, will monitor long-term temperature stability of the subsurface. Temperature in the shallow subsurface would fluctuate with the diurnal, annual, and precession cycles of the Moon. These thermal waves penetrate to different depths into the regolith. Longterm monitoring of the subsurface temperature would allow us to accurately delineate these cyclic signals and separate them from the signal associated with the outward flow of the Moon s endogenic heat. Further, temperature toward bottom of the 3-m hole should be fairly stable after the heat generated during the excavation dissipates into the surrounding soil. The geothermal gradient may be determined reliably from temperature measurements at the RTDs near the bottom. In order to minimize the heat conduction along the stem from affecting the geothermal gradient measurements, we plan to use low-conductive materials for the stem and develop a mechanism to achieve close coupling between the RTDs and the wall of the excavated hole.

Nagihara, Seiichi↗
NASA NTRS1993

The thumbprint terrain: What will Mars Observer tell us?

Some of the more puzzling features seen on Mars are those known as curvilinear features or 'thumb print' terrain, types of patterned ground found in the northern plains. The thumbprint terrain, named for its resemblance to the lines of a human thumbprint is found on what appears to be level, relatively low-lying ground near the crustal dichotomy boundary. It is often found near the termini of large channels. There are three types of thumbprint terrain in the classification of Rossbacher and Judson. The first consists of ridges and depressions about a kilometer wide, separated by a few kilometers, and with an apparent relief of 40-100 m. The second consists of steep-sided, flat-floored depressions about a kilometer wide, separated by a few kilometers, and with an apparent relief of 10-290 m. The third consists only of albedo markings, with a typical scale of about a kilometer. Many models were proposed to explain the origin of this terrain. It was suggested that it was caused by lava flows, removal of debris mantles, glaciers, or karst. However, the most popular models at present involve the action of subsurface ice to form such thermokarst features as striped ground, solifluction lobes, and/or linear, ice-cored ridges. There are several instruments on the Mars Observer spacecraft that will be able to provide us with information useful in distinguishing between these models.

Schaefer, M. W.↗
DOE OSTI2021

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗
NASA NTRS1991

Strain accommodation beneath structures on Mars

A recent review of tectonic features on Mars shows that most of their subsurface structures can be confidently extended only a few kilometers deep (exceptions are rifts, in which bounding normal faults penetrate the entire brittle lithosphere, with ductile flow at deeper levels). Nevertheless, a variety of estimates of elastic lithosphere thickness and application of accepted failure criteria under likely conditions on Mars suggest a brittle lithosphere that is many tens of kilometers thick. This raises the question of how the strain (extension or shortening) accommodated by grabens and wrinkle ridges within the upper few kilometers is being accommodated at deeper levels in the lithosphere. Herein, the nonrift tectonic features present on Mars are briefly reviewed, along with their likely subsurface structures, and some inferences and implications are presented for behavior of the deeper lithosphere.

Golombek, M. P.↗
DOE OSTI2022

Experimental study of fracture slip due to stress perturbation in fractured geo-resources

Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.

58 GEOSCIENCES↗
DOE OSTI2021

Predicting Cement-Sheath Integrity with Consideration of Initial State of Stress and Thermoporoelastic Effects

In conventional wellbore-integrity analysis, the cement sheath’s initial state of stress and transient thermoporoelastic effects are often neglected. However, the initial state of stress is prerequisite information for accurately predicting the safe operating conditions that prevent a cemented well from being damaged. In addition, transient thermoporoelastic effects can have a profound effect on when damage will occur. Here in this paper, we propose a model that includes these effects to predict the safe operating pressures and temperatures that will prevent cement-sheath failure. For the initial state of stress, we proposed an empirical model using measurements. Subsequent stress changes are evaluated by a fully coupled transient thermoporoelastic model to analyze the mechanical behavior of the cement sheath. We predict the safe operating envelope (SOE) for shear, tensile, and debonding cement-sheath failures caused by pressure and temperature perturbations after the cement sets. Our model predicts that pore pressure is a key factor for cement failure, especially for rapid temperature changes. If the formation is low permeability, the transient pore pressures are amplified, causing the risk of damage to increase. Compared with conventional thermoelastic models, the thermoporoelastic model predicts a smaller SOE when heating the internal casing fluid and a larger envelope when cooling the internal casing fluid. Finally, the heating rate was considered with respect to field applications. The heating rate was also considered, and slower heating/cooling rates can prevent damage to the cement sheath. Finally, the thermoporoelastic model was applied to explain several laboratory and field experiments and achieved good matches.

02 PETROLEUM↗
DOE OSTI2022

Colloid-Facilitated Actinide Transport in a Cementitious-Impacted SRS Groundwater

The objective of this report is to provide a literature review relevant to the question of whether colloid facilitated transport of plutonium is enhanced in a cementitious groundwater system on the Savannah River Site (SRS). Contaminant transport is commonly described as taking place in a system with a mobile aqueous phase and an immobile solid phase. There has been an increasing awareness of a third phase, a mobile solid phase, also referred to as a mobile colloidal phase. Mobile colloids consist of organic and/or inorganic submicron-particles that move with groundwater flow. When radionuclides are associated with the mobile colloids, the net effect is that radionuclides can move faster through the subsurface system than would be predicted by transport models that do not include mobile colloids. It is important to distinguish between subsurface colloids and subsurface mobile colloids. The subsurface environment includes an enormous reservoir of colloids, but only a tiny fraction, if any, are mobile. Mobile colloid formation is commonly described as involving a three-step process: genesis, stabilization, and transport. It was concluded that there is a strong likelihood that submicron colloids of plutonium exist near the source term, (i.e., the genesis step is likely completed). While such particles have not been directly detected, it is highly likely that plutonium could either attach to submicron particles in the sediment or that submicron plutonium fragments exist at the source. However, the tendency for these plutonium colloids to move in the SRS subsurface is extremely low, especially in engineered cementitious environments where the ionic strength of the solution and the elevated concentrations of divalent cations greatly curtails colloid suspension stability (i.e., the stabilization step is likely not completed). Together these data strongly indicate that plutonium-bearing colloids would likely exist in the near field of a tank closure facility, but the colloids would be unlikely to be mobile, thereby providing a vector for enhanced transport. Under cementitious leachate impacted groundwater conditions, it is reasonable to assume plutonium transport occurs primarily as a two-phase system, a mobile aqueous phase and an immobile solid phase.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1975

Origin of Martian channels - Clathrates and water

Criticism is directed at the suggestion that Martian channels may have been eroded by liquid water produced by the depressurization of CO2 hydrate. The release of pre-existing subsurface liquid water, such as that trapped under a permafrost layer, by meteorite impact or tectonic activity could produce sufficient flow and would not require heat transfer. The presence of water in a CO2 hydrate is shown to be detrimental to its release from an underground reservoir.

Peale, S. J.↗
NASA NTRS1990

Sunspot dynamics

The goal of this research was the understanding of the various oscillatory, transient, and quasi-steady motions in sunspots and the basic structure of a sunspot. The research involved both theoretical modeling (based on thermohydrodynamic theory) and observations of dynamical phenomena in sunspots. The principal topics of the research were sunspot seismology (the interaction of solar p-modes with a sunspot as a probe of the subsurface structure of a sunspot); three minute umbral oscillations and their relation to the structure of the umbral atmosphere; siphon flows in isolated magnetic flux tubes and their relation to the photospheric Evershed flow and to intense magnetic elements outside of sunspots; and more general theoretical work on magneto-atmospheric waves. Here, a summary of results is given.

Thomas, John H.↗
NASA NTRS2022

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗