Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sim”

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 721 records · Page 40

Picometer Level Modeling of a Shared Vertex Double Corner Cube in the Space Interferometry Mission Kite Testbed

The Space Interferometry Mission (SIM) is a microarcsecond interferometric space telescope that requires picometer level precision measurements of its truss and interferometer baselines. Single-gauge metrology errors due to non-ideal physical characteristics of corner cubes reduce the angular measurement capability of the science instrument. Specifically, the non-common vertex error (NCVE) of a shared vertex, double corner cube introduces micrometer level single-gauge errors in addition to errors due to dihedral angles and reflection phase shifts. A modified SIM Kite Testbed containing an articulating double corner cube is modeled and the results are compared to the experimental testbed data. The results confirm modeling capability and viability of calibration techniques.

metrology↗

Space Interferometry Mission Science Overview

Conclusions: SIM is on track for a launch in 2010. Project 'Phase C' start in October 2005. SIM is a flexibly-pointed instrument Astrometric capability (wide angle) a) 4 pas parallax precision on targets down to V = 20 = Stellar and Galactic astrophysics; and b) 1 pas single measurement accuracy. Sensitive to terrestrial-mass planets around the nearest stars New proposal opportunities are: a) AO-2 for new (large) Key projects - 2005; and b) General Observer program, for smaller proposals - 2006.

planet searching↗

Mask Design for the Space Interferometry Mission Internal Metrology

This slide presentation reviews the mask design used for the internal metrology of the Space Interferometry Mission (SIM). Included is information about the project, the method of measurements with SIM, the internal metrology, numerical model of internal metrology, wavefront examples, performance metrics, and mask design

interferometry↗

Design and Performance Evaluation of a Resistive Control Using a Hydraulic PTO System for the TALOS Wave Energy Converter

This study is focused on developing a numerical model to evaluate the performance of a hydraulic PTO system for the TALOS Wave Energy Converter. The WEC device is described and the architecture of the hydraulic PTO system is presented with detail. The WEC is modeled using WEC-Sim, and the PTO is modeled using the Simscape Fluids library from Simulink. The hydraulic PTO is based on a constant pressure configuration that is suitable for WEC passive control. The hydraulic system is composed by a set of rectifying valves and two hydraulic accumulators that reduce the stiffness of the system and also serve as energy storage devices. One of the advantages of this hydraulic PTO architecture is the possibility of controlling the electric generator to operate around the optimal efficiency operating point. The main components of the hydraulic PTO are off-the-shelf devices that are commercially available, which will facility a future deployment of the designed system. The design variables used for this study are the accumulator size, the maximum pressure in the accumulators, the hydraulic motor maximum displacement, and the shaft speed in the electric generator. The performance of the system is evaluated individually, using sinusoidal inputs that replicates regular wave conditions. In addition to this, the numerical model of the PTO is coupled to a WEC-Sim simulation of the TALOS Wave Energy Converter with six PTOs to generate a wave-to-wire model. The main objective of this work is to present a comprehensive design methodology that could serve as a guideline for future research efforts focused on implementing control algorithms on multi degree of freedom WECs.

hydraulic systems↗

Moments of axial-vector GPD from lattice QCD: quark helicity, orbital angular momentum, and spin-orbit correlation

In this work, we present a lattice QCD calculation of the Mellin moments of the twist-2 axial-vector generalized parton distribution (GPD), $\overset{\sim }{H}\left(x,\xi, t\right)$ , at zero skewness, ξ, with multiple values of the momentum transfer, t. Our analysis employs the short-distance factorization framework on ratio-scheme renormalized quasi-GPD matrix elements. The calculations are based on an N f = 2 + 1 + 1 twisted mass fermions ensemble with clover improvement, a lattice spacing of a = 0.093 fm, and a pion mass of m π = 260 MeV. We consider both the iso-vector and iso-scalar cases, utilizing next-to-leading-order perturbative matching while omitting the disconnected contributions and gluon mixing in the iso-scalar case. For the first time, we determine the Mellin moments of $\overset{\sim }{H}$ up to the fifth order. From these moments, we discuss the quark helicity and orbital angular momentum contributions to the nucleon spin, as well as the spin-orbit correlations of the quarks. Additionally, we perform a Fourier transform over the momentum transfer, which allows us to explore the spin structure in the impact-parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Search for supersymmetry in final states with missing transverse momentum and charm-tagged jets using 139 fb−1 of proton-proton collisions at s = 13 TeV with the ATLAS detector

