Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tidal Prediction”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Compact X-ray binaries in hierarchical triples. I - Tidal angular momentum loss and GX 17 + 2

A general formalism is developed for the enhanced mass transfer rate expected in a close binary with a (distant) third companion star. Such a hierarchical triple model is applied to the X-ray burster and QPO source GX 17 + 2, consisting of an inner mass-transferring binary comprising a main-sequence secondary and an accreting neutron star, and a more distant companion. The model is shown to account for the anomalously high mass transfer rate for this burster and other systems with short orbital periods. A G star, which does not appear to be the binary companion of the neutron star but is consistent with the sub-arc second radio error box for GX 17 + 2, may be the triple companion. A predicted velocity curve for the triple companion is presented.

Bailyn, Charles D.↗
DOE OSTI2025

The GD-1 Stellar Stream Perturber as a Core-collapsed Self-interacting Dark Matter Halo

The GD-1 stellar stream exhibits spur and gap structures that may result from a close encounter with a dense substructure. When interpreted as a dark matter subhalo, the perturber is denser than predicted in the standard cold dark matter (CDM) model. In self-interacting dark matter (SIDM), however, a halo could evolve into a phase of gravothermal collapse, resulting in a higher central density than its CDM counterpart. We conduct high-resolution controlled N-body simulations to show that a collapsed SIDM halo could account for the GD-1 perturber's high density. We model a progenitor halo with a mass of 3 × 10 8 M ⊙ , motivated by a cosmological simulation of a Milky Way analog, and evolve it in the Milky Way's tidal field. For a cross section per mass of σ/m ≈ 30–100 cm 2 g −1 at ${V}_{{\rm{\max }}}\unicode{x0007E}10\,{\rm{km}}\,{{\rm{s}}}^{-1}$, the enclosed mass of the SIDM halo within the inner 10 pc can be increased by more than 1 order of magnitude compared to its CDM counterpart, leading to a good agreement with the properties of the GD-1 perturber. Our findings indicate that stellar streams provide a novel probe into the self-interacting nature of dark matter.

dark matter↗
DOE OSTI2026

Revisiting the validity of eddy viscosity models for predicting airflow over water waves

In this study, we revisit the validity of eddy viscosity models for predicting wave-induced airflow disturbances over ocean surface waves. We first derive a turbulence curvilinear model for the phase-averaged Navier–Stokes equations, extending the work of Cao, Deng & Shen (2020 J. Fluid Mech. 901, A27), by incorporating turbulence stress terms previously neglected in the linearised viscous curvilinear model. To verify our formulation, we perform a priori tests by numerically solving the model using mean wind and turbulence stress profiles from large-eddy simulations (LES) of airflow over waves across various wave ages. Results show that including turbulence stress terms improves wave-induced airflow predictions compared with the previous viscous curvilinear model. We further show that using a standard mixing-length eddy viscosity yields inaccurate predictions at certain wave ages, as it fails to capture wave-induced turbulence, which fundamentally differs from mean shear-driven turbulence. The LES data show that accurate representations of wave-induced stresses require a complex-valued eddy viscosity. The maximum magnitude of this eddy viscosity scales as ∼𝑢 𝜏 ⁢𝜁 𝑖𝑛𝑛𝑒𝑟 , where 𝑢 𝜏 is the friction velocity and 𝜁 𝑖𝑛𝑛𝑒𝑟 is the inner-layer thickness, the height at which the eddy-turnover time matches the wave advection time scale. This scaling aligns with the prediction by Belcher & Hunt (1993 J. Fluid Mech. 251, 109–148). Overall, the findings demonstrate that traditional eddy viscosity models are inadequate for capturing wave-induced turbulence. More sophisticated turbulence models are essential for the accurate prediction of airflow disturbances and form drag in wind–wave interaction models.

16 TIDAL AND WAVE POWER↗
NASA NTRS1996

Dynamics of the Impacts

