Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Using Radiative Convective Equilibrium to Explore Clouds and Climate in the Community Atmosphere Model

Abstract Characteristics of, and fundamental differences between, the radiative‐convective equilibrium (RCE) climate states following the Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP) protocols in the Community Atmosphere Model version 5 (CAM5) and version 6 (CAM6) are presented. This paper explores the characteristics of clouds, moisture, precipitation and circulation in the RCE state, as well as the tropical response to surface warming, in CAM5 and CAM6 with different parameterizations. Overall, CAM5 simulates higher precipitation rates that result in larger global average precipitation, despite lower outgoing longwave radiation compared to CAM6. Differences in the structure of clouds, particularly the amount and vertical location of cloud liquid, exist between the CAM versions and can, in part, be related to distinct representations of shallow convection and boundary layer processes. Both CAM5 and CAM6 simulate similar peaks in cloud fraction, relative humidity, and cloud ice, linked to the usage of a similar deep convection parameterization. These anvil clouds rise and decrease in extent in response to surface warming. More generally, extreme precipitation, aggregation of convection, and climate sensitivity increase with warming in both CAM5 and CAM6. This analysis provides a benchmark for future studies that explore clouds, convection, and climate in CAM with the RCEMIP protocols now available in the Community Earth System Model. These results are discussed within the context of realistic climate simulations using CAM5 and CAM6, highlighting the usefulness of a hierarchical modeling approach to understanding model and parameterization sensitivities to inform model development efforts.

Reed, Kevin A.↗
DOE OSTI2023

Soil Respiration Responses to Throughfall Exclusion Are Decoupled From Changes in Soil Moisture for Four Tropical Forests, Suggesting Processes for Ecosystem Models

Abstract Climatic drying is predicted for many tropical forests yet models remain poorly parameterized for these ecosystems, hampering predictions of forest‐climate interactions. We applied an integrated model–experiment approach, parameterizing an ecosystem model with tropical forest observational data and comparing model predictions to a field drying manipulation. We hypothesized that drying suppresses soil CO 2 fluxes (i.e., respiration) in relatively dry tropical forests but increases CO 2 fluxes in wetter tropical forests by alleviating anaerobiosis. We measured soil CO 2 fluxes during wet‐dry cycles from 2015 to 2022 in four Panamanian forests that vary in rainfall and soil fertility. Measured soil CO 2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in lower soil moisture optima for respiration than previously modeled. We then parameterized the model using field data and the new moisture‐respiration response functions. The updated model predicted increased soil CO 2 fluxes with drying in wetter and fertile forests and suppressed fluxes in drier, infertile forests. In contrast to model predictions, a chronic throughfall exclusion experiment initially suppressed soil respiration across forests, with sustained suppression for four years in the wettest forest only (−28% ± 4% during the dry season). In the fertile forest, drying eventually elevated CO 2 fluxes over this period (+75% ± 28% during the late wet season). The unexpected negative drying effect in the wettest, infertile forest could have resulted from reduced vertical flushing of nutrients into soils. Including hydro‐nutrient interactions in ecosystem models could improve predictions of tropical forest‐climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Aircraft Observations of Turbulence in Cloudy and Cloud-Free Boundary Layers Over the Western North Atlantic Ocean From ACTIVATE and Implications for the Earth System Model Evaluation and Development

This study examines boundary layer turbulence derived from high temporal resolution meteorological measurements from 40 research flights over the western North Atlantic Ocean during the 2020 deployments of ACTIVATE. Frequency distributions of various turbulent quantities reveal stronger turbulence during the winter deployment than in summer and for cloud‐topped than in cloud‐free boundary layers during the summer deployment. Maximum turbulence kinetic energy (TKE) is most often within cloud from observations in winter and summer, whereas it is mostly below cloud in both seasons by a global model turbulence parameterization. Further, bivariate frequency distributions are consistent with the bivariate Gaussian probability distribution functions assumed for the closure of higher‐order turbulence/shallow convection parameterizations used by some global models. Turbulence simulated by the Community Atmosphere Model version 6 and the Energy Exascale Earth System Model Atmosphere Model version 2 using such parameterizations is not as strong as observed, with more TKE going into vertical wind perturbations rather than into zonal wind perturbations as observed, suggesting that the treatment of turbulence in Earth system models still needs to be further improved.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assessing Clouds Using Satellite Observations Through Three Generations of Global Atmosphere Models