The paper presents a search for supersymmetric particles produced in proton-proton collisions at s$$ \sqrt{s} $$ = 13 TeV and decaying into final states with missing transverse momentum and jets originating from charm quarks. The data were taken with the ATLAS detector at the Large Hadron Collider at CERN from 2015 to 2018 and correspond to an integrated luminosity of 139 fb−1. No significant excess of events over the expected Standard Model background expectation is observed in optimized signal regions, and limits are set on the production cross-sections of the supersymmetric particles. Pair production of charm squarks or top squarks, each decaying into a charm quark and the lightest supersymmetric particle χ~10$$ {\overset{\sim }{\chi}}_1^0 $$, is excluded at 95% confidence level for squarks with masses up to 900 GeV for scenarios where the mass of χ~10$$ {\overset{\sim }{\chi}}_1^0 $$ is below 50 GeV. Additionally, the production of leptoquarks with masses up to 900 GeV is excluded for the scenario where up-type leptoquarks decay into a charm quark and a neutrino. Model-independent limits on cross-sections and event yields for processes beyond the Standard Model are also reported.

Experimental Nuclear Physics↗

Search for supersymmetry in final states with missing transverse momentum and multiple b -jets in proton-proton collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

A search for supersymmetry involving the pair production of gluinos decaying via third-generation squarks into the lightest neutralino ($\sim\atop{\chi}$ $0\atop{1}$) is reported. It uses LHC proton-proton collision data at a centre-of-mass energy s√=13 TeV with an integrated luminosity of 36.1 fb -1 collected with the ATLAS detector in 2015 and 2016. The search is performed in events containing large missing transverse momentum and several energetic jets, at least three of which must be identified as originating from b-quarks. To increase the sensitivity, the sample is divided into subsamples based on the presence or absence of electrons or muons. No excess is found above the predicted background. For $\sim\atop{\chi}$ $0\atop{1}$ masses below approximately 300 GeV, gluino masses of less than 1.97 (1.92) TeV are excluded at 95% confidence level in simplified models involving the pair production of gluinos that decay via top (bottom) squarks. An interpretation of the limits in terms of the branching ratios of the gluinos into third-generation squarks is also provided. Finally, these results improve upon the exclusion limits obtained with the 3.2 fb -1 of data collected in 2015.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Light (and darkness) from a light hidden Higgs

Here, we examine light diphoton signals from extended Higgs sectors possessing (approximate) fermiophobia with Standard Model (SM) fermions as well as custodial symmetry. This class of Higgs sectors can be realized in various beyond the SM scenarios and is able to evade many experimental limits, even at light masses, which are otherwise strongly constraining. Below the $WW$ threshold, the most robust probes of the neutral component are di and multi-photon searches. Utilizing the dominant Drell-Yan Higgs pair production mechanism and combining it with updated LHC diphoton data, we derive robust upper bounds on the allowed branching ratio for masses between $45 - 160$ GeV. Furthermore, masses $\lesssim 110$ GeV are ruled out if the coupling to photons is dominated by $W$ boson loops. We then examine two simple ways to evade these bounds via cancellations between different loop contributions or by introducing decays into an invisible sector. This also opens up the possibility of future LHC diphoton signals from a light hidden Higgs sector. As explicit realizations, we consider the Georgi-Machacek (GM) and Supersymmetric GM (SGM) models which contain custodial (degenerate) Higgs bosons with suppressed couplings to SM fermions and, in the SGM model, a (neutralino) LSP. We also breifly examine the recent $\sim 3\sigma$ CMS diphoton excess at $\sim 95$ GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Light signals from a lighter Higgs

Here, with the Higgs search program already quite mature, there is the exciting possibility of discovering a new particle with rates near that of the SM Higgs. We consider models with a signal in $\gamma \gamma$ below the SM Higgs mass, taking the recent $2.9\, \sigma$ (local) CMS excess at 95 GeV as a target. We discuss singlet models with additional vectorlike matter, but argue that a Type-I two Higgs doublet model provides a more economical scenario. In such a setup, going into regions of moderate-to-strong fermiophobia, the enhanced $\gamma \gamma$ branching ratio allows signals from $VH$+VBF production to yield $\sigma \times BR_{\gamma\gamma} $ comparable to total SM rates. Light $H$ production can be dominated via rare top decays $t \rightarrow b H^+ \rightarrow b W^{*} H$, which provides an alternate explanation of the excess. We consider this in the context of other Higgs anomalies, namely the LEP Higgs excess near the same mass, and excesses in $t\bar{t}h$ searches at Tevatron and LHC. We find that with $140\, \mathrm{GeV} < m_{H^+} < 160\, \mathrm{GeV}$, $\tan \beta \sim 5$ and a coupling to gauge bosons of $\sin^2 \delta \sim 0.1$, such a scenario can produce all the excesses simultanously, where $tth$ arise from contamination from the rare top decays, as previously proposed. An implication of the Type-I scenario is that any $\gamma \gamma$ excess should be associated with additional elements that could reduce background, including $b$-jets, forward jets or signs of vector boson production.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Orientation reversal and the Chern-Simons natural boundary

