Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pattern separation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Unified and optimal frame choice for generalized parton distributions

Reconstructing the internal three-dimensional quark and gluon structures of hadrons through generalized parton distributions (GPDs) from hard exclusive scattering processes is one of the most challenging tasks in nuclear and particle physics. In this paper, we introduce a new optimized reference frame that, for the first time, enables a unified view of all the reactions sensitive to GPDs and facilitates the interpretation of a variety of phase-space patterns that were previously hardly accessible and interpretable. Similarly to how the heliocentric description advanced our understanding of the solar system and gravitation, our new frame centers around a quasireal state, allows for a consistent separation of physical scales, and reveals a novel quantum interference mechanism. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experiment and model for a Stokes layer in a strongly coupled dusty plasma

A Stokes layer, which is a flow pattern that arises in a viscous fluid adjacent to an oscillatory boundary, was observed in an experiment using a two-dimensional strongly coupled dusty plasma. Liquid conditions were maintained using laser heating, while a separate laser manipulation applied an oscillatory shear that was localized and sinusoidal. The evolution of the resulting flow was analyzed using space-time diagrams. These figures provide an intuitive visualization of a Stokes layer, including features such as the depth of penetration and wavelength. Another feature, the characteristic speed for the penetration of the oscillatory flow, also appears prominently in space-time diagrams. To model the experiment, the Maxwell-fluid model of a Stokes layer was generalized to describe a two-phase liquid. In our experiment, the phases were gas and dust, where the dust cloud was viscoelastic due to strong Coulomb coupling. Furthermore, the model is found to agree with the experiment, in the appearance of the space-time diagrams, and in the values of the characteristic speed, depth of penetration, and wavelength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient loading of reduced data ensembles produced at ORNL SNS/HFIR neutron time-of-flight facilities

We present algorithmic improvements to the loading operations of certain reduced data ensembles produced from neutron scattering experiments at Oak Ridge National Laboratory (ORNL) facilities. Ensembles from multiple measurements are required to cover a wide range of the phase space of a sample material of interest. They are stored using the standard NeXus schema on individual HDF5 files. This makes it a scalability challenge, as the number of experiments stored increases in a single ensemble file. The present work follows up on our previous efforts on data management algorithms, to address identified input output (I/O) bottlenecks in Mantid, an open-source data analysis framework used across several neutron science facilities around the world. We reuse an in-memory binary-tree metadata index that resembles data access patterns, to provide a scalable search and extraction mechanism. In addition, several memory operations are refactored and optimized for the current common use cases, ranging most frequently from 10 to 180, and up to 360 separate measurement configurations. Results from this work show consistent speed ups in wall-clock time on the Mantid LoadMD routine, ranging from 19% to 23% on average, on ORNL production computing systems. The latter depends on the complexity of the targeted instrument-specific data and the system I/O and compute variability for the shared computational resources available to users of ORNL’s Spallation Neutron Source (SNS) and the High Flux Isotope Reactor (HFIR) instruments. Nevertheless, we continue to highlight the need for more research to address reduction challenges as experimental data volumes, user time and processing costs increase.

Godoy, William↗

Additive Manufacturing of Coreless Flyback Transformers Using Aerosol Jet Printing and Electrochemical Deposition

An additive manufacturing approach combining aerosol jet printing (AJP) and electrodeposition opens a new pathway to the production of lightweight coreless flyback transformer devices for power electronics. AJP of seed layers with resolution on the order of 30 μm is combined with electrodeposition of Cu and Ni for decreased resistance. This combined approach addresses known shortcomings of AJP and electrodeposition. Nanoparticle inks used in AJP of metals have low conductivity versus bulk materials due to their high grain boundary resistance. There is a lack of readily available high-resolution patterning techniques for electrodeposition outside of expensive clean-room-based lithography techniques. Combining these two techniques enables the patterning of high-resolution, high-conductivity components. Here in this manuscript, we report on the construction of coreless flyback transformers consisting of two-layer primary and two-layer secondary spiral inductors separated by layers of a printed UV-curable dielectric. An input voltage of 17 V at 400 kHz was amplified to an output of 1250 V corresponding to a gain of 73.5. COMSOL modeling at the individual inductor level and at the transformer level was used to compare expected inductance, equivalent series resistance (ESR), and coupling with experimentally measured values.