The SL9 impacts are best known by their plumes, several of which were imaged towering 3000 km above Jupiter's limb. The heat released when they fen back on the atmosphere produced the famous infrared main events. The reentry shock was observed directly in CO emission and indirectly through thermally glowing dust. The relatively low color temperature of the latter implies that most of the energy released on reentry was radiated. If so, the R impact released some 3 x 10(exp 26) ergs, and the larger L, K, and G impacts as much as 3 x 10(exp 27) ergs. These estimates agree well with impact energies deduced from tidal disruption computations. The duration of the infrared events measures time aloft and hence provides a second measure of plume height. These indicate that plume height was only weakly dependent on impact energy. Evidently all the plumes were launched at roughly 10-13 km/s. Using a semi-analytic model for the deceleration, disintegration, and destruction of intruding bodies by an ever-vigilant atmosphere, we find that similar plume heights is a direct consequence of smaller impactors exploding at higher altitudes, in such a way that the different explosions were geometrically similar. We then compare the predictions of our model to the Venerian cratering record, which provides a good statistical test of impact deceleration by a thick atmosphere. Chemistry should have provided an independent measure of explosion altitude: abundant shock-generated CS, CS2 and HCN indicates a source in dry jovian air, above the putative water clouds. However, the Galileo Probe results seem to imply that we should expect no more. Observed water and S2 are consistent with a somewhat oxidized gas (presumably the comet itself), but the absence of SO2 and CO2 shows that conditions were neither too oxidizing nor the shocks too hot.

Zahnle, Kevin↗
DOE OSTI2024

TEAMER - Extreme Events Modeling for the MARMOK-OWC Wave Energy Converter

Through the TEAMER program, Sandia National Laboratories (SNL) collaborated with IDOM Incorporated to study their MARMOK-Oscillating Water Column (MARMOK-OWC) wave energy conversion device. The study yielded a quantitative understanding of hydrodynamic pressures on the oscillating water column (OWC) device surfaces, the mooring tensions, and the dynamic performance of the device under extreme ocean wave conditions. This project utilized a comprehensive multi-phase Navier-Stokes flow solver with an overset body-fit mesh to predict fluid velocities and hydrodynamic forces on the MARMOK-OWC device. Computational Fluid Dynamics (CFD) analysis were conducted using OpenFOAM. This data includes the OpenFOAM cases (setup and data) to run the extreme events developed during the project. This project is part of the TEAMER RFTS 4 (request for technical support) program.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Wave

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production using a single, simulated wave energy conversion device. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel Hydrogen. While the unit supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. For the wave energy, NLR used a wave energy converter model from PacWave. These devices can be equipped with accumulators and pressure relief values to smooth the power output by storing and releasing hydraulic energy. Using a peak power output of 10 MW, the model created two 25-minute profiles: one with and one without the accumulators and pressure relief valves. To down select the profile data from the native resolution of 20 Hz to 1 Hz, NLR took the mean of every 20 data points. NLR experimented with two simulated wave energy power plants: one that peaks at 10 MW, and one that peaks at 5 MW. These profiles were scaled for the physical 1.25 MW electrolyzer by multiplying the original profiles by one eighth and one quarter, respectively. The first profile matches the capacity rating of eight of the 1.25 MW electrolyzers, while the second matches four electrolyzers. Finally, NLR experimented with two settings for the electrolyzer power supply minimum and maximum current ramp rates (gain and slew): 200 and 400 amperes per second. The simulated profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wave electrolysis experiment and is formatted as follows: {technology}-{accumulator?}_{number of 1.25 MW electrolyzers connected}-{electrolyzer ramp rate in amperes/second} For instance, “wavePacWave-Noacc_4-400.zip” represents the 25 minute-long experiment using the PacWave’s wave energy converter model, equipped with no accumulator, connected to four 1.25-MW electrolyzers with their power supplies set to a maximum current ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. An experiment, labeled “characterization_200.zip”, demonstrates the MC250 electrolyzer steady-state response with 30 minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all wave profiles combined into one dataset labeled "combined_wave_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis.

08 HYDROGEN↗
NASA NTRS1997