We show that the fundamental property of preservation of relations, underlying resurgent analysis, provides a new perspective on crossing a natural boundary, an important general problem in theoretical and mathematical physics. This reveals a deeper rigidity aspect of resurgence in a quantum field theory path integral. The physical context here is the non-perturbative completion of complex Chern-Simons theory that associates to a 3-manifold a collection of q-series invariants labeled by Spinc structures, for which crossing the natural boundary corresponds to orientation reversal of the 3-manifold. Our new resurgent perspective leads to a practical numerical algorithm that generates q-series which are dual to unary q-series composed of false theta functions. Until recently, these duals were only known in a limited number of cases, essentially based on Ramanujan’s mock theta functions, and the common belief was that the duals might not even exist in the general case. Resurgence analysis identifies as primary objects Mordell integrals: up to changes of variables, they are Laplace transforms of resurgent functions. Their unique Borel summed transseries decomposition on either side of the Stokes line is simply the unique decomposition into real and imaginary parts. In turn, the latter are combinations of unary q-series in terms of q and its modular counterpart $\overset{\sim }{q}$ , and are resurgent by construction. The Mordell integral is analytic across the natural boundary of the q and $\overset{\sim }{q}$ series, and uniqueness of a similar decomposition which preserves algebraic relations on the other side of the boundary defines the unique boundary crossing of the q series. We demonstrate that this continuation can be efficiently implemented numerically. In the cases where unique mock modular identities are known, they are found by this numerical procedure, but the procedure can go well beyond the known list of identities. A particularly interesting feature of the resurgent approach is that it reveals new aspects, and is very different from other known approaches based on indefinite theta series, Appell-Lerch sums, and representation theory of logarithmic vertex operator algebras.

Chern-Simons theories↗

Retained Austenite Attributes and Mechanical Performance of Different Compositions of TRIP Steel Alloys

The main objective of this work is to correlate the different attributes of retained austenite phase in TRIP steel structure with its mechanical properties. Three different chemical compositions were used to produce TRIP steel alloys. The first alloy contains only silicon as a carbide precipitation inhibitor; in the second alloy, silicon is partially replaced with aluminum, while in the third alloy phosphorus is added to the silicon–aluminum alloy. These compositions were proposed in order to avoid the negative effect of silicon during galvanizing. These three steel alloys were cast, homogenized, hot-forged and then subjected to intercritical annealing treatment followed by cooling and isothermal treatment in salt bath with subsequent air-cooling. Microstructural characterization, using optical microscope, SEM and XRD, along with tensile testing was carried out. Retained austenite attributes include volume fraction, carbon content, total carbon content, retained austenite transformable index and its stability coefficient against strain-induced martensite (SIM) transformation. The effect of these attributes on the mechanical properties of the studied TRIP steel alloys was evaluated using tensile test date, i.e., ultimate tensile strength (UTS), ductility, yield strength to UTS ratio, performance index and strain-hardening exponent. The variation of strain-hardening rate and strain-hardening exponent with strain was also considered. The failure behavior of the studied steel alloys was investigated using SEM on tensile-tested specimens. The results show that Si/Al-alloy has the lowest retained austenite transformable index, lowest stability coefficient against SIM transformation, highest performance index and more stable and uniform variation in strain-hardening rate and strain-hardening exponent with strain. These properties satisfy the industrial constrains of the sheets used in automotive industry relating to better mechanical properties and galvanizability.

36 MATERIALS SCIENCE↗

The MBARI-WEC: a power source for ocean sensing

Abstract Interest in wave energy converters to provide autonomous power to various ocean-bound systems, such as autonomous underwater vehicles, sensor systems, and even aquaculture farms, has grown in recent years. The Monterey Bay Aquarium Research Institute has developed and deployed a small two-body point absorber wave energy device suitable to such needs. This paper provides a description of the system to support future open-source access to the device and further the general development of similar wave energy systems. Additionally, to support future control design and system modification efforts, a set of hydrodynamic models are presented and cross-compared. To test the viability of using a linear frequency-domain admittance model for controller tuning, the linear model is compared against four WEC-Sim models of increasing complexity. The linear frequency-domain model is found to be generally adequate for capturing system dynamics, as the model agreement is good and the degree of nonlinearity introduced in the WEC-Sim models is generally less than 2.5%.

