SEARCH · Search NASA

Results for “sensitivity gradient”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Toward accurate measurement of electromagnetic field by retrieving and refining the center position of non-uniform diffraction disks in Lorentz 4D-STEM

Recent advancement in scanning transmission electron microscopy (STEM) allows the use of 4D-STEM, a technique that captures an electron diffraction pattern at each scan point in STEM, to measure electrostatic and magnetic potential and field in materials. However, accurate measurement, separation of the magnetic and electric signals, and removal of artifacts remain challenging, especially in the presence of complex non-uniform diffraction contrast within the disks. In this work, based on dynamic simulations of 4D-STEM patterns built upon superstructures consisting of millions of atoms to account for different sample thickness and edge geometries, we show how the shape and intensity distribution of the central disk are affected by multiple scattering. We propose a robust refinement procedure through iteration of the spin-sensitive peak position of the disk-center in the circular Hough transform filtered images from experimental Lorentz 4D-STEM dataset after minimizing the possible artifacts, such as those due to the change of thickness, dynamic scattering, and scanning process. We verify that caution must be taken as in practice the rigid-disk-shift model used to reconstruct induction maps can easily break down due to disk-protrusion when there exists a nonconstant phase gradient or thickness within the width of the probe. Through quantitative analysis and comparing experiment with calculation the effect of the non-spin-related intensity distribution inside the disk as well as that causes the disk shift due to the intensity-protrusion can be removed, and high-quality magnetic field mapping is possible.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1999

Climate Change and the Brewer-Dobson Circulation

Using GCM simulations of the troposphere-stratosphere circulation in different climate regimes, the sensitivity of the Brewer-Dobson Circulation to climate change will be discussed. Included will be simulations of the stratosphere with increasing/doubled C02 as done with different versions of the GISS Global Climate Middle Atmosphere Model, with and without ozone response. To put the results in a broader perspective, stratospheric circulation changes from several paleoclimate simulations will also be presented, including the Last Glacial Maximum and the Paleocene, each with several different (plausible) latitudinal sea surface temperature gradients. Hence the climate changes investigated will include strong variations in baroclinic and topographic tropospheric planetary wave and gravity wave forcing, as well as global mean temperature and atmospheric C02 levels. Results will be analyzed in terms of wave driving via E-P flux convergences and gravity wave effects, and will be shown to be strongly dependent on wave propagation characteristics.

Rind, David↗

NASA NTRS2010

Mission Simulation of Space Lidar Measurements for Seasonal and Regional CO2 Variations

