Search NASA⌕ Search

SEARCH · Search NASA

Results for “column generation”

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 91 records · Page 5

Observing System Experiments for Evaluating the Impact of Satellite Sea Surface Salinity on Seasonal Predictions from the GMAO S2S System

We assess the impact of satellite sea surface salinity (SSS) observations on dynamical ENSO (El Nino-Southern Oscillation) forecasts. Assimilation of SSS improves the mixed layer depth (MLD) and modulates the Kelvin waves associated with ENSO. In column 2, the initialization differences between experiments that assimilate SSS minus those withholding SSS assimilation are presented. Column 3 shows examples of forecasts generated for the different phases of ENSO assimilating the different satellite SSS. In general, for all phases of ENSO, SSS assimilation improves forecasts. The far right column compares ensemble means for assimilation of individual and combined SMOS, Aquarius, SMAP (Soil Moisture and Ocean Salinity, NASA Aquarius, Soil Moisture Active Passive) SSS forecasts. Finally, the latest forecasts are presented comparing assimilation versus no- assimilation of satellite SSS for single forecasts over the last year.

Hackert, E.↗

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the NASA/GMAO Seasonal Forecast System

We assess the impact of satellite sea surface salinity (SSS) observations on dynamical ENSO forecasts. Assimilation of SSS improves the mixed layer depth (MLD) and modulates the Kelvin waves associated with ENSO. In column 2, the initialization differences between experiments that assimilate SSS minus those withholding SSS assimilation are presented. Column 3 shows examples of forecasts generated for the different phases of ENSO assimilating the different satellite SSS. In general, for all phases of ENSO, SSS assimilation improves forecasts. The far right column compares ensemble means for assimilation of individual and combined SMOS, Aquarius, and SMAP SSS forecasts. Finally, the latest forecasts are presented comparing assimilation versus no- assimilation of satellite SSS for single forecasts over the last year.

Hackert, E.↗

A comparison of coupling efficiencies for a Stix coil and an m equals 1 coil

This theoretical and experimental study compares the ion-cyclotron wave generating characteristics of a Stix coil (which generates waves with azimuthal mode number m = 0) with those of a coil which produces primarily m = + or -1 ion-cyclotron modes. The theoretical work of J.E. Hipp, which predicted very good coupling for the m = 1 coil, was extended to determine the scaling laws for plasma column radius and coil wavelength. Experimentally, an m = 1 coil and an m = 0 coil were used to generate ion-cyclotron waves on a beam generated plasma column with electron density = 10 to the 12th power/cu cm. Coupling resonances with peak efficiencies of approximately 40 to 50 percent were measured for both coils in low power (approximately 10k W) experiments. For equal power transfer to the plasma, the m = 0 coil voltage was more than a factor of two greater than that for the m = 1 coil.

Sigman, D. R.↗

Stiffness coupling application to modal synthesis program, users guide

A FORTRAN IV computer program used to perform modal synthesis of structures by stiffness coupling, using the dynamic transformation method is described. The program was named SCAMP (Stiffness Coupling Approach Modal-Synthesis Program). The program begins with the entry of a substructure's physical mode shapes and eigenvalues or a substructure's mass and stiffness matrix. If the mass and stiffness matrices are entered, the eigen problem for the individual substructure is solved. Provisions are included for a maximum of 20 substructures which are coupled by stiffness matrix springs. Each substructure has a number degrees of freedom (DOF), except that for DOF greater than 100; vector sets having maximum row and column size of 100 were generated prior to entering SCAMP. The substructures are then coupled together via coupling springs, and the dynamic transformation is used to reduce the size of the eigen problem.

Kuhar, E. J.↗

On the thermal structure of Triton's thermosphere

The analysis of the Voyager 2 Ultraviolet Spectrometer solar occultation data obtained at Triton is consistent with a spherically symmetric, isothermal thermosphere above 400 km at T(infinity) = 96 K. A detailed calculation of energy loss processes in a pure N2 atmosphere indicates that solar heating, with calculated T(infinity) = 70 K, is insufficient to account for the inferred T(infinity) = 96 K. The magnetosphere must deposit twice as much power as the sun to heat the thermosphere to 96 K and generate the observed N2 tangential column densities above 450 km. The thermal escape of H and N atoms and the downward diffusion of N atoms to recombine below 130 km results in local ionospheric heating efficiency of 24 percent. An upper limit on the tropopause CO mixing ratio of 2 x 10 exp -4 is inferred in the absence of aerosol heating to balance its efficient cooling by LTE rotational line emission.