16 TIDAL AND WAVE POWER↗

Hydrodynamic characterization of the coastal pioneer array ocean observing system

Ocean observation buoys require relatively small amounts of power, yet traditionally necessitate costly resupply trips for battery replacement. With the offshore location of the buoys and small power requirements, wave energy may be an effective solution for providing consistent and reliable power to support the buoy instrumentation. The US National Science Foundation Ocean Observatories Initiative (OOI) includes arrays of point absorber-like buoy systems used for ocean observation that have been deployed at multiple locations including the Southern Mid-Atlantic Bight. A study is currently underway to design a pitch resonator wave energy converter to supplement existing renewable energy generation for powering observation instrumentation. This paper details field measurements from surface moorings of the OOI Coastal Pioneer Array, which informs the subsequent development of a numerical model for the moored observation system. The model is developed in Wave Energy Converter Simulator (WEC-Sim), which leverages the Simscape multibody solver within the MATLAB/Simulink framework and linear potential flow theory to simulate the hydrodynamic interactions and multibody dynamics in 6 degrees of freedom. Multiple tuning variables are considered to produce a model for the system that matches well with empirical data (about 8% error). In conclusion, the WEC-Sim model will serve as a platform for integrating the pitch resonator wave energy converter concept and deployment preparation (detailed design including power take-off and control systems, response evaluation, etc.).

hydrodynamic modeling↗

Improved natural convection heat transfer correlations for reactor cavity cooling systems of high-temperature gas-cooled reactors: From computational fluid dynamics to Pronghorn

The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.

42 ENGINEERING↗

Effect of the SiO 2 interlayer properties with solid-source hydrogenation on passivated contact performance and surface passivation

We investigate how SiO x oxide interlayers prepared by different techniques (chemical, thermal) in combination with hydrogen released from an ALD Al 2 O 3 source layer govern passivation in 1) passivated contacts based on doped poly-Si layers and tunneling SiO 2 , and 2) wafer surface passivation by Al 2 O 3 . Profiles of O and H in these structures with engineered, buried SiOx interlayers were measured by Time-of-Flight SIMS (TOF-SIMS) at nanometer resolution. Passivated contacts perform best with thermally oxidized SiOx, while chemical SiO x causes poly-Si film blistering and performance degradation. ALD Al 2 O 3 acts as passivating H source, significantly improving B-doped and intrinsic poly-Si contacts for IBC cells. Fast-diffusing hydrogen from the Al 2 O 3 source layer appears to penetrate Si wafer thickness, improving the passivation of structures at the opposite side. In contrast to the passivated contacts, chemical SiOx interlayer promotes wafer surface passivation by ALD Al 2 O 3 , while similarly thin thermal SiO 2 suppresses passivation and built-in charge.

14 SOLAR ENERGY↗

Coordinated U–Pb geochronology, trace element, Ti-in-zircon thermometry and microstructural analysis of Apollo zircons

Here, we present the results of a coordinated SIMS U–Pb, trace element, Ti-in-zircon thermometry, and microstructural study of 155 lunar zircons separated from Apollo 14, 15, and 17 breccia and soil samples that help resolve discrepancies between the zircon data, the lunar whole rock history and lunar magma ocean crystallization models. The majority of lunar grains are detrital fragments, some nearly 1 mm in length, of large parent crystals suggesting that they crystallized in highly enriched KREEP magmas. The zircon age distributions for all three landing sites exhibit an abundance of ages at ~4.33 Ga, however they differ in that only Apollo 14 samples have a population of zircons with ages between 4.1 and 3.9 Ga. These younger grains comprise only 10% of all dated lunar zircons and are usually small and highly shocked making them more susceptible to Pb-loss. These observations suggest that the majority of zircons crystallized before 4.1 Ga and that KREEP magmatism had predominantly ceased by this time. We also observed that trace element analyses are easily affected by contributions from inclusions (typically injected impact melt) within SIMS analyses spots. After filtering for these effects, rare-earth element (REE) abundances of pristine zircon are consistent with one pattern characterized by a negative Eu anomaly and no positive Ce anomaly, implying that the zircons formed in a reducing environment. This inference is consistent with crystallization temperatures based on measured Ti concentrations and new estimates of oxide activities which imply temperatures ranging between 958 ± 57 and 1321 ± 100 °C, suggesting that zircon parent magmas were anhydrous. Together, the lunar zircon ages and trace elements are consistent with a ≤300 My duration of KREEP magmatism under anhydrous, reducing conditions. We also report two granular texture zircons that contain baddeleyite cores, which both yield 207 Pb– 206 Pb ages of 4.33 Ga. These grains are our best constraints on impact ages within our sample population, and suggest at least one large impact is contemporaneous with the most common time of magmatic zircon formation on the Moon’s crust visited by the Apollo missions.

