Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutral model”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Monte Carlo Modeling of Sodium in Mercury's Exosphere During the First Two MESSENGER Flybys

We present a Monte Carlo model of the distribution of neutral sodium in Mercury's exosphere and tail using data from the Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft during the first two flybys of the planet in January and September 2008. We show that the dominant source mechanism for ejecting sodium from the surface is photon-stimulated desorption (PSD) and that the desorption rate is limited by the diffusion rate of sodium from the interior of grains in the regolith to the topmost few monolayers where PSD is effective. In the absence of ion precipitation, we find that the sodium source rate is limited to approximately 10(exp 6) - 10(exp 7) per square centimeter per second, depending on the sticking efficiency of exospheric sodium that returns to the surface. The diffusion rate must be at least a factor of 5 higher in regions of ion precipitation to explain the MASCS observations during the second MESSENGER f1yby. We estimate that impact vaporization of micrometeoroids may provide up to 15% of the total sodium source rate in the regions observed. Although sputtering by precipitating ions was found not to be a significant source of sodium during the MESSENGER flybys, ion precipitation is responsible for increasing the source rate at high latitudes through ion-enhanced diffusion.

Burger, Matthew H.↗
NASA NTRS2013

A Global Atmospheric Model of Meteoric Iron

The first global model of meteoric iron in the atmosphere (WACCM-Fe) has been developed by combining three components: the Whole Atmosphere Community Climate Model (WACCM), a description of the neutral and ion-molecule chemistry of iron in the mesosphere and lower thermosphere (MLT), and a treatment of the injection of meteoric constituents into the atmosphere. The iron chemistry treats seven neutral and four ionized iron containing species with 30 neutral and ion-molecule reactions. The meteoric input function (MIF), which describes the injection of Fe as a function of height, latitude, and day, is precalculated from an astronomical model coupled to a chemical meteoric ablation model (CABMOD). This newly developed WACCM-Fe model has been evaluated against a number of available ground-based lidar observations and performs well in simulating the mesospheric atomic Fe layer. The model reproduces the strong positive correlation of temperature and Fe density around the Fe layer peak and the large anticorrelation around 100 km. The diurnal tide has a significant effect in the middle of the layer, and the model also captures well the observed seasonal variations. However, the model overestimates the peak Fe+ concentration compared with the limited rocket-borne mass spectrometer data available, although good agreement on the ion layer underside can be obtained by adjusting the rate coefficients for dissociative recombination of Fe-molecular ions with electrons. Sensitivity experiments with the same chemistry in a 1-D model are used to highlight significant remaining uncertainties in reaction rate coefficients, and to explore the dependence of the total Fe abundance on the MIF and rate of vertical transport.

meteoric inputfunction↗
NASA NTRS2016

Solar Wind Charge Exchange Contribution To The ROSAT Sky Survey Maps

DXL (Diffuse X-ray emission from the Local Galaxy) is a sounding rocket mission designed to estimate the contribution of solar wind charge eXchange (SWCX) to the diffuse X-ray background and to help determine the properties of the Local Hot Bubble. The detectors are large area thin-window proportional counters with a spectral response that is similar to that of the PSPC (Position Sensitive Proportional Counters) used in the ROSAT All Sky Survey (RASS). A direct comparison of DXL and RASS data for the same part of the sky viewed from quite different vantage points in the solar system, and the assumption of approximate isotropy for the solar wind, allowed us to quantify the SWCX contribution to all six RASS bands (R1-R7, excluding R3). We find that the SWCX contribution at l = 140 degrees, b = 0 degrees, where the DXL path crosses the Galactic plane, is 33 percent plus or minus 6 percent (statistical) plus or minus 12 percent (systematic) for R1, 44 percent plus or minus 6 percent plus or minus 5 percent for R2, 18 percent plus or minus 12 percent plus or minus 11 percent for R4, 14 percent plus or minus 11 percent plus or minus 9 percent for R5, and negligible for the R6 and R7 bands. Reliable models for the distribution of neutral H and He in the solar system permit estimation of the contribution of interplanetary SWCX emission over the the whole sky and correction of the RASS maps. We find that the average SWCX contribution in the whole sky is 26 percent plus or minus 6 percent plus or minus 13 percent for R1, 30 percent plus or minus 4 percent plus or minus 4 percent for R2, 8 percent plus or minus 5 percent plus or minus 5 percent for R4, 6 percent plus or minus 4 percent plus or minus 4 percent for R5, and negligible for R6 and R7.