Paleobiology of the Mesoproterozoic-Neoproterozoic transition: the Sukhaya Tunguska Formation, Turukhansk Uplift, Siberia

Silicified carbonates of the latest Mesoproterozoic Sukhaya Tunguska Formation, northwestern Siberia, contain abundant and diverse permineralized microfossils. Peritidal environments are dominated by microbial mats built by filamentous cyanobacteria comparable to modern species of Lyngbya and Phormidium. In subtidal to lower intertidal settings, mat-dwelling microbenthos and possible coastal microplankton are abundant. In contrast, densely woven mat populations with few associated taxa characterize more restricted parts of tidal flats; the preservation of vertically oriented sheath bundles and primary fenestrae indicates that in these mats carbonate cementation was commonly penecontemporaneous with mat growth. Eoentophysalis mats are limited to restricted environments where microlaminated carbonate precipitates formed on or just beneath the sediment surface. Most microbenthic populations are cyanobacterial, although eukaryotic microfossils may occur among the simple spheroidal cells interpreted as coastal plankton. Protists are more securely represented by large (up to 320 micrometers in diameter) but poorly preserved acritarchs in basinal facies. The Sukhaya Tunguska assemblage contains 27 species in 18 genera. By virtue of their stratigraphic longevity and their close and predictable association with specific paleoenvironmental conditions, including substrates, Proterozoic cyanobacteria support a model of bacterial evolution in which populations adapt rapidly to novel environments and, thereafter, resist competitive replacement. The resulting evolutionary pattern is one of accumulation and stasis rather than the turnover and replacement characteristic of Phanerozoic plants and animals.

NASA Discipline Exobiology↗
DOE OSTI2025

ELM–Wet: Inclusion of a Wet–Landunit With Sub–Grid Representation of Eco–Hydrological Patches and Hydrological Forcing Improves Methane Emission Estimations in the E3SM Land Model (ELM)

Wetlands are the largest emitters of biogenic methane (CH 4 ) and represent the highest source of uncertainty in global CH 4 budgets. Here, we aim to improve the realism of wetland representation in the U.S. Department of Energy's Exascale Earth System Model land surface model, ELM, thereby reducing uncertainty of CH 4 flux predictions. We develop an updated version, ELM-Wet, where we activate a separate landunit for wetlands that handles multiple wetland-specific eco-hydrological patch functional types. We introduce more realistic hydrological forcing through prescribing site-level constraints on surface water elevation, which allows resolving different sustained inundation depth for different patches, and if data exists, prescribing inundation depth. We modified the calculation of aerenchyma transport diffusivity based on observed conductance per leaf area for different vegetation types. We use Bayesian Optimization to parameterize CO 2 and CH 4 fluxes in the developed wet-landunit. Site-level simulations of a coastal non-tidal freshwater wetland in Louisiana were performed with the updated model. Eddy covariance observations of CO 2 and CH 4 fluxes from 2012 to 2013 were used to train the model and data from 2021 were used for validation. Patch-specific chamber flux observations and observations of CH 4 concentration profiles in the soil porewater from 2021 were used for evaluation of the model performance. Our results show that ELM-Wet reduces the model's CH 4 emission root mean squared error by up to 33% and is able to represent inter-daily CO 2 and CH 4 flux variability across the wetland's eco-hydrological patches, including during periods of extreme dry or wet conditions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Venus' Free Obliquity

The predicted orientation of Venus' rotation axis relative to its orbit ran be uniquely determined given knowledge of its J2 gravity coefficient and polar moment of inertia C if its free obliquity is fully damped. This assumption seems warranted given the dominant damping mechanism: turbulent fluid friction at a core mantle boundary (CMB). This skin friction results from differential obliquity of mantle and core spin axes, and the associated damping rate could be as short as 1/10(exp 6) year. However, the observed pole orientation indicates a free obliquity amplitude epsilon approx. = 2.1 deg compared with a nominal forced amplitude of 0.5 deg. There are two plausible explanations. The most likely is that the observed obliquity is a tidally evolved end state in which core friction, modulated by CMB ellipticity and core obliquity amplitude, counterbalances solid and atmospheric tidal torques. This concept is similar to the explanation for the retrograde spin omega as an end state in which solid and atmospheric thermal tidal torques balance at the present spin rate because of the omega(sup -1) dependence of the axial thermal torque. Weaker effects such as solid and thermal tides can then compete with core friction and for plausible models, their sum tends to increase free obliquity. The obliquity balance is controlled by the nonlinear (and nearly quadratic) dependence of the CMB turbulent "skin friction" torque on obliquity.Core ellipticity also has a profound effect on tidal evolution of Venus obliquity.