Stevens, Michael H.↗

Time-Dependent Photodissociation Regions

We present theoretical models of the time-dependent thermal and chemical structure of molecular gas suddenly exposed to far-ultraviolet (FUV) (6 eV less than hv less than 13.6 eV) radiation fields and the consequent time- dependent infrared emission of the gas. We focus on the response of molecular hydrogen for cloud densities ranging from n = 10(exp 3) to 10(exp 6)/cu cm and FUV fluxes G(sub 0) = 10(exp 3)-10(exp 6) times the local FUV interstellar flux. For G(sub 0)/n greater than 10(exp -2) cu cm, the emergent H(sub 2) vibrational line intensities are initially larger than the final equilibrium values. The H(sub 2) lines are excited by FUV fluorescence and by collisional excitation in warm gas. Most of the H(sub 2) intensity is generated at a characteristic hydrogen column density of N approximately 10(exp 21)/sq cm, which corresponds to an FUV optical depth of unity caused by dust opacity. The time dependence of the H(sub 2) intensities arises because the initial abundances of H(sub 2) at these depths is much higher than the equilibrium values, so that H(sub 2) initially competes more effectively with dust in absorbing FUV photons. Considerable column densities of warm (T approximately 1000) K H(sub 2) gas can be produced by the FUV pumping of H(sub 2) vibrational levels followed by collisional de-excitation, which transfers the energy to heat. In dense (n greater than or approximately 10(exp 5)/cu cm) gas exposed to high (G(sub 0) greater than or approximately 10(exp 4)) fluxes, this warm gas produces a 2-1 S(1)/1-0 S(l) H(sub 2) line ratio of approximately 0.1, which mimics the ratio found in shocked gas. In lower density regions, the FUV pumping produces a pure-fluorescent ratio of approximately 0.5. We also present calculations of the time dependence of the atomic hydrogen column densities and of the intensities of 0 I 6300 A, S II 6730 A, Fe II 1.64 microns, and rotational OH and H20 emission. Potential applications include star-forming regions, clouds near active galactic nuclei, and planetary nebulae. We apply our models to five planetary nebulae and conclude that only BD +30deg3639 shows evidence of enhanced H(sub 2) emission due to (high) nonequilibrium H(sub 2) abundances.

Hollenbach, David↗

The Characterization of a DC Glow Discharge in a Tripped Hypersonic Boundary Layer

A direct current glow discharge plasma was generated on a canonical 2.75° half-angle wedge test article. Experiments were conducted in the Texas A&M University Actively Controlled Expansion tunnel at 𝑀 = 5.7 and 𝑅𝑒 = 6 × 10 6 /m. The effect of Reynolds number, polarity, and upstream perturbation were all independently studied via imaging, power measurements, and optical emission spectroscopy (OES). These measurements are the first to be conducted in a tripped hypersonic boundary layer and allow a deeper understanding of the electrohydrodynamics. The plasma, formally characterized as a normal glow discharge, had negligible Joule and cathode heating effects on the boundary layer. Analogously, the installation of trips upstream of the electrodes changed only the appearance of the plasma, creating periodic streaks corresponding to the wakes and vortices produced by each trip but not altering the power across the electrodes. The plasma power peaked at 47 W, which OES revealed was sufficient to generate NO in the positive column and negative glow.

Casey J. Broslawski↗

The Characterization of a DC Glow Discharge in a Tripped Hypersonic Boundary Layer

A direct current glow discharge plasma was generated on a canonical 2.75° half-angle wedge test article. Experiments were conducted in the Texas A&M University Actively Controlled Expansion tunnel at 𝑀 = 5.7 and 𝑅𝑒 = 6 × 10 6 /m. The effect of Reynolds number, polarity, and upstream perturbation were all independently studied via imaging, power measurements, and optical emission spectroscopy (OES). These measurements are the first to be conducted in a tripped hypersonic boundary layer and allow a deeper understanding of the electrohydrodynamics. The plasma, formally characterized as a normal glow discharge, had negligible Joule and cathode heating effects on the boundary layer. Analogously, the installation of trips upstream of the electrodes changed only the appearance of the plasma, creating periodic streaks corresponding to the wakes and vortices produced by each trip but not altering the power across the electrodes. The plasma power peaked at 47 W, which OES revealed was sufficient to generate NO in the positive column and negative glow.