Abstract Clouds are parameterized in climate models using quantities on the model grid‐scale to approximate the cloud cover and impact on radiation. Because of the complexity of processes involved with clouds, these parameterizations are one of the key challenges in climate modeling. Differences in parameterizations of clouds are among the main contributors to the spread in climate sensitivity across models. In this work, the clouds in three generations of an atmosphere model lineage are evaluated against satellite observations. Satellite simulators are used within the model to provide an appropriate comparison with individual satellite products. In some respects, especially the top‐of‐atmosphere cloud radiative effect, the models show generational improvements. The most recent generation, represented by two distinct branches of development, exhibits some regional regressions in the cloud representation; in particular the southern ocean shows a positive bias in cloud cover. The two branches of model development show how choices during model development, both structural and parametric, lead to different cloud climatologies. Several evaluation strategies are used to quantify the spatial errors in terms of the large‐scale circulation and the cloud structure. The Earth mover's distance is proposed as a useful error metric for the passive satellite data products that provide cloud‐top pressure‐optical depth histograms. The cloud errors identified here may contribute to the high climate sensitivity in the Community Earth System Model, version 2 and in the Energy Exascale Earth System Model, version 1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed‐Phase Arctic Stratus Cloud Layer

Abstract Mixed‐phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long‐lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice‐nucleating particle (INP) reservoir dynamics governing immersion‐mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol‐cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time‐independent (singular) number‐ and surface area‐based descriptions and a time‐dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT‐based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Representing Anthropogenic Dust From Agricultural Sources in E3SMv1: Implementation, Evaluation, and Assessment of the Radiative Forcing

Dust emissions related to anthropogenic activities (i.e., anthropogenic dust) is not represented in most global climate models and its radiative impact remains unassessed. In this study, we develop a new and physically based method to parameterize anthropogenic dust emission from agricultural sources (i.e., agricultural dust, or AD) based on the DOE's Energy Exascale Earth System Model version 1 (E3SMv1). This method relates AD emission to the crop land use fraction in the E3SMv1 land component. Major AD sources simulated by our parameterization include those over Central America, the Sahel, North India, and North China. The annual averaged AD emission with diameter less than 10 μm is 567 Tg yr −1 in present-day (year 2000), which contributes to 13.3% of total dust emission. Model evaluation against satellite and ground-based observations shows that the new parameterization can represent AD emissions and global dust cycle reasonably well. We find that the total dust emission increases by 13.1% (495 Tg yr −1 ) from 1850 to 2000 due to the cropland land use fraction changes without considering the impact from climate changes. This induces a net dust direct effective radiative forcing of −0.041 W m −2 at top of the atmosphere. This dust-induced cooling exceeds 10% of the total anthropogenic aerosol direct effective radiative forcing from 1750 to 2014 estimated by the Intergovernmental Panel on Climate Change Sixth Assessment Report. Our findings indicate an important role of anthropogenic dust in global climate change, which should be included in future climate change assessments.

agricultural dust↗
DOE OSTI2024

Improving River Routing Using a Differentiable Muskingum‐Cunge Model and Physics‐Informed Machine Learning

Recently, rainfall-runoff simulations in small headwater basins have been improved by methodological advances such as deep neural networks (NNs) and hybrid physics-NN models—particularly, a genre called differentiable modeling that intermingles NNs with physics to learn relationships between variables. However, hydrologic routing simulations, necessary for simulating floods in stem rivers downstream of large heterogeneous basins, had not yet benefited from these advances and it was unclear if the routing process could be improved via coupled NNs. We present a novel differentiable routing method (δMC-Juniata-hydroDL2) that mimics the classical Muskingum-Cunge routing model over a river network but embeds an NN to infer parameterizations for Manning's roughness (n) and channel geometries from raw reach-scale attributes like catchment areas and sinuosity. The NN was trained solely on downstream hydrographs. Synthetic experiments show that while the channel geometry parameter was unidentifiable, n can be identified with moderate precision. With real-world data, the trained differentiable routing model produced more accurate long-term routing results for both the training gage and untrained inner gages for larger subbasins (>2,000 km2) than either a machine learning model assuming homogeneity, or simply using the sum of runoff from subbasins. The n parameterization trained on short periods gave high performance in other periods, despite significant errors in runoff inputs. The learned n pattern was consistent with literature expectations, demonstrating the framework's potential for knowledge discovery, but the absolute values can vary depending on training periods. The trained n parameterization can be coupled with traditional models to improve national-scale hydrologic flood simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of a New Approach for Entrainment and Detrainment Rate Estimation