sxcx↗
NASA NTRS2012

Kuang's Semi-Classical Formalism for Electron Capture Cross-Sections in Ion-Ion Collisions at Approximately to MeV/amu: Application to ENA Modeling

Recent discovery by STEREO A/B of energetic neutral hydrogen is spurring an interest and need for reliable estimates of electron capture cross sections at few MeVs per nucleon as well as for multi-electron ions. Required accuracy in such estimates necessitates detailed and involved quantum-mechanical calculations or expensive numerical simulations. For ENA modeling and similar purposes, a semi-classical approach offers a middle-ground approach. Kuang's semiclassical formalism to calculate electron-capture cross sections for single and multi-electron ions is an elegant and efficient method, but has so far been applied to limited and specific laboratory measurements and at somewhat lower energies. Our goals are to test and extend Kuang s method to all ion-atom and ion-ion collisions relevant to ENA modeling, including multi-electron ions and for K-shell to K-shell transitions.

Barghouty, A. F.↗
DOE OSTI2022

Density-matrix approach for sequential dissociative double ionization of molecules

Intense femtosecond infrared (IR) laser pulses have been used in recent years to study the breakup dynamics of molecules. However, observables such as kinetic energy release and branching ratios of molecular fragments are often difficult to predict by theory, making information of the molecular dynamics difficult to retrieve. In this work, we develop a simple model for sequential double ionization of molecules based on a density-matrix approach. The model describes tunneling ionization of the neutral and the ion as well as laser couplings between the ionic states simultaneously. Population of different doubly charged states can be obtained at a low computational cost. Furthermore, we applied our model to N 2 and obtained a good agreement on the kinetic energy release spectrum with a previous experiment. This theoretical development could open up the opportunities to use intense short IR laser pulses with coincidence measurement to probe molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1986

A method for extracting meridonal winds from ionosonde measurements by using ionospheric models

There has been great progress in modelling and measuring the dynamics of the neutral upper atmosphere in recent years. However, future progress will depend on the availability of global measurements of neutral winds. Attention is drawn to a relatively cheap means of supplementing the data base of neutral winds provided by radar and optical measurements with data obtained by ionosondes. Rishbeth in his review of F-region dynamics, derived the relationship between the height of the F2 layer and the component of the neutral wind parallel to the magnetic field of the Earth. The sensitivity is examined of the height and density of the F2 layer over Boulder, Colorado on 30 July 1982 to changes in meridional wind speed using a comprehensive interhemispheric numerical model that solves the continuity and momentum equations for H+ and O+, the energy equations for Te and Ti, and the 2-stream photoelectron equations to obtain electron heating rates. For the neutral atmosphere temperature and densities, the mass spectrometer incoherent scatter radar model of Hedin was used.

Richards, P. G.↗
NASA NTRS1994

Magnetospheric space plasma investigations

A time dependent semi-kinetic model that includes self collisions and ion-neutral collisions and chemistry was developed. Light ion outflow in the polar cap transition region was modeled and compared with data results. A model study of wave heating of O+ ions in the topside transition region was carried out using a code which does local calculations that include ion-neutral and Coulomb self collisions as well as production and loss of O+. Another project is a statistical study of hydrogen spin curve characteristics in the polar cap. A statistical study of the latitudinal distribution of core plasmas along the L=4.6 field line using DE-1/RIMS data was completed. A short paper on dual spacecraft estimates of ion temperature profiles and heat flows in the plasmasphere ionosphere system was prepared. An automated processing code was used to process RIMS data from 1981 to 1984.