36 MATERIALS SCIENCE↗

Quantum Mechanical Simulations of Complex Nanostructures for Photovoltaic Applications

A quantitative understanding of the electronic excitations in nanostructures, especially complex nanostructures, is crucial for making new-generation photovoltaic (PV) cells based on nanotechnology, which have high efficiency and low cost. Yet current quantum mechanical simulation methods are either computationally too expensive or not accurate and reliable enough, hindering the rational design of the nanoscale PV cells. The PI seeks to develop new methodologies to overcome the challenges in this very difficult and long-lasting problem, pushing the field forward so that electronic excitations can be accurately predicted for systems involving thousands of atoms. The primary objective of this project is to develop new approaches for electronic excitation calculations that are more accurate than traditional density functional theory (DFT) and are applicable to systems larger than what current beyond-DFT methods can treat. In this proposal, the PI will first address the excited-state problem within the DFT framework to obtain quasiparticle energies from both Kohn-Sham (KS) eigenvalues and orbitals; and the electron-hole binding energy will be computed based on screened Coulomb interaction of corresponding DFT orbitals. The accuracy of these approaches will be examined against many-body methods of GW/BSE and quantum Monte Carlo (QMC). The PI will also work on improving the accuracy and efficiency of the GW/BSE and QMC methods in electronic excitation computations by using better KS orbitals obtained from orbital-dependent DFT as inputs. Then an extended QMC database of ground- and excited-state properties will be developed, and this will be spot checked and supplemented with data from GW/BSE calculations. The investigation will subsequently focus on the development of an improved exchange-correlation (XC) density functional beyond the current generalized gradient approximation (GGA) level of parameterization, with parameters fitted to the QMC database. This will allow the ground-state properties of focus systems to be more precisely predicted using DFT. These new developments will then be applied to investigate a chosen set of complex nanostructures that have great potential for opening new routes in designing materials with improved transport, electronic, and optical properties for PV and other optoelectronic usages: (1) Hybrid interfaces between materials with distinct electronic and optical properties, such as organic molecules (conjugated polymers, e.g. P3HT) and inorganic semiconducting materials (Si and ZnO). Complicated interface structures, including interface bonding configurations, compositional and geometrical blending patterns, interfacial defects, and various sizes and shapes of inorganic nanomaterials, will be considered for the purpose of understanding the working mechanisms of present organic/nano PV systems and designing optimum interface structures for fast charge separation and injection. (2) Complex-structured semiconducting nanomaterials that could induce charge separation without pn- or hetero-junctions. The new methodology will allow the PI to investigate the performance of realistic semiconducting nanomaterials of internal (impurities, defects, etc.) and external (uneven surface, mechanical twisting and bending, surface chemistry, etc.) complexities on optical absorption and charge transport against charge trapping and recombination. Of particular interest is whether such structural complexity in a single material could even be beneficial for PV usage, for example, charge separation through morphology control. Successful completion of the proposed DFT methodology would have a far-reaching impact on our ability to study and exploit the nature of electronic excitations in complex materials, advancing the design of next-generation electronic and optoelectronic devices in all facets of renewable energy conversion and storage, including photovoltaics, thermoelectricity, photochemistry, etc.

14 SOLAR ENERGY↗

Performance Validation of Pulsed Thermal Imaging System for In-Service Applications

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of complex topology nuclear reactor parts from high-strength corrosion resistance alloys, such as stainless steel and Inconel. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which has the capability of melting metallic powder and net shaping the structures with relatively high precision. Some of the challenges with using LPBF method for nuclear manufacturing include the possibility of introducing pores into metallic structures. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures in high temperature nuclear reactor environment. Currently, there exist limited capabilities to evaluate actual AM structures non-destructively. Pulsed Thermography Imaging (PTI) provides a capability for non-destructive evaluation (NDE) of subsurface defects in arbitrary size structures. The PTI method is based on recording material surface temperature transients with infrared (IR) camera following thermal pulse delivered on material surface with flash light. The PTI method has advantages for NDE of actual AM structures because the method involves one-sided non-contact measurements and fast processing of large sample areas captured in one image. Following initial qualification of an AM component for deployment in a nuclear reactor, a PTI system can also be used for in-service nondestructive evaluation (NDE) applications. In this report, we describe recent progress in enhancing PTI capabilities in detecting microscopic defects in metallic specimens. SS316 and IN718 specimens were developed with a pattern of subsurface calibrated flat bottom hole (FBH) defects with diameters from 500µm to 200µm. FBH’s were created with EDM (electron discharge machining) drill. PTI imaging data was processed Spatial Temporal Denoised Thermal Source Separation (STDTSS) unsupervised machine learning (ML) algorithm. We show that defects as small as 200µm in SS316 and IN718 can be detected with STDTSS algorithm. To the best of our knowledge, these are the smallest detected defects which are reported in literature.