Results of mission simulation studies are presented for a laser-based atmospheric [82 sounder. The simulations are based on real-time carbon cycle process modeling and data analysis. The mission concept corresponds to the Active Sensing of [82 over Nights, Days, and Seasons (ASCENDS) recommended by the US National Academy of Sciences Decadal Survey of Earth Science and Applications from Space. One prerequisite for meaningful quantitative sensor evaluation is realistic CO2 process modeling across a wide range of scales, i.e., does the model have representative spatial and temporal gradients? Examples of model comparison with data will be shown. Another requirement is a relatively complete description of the atmospheric and surface state, which we have obtained from meteorological data assimilation and satellite measurements from MODIS and [ALIPS0. We use radiative transfer model calculations, an instrument model with representative errors ' and a simple retrieval approach to complete the cycle from "nature" run to "pseudo-data" CO2, Several mission and instrument configuration options are examined/ and the sensitivity to key design variables is shown. We use the simulation framework to demonstrate that within reasonable technological assumptions for the system performance, relatively high measurement precision can be obtained, but errors depend strongly on environmental conditions as well as instrument specifications. Examples are also shown of how the resulting pseudo - measurements might be used to address key carbon cycle science questions.

Kawa, Stephan↗

DOE OSTI2024

Offshore application of landslide susceptibility mapping using gradient-boosted decision trees: a Gulf of Mexico case study

Abstract Among natural hazards occurring offshore, submarine landslides pose a significant risk to offshore infrastructure installations attached to the seafloor. With the offshore being important for current and future energy production, there is a need to anticipate where future landslide events are likely to occur to support planning and development projects. Using the northern Gulf of Mexico (GoM) as a case study, this paper performs Landslide Susceptibility Mapping (LSM) using a gradient-boosted decision tree (GBDT) model to characterize the spatial patterns of submarine landslide probability over the United States Exclusive Economic Zone (EEZ) where water depths are greater than 120 m. With known spatial extents of historic submarine landslides and a Geographic Information System (GIS) database of known topographical, geomorphological, geological, and geochemical factors, the resulting model was capable of accurately forecasting potential locations of sediment instability. Results of a permutation modelling approach indicated that LSM accuracy is sensitive to the number of unique training locations with model accuracy becoming more stable as the number of training regions was increased. The influence that each input feature had on predicting landslide susceptibility was evaluated using the SHapely Additive exPlanations (SHAP) feature attribution method. Areas of high and very high susceptibility were associated with steep terrain including salt basins and escarpments. This case study serves as an initial assessment of the machine learning (ML) capabilities for producing accurate submarine landslide susceptibility maps given the current state of available natural hazard-related datasets and conveys both successes and limitations.

Dyer, Alec S. (ORCID:0000000219813904)↗

DOE OSTI2025

Prediction of electric and magnetic fields from spectral data using machine learning algorithms for Doppler-free saturation spectroscopy diagnostics

The prediction of electric and magnetic field amplitudes from atomic spectral data is critical for plasma control in fusion devices such as tokamaks. Conventional approaches that rely on physics-based models are computationally expensive and unsuitable for real-time applications. In this work, we develop and benchmark three machine learning algorithms—simulation-based inference (SBI), fully connected neural networks (FCNN), and histogram-based gradient boosting regression (GBR-Hist)—to infer field intensities directly from Doppler-free saturation spectroscopy (DFSS) spectra. Synthetic datasets of spectra were generated using the EZSSS code and evaluated both with and without added Poisson noise to mimic experimental conditions. We find that SBI achieves the highest accuracy and robustness, FCNN provides a strong balance of accuracy and computational efficiency for real-time applications, and GBR-Hist offers the fastest inference but is more sensitive to noise. Furthermore, these results demonstrate the potential of machine learning to accelerate DFSS analysis and enhance its utility for plasma diagnostics and control.

Doppler-free saturation spectroscopy↗

DOE OSTI2016

Warming and provenance limit tree recruitment across and beyond the elevation range of subalpine forest

Climate niche models project that subalpine forest ranges will extend upslope with climate warming. These projections assume that the climate suitable for adult trees will be adequate for forest regeneration, ignoring climate requirements for seedling recruitment, a potential demographic bottleneck. Moreover, local genetic adaptation is expected to facilitate range expansion, with tree populations at the upper forest edge providing the seed best adapted to the alpine. Here, we test these expectations using a novel combination of common gardens, seeded with two widely distributed subalpine conifers, and climate manipulations replicated at three elevations. Infrared heaters raised temperatures in heated plots, but raised temperatures more in the forest than at or above treeline because strong winds at high elevation reduced heating efficiency. Watering increased season-average soil moisture similarly across sites. Contrary to expectations, warming reduced Engelmann spruce recruitment at and above treeline, as well as in the forest. Warming reduced limber pine first-year recruitment in the forest, but had no net effect on fourth-year recruitment at any site. Watering during the snow-free season alleviated some negative effects of warming, indicating that warming exacerbated water limitations. Contrary to expectations of local adaptation, low-elevation seeds of both species initially recruited more strongly than high-elevation seeds across the elevation gradient, although the low-provenance advantage diminished by the fourth year for Engelmann spruce, likely due to small sample sizes. High- and low-elevation provenances responded similarly to warming across sites for Engelmann spruce, but differently for limber pine. In the context of increasing tree mortality, lower recruitment at all elevations with warming, combined with lower quality, high-provenance seed being most available for colonizing the alpine, portends range contraction for Engelmann spruce. The lower sensitivity of limber pine to warming indicates a potential for this species to become more important in subalpine forest communities in the coming centuries.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1993

Automated Finite Element Modeling of Wing Structures for Shape Optimization

The displacement formulation of the finite element method is the most general and most widely used technique for structural analysis of airplane configurations. Modem structural synthesis techniques based on the finite element method have reached a certain maturity in recent years, and large airplane structures can now be optimized with respect to sizing type design variables for many load cases subject to a rich variety of constraints including stress, buckling, frequency, stiffness and aeroelastic constraints (Refs. 1-3). These structural synthesis capabilities use gradient based nonlinear programming techniques to search for improved designs. For these techniques to be practical a major improvement was required in computational cost of finite element analyses (needed repeatedly in the optimization process). Thus, associated with the progress in structural optimization, a new perspective of structural analysis has emerged, namely, structural analysis specialized for design optimization application, or.what is known as "design oriented structural analysis" (Ref. 4). This discipline includes approximation concepts and methods for obtaining behavior sensitivity information (Ref. 1), all needed to make the optimization of large structural systems (modeled by thousands of degrees of freedom and thousands of design variables) practical and cost effective.

Harvey, Michael Stephen↗

DOE OSTI2015

Systematic evaluation of photodetector performance for plastic scintillation dosimetry

Purpose: The authors’ objective was to systematically assess the performance of seven photodetectors used in plastic scintillation dosimetry. The authors also propose some guidelines for selecting an appropriate detector for a specific application. Methods: The plastic scintillation detector (PSD) consisted of a 1-mm diameter, 10-mm long plastic scintillation fiber (BCF-60), which was optically coupled to a clear 10-m long optical fiber of the same diameter. A light-tight plastic sheath covered both fibers and the scintillator end was sealed. The clear fiber end was connected to one of the following photodetectors: two polychromatic cameras (one with an optical lens and one with a fiber optic taper replacing the lens), a monochromatic camera with an optical lens, a PIN photodiode, an avalanche photodiode (APD), or a photomultiplier tube (PMT). A commercially available W1 PSD was also included in the study, but it relied on its own fiber and scintillator. Each PSD was exposed to both low-energy beams (120, 180, and 220 kVp) from an orthovoltage unit and high-energy beams (6 and 23 MV) from a linear accelerator. Various dose rates were tested to identify the operating range and accuracy of each photodetector. Results: For all photodetectors, the relative uncertainty was less than 5% for dose rates higher than 3 mGy/s. The cameras allowed multiple probes to be used simultaneously, but they are less sensitive to low-light signals. The PIN, APD, and PMT had higher sensitivity, making them more suitable for low dose rate and out-of-field dose monitoring. The relative uncertainty of the PMT was less than 1% at the lowest dose rate achieved (0.10 mGy/s), suggesting that it was optimal for use in live dosimetry. Conclusions: For dose rates higher than 3 mGy/s, the PIN diode is the most effective photodetector in terms of performance/cost ratio. For lower dose rates, such as those seen in interventional radiology or high-gradient radiotherapy, PMTs are the optimal choice.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Role of Maritime Continent Land Convection on the Mean State and MJO Propagation

The Maritime Continent (MC) region is known as a ‘barrier’ in the life-cycle of the Madden-Julian oscillation (MJO). During boreal winter, the MJO detours the equatorial MC land region southward and propagates through the oceanic region. Also, about half of the MJO events that initiate over the Indian Ocean cease around the MC. The mechanism through which the MC affects MJO propagation, however, has remained unanswered. The current study investigates the MJO-MC interaction with a particular focus on the role of MC land convection. Using a global climate model that simulates both mean climate and MJO realistically, we performed two sensitivity experiments in which updraft plume radius is set to its maximum and minimum value only in the MC land grid points, making convective top deeper and shallower, respectively. In this work, our results show that MC land convection plays a key role in shaping the 3-D climatological moisture distribution around the MC through its local and non-local effects. Shallower and weaker MC land convection results in a steepening of the vertical and meridional mean moisture gradient over the MC region. The opposite is the case when MC land convection becomes deeper and stronger. The MJO’s eastward propagation is enhanced (suppressed) with the steeper (lower) mean moisture gradient. The moist static energy (MSE) budget of the MJO reveals the vertical and meridional advection of the mean MSE by MJO wind anomalies as the key processes that are responsible for the changes in MJO propagation characteristics. Our results pinpoint the critical role of the background moisture gradient on MJO propagation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Laser and Electron Beam Technology for Parity Violating Electron Scattering Measurements

Parity violating electron scattering (PVES) is a precision tool used in a broad program of experiments which include studying the structure of protons and nuclei and searching for new Beyond the Standard Model (BSM) physics. In parity violation electron scattering experiments, a longitudinally polarized electron beam is incident on an unpolarized target. The sign of the longitudinal polarization is changed (making a parity transformation), and the fractional rate difference between right and left helicity states, dσ R and dσ L , is measured. An interference between the electromagnetic and weak amplitudes, M γ and M Z respectively, gives rise to a parity violating asymmetry A PV = $\frac{σ_{R}-σ_{L}}{σ_{R}-σ_{L}}$, defined by the fractional rate difference between right and left helicity states. A PV can be measured to extraordinarily high precision and is proportional to the ratio of the weak and electromagnetic amplitudes. This thesis highlights three such PVES experiments as well as a new polarized beam source which will improve their precision. The PVES experiments are PREX-II, MOLLER, and a 12 C transverse asymmetry measurement at Mainz. Chapters 1-3 of this thesis cover introductions to MOLLER and PREX-II, Chapter 4 contains the results of the 12 C transverse asymmetry measurement at Mainz, and Chapters 5-6 contain the original work performed by this author. PREX-II, the Lead Radius Experiment, measures the weak skin of the 208 Pb nucleus, provides a clean measurement of R N , the RMS radius of neutrons in a heavy nucleus, and constrains the equation of state (EOS) of highly dense matter which is important for describing neutron star structure, heavy ion collisions, and atomic parity violation experiments. The 12 C transverse asymmetry measurement at Mainz, and the future Mainz program of measuring the 208 Pb transverse asymmetry as well, addresses an important systematic correction for PREX and other PVES experiments preformed with heavy-nuclei. MOLLER (Measurement Of Lepton Lepton Electroweak Reaction) is an extremely precise PVES experiment searching for new neutral currents in electron-electron scattering. It’s called MOLLER because it will measure A PV in the Moller ee scattering process, and will infer the weak charge of the electron $Q^e_W$ to extremely high precision. MOLLER is sensitive to new Beyond the Standard Model physics at MeV and multi-TeV scales and will serve as an indirect complementary measurement to direct searches at high energy colliders. To achieve high precision measurements on A PV for MOLLER, we have developed an innovative Rubidium Titanyl Phosphate (RTP) Pockels cell in the polarized source to satisfy both statistical and systematic requirements as regards the electron beam produced. This new ultra-fast RTP cell design uses electric field gradients to provide unprecedented control over helicity correlated beam asymmetries and has been demonstrated to be capable of producing precisely controlled polarized electron beam at Jefferson Laboratory, controlling beam steering down to the nm-level. The precision reached with the RTP cell offers sufficient control over and minimization of helicity correlated beam asymmetries to perform PREX II. The RTP Pockels cell system will provide fast flipping and suitable control helicity correlated beam asymmetries and parity quality beam for the future MOLLER experiment, providing an unprecedented precision on the electron weak charge and electroweak mixing angle.

Palatchi, Caryn↗

DOE OSTI2026

Refining a Novel Process Monitoring Method to Safeguard Continuously Cycling Designs Using Isotopic Ratios

In advanced reactor (AR) designs, a common feature is continuous chemical processing and circulation of the nuclear material. This work bridges a significant measurement gap in safeguarding reactors with circulating fuel or continuous refueling by leveraging and building on the isotope ratio method first developed by our team under an FY21 Advanced Reactors International Safeguards Engagement (ARISE) project (Uribe et al. 2021). In circulating fuel designs, the radioisotope inventory changes from traditional effects (e.g., radioactive decay, fission) but also includes material transport due to pressure and temperature gradients. Such designs may also require regular or continuous additions or removals during operation, which significantly increases the rate of inventory change compared to traditional pressurized water reactor (PWR)s. Thus, directly tracking the nuclear inventory is ineffective since the isotopes are continuously added and removed. The isotope ratio method instead focuses on detecting changes to the input and output flows of radioisotopes. Previous work showed that for well-chosen pairs of isotopes, the isotopic ratio provides a sensitive and lasting indicator of deviation from normal conditions (e.g., startup, shutdown, diversion). The isotope ratio method is a process monitoring method with potential for application in for forward-looking approaches to International Atomic Energy Agency (IAEA) safeguards. The original process monitoring method was developed for a specific case—the decay tank of a thorium-fueled molten salt breeder reactor. In this expanded work, we explored other types of reactors and processes with nonstationary (e.g., flowing) nuclear material, which are difficult to safeguard with traditional methods because of the transient nature of the systems. The goal of the present work is to generalize the isotope ratio method for use in processes with continuously flowing nuclear material. All continuous processes have an average time for isotopes to be replaced in the system. The isotope ratio method works by choosing radioisotopes with half-lives both above and below the processing time. The present work seeks to explore which isotopes are suitable for the method by simulating the nuclear inventory, radioisotope emissions, and detector responses for several classes of advanced reactors. While the method can in principle be applied to other processes (e.g., enrichment or reprocessing facilities), the present work limits scope to ARs with continuously flowing fuel. Section 2 details the mathematics supporting the isotope ratio method, and Section 3 introduces the representative ARs selected for this work. Section 4 discusses how each reactor was analyzed, and Section 5 showcases the results for each representative reactor. Finally, Section 6 provides concluding remarks and suggests pathways for further analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DOE OSTI2017

Describing rainfall in northern Australia using multiple climate indices

Savanna landscapes are globally extensive and highly sensitive to climate change, yet the physical processes and climate phenomena which affect them remain poorly understood and therefore poorly represented in climate models. Both human populations and natural ecosystems are highly susceptible to precipitation variation in these regions due to the effects on water and food availability and atmosphere–biosphere energy fluxes. Here we quantify the relationship between climate phenomena and historical rainfall variability in Australian savannas and, in particular, how these relationships changed across a strong rainfall gradient, namely the North Australian Tropical Transect (NATT). Climate phenomena were described by 16 relevant climate indices and correlated against precipitation from 1900 to 2010 to determine the relative importance of each climate index on seasonal, annual and decadal timescales. Precipitation trends, climate index trends and wet season characteristics have also been investigated using linear statistical methods. In general, climate index–rainfall correlations were stronger in the north of the NATT where annual rainfall variability was lower and a high proportion of rainfall fell during the wet season. This is consistent with a decreased influence of the Indian–Australian monsoon from the north to the south. Seasonal variation was most strongly correlated with the Australian Monsoon Index, whereas yearly variability was related to a greater number of climate indices, predominately the Tasman Sea and Indonesian sea surface temperature indices (both of which experienced a linear increase over the duration of the study) and the El Niño–Southern Oscillation indices. These findings highlight the importance of understanding the climatic processes driving variability and, subsequently, the importance of understanding the relationships between rainfall and climatic phenomena in the Northern Territory in order to project future rainfall patterns in the region.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

Automated path planning of the Payload Inspection and Processing System

The Payload Changeout Room Inspection and Processing System (PIPS) is a highly redundant manipulator intended for performing tasks in the crowded and sensitive environment of the Space Shuttle Orbiter payload bay. Its dexterity will be exploited to maneuver the end effector in a workspace populated with obstacles. A method is described by which the end effector of a highly redundant manipulator is directed toward a target via a Lyapunov stability function. A cost function is constructed which represents the distance from the manipulator links to obstacles. Obstacles are avoided by causing the vector of joint parameters to move orthogonally to the gradient of the workspace cost function. A C language program implements the algorithm to generate a joint history. The resulting motion is graphically displayed using the Interactive Graphical Robot Instruction Program (IGRIP) produced by Deneb Robotics. The graphical simulation has the potential to be a useful tool in path planning for the PIPS in the Shuttle Payload Bay environment.

Byers, Robert M.↗

NASA NTRS1999

Quantitative Species Measurements in Microgravity Combustion Flames using Near-Infrared Diode Lasers

Understanding the physical phenomena controlling the ignition and spread of flames in microgravity has importance for space safety as well as for characterizing dynamical and chemical combustion processes which are normally masked by buoyancy and other gravity-related effects. Unfortunately, combustion is highly complicated by fluid mechanical and chemical kinetic processes, requiring the use of numerical modeling to compare with carefully designed experiments. More sophisticated diagnostic methods are needed to provide the kind of quantitative data necessary to characterize the properties of microgravity combustion as well as provide accurate feedback to improve the predictive capabilities of the models. Diode lasers are a natural choice for use under the severe conditions of low gravity experiments. Reliable, simple solid state operation at low power satisfies the operational restrictions imposed by drop towers, aircraft and space-based studies. Modulation wavelength absorption spectroscopy (WMS) provides a means to make highly sensitive and quantitative measurements of local gas concentration and, in certain cases, temperature. With near-infrared diode lasers, detection of virtually all major combustion species with extremely rapid response time is possible in an inexpensive package. Advancements in near-infrared diode laser fabrication technology and concurrent development of optical fibers for these lasers led to their use in drop towers. Since near-infrared absorption line strengths for overtone and combination vibrational transitions are weaker than the mid-infrared fundamental bands, WMS techniques are applied to increase detection sensitivity and allow measurement of the major combustion gases. In the first microgravity species measurement, Silver et al. mounted a fiber-coupled laser at the top of the NASA 2.2-sec drop tower and piped the light through a single-mode fiber to the drop rig. A fiber splitter divided the light into eight channels that directed the laser beam across a methane or propane diffusion jet flame. The light beams were recaptured by a set of gradient index lenses, coupled back into separate fiber optic lines, and transmitted back to detectors and electronics in the instrument package. In these experiments a 6-mm od fiber cable (containing the nine optical fibers) fell with the drop rig. Using separate detection and demodulation channels, spatial and temporal (up to 20 Hz) maps of water vapor and methane concentrations were obtained at differing heights in the flames. While this apparatus was useful from a demonstration standpoint, several drawbacks needed attention before useful scientific measurements could be obtained. First, eight lines of sight are somewhat insufficient for detailing the spatial profiles of the gas. Second, multiple detection channels operating in parallel are both expensive and present a challenge for accurate calibration. As a result, a newer scanning system was developed in our first contract under this program. The primary characteristic of this system is that it contains a single detection channel and achieves "continuous" spatial resolution by scanning the laser beam across the flame region, then directing this beam onto a single detector. Thus spatial measurements are converted to a temporal series of data. The true spatial resolution is limited only by the beam diameter and width of the sweep. In these experiments the beam is focused to about 1-mm diameter and scans across a region up to 4-cm wide.

Silver, Joel A.↗

NASA NTRS1998

3-D Modeling of Double-Diffusive Convection During Directional Solidification of a Non-Dilute Alloy with Application to the HgCdTe Growth Under Microgravity Conditions

A numerical calculation for a non-dilute alloy solidification was performed using the FIDAP finite element code. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. Estimates of convection contributions for various gravity conditions was performed parametrically. For gravity levels higher then 1 0 -7 of earth's gravity it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by factor of 50% for precise orientation of the ampoule in the microgravity environment. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D modeling are compared with previous 2-D finding. A video film featuring melt convection will be presented.

Bune, Andris V.↗

NASA NTRS2012

Effects of Mesh Irregularities on Accuracy of Finite-Volume Discretization Schemes

The effects of mesh irregularities on accuracy of unstructured node-centered finite-volume discretizations are considered. The focus is on an edge-based approach that uses unweighted least-squares gradient reconstruction with a quadratic fit. For inviscid fluxes, the discretization is nominally third order accurate on general triangular meshes. For viscous fluxes, the scheme is an average-least-squares formulation that is nominally second order accurate and contrasted with a common Green-Gauss discretization scheme. Gradient errors, truncation errors, and discretization errors are separately studied according to a previously introduced comprehensive methodology. The methodology considers three classes of grids: isotropic grids in a rectangular geometry, anisotropic grids typical of adapted grids, and anisotropic grids over a curved surface typical of advancing layer grids. The meshes within the classes range from regular to extremely irregular including meshes with random perturbation of nodes. Recommendations are made concerning the discretization schemes that are expected to be least sensitive to mesh irregularities in applications to turbulent flows in complex geometries.

Diskin, Boris↗

NASA NTRS1985

Effects of dynamic heat fluxes on model climate sensitivity Meridional sensible and latent heat fluxes

The high- and low-latitude radiative-dynamic (HLRD) climatic model of Wang et al. (1984) was used to study the effect of meridional heat (MH) fluxes on climate changes caused by increases of CO2 abundance and solar constant variations. However, the empirical MH parameterization of the HLRD model was replaced by physically based parameterization, which gives separate meridional sensible and latent heat fluxes and provides a complete representation of the dependence of the flux on the mean temperature field. Both parameterization methods yielded about the same changes in global mean surface temperature and ice line, and both produced only small changes in meridional temperature gradient, although the latter were even smaller with the physically based parameterizations. At any latitude, the hemispheric mean surface temperature, rather than MH fluxes, dominates the surface temperature changes.

Gutowski, W. J., Jr.↗

DOE OSTI2018

The KPP Boundary Layer Scheme for the Ocean: Revisiting Its Formulation and Benchmarking One‐Dimensional Simulations Relative to LES

Abstract We evaluate the Community ocean Vertical Mixing project version of the K‐profile parameterization (KPP) for modeling upper ocean turbulent mixing. For this purpose, one‐dimensional KPP simulations are compared across a suite of oceanographically relevant regimes against horizontally averaged large eddy simulations (LESs). We find the standard configuration of KPP consistent with LES across many forcing regimes, supporting its physical basis. Our evaluation also motivates recommendations for KPP best practices within ocean circulation models and identifies areas where further research is warranted. The original treatment of KPP recommends the matching of interior diffusivities and their gradients to the KPP‐predicted values computed in the ocean surface boundary layer (OSBL). However, we find that difficulties in representing derivatives of rapidly changing diffusivities near the base of the OSBL can lead to loss of simulation fidelity. To mitigate this difficulty, we propose and evaluate two computationally simpler approaches: (1) match to the internal predicted diffusivity alone and (2) set the KPP diffusivity to 0 at the OSBL base. We find the KPP entrainment buoyancy flux to be sensitive to vertical grid resolution and details of how to diagnose the KPP boundary layer depth. We modify the KPP turbulent shear velocity parameterization to reduce resolution dependence. Additionally, an examination of LES vertical turbulent scalar flux budgets shows that the KPP‐parameterized nonlocal tracer flux is incomplete due to the assumption that it solely redistributes the surface tracer flux. This result motivates further studies of the nonlocal flux parameterization.

54 ENVIRONMENTAL SCIENCES↗