Casey Broslawski↗

High accuracy in-situ direct gas analysis of Li-ion batteries

Cycling of lithium-ion batteries containing Ni-rich NMC cathodes at high voltage involves intense gas generation. From the safety standpoint, it is critical to understand how different gas species respond to changes of upper cut-off voltages. In this manuscript, we introduce a novel experimental set up for real-time analysis of gas generation in prismatic pouch cells. In a typical experiment, a lithium-ion pouch cell is directly connected to a quadruple mass spectrometer by glass capillary. Pressure difference helps move the generated gases to the mass spectrometer column for full analysis. The gaseous species are probed during both formation cycle and aging cycles in Li-ion pouch cells comprising NMC811 cathode and graphite anode. The gases are generated upon the creation of the solid electrolyte interphase on graphite during the first formation charge. Besides ethylene gas, C 2 H 4 , whose occurrence is related to the decomposition of ethylene carbonate, a minor solvent decomposition happened at a lower voltage is demonstrated to be associated with residue water on the electrode. For aging cycles, three upper cut-off voltages have been selected to investigate the gas evolution behavior upon charging to 4.2, 4.4 and 4.6 V. Higher upper cut-off voltages do not affect the amount of O 2 measured as a part of the generated oxygen could have reacted with electrolytes. Quite similar and low amounts of O 2 remain in the cells upon charging to the selected cut-off voltages. On the other hand, the generation of CO 2 has been found to be very sensitive to upper cut-off voltage. During similar aging cycles, 6167 nmol $g_{NMC}^{-1}$ of CO 2 is generated at 4.6 V, versus 1650 nmol $g_{NMC}^{-1}$ of CO 2 at 4.4 V, versus 91 nmol $g_{NMC}^{-1}$ of CO 2 at 4.2 V.

25 ENERGY STORAGE↗

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗

GOES weather satellite observations and measurements of the May 18, 1980, Mount St. Helens eruption

We demonstrate the use of Geostationary Operational Environmental Satellite (GOES) images of the May 18, 1980, Mount St. Helens volcanic plume in providing details of the dynamics and changing character of this major explosive eruption. Visible and thermal infrared (IR) data from a sequence of images at 30-min intervals from 0850 to 1720 Local Time (LT) give information on dispersal and plume top temperature. Initial visible and IR images at 0850 show the top of a spreading co-ignimbrite-like umbrella plume and an overshooting column emerging from it, both rising off the ground-hugging pyroclastic gravity flow generated by the opening directed blast. The overshooting column had a minimum temperature significantly colder than local ambient atmosphere, indicating substantial undercooling, and a maximum altitude of 31 +/- 2 km at 0920. This large plume system then formed a high-velocity, radially spreading, gravitationally driven current before becoming advected in the wind field at an average downwind velocity of 29 m/s. Reflectance values from visible GOES data change from lower to higher during periods of transition from darker toned Plinian to lighter toned co-ignimbrite plumes indicating that in this case satellite data resolved changes in eruptive style from plumes with a coarser to a finer dominant particle size.

Holasek, R. E.↗

Creation of Synthetic Electric Grids (SPP/MISO) Supporting PERFORM (Final Report)

Over the course of the project, two “realistic but not real” synthetic transmission-level grid models over the SPP-MISO and ERCOT footprints were created to provide more realistic data and increase the reliability and resiliency of the grids under a variety of scenarios. The synthetic ERCOT transmission grid is compatible with the distribution grid developed in collaboration with NREL. All generators are based on the EIA 860 data and a column with EIA plant code and Gen ID is added to generators of both grids so that they can be easily mapped. The improvements are also made to electric grids including N-1 contingencies with some remedial actions, improving the transmission lines to avoid lines in lakes, including an HVDC line to the SPP-MISO case, providing several generator parameters and their temporal constraints that were not included in EIA 860 form, generators’ cost curves, load offer curves, adding phase shifters and tap changers with impedance correction tables, adding reactive power control and partitioning the grids into active and reactive reserve zones and determine different types of the required reserve for each zone. Hourly load time series at the bus level were generated to create scenarios for solving power flow in different loading conditions. Weather measurement information and the models of renewable generators are used to directly include the impact of weather on the grids. Based on a variety of load and weather conditions the grids are improved to accommodate different conditions. The ERCOT 7k-bus grids were also modeled for the year 2030 with predicted improvements in renewable resources. The renewable generation model was also improved with historic weather data included. The impact of electric vehicles on the ERCOT grid is also modeled.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Recovery of high specific activity molybdenum-99 from accelerator-induced fission on low-enriched uranium for technetium-99m generators