Yoder, Charles F.↗
DOE OSTI2026

A Flexible Quasi-Static Mooring Design Optimization Method for Floating Structures

This paper presents a flexible and efficient design method for optimizing the mooring systems of floating structures. Mooring system optimization is challenging because of the strong nonlinearity of mooring system behavior and the many technical constraints that must be satisfied. Furthermore, different mooring configurations can have very different design spaces. While some successful examples of mooring design optimization exist in the literature, developing an optimization approach that can work across various mooring design problems is a larger challenge. We present such a method based on a flexible parameterization that allows a wide variety of mooring designs to be described by a list of variables, a quasi-static mooring model that provides efficient evaluation of a mooring design without directly considering mooring system dynamics, and an optimization framework that generates, evaluates, and adjusts the mooring design while considering user-specified constraints such as offset limits, strength safety factors, and seabed contact limits. We demonstrate the design optimization framework on four mooring design problems, each for a different type of mooring system. We compare the use of different design modes to simplify the optimization problem, showing that they can reduce the computation time by up to 75%. We also compare different optimization algorithms and find that the resulting computational speed can vary by up to 51 times. We perform a sensitivity study on one design and find that the local sensitivity of anchoring radius to water depth has a positive correlation of 0.29, but the global sensitivity shows large nonlinearities. Lastly, we perform a coupled dynamic analysis on one of the optimized designs and find that the predicted mean platform motions and mooring line tensions are within 1% of dynamic results and the extreme motions and tensions are within 14%. Lastly, we show that a DEA-Chain-Polyester mooring configuration is cost-optimal for the given design problem of the demonstrations, which aligns with general industry practice.

16 TIDAL AND WAVE POWER↗
NASA NTRS2007

Fragmentation hierarchy of bright sungrazing comets and the birth and orbital evolution of the Kreutz System. II. the case for cascading fragmentation

We examine the process of cascading fragmentation for the Kreutz sungrazer system to continue our exploration of its birth, orbital evolution, and temporal clumping. We modify and broaden the two-superfragment model from Paper I to include clusters of ~30 bright comets spanning four centuries and 1000 SOHO sungrazers from 1996 to 2006. The spectacular parent sungrazer X/1106 C1 is assumed to have tidally split shortly after perihelion into a train of major protofragments immersed in a cloud of particulate debris, which at larger heliocentric distances were breaking up nontidally over and over again. We describe potential evolutionary paths for the Kreutz system by linking X/1106 C1 in subgroup I-type orbit with the comet of February 423 in one scenario or with the comet of February 467 in another. The latter scenario accounts for sungrazer clusters in as early as the 16th century, suggests that the progenitor object may have been observed as the comet of 214 BCE, is quite consistent with the orbital distribution of the SOHO sungrazers that sample the central filament of the Kreutz system between the clusters of major sungrazers, and predicts future clusters until ~2120. Comet X/1106 C1 and the common parent of C/1882 R1 and C/1965 S1 were two first-generation fragments of the progenitor that split nontidally on the way to its 5th century perihelion, reminiscent of the superfragments in Paper I. We provide computational tools needed for solving the problem of the Kreutz system’s orbital evolution, but no unique scenarios are presented for the individual comets. Another cluster of bright sungrazers is expected to arrive in the coming decades, its earliest member possibly just several years from now.

Chodas, Paul W.↗
NASA NTRS2000

The Mathematics of Navigating the Solar System