42 ENGINEERING↗

Search for Electron Neutrinos in Multiple Topologies with the MicroBooNE Experiment

This note presents the status of the measurement of electron neutrinos from the Fermilab Booster Neutrino Beamline (BNB) with the MicroBooNE experiment. The analysis is aimed at investigating the nature of the low energy excess of electromagnetic activity observed by the MiniBooNE experiment. The $ν_e$ event selection relies on topological and calorimetric information to characterise particles produced in these interactions, leveraging the Pandora multi-algorithm reconstruction framework as well as custom particle identification and pattern recognition tools. Results presented in this note use 5.88 × 10 20 protons on target of data collected between 2015 and 2018. These include the performance of particle identification tools for $\mu/p$ and $e/\gamma$ separation, along with electromagnetic shower calibration with the largest sample of $\pi^0$ events measured on argon. Progress towards the completion of this analysis is shown through the measurement of high-energy charged-current $ν_e$ interactions and their kinematic distributions using both inclusive (1eX) and exclusive (1$e$0$p$0$\pi$ and 1$e$N$p$0$\pi$) channels. Measurements of charged-current $ν_µ$ interactions aimed at constraining flux and cross-section systematic uncertainties are also shown. We present the analysis’ preliminary sensitivity of an electron-like signal hypothesis to the MiniBooNE excess which includes flux, cross-section, and detector systematic uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A visual study of radial inward choked flow of liquid nitrogen.

Data and high speed movies were acquired on pressurized subcooled liquid nitrogen flowing radially inward through a 0.0076 cm gap. The stagnation pressure ranged from 0.7 to 4 MN/sq m. Steady radial inward choked flow appears equivalent to steady choked flow through axisymmetric nozzles. Transient choked flows through the radial gap are not uniform and the discharge pattern appears as nonuniform impinging jets. The critical mass flow rate data for the transient case appear different from those for the steady case. On the mass flow rate vs pressure map, the slope and separation of the isotherms appear to be less for transient than for steady radial choked flow.

Hendricks, R. C.↗

Scale model studies for establishing the flow patterns of a low-speed tunnel

Experiments were conducted in a model tunnel scaled down in a ratio of 1:24 from the full size prototype V/STOL tunnel at NASA Langley Research Center. The flow characteristics around the tunnel were studied and the location and causes of local flow separation and local recirculation were documented. The study of the interaction between the various components from which the tunnel is built, an experimental setup similar to the sequence of tunnel circuit components was constructed which enabled the various components to be tested either individually or in combination. The components were tested by the simple technique of blowing air through them, then measuring the velocity distribution at relevant sections.

Barna, P. S.↗

Analysis of Martian terrains using optical power spectra

Planetary geological studies are almost entirely based on the analysis of orbital imagery. In the case of Mars, optical power spectra are providing the photogeologist with an additional aid in his task of classification and characterization of diverse terrains. Statistical pattern recognition techniques using optical power spectral data may be especially valuable in subdividing terrain units with characteristics that are only subtly different and in correlation of isolated patches of similar materials that are widely separated on the planet's surface.

Wolfe, R. W.↗

Simplified Learning Scheme For Analog Neural Network

Synaptic connections adjusted one at a time in small increments. Simplified gradient-descent learning scheme for electronic neural-network processor less efficient than better-known back-propagation scheme, but offers two advantages: easily implemented in circuitry because data-access circuitry separated from learning circuitry; and independence of data-access circuitry makes possible to implement feedforward as well as feedback networks, including those of multiple-attractor type. Important in such applications as recognition of patterns.

Eberhardt, Silvio P.↗

Global Geodetic Signatures of the Antarctic Ice Sheet