Entrainment and detrainment rates (ε and δ) constitute the most critical free parameters in mass flux schemes commonly employed for cumulus parameterizations. Recently, Zhu et al. (2021) introduced a new approach that utilizes aircraft observations to simultaneously estimate ε and δ for cumulus clouds, overcoming the limitation of other observation-based approaches that solely yield ε without offering insights into δ. This study aims to comprehensively evaluate the reliability of this new approach. First, evaluation using an Explicit Mixing Parcel Model demonstrates the capability of the new approach to back-calculate predetermined ε and δ based on the physical properties before and after the entrainment mixing. Second, evaluation using large-eddy simulations illustrates that the new approach yields consistent ε and δ profiles compared to the traditional approach. Sensitivity tests indicate a weak sensitivity of the estimated δ with the new approach to the entrained air source. A decrease in the proportion of cloudy air in the assumed detrained air leads to a reduction in the estimated δ, while ε remains unaffected. Finally, the most appropriate assumptions for entrained and detrained air are discussed. Estimating ε for cumulus parameterizations involves acquiring ambient air more than 500 m away from the cloud edge as entrained air. Due to implicit mean field approximations in the traditional approach, determining the optimal assumption for detrained air properties proves challenging. Finally, this study confirms the reliability of the new approach in estimating ε and δ, providing confidence in its application to extensive observational data and advancement in parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mechanisms in Regulating the Quasi‐Biennial Oscillation in Exascale Earth System Model Version 2

Abstract The update to the deep convection parameterization in the energy exascale Earth system model version 2 (E3SMv2) makes the simulated quasi‐biennial oscillation (QBO) have a shorter period than without the update when other tunable parameters the same since the update makes convection more intense but less frequent while leaving the time‐mean convective heating tendency almost unchanged. Amplitudes of momentum fluxes of parameterized gravity waves (GWs) are intensified since they are determined by the square of convective heating tendency. In addition, stronger planetary waves are simulated in E3SMv2 with the convection scheme update, partially due to the tropospheric precipitation change. Furthermore, there is evidence that planetary waves are intensified in the stratosphere due to the dissipation of enhanced parameterized GWs. These factors are found to be responsible for modulating the QBO simulated in E3SMv2.