Comfort, Richard H.↗
DOE OSTI2021

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A three-dimensional model of terrain-induced updrafts for movement ecology studies

Abstract Background Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. Methods In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. Results Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( $$\sigma =0.28$$ σ = 0.28 m/s), compared to 0.85 m/s ( $$\sigma =0.58$$ σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( $$\sigma =0.25$$ σ = 0.25 m/s), compared to baseline 0.54 m/s ( $$\sigma =0.45$$ σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles ( Aquila chrysaetos ) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. Conclusions We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

17 WIND ENERGY↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Far-Infrared Water Emissions from Magnetohydrodynamic Shock Waves

Nondissociative, magnetohydrodynamic, C-type shock waves are expected to be a prodigious source of far-infrared water emissions in dense interstellar regions. We have constructed a model to calculate the farinfrared H20 line spectra that emerge from such shocks. Using the best estimates currently available for the radiative cooling rate and the degree of ion-neutral coupling within the shocked gas, we modeled the temperature structure of MHD shocks using standard methods in which the charged and neutral particles are treated separately as two weakly coupled, interpenetrating fluids. Then we solved the equations of statistical equilibrium to find the populations of the lowest 179 and 170 rotational states of ortho- and para-H2O We have completed an extensive parameter study to determine the emergent H2O line luminosities as a function of preshock density in the range n(H2) equals 10(exp 4) - 10(sup 6.5)/cc and shock velocity in the range upsilon(sub s) = 5 - 40 km/ s. We find that numerous rotational transitions of water are potentially observable using the Infrared Space Observatory and the Submillimeter Wave Astronomy Satellite and may be used as diagnostics of the shocked gas. We have also computed the rotational and ro-vibrational emissions expected from H2, CO, and OH, and we discuss how complementary observations of such emissions may be used to further constrain the shock conditions. In common with previous studies, we come close to matching the observed H2, and high-J CO emissions from the Orion-KL star-forming region on the basis of a single shock model. We present our predictions for the strengths of H2O line emission from the Orion shock, and we show how our results may be scaled to other regions where molecular shocks are likely to be present.

Kaufman, Michael J.↗
DOE OSTI2023

Influence of strong Coulomb coupling on diffusion in atmospheric pressure plasmas

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion–ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions (x i ≲ 10 —6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions (10 —6 ≲ x i ≲ 10 —2 ) there is a transition from ion-neutral to ion–ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion–ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion–ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions (x i ≳ 10 —2 ), strongly correlated ion–ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. Importantly, the results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

fast neutral gas heating↗
NASA NTRS1998

Io's Plasma Environment During the Galileo Flyby: Global Three-Dimensional MHD Modeling with Adaptive Mesh Refinement

The first results for applying a three-dimensional multimedia ideal MHD model for the mass-loaded flow of Jupiter's corotating magnetospheric plasma past Io are presented. The model is able to consider simultaneously physically realistic conditions for ion mass loading, ion-neutral drag, and intrinsic magnetic field in a full global calculation without imposing artificial dissipation. Io is modeled with an extended neutral atmosphere which loads the corotating plasma torus flow with mass, momentum, and energy. The governing equations are solved using adaptive mesh refinement on an unstructured Cartesian grid using an upwind scheme for AHMED. For the work described in this paper we explored a range of models without an intrinsic magnetic field for Io. We compare our results with particle and field measurements made during the December 7, 1995, flyby of to, as published by the Galileo Orbiter experiment teams. For two extreme cases of lower boundary conditions at Io, our model can quantitatively explain the variation of density along the spacecraft trajectory and can reproduce the general appearance of the variations of magnetic field and ion pressure and temperature. The net fresh ion mass-loading rates are in the range of approximately 300-650 kg/s, and equivalent charge exchange mass-loading rates are in the range approximately 540-1150 kg/s in the vicinity of Io.

Combi, M. R.↗
DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗
NASA NTRS1986

Two-dimensional quasi-neutral description of particles and fields above discrete auroral arcs

Models are presented for particle distributions, electric fields and currents in an adiabatic treatment of auroral electrostatic potential distributions in order to describe the quiet-time evening auroral arcs featuring both upward and return currents. The models are consistent with current continuity and charge balance requirements for particle populations controlled by adiabatic invariants and quasi-neutrality in the magnetosphere. The effective energy of the cool electron population is demonstrated to have a significant effect on the latitudinal breadth of the auroral electrostatic potential structure and the extent of the penetration of the accelerating potential into the ionosphere. Another finding is that the energy of any parallel potential drop in the lowest few thousand kilometers of the field line is of the same order of magnitude as the thermal energy of the cool electrons. Additional predictions include density cavities along field lines that support large potential drops, and density enhancements along field lines at the edge of an inverted V with a small potential drop.

Newman, A. L.↗
NASA NTRS1989

PIC code modeling of spacecraft charging potential during electron beam injection into a background of neutral gas and plasma, part 1

Injections of nonrelativistic electron beams from an isolated equipotential conductor into a uniform background of plasma and neutral gas were simulated using a 2-D electrostatic particle code. The ionization effects on spacecraft charging are examined by including interactions of electrons with neutral gas. The simulations show that the conductor charging potential decreases with increasing neutral background density due to the production of secondary electrons near the conductor surface. In the spacecraft wake, the background electrons accelerated towards the charged spacecraft produce an enhancement of secondary electrons and ions. Simulations run for longer times indicate that the spacecraft potential is further reduced and short wavelength beam-plasma oscillations appear. The results are applied to explain the spacecraft charging potential measured during the SEPAC experiments from Spacelab 1.

Koga, J. K.↗
NASA NTRS2004

Remote Sensing of Icy Galilean Moon Surface and Atmospheric Composition Using Low Energy (1 eV-4 keV) Neutral Atom Imaging

We describe a low energy neutral atom imager suitable for composition measurements Europa and other icy Galilean moons in the Jovian magnetosphere. This instrument employs conversion surface technology and is sensitive to either neutrals converted to negative ions, neutrals converted to positive ions and the positive ions themselves depending on the power supply. On a mission such as the Jupiter Icy Moons Orbiter (JIMO), two back-to-back sensors would be flown with separate power supplies fitted to the neutral atom and iodneutral atom sides. This will allow both remote imaging of 1 eV < E < 4 keV neutrals from icy moon surfaces and atmospheres, and in situ measurements of ions at similar energies in the moon ionospheres and Jovian magnetospheric plasma. The instrument provides composition measurements of the neutrals and ions that enter the spectrometer with a mass resolution dependent on the time-of-flight subsystem and capable of resolving molecules. The lower energy neutrals, up to tens of eV, arise from atoms and molecules sputtered off the moon surfaces and out of the moon atmospheres by impacts of more energetic (keV to MeV) ions from the magnetosphere. Direct Simulation Monte Carlo (DSMC) models are used to convert measured neutral abundances to compositional distributions of primary and trace species in the sputtered surfaces and atmospheres. The escaping neutrals can also be detected as ions after photo- or plasma-ionization and pickup. Higher energy, keV neutrals come from charge exchange of magnetospheric ions in the moon atmospheres and provide information on atmospheric structure. At the jovicentric orbits of the icy moons the presence of toroidal gas clouds, as detected at Europa's orbit, provide M e r opportunities to analyze both the composition of neutrals and ions originating from the moon surfaces, and the characteristics of magnetospheric ions interacting with neutral cloud material. Charge exchange of low energy ions near the moons, and directional distributions of the resultant neutrals, allow indirect global mapping of magnetic field structures around the moons. Temporal variation of the magnetic structures can be linked to induced magnetic fields associated with subsurface oceans.

Collier, M. R.↗