Abstract A new process was developed to recover high specific activity (no carrier added) 99 Mo from electron-accelerator irradiated U 3 O 8 or uranyl sulfate targets. The process leverages a novel solvent extraction scheme to recover Mo using di(2-ethylhexyl) phosphoric acid following uranium and transuranics removal with tri- n -butyl phosphate. An anion-exchange concentration column step provides a final purification, generating pure 99 Mo intended for making 99 Mo/ 99m Tc generators. The process was demonstrated with irradiated uranium targets resulting in more than 95% 99 Mo recovery and without presence of fission products or actinides in the product.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

CALIPSO’s New ODCOD Product Compared to Langley’s HSRL and HSRL-2

CALIPSO’s Version 4.51 Level 2 data release introduces an all-new group of science data sets containing estimates of total column two-way transmittances and effective optical depths derived from CALIOP ocean surface backscatter measurements and MERRA-2 reanalysis wind speed data. This new data product is known as the Ocean Derived Column Optical Depth product (ODCOD). These estimates use CALIOP’s ocean surface lidar return to estimate the total column optical depth overhead. This new retrieval will be provided for the entire mission, day and night, at 532nm, and are reported at single shot, 1km, and 5km resolutions for all profiles in which a valid lidar ocean surface return is detected. This poster will present a collection of HSRL and HSRL-2 CALIPSO under flights and compare the HSRL measured aerosol optical depth with the estimated total column optical depth of ODCOD. The ODCOD algorithm uses CALIOP’s ocean surface lidar return to estimate the effective total column optical depth overhead; i.e., the single scattering optical depth reduced by the contributions of multiple scattering to the ocean surface signal. Because the retrieval generates a single estimate for the entire column based solely on the surface backscatter, there is currently no way to correct for multiple scattering effects from different kinds of features within the column. This new retrieval will be provided for the entire mission, day and night, at 532nm, and is reported at single shot, 1km, and 5km resolutions for all profiles in which a valid lidar ocean surface return is detected. This poster presents a collection of HSRL and HSRL-2 CALIPSO under-flight case studies that compare the HSRL measured aerosol optical depth with the estimated total column optical depth derived by ODCOD.

R. Ryan↗

The Earth Model Column Collaboratory (EMC2) v1.1: An Open-Source Ground-Based Lidar and Radar Instrument Simulator and Subcolumn Generator for Large-Scale Models

Climate models are essential for our comprehensive understanding of Earth's atmosphere and can provide critical insights on future changes decades ahead. Because of these critical roles, today's climate models are continuously being developed and evaluated using constraining observations and measurements obtained by satellites, airborne, and ground-based instruments. Instrument simulators can provide a bridge between the measured or retrieved quantities and their sampling in models and field observations while considering instrument sensitivity limitations. Here we present the Earth Model Column Collaboratory (EMC2), an open-source ground-based lidar and radar instrument simulator and subcolumn generator, specifically designed for large-scale models, in particular climate models, but also applicable to high-resolution model output. EMC2 provides a flexible framework enabling direct comparison of model output with ground-based observations, including generation of subcolumns that may statistically represent finer model spatial resolutions. In addition, EMC2 emulates ground-based (and air- or space-borne) measurements while remaining faithful to large-scale models' physical assumptions implemented in their cloud or radiation schemes. The simulator uses either single particle or bulk particle size distribution lookup tables, depending on the selected scheme approach, to perform the forward calculations. To facilitate model evaluation, EMC2 also includes three hydrometeor classification methods, namely, radar- and sounding-based cloud and precipitation detection and classification, lidar-based phase classification, and a Cloud Feedback Model Intercomparison Project Observational Simulator Package (COSP) lidar simulator emulator. The software is written in Python, is easy to use, and can be straightforwardly customized for different models, radars, and lidars. Following the description of the logic, functionality, features, and software structure of EMC2, we present a case study of highly supercooled mixed-phase cloud based on measurements from the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE). We compare observations with the application of EMC2 to outputs from four configurations of the NASA Goddard Institute for Space Studies (GISS) climate model (ModelE3) in single-column model (SCM) mode and from a large-eddy simulation (LES) model. We show that two of the four ModelE3 configurations can form and maintain highly supercooled precipitating cloud for several hours, consistent with observations and LES. While our focus is on one of these ModelE3 configurations, which performed slightly better in this case study, both of these configurations and the LES results post-processed with EMC2 generally provide reasonable agreement with observed lidar and radar variables. As briefly demonstrated here, EMC2 can provide a lightweight and flexible framework for comparing the results of both large-scale and high-resolution models directly with observations, with relatively little overhead and multiple options for achieving consistency with model microphysical or radiation scheme physics.