Li, Yuanpu [Climate and Global Dynamics Division U↗
DOE OSTI2025

Ship‐Based Lidar Evaluation of Southern Ocean Low Clouds in the Storm‐Resolving General Circulation Model ICON and the ERA5 and MERRA‐2 Reanalyses

Global storm resolving models (GSRMs) represent the next generation of global climate models. One of them is a 5-km Icosahedral Nonhydrostatic Weather and Climate Model (ICON). Its high resolution means that parameterizations of convection and clouds, including subgrid-scale clouds, are omitted, relying on explicit simulation but necessarily utilizing microphysics and turbulence parameterizations. Standard-resolution (10–100 km) models, which use convection and cloud parameterizations, have substantial cloud biases over the Southern Ocean (SO), adversely affecting radiation and sea surface temperature. The SO is dominated by low clouds, which cannot be observed accurately from space due to overlapping clouds, attenuation, and ground clutter. We evaluated SO clouds in ICON and the ERA5 and MERRA-2 reanalyzes using approximately 2400 days of lidar observations and 2300 radiosonde profiles from 31 voyages and a Macquarie Island station during 2010–2021, compared to the model and reanalyzes using a ground-based lidar simulator. We found that ICON and the reanalyzes underestimate the total cloud fraction by about 10% and 20%, respectively. ICON and ERA5 overestimate the cloud occurrence peak at about 500 m, associated with underestimated lower tropospheric stability and overestimated lifting condensation level. The reanalyzes strongly underestimate fog and very low-level clouds, and MERRA-2 underestimates cloud occurrence at almost all heights. Outgoing shortwave radiation is overestimated in MERRA-2, implying a “too few, too bright” cloud problem. SO cloud and fog biases are a substantial issue in the analyzed model and reanalyzes and result in shortwave and longwave radiation biases.

Kuma, Peter [Stockholm Univ. (Sweden); Univ. of Ca↗
DOE OSTI2025

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape

Arctic warming is altering vegetation and carbon dynamics with global implications, yet Earth System Model (ESM) predictions in the Arctic remain highly uncertain, in part due to historically limited data for model parameterization and validation. As such, ESMs typically represent Arctic ecosystems in an oversimplified manner. Recently, nine plant functional types (PFTs) designed to realistically represent tundra vegetation were integrated into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) and parameterized using plot-scale observations from a single site. Additional evaluation was needed to determine their transferability across the Arctic. Here, in this study, we evaluated whether refined representation of tundra vegetation improved model accuracy by conducting spatially explicit 100 × 100 m resolution ELM simulations on Alaska's Seward Peninsula. Simulations with the default two-PFT configuration and with the nine Arctic-specific PFTs were benchmarked against observations of net ecosystem exchange, gross primary production, and aboveground biomass from multiple data streams including an eddy covariance flux tower, flux chambers, and aircraft and unoccupied aerial system hyperspectral remote sensing. Evaluation revealed that Arctic-specific PFT simulations produced more realistic landscape-level carbon exchanges, and better captured observed heterogeneity in biomass and productivity, explaining 60%–70% of spatial variance (R 2 = 0.6–0.7) compared to just 12%–18% (R 2 = 0.12–0.18) with the default configuration. However, the refined model failed to reproduce observed aboveground biomass for highly productive alder-willow communities, requiring further evaluation of carbon allocation parameterizations for tall shrubs that are increasingly expanding across tundra landscapes. Our results demonstrate that enhanced representation of vegetation heterogeneity boosts predictive understanding of tundra carbon dynamics, facilitating regional to pan-Arctic model and remote-sensing scaling.

Murphy, Bailey A. [Oak Ridge National Laboratory (↗
DOE OSTI2018

Anthropogenic influences on major tropical cyclone events

There is no consensus on whether climate change has yet affected the statistics of tropical cyclones, owing to their large natural variability and the limited period of consistent observations. In addition, projections of future tropical cyclone activity are uncertain, because they often rely on coarse-resolution climate models that parameterize convection and hence have difficulty in directly representing tropical cyclones. Here we used convection-permitting regional climate model simulations to investigate whether and how recent destructive tropical cyclones would change if these events had occurred in pre-industrial and in future climates. We found that, relative to pre-industrial conditions, climate change so far has enhanced the average and extreme rainfall of hurricanes Katrina, Irma and Maria, but did not change tropical cyclone wind-speed intensity. In addition, future anthropogenic warming would robustly increase the wind speed and rainfall of 11 of 13 intense tropical cyclones (of 15 events sampled globally). Additional regional climate model simulations suggest that convective parameterization introduces minimal uncertainty into the sign of projected changes in tropical cyclone intensity and rainfall, which allows us to have confidence in projections from global models with parameterized convection and resolution fine enough to include tropical cyclones.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

VR-SCOSMO: A smooth conductor-like screening model with charge-dependent radii for modeling chemical reactions

To better represent the solvation effects observed along reaction pathways, and of ionic species in general, a charge-dependent variable-radii smooth conductor-like screening model (VR-SCOSMO) is developed. This model is implemented and parameterized with a third order density-functional tight binding quantum model, DFTB3/3OB-OPhyd, a quantum method which was developed for organic and biological compounds, utilizing a specific parameterization for phosphate hydrolysis reactions. Unlike most other applications with the DFTB3/3OB model, an auxiliary set of atomic multipoles is constructed from the underlying DFTB3 density matrix which is used to interact the solute with the solvent response surface. The resulting method is variational, produces smooth energies, and has analytic gradients. As a baseline, a conventional SCOSMO model with fixed radii is also parameterized. The SCOSMO and VR-SCOSMO models shown have comparable accuracy in reproducing neutral-molecule absolute solvation free energies; however, the VR-SCOSMO model is shown to reduce the mean unsigned errors (MUEs) of ionic compounds by half (about 2-3 kcal/mol). The VR-SCOSMO model presents similar accuracy as a charge-dependent Poisson-Boltzmann model introduced by Hou et al. [J. Chem. Theory Comput. 6, 2303 (2010)]. VR-SCOSMO is then used to examine the hydrolysis of trimethylphosphate and seven other phosphoryl transesterification reactions with different leaving groups. Two-dimensional energy landscapes are constructed for these reactions and calculated barriers are compared to those obtained from ab initio polarizable continuum calculations and experiment. Results of the VR-SCOSMO model are in good agreement in both cases, capturing the rate-limiting reaction barrier and the nature of the transition state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Transferable density functional tight binding for carbon, hydrogen, nitrogen, and oxygen: Application to shock compression

A new parameterization for density functional tight binding (DFTB) theory, lanl31, has been developed for molecules containing carbon, hydrogen, nitrogen, and oxygen. Optimal values for the Hubbard Us, on-site energies, and the radial dependences of the bond integrals and repulsive potentials were determined by numerical optimization using simulated annealing to a modest database of ab initio-calculated atomization energies and interatomic forces. The transferability of the optimized DFTB parameterization has been assessed using the CHNO subset of the QM-9 database [R. Ramakrishnan et al., Sci. Data, 1, 140022 (2014)]. These analyses showed that the errors in the atomization energies and interatomic forces predicted by our model are small and in the vicinity of the di erences between density functional theory calculations with di erent basis sets and exchange-correlation functionals. Good correlations between the molecular dipole moments and HOMO-LUMO gaps predicted by lanl31 and the QM-9 data set are also found. Furthermore, the errors in the atomization energies and forces derived from lanl31 are signi cantly smaller than those obtained from the ReaxFF-lg reactive force eld for organic materials [L. Liu et al., J. Phys. Chem. A, 115, 11016 (2011)]. The lanl31 DFTB parameterization for C, H, N, and O has been applied the to the molecular dynamics simulation of the principal Hugoniot of liquid nitromethane, liquid benzene, liquid nitrogen, pentaerythritol tetranitrate, trinitrotoluene, and cyclotetramethylene tetranitramine. The computed and measured Hugoniot loci are in excellent agreement with experiment and we discuss the sensitivity of the loci to the underestimated shock heating that is a characteristic of classical molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Machine Learned Hückel Theory: Interfacing Physics and Deep Neural Networks

The Hückel Hamiltonian is an incredibly simple tight-binding model known for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these empirical parameters with machine-learned dynamic values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability, while the deep neural network parameterization is smooth and accurate and reproduces insightful features of the original empirical parameterization. Altogether, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Transition from ITG to MTM linear instabilities near pedestals of high density plasmas

Here, investigation of linear gyrokinetic ion-scale modes ( k θ ρ s = 0.3) finds that a transition from ion temperature gradient to microtearing mode (MTM) dominance occurs as the density is increased near the pedestal region of a parameterized DIII-D sized tokamak. H-modes profile densities, temperatures, and equilibria are parameterized utilizing the OMFIT PRO_create module. With these profiles, linear gyrokinetic ion-scale instabilities are predicted with CGYRO. This transition (nMTM) has a weak dependence on radial location in the region near the top of the pedestal ( ρ = 0.7 – 0.9), which allows simulating single radii to examine the approximate scaling of nMTM with global parameters. The critical nMTM is found to scale with plasma current. Additionally, increasing the minor radius by decreasing the aspect ratio and increasing the major radius are found to reduce nMTM. However, any relationship between nMTM and density limit physics remains unclear as nMTM increases relative to the Greenwald density with larger minor radius and with larger magnetic field, suggesting that the transport due to MTM may be less important for a reactor. Additionally, nMTM is sensitive to the pedestal temperature, the local electron and ion gradients, the ratio of ion to electron temperature T i / T e, and the current profile. MTMs are predicted to be the dominant instability in the core at similar Greenwald fractions for DIII-D, NSTX, and NSTX-U H-mode experiments, supporting the results of the parameterized study. Additionally, MTMs continue to be the dominant linear instability in a DIII-D L-mode after an H–L transition as the plasma approaches a density limit disruption despite the large change in plasma profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Speeding up low-mass planetary microlensing simulations and modeling: The caustic region of influence

Extensive simulations of planetary microlensing are necessary both before and after a survey is conducted: before to design and optimize the survey and after to understand its detection efficiency. The major bottleneck in such computations is the computation of light curves. However, for low-mass planets, most of these computations are wasteful, as most light curves do not contain detectable planetary signatures. In this paper, I develop a parameterization of the binary microlens that is conducive to avoiding light curve computations. I empirically find analytic expressions describing the limits of the parameter space that contain the vast majority of low-mass planet detections. Through a large-scale simulation, I measure the (in)completeness of the parameterization and the speed-up it is possible to achieve. For Earth-mass planets in a wide range of orbits, it is possible to speed up simulations by a factor of ∼30-125 (depending on the survey's annual duty-cycle) at the cost of missing ∼1% of detections (which is actually a smaller loss than for the arbitrary parameter limits typically applied in microlensing simulations). The benefits of the parameterization probably outweigh the costs for planets below 100 M{sub ⊕}. For planets at the sensitivity limit of AFTA-WFIRST, simulation speed-ups of a factor ∼1000 or more are possible.

79 ASTRONOMY AND ASTROPHYSICS↗