In navigating spacecraft throughout the solar system, the space navigator relies on three academic disciplines - optimization, estimation, and control - that work on mathematical models of the real world. Thus, the navigator determines the flight path that will consume propellant and other resources in an efficient manner, determines where the craft is and predicts where it will go, and transfers it onto the optimal trajectory that meets operational and mission constraints. Mission requirements, for example, demand that observational measurements be made with sufficient precision that relativity must be modeled in collecting and fitting (the estimation process) the data, and propagating the trajectory. Thousands of parameters are now determined in near real-time to model the gravitational forces acting on a spacecraft in the vicinity of an irregularly shaped body. Completing these tasks requires mathematical models, analyses, and processing techniques. Newton, Gauss, Lambert, Legendre, and others are justly famous for their contributions to the mathematics of these tasks. More recently, graduate students participated in research to update the gravity model of the Saturnian system, including higher order gravity harmonics, tidal effects, and the influence of the rings. This investigation was conducted for the Cassini project to incorporate new trajectory modeling features in the navigation software. The resulting trajectory model will be used in navigating the 4-year tour of the Saturnian satellites. Also, undergraduate students are determining the ephemerides (locations versus time) of asteroids that will be used as reference objects in navigating the New Millennium's Deep Space 1 spacecraft autonomously.

Hintz, Gerald↗
DOE OSTI2024

Analysis on Evaluations of Monterey Bay Aquarium Research Institute’s Wave Energy Converter’s Field Data Using WEC-Sim and Gazebo: A Simulation Tool Comparison

Although many studies have validated wave energy converter (WEC) numerical models against scaled prototype experimental data, there remains a notable lack of validation using data from full-scale deployed WECs. This paper compares two numerical models of Monterey Bay Aquarium Research Institute’s Wave Energy Converter (MBARI-WEC), a two-body point absorber with an electro-hydraulic power take-off system (PTO). The models are implemented in WEC-Sim/Simscape and Gazebo Simulator. A statistical analysis of the models was performed, and field results were obtained to compare the models’ accuracy in predicting the RMS piston velocity, RMS motor speed, and mean electric power compared to field data for 56 observations across varying sea states. The Gazebo model demonstrated a closer agreement across all three parameters for a majority of the observations. When compared to the field data, the Gazebo and WEC-Sim models exhibited average mean electric power overestimations of 13% and 22%, respectively.

16 TIDAL AND WAVE POWER↗
NASA NTRS1977

Lagrangian drifter design for the determination of surface currents by remote sensing

In estuaries, the interaction of wind, tidal current, and mixing of fresh and saline water produces a variable depth profile of current, with foam lines and convergence zones between water types. Careful measurement of surface currents via Lagrangian drifters requires a drifter design appropriate to both the depth of current to be measured and the tide and wind conditions of interest. The use of remote sensing to track drifters contributes additional constraints on drifter design. Several designs of biodegradable drifters which emit uranine dye plumes, resolvable in aerial imagery to 1:60,000 scale, were tested for wind drag in field conditions against data from calibrated current meters. A 20 cm-vaned wooden drifter and a window shade drifter set to 1.5 m depth had negligible wind drag in winds to 8 m/sec. Prediction of oil slick trajectories using surface current data and a wind factor should be approached cautiously, as surface current data may be wind-contaminated, while the usual 3.5% wind factor is appropriate only for currents measured at depth.

Gordon, H. H.↗
NASA NTRS2013

A Storm Surge and Inundation Model of the Back River Watershed at NASA Langley Research Center

This report on a Virginia Institute for Marine Science project demonstrates that the sub-grid modeling technology (now as part of Chesapeake Bay Inundation Prediction System, CIPS) can incorporate high-resolution Lidar measurements provided by NASA Langley Research Center into the sub-grid model framework to resolve detailed topographic features for use as a hydrological transport model for run-off simulations within NASA Langley and Langley Air Force Base. The rainfall over land accumulates in the ditches/channels resolved via the model sub-grid was tested to simulate the run-off induced by heavy precipitation. Possessing both the capabilities for storm surge and run-off simulations, the CIPS model was then applied to simulate real storm events starting with Hurricane Isabel in 2003. It will be shown that the model can generate highly accurate on-land inundation maps as demonstrated by excellent comparison of the Langley tidal gauge time series data (CAPABLE.larc.nasa.gov) and spatial patterns of real storm wrack line measurements with the model results simulated during Hurricanes Isabel (2003), Irene (2011), and a 2009 Nor'easter. With confidence built upon the model's performance, sea level rise scenarios from the ICCP (International Climate Change Partnership) were also included in the model scenario runs to simulate future inundation cases.