Four scenarios of present day Antarctic ice sheet mass change are developed from comprehensive reviews of the available glaciological and oceanographic evidence. The gridded scenarios predict widely varying contributions to secular sea level change xi ranging from -1.1 to 0.45 mm/yr, and predict polar motion m and time-varying low-degree gravitational coefficients J(sub l) that differ significantly from earlier estimates. A reasonably linear relationship between the rate of sea level change from Antarctica xi(sub A) and the predicted Antarctic J(sub l) is found for the four scenarios. This linearity permits a series of forward models to be constructed that incorporate the effects of ice mass changes in Antarctica, Greenland, and distributed smaller glaciers, as well as postglacial rebound (assuming the ICE-3G deglaciation history), with the goal of obtaining optimum reconciliation between observed constraints on J(sub l) and sea level rise xi. Numerous viable combinations of lower mantle viscosity and hydrologic sources are found that safely "observed" in the range of 1 to 2-2.5 mm/yr and observed J(sub l) for degrees 2, 3, and 4. In contrast, rates of global sea level rise above 2.5 mm/yr are inconsistent with available J(sub l) observations. The successful composite models feature a pair of lower mantle viscosity solutions arising from the sensitivity of J(sub l) to glacial rebound. The paired values are well separated at xi = 1 mm/yr, but move closer together as xi is increased, and, in fact, merge around xi = 2 - 2.5 mm/yr, revealing an intimate relation between xi and preferred lower mantle viscosity. This general pattern is quite robust and persists for different J(sub l) solutions, for variations in source assumptions, and for different styles of lower mantle viscosity stratification. Tighter J(sub l) constraints for l greater than 2 may allow some viscosity stratification schemes and source assumptions to be excluded in the future. For a given total observed xi the sea level rise from Antarctica xi(sub A) is tightly constrained and ranges from 0 to + 1 mr/yr (corresponding to an ablating ice sheet) as estimates of xi are raised from 1 to 2.5 mm/yr. However, when the degree 3 zonal harmonic constraint is removed, the solutions show little sensitivity to Antarctic mass balance, emphasizing the need for a well determined odd-degree secular zonal harmonic for determining polar ice mass balance.

James, Thomas S.↗

H-alpha Proxies for EIT Crinkles: Further Evidence for Pre-Flare "Breakout"-Type Activity in an Ejective Solar Eruption

We present Halpha observations from Big Bear Solar Observatory of an eruptive flare in NOAA AR 8210, occurring near 22:30 UT on 1998 May 1. Previously, using the EUV Imaging Telescope (EIT) on the SOHO spacecraft, we found that a pattern of transient, localized brightenings, which we call "EIT crinkles," appears in the neighborhood of the eruption near the time of flare onset. These EIT crinkles occur at a location in the active region well separated from the sheared core magnetic fields, which is where the most intense features of the eruption are concentrated. We also previously found that high-cadence images from the Soft X-ray Telescope (SXT) on Yohkoh indicate that soft X-ray intensity enhancements in the core begin after the start of the EIT crinkles. With the Halpha data, we find remote flare brightening counterparts to the EIT crinkles. Lightcurves as functions of time of various areas of the active region show that several of the remote flare brightenings undergo intensity increases prior to onset of principle brightenings in the core region, consistent with our earlier findings from EIT and SXT data. These timing relationships are consistent with the eruption onset mechanism known as the breakout model, introduced by Antiochos and colleagues, which proposes that eruptions begin with reconnection at a magnetic null high above the core region. Our observations are also consistent with other proposed mechanisms which do not involve early reconnection in the core region. As a corollary, our observations are not consistent with the so-called tether cutting models, which say that the eruption begins with reconnection in the core. The Halpha data further show that a filament in the core region becomes activated near the time of EIT crinkle onset, but little if any of the filament actually erupts, despite the presence of a halo Coronal Mass Ejection (CME) associated with this event.

Sterling, Alphonse C.↗

H(alpha) Proxies for EIT Crinkles: Further Evidence for Preflare "Breakout"-Type Activity in an Ejective Solar Eruption