58 GEOSCIENCES↗

Toward an Understanding of Surface Layer Formation, Growth, and Transformation at the Glass–Fluid Interface

Silicate glass is a metastable and durable solid that has application to a number of energy and environmental challenges (e.g., microelectronics, fiber optics, and nuclear waste storage). If allowed to react with water over time silicate glass develops an altered layer at the solid-fluid interface. In this study, we used borosilicate glass (LAWB45) as model material to develop a robust understanding of altered layer formation (i.e., amorphous hydrated surface layer and crystalline reaction products). Experiments were conducted at high surface area-to-volume ratio (~200,000 m-1) and 90 °C in the pressurized unsaturated flow (PUF) apparatus for 1.5-years to facilitate the formation of thick altered layers and allow for the effluent solution chemistry to be monitored continuously. A variety of microscopy techniques were used to characterized reacted grains and suggest the average altered layer thickness is 13.2 ± 8.3 μm with the hydrated and clay layer representing 74.8% and 25.2% of the total altered layer, respectively. This estimate is within the experimental error of the value estimated from the B release rate data (~10 ±1 μm/yr) over the 1.5-year duration. PeakForce® quantitative nanomechanical mapping results suggest the hydrated layer has a modulus that ranges between ~20 to 40 GPa, which is in the range of porous silica that contains from ~20 to ~50% porosity, yet significantly lower than dense silica (~70 to 80 GPa). Scanning transmission electron microscopy (STEM) images confirm the presence of pores and an analysis of a higher resolution image of a region provides a qualitative estimate of ≥ 22% porosity in this layer with variations in the hydrated layer in void volume with increasing distance from the unaltered glass. Chemical composition analyses, based on a combination of time-of-flight secondary-ion mass spectrometry (ToF-SIMS) and scanning electron microscopy with X-ray energy dispersive spectroscopy (EDS) and STEM-EDS, clearly show that the altered layer is mainly composed of Al, H, Si, and O with the clay layer being enriched in Li, Zn, Fe, and Mg. The amorphous hydrated layer is enriched in Ca, H, and Zr with a minor amount of K. Furthermore, ToF-SIMS results also suggest the B profile is anti-correlated with the H profile in the hydrated layer. Our selected-area electron diffraction results suggest the structure of the hydrated layer closely resembles opal-AG (amorphous gel-like) with an average crystallite size of ~0.7 nm which is smaller than the critical nucleus for silica nanoparticles (i.e., 1.4 to 3 nm). These results suggest the hydrated layer is more consistent with a polymeric gel rather than a colloidal gel and is comprised of molecular units (<1 nm in size) that result from the difficult to hydrolyze bonds, such as Si—O—Zr units, during the glass corrosion process. The size of individual particles or molecular units is a function of formation conditions (e.g., pH, ionic strength, nano-confinement, solute composition) in the hydrated layer.

36 MATERIALS SCIENCE↗

Constraining the magnitude of the carbon isotope excursion during the Paleocene-Eocene thermal maximum using larger benthic foraminifera

The Paleocene-Eocene thermal maximum (PETM) was an extraordinary pulse of global warming approximately 56 Ma. This warmingevent is associated with an addition of large amounts of 13C-depleted carbon into the atmosphere-ocean system, but the magnitude of the negative carbon isotope excursion (CIE)signaling the PETMand often used to estimate mass of the released carbonis still debated. Here we gauge the CIE magnitude usingsecondary ion mass spectrometry (SIMS) to perform in situ δ 13 C measurements within individual larger benthic foraminifera preserved in a tropical shallow marine limestone section in south Tibet. This SIMS-based δ 13 C record yields a negative CIE (Δ~7‰) comparable in magnitude to that registered by terrestrial records but larger than the ~4‰CIE returned by surface-dwelling planktonic foraminiferain deep-sea records. We therefore posit that the CIE magnitude was~7‰ in the surface ocean, and that theprevious~4‰ estimates for the CIE areattenuated by incomplete preservation and diagenetic overprinting. The large CIE magnitude in the surface ocean (~7‰) and the evidently smaller CIE magnitude registeredby deep-sea benthic foraminifera(~2-3.5‰) arelikely a function of rateand duration of carbon release as well as ocean stratification during the PETM.

58 GEOSCIENCES↗