Loftis, Jon Derek↗
NASA NTRS2020

Candidate Electromagnetic Counterpart to the Binary Black Hole Merger Gravitational-Wave Event S190521g

We report the first plausible optical electromagnetic counterpart to a (candidate) binary black hole merger. Detected by the Zwicky Transient Facility, the electromagnetic flare is consistent with expectations for a kicked binary black hole merger in the accretion disk of an active galactic nucleus [B. McKernan, K. E. S. Ford, I. Bartoset al., Astrophys. J. Lett.884, L50 (2019)] and is unlikely [<O(0.01%))] due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernova and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing ora tidal disruption event, are ruled out or constrained to be <O(0.1%). If the flare is associated withS190521g, we find plausible values of total mass M(BBH) ∼ 100 Mꙩ, kick velocity v(k) ∼ 200 km/s at θ ∼ 60° in a disk with aspect ratio H/a ∼ 0.01(i.e., disk height H at radius a) and gas density ρ ∼ 10^(−10)g/cu.cm. The merger could have occurred at a disk migration trap (a ∼ 700 r(g); r(g) ≡ GM(SMBH)/sq.c, where M(SMBH) is the mass of the active galactic nucleus supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a reencountering with the disk in ∼ 1.6 yr(M(SMBH)/10^(8) Mꙩ)(a/10^(3)r(g))^(3/2).

M. J. Graham↗
NASA NTRS1974

Starspots on flare stars

Sizes of starspots on flare stars can be derived from the author's convection-cell hypothesis. The sizes are in fair agreement with those observed on YY Gem, CC Eri, and BY Dra by Bopp and Evans (1973). The hypothesis predicts that periodic brightness variations due to starspots are restricted to stars brighter than a critical absolute visual magnitude. A convective model of a starspot on YY Gem has been computed, assuming that the missing flux is in the form of Alfven waves. It is found that the surface field must exceed 10,000 G, and is probably less than about 30,000 G. With a surface field of 20,000 G, the effective temperature of the spot is in the range from 1590 to 1890 K, depending on the field gradient. These figures are to be compared with an effective temperature of 2000 K estimated from observations by Bopp and Evans. Efficient dynamo action is shown to be a possible mechanism for generating such large surface fields. There is a possibility that tidal effects may influence starspot formation.

Mullan, D. J.↗
NASA NTRS2010

Ocean Turbulence, III: New GISS Vertical Mixing Scheme

We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Gamma-alpha (alpha stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation Epsilon. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Gamma-alpha s are functions of both Ri and Rq (Richardson number and density ratio representing double diffusion, DD); the Gamma-alpha are different for heat, salt and momentum; in the case of heat, the traditional value Gamma-h = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Gamma-h can be three times as large, a result that we reproduce. The salt Gamma-s is given in terms of Gamma-h. The momentum Gamma-m has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Gamma-m(Ri,R-rho). Having tested Gamma-h, we then test the momentum Gamma-m by showing that the turbulent Prandtl number Gamma-m/Gamma-h vs. Ri reproduces the available data quite well. As for the dissipation epsilon, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean;s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable epsilon/N(exp 2) from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Gamma-alpha = 0.2 with the RSM result that brings into the problem the Ri, R-rho dependence of the Gamma-alpha; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Gamma-h(Ri,Rq). The RSM model reproduces well the NATRE data. Bimodal epsilon-distribution. NATRE data show that epsilon (Ri < 1) approximately equals 10epsilon(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri much greater than 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.

Canuto, V. M.↗