We present H(alpha) observations from Big Bear Solar Observatory of an eruptive flare in NOAA Active Region 8210, occurring near 22:30 UT on 1998 May 1. Previously, using the Extreme Ultraviolet Imaging Telescope (EIT) on the Solar and Heliospheric Observatory (SOHO) spacecraft, we found that a pattern of transient, localized brightenings, which we call 'EIT crinkles,' appears in the neighborhood of the eruption near the time of flare onset. These EIT crinkles occur at a location in the active region well separated from the sheared core magnetic fields, which is where the most intense features of the eruption are concentrated. We also previously found that high-cadence images from the Soft X-ray Telescope (SXT) on Yohkoh indicate that soft X-ray intensity enhancements in the core begin after the start of the EIT crinkles. With the H(alpha) data, we find remote flare brightening counterparts to the EIT crinkles. Light curves as functions of time of various areas of the active region show that several of the remote flare brightenings undergo intensity increases prior to the onset of principal brightenings in the core region, consistent with our earlier findings from EIT and SXT data. These timing relationships are consistent with the eruption onset mechanism known as the breakout model, introduced by Antiochos and colleagues, which proposes that eruptions begin with reconnection at a magnetic null high above the core region. Our observations are also consistent with other proposed mechanisms that do not involve early reconnection in the core region. As a corollary, our observations are not consistent with the so-called tether-cutting models, which say that the eruption begins with reconnection in the core. The H(alpha) data further show that a filament in the core region becomes activated near the time of EIT crinkle onset, but little if any of the filament actually erupts, despite the presence of a halo coronal mass ejection (CME) associated with this event.

Sterling, Alphonse C.↗

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗

The Mars Science Laboratory Curiosity Rover Mastcam Instruments: Preflight and In-Flight Calibration, Validation, and Data Archiving

The NASA Curiosity rover Mast Camera (Mastcam) system is a pair of fixed-focal length, multispectral, color CCD imagers mounted approximately 2 m above the surface on the rover's remote sensing mast, along with associated electronics and an onboard calibration target. The left Mastcam (M-34) has a 34 mm focal length, an instantaneous field of view (IFOV) of 0.22 mrad, and a FOV of 20 deg × 15 deg over the full 1648 × 1200 pixel span of its Kodak KAI-2020 CCD. The right Mastcam (M-100) has a 100 mm focal length, an IFOV of 0.074 mrad, and a FOV of 6.8 deg × 5.1 deg using the same detector. The cameras are separated by 24.2 cm on the mast, allowing stereo images to be obtained at the resolution of the M-34 camera. Each camera has an eight-position filter wheel, enabling it to take Bayer pattern red, green, and blue (RGB) 'true color' images, multispectral images in nine additional bands spanning approximately 400-1100 nm, and images of the Sun in two colors through neutral density-coated filters. An associated Digital Electronics Assembly provides command and data interfaces to the rover, 8 Gb of image storage per camera, 11 bit to 8 bit companding, JPEG compression, and acquisition of high-definition video. Here we describe the preflight and in-flight calibration of Mastcam images, the ways that they are being archived in the NASA Planetary Data System, and the ways that calibration refinements are being developed as the investigation progresses on Mars. We also provide some examples of data sets and analyses that help to validate the accuracy and precision of the calibration.

MSL Mastcam; quantitative multispectral↗

Origami-Inspired Optical Shield for a Starshade Inner Disk Testbed: Design, Fabrication, and Analysis

In 2019, a 10 m-diameter starshade inner disk test article was assembled; this test article demonstrated deployment accuracy sufficient for starshade mission concepts. Here, we describe the design, fabrication, and computational structural analysis of a key inner disk component realized for this effort: the origami-folded optical shield. The optical shield is a 10 m-diameter lightweight cover that makes the inner disk opaque. It stows within a 2.3 m-diameter volume by using an origami-inspired wrapping pattern, and is deployed passively by the inner disk perimeter truss. The fold pattern was generated by a custom algorithm based one existing generative design approaches. This prototype demonstrated critical functions: stowage in a compact volume, static equilibrium when stowed and deployed, and low strain when stowed. Four separate computational structural analysis models of the optical shield were de- veloped, at varying levels of fidelity and using a variety of software solutions: Abaqus, RAPID, and ADAMS. These models were intended to pathfind approaches for modeling the stowage and deployment of the inner disk, and to demonstrate that the optical shield is amenable existing structural modeling approaches. These models were found to adequately capture pertinent stowage and deployment behavior of the optical shield prototype.

Hoffman, Laura↗

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↗