Earth Model Column Collaboratory↗

Iterated Gauss-Seidel GMRES

The GMRES algorithm of Saad and Schultz [SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869] is an iterative method for approximately solving linear systems Ax = b, with initial guess x0 and residual r0 = b Ax0. The algorithm employs the Arnoldi process to generate the Krylov basis vectors (the columns of Vk ). It is well known that this process can be viewed as a QR factorization of the matrix Bk = [r0, AVk] at each iteration. Despite an O (..epsilon..)..kappa.. (Bk ) loss of orthogonality, for unit roundoff ..epsilon..and condition number ..kappa.. , the modified Gram-Schmidt formulation was shown to be backward stable in the seminal paper by Paige et al. [SIAM J. Matrix Anal.Appl., 28 (2006), pp. 264-284]. We present an iterated Gauss-Seidel formulation of the GMRES algorithm (IGS-GMRES) based on the ideas of Ruhe [Linear Algebra Appl., 52 (1983), pp. 591-601] and Swirydowicz et al. [Numer. Linear Algebra Appl., 28 (2020), pp. 1-20]. IGS-GMRES maintains orthogonality to the level O (..epsilon..)..kappa.. (Bk ) or O (..epsilon..), depending on the choice of one or two iterations; for two Gauss-Seidel iterations, the computed Krylov basis vectors remain orthogonal to working accuracy and the smallest singular value of Vk remains close to one. The resulting GMRES method is thus backward stable. We show that IGS-GMRES can be implemented with only a single synchronization point per iteration, making it relevant to large-scale parallel computing environments. We also demonstrate that, unlike MGS-GMRES, in IGS-GMRES the relative Arnoldi residual corresponding to the computed approximate solution no longer stagnates above machine precision even for highly nonnormal systems.

Arnoldi-QR↗

Results of induced atmosphere measurements from the Apollo program

Experiments on Apollo missions 15, 16, and 17 were utilized in an attempt to learn about the induced environment in the vicinity of manned spacecraft. Photographic sequences were examined to obtain scattered light data from the spacecraft-generated particulates during quiescence periods and after liquid dumps. The results allowed estimates of the obscuration factor and the clearing times after dumps. It was found that the clearing times were substantially longer than anticipated. The mass spectrometer detected a high molecular flux in lunar orbit which was induced by the spacecraft. It is shown that this is most likely caused by small ice crystals being continually produced in lunar orbit. Other data from the ultraviolet spectrometer and the stellar camera are also analyzed, and estimated values or upper limits are placed on the total scattering background, the size and number of particles generated, the velocity range, and the column density.

Naumann, R. J.↗

Banded surface flow maintained by convection in a model of the rapidly rotating giant planets

In three-dimensional numerical simulations of a rapidly rotating Boussinesq fluid shell, thermally driven convection in the form of columns parallel to the rotation axis generates an alternately directed mean zonal flow with a cylindrical structure. The mean structure at the outer spherical surface consists of a broad eastward flow at the equator and alternating bands of westward and eastward flows at higher latitudes in both hemispheres. The banded structure persists even though the underlying convective motions are time-dependent. These results, although still far from the actual motions seen on Jupiter and Saturn, provide support for theoretical suggestions that thermal convection can account for the remarkable banded flow structures on these planets.

Sun, Zi-Ping↗