Search NASA⌕ Search

Engineering topics

Wan, Hui

Publications and source records attributed to Wan, Hui.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quantifying structural errors in cloud condensation nuclei activity from reduced representation of aerosol size distributions

Aerosol effects on clouds and radiation are the dominant contribution to uncertainty in radiative forcing relative to the pre-industrial atmosphere. While previous studies have assessed the impact of parametric uncertainty on modeled forcing, structural errors from the numerical representation of particle distributions have not been well quantified. Here we present a framework for quantifying error in aerosol size distributions and cloud condensation nuclei activity, which we apply to the widely used 4-mode version of the Modal Aerosol Module (MAM4). Box model predictions from the MAM4 are evaluated against the Particle Monte Carlo Model for Simulating Aerosol Interactions and Chemistry (PartMC-MOSAIC), a benchmark model that tracks the evolution of individual particles. We show that size distributions simulated by MAM4 diverge from those simulated by PartMC-MOSAIC after only a few hours of aging by condensation and coagulation in polluted conditions, which leads to large errors in modeled cloud condensation nuclei concentrations. We find that differences between MAM4 and PartMC-MOSAIC are largest under polluted conditions, where the size distribution evolves rapidly though aging by condensation of semi-volatile substances and coagulation among particles. These findings suggest that structural error in modeled aerosol properties contributes to the large inter-model variability in aerosol radiative forcing.

Fierce, Laura M.↗
DOE OSTI2023

E3SMv1 simulation output from sensitivity study on aerosol process coupling

This dataset contains annual mean and 6-hourly output data from a set of E3SMv1 simulations that demonstrates the impact of numerical process coupling on the simulated dust life cycle. The simulations are documented and analyzed in the attached manuscript (pdf file). Publishing the manuscript in a peer-reviewed journal will require making the model output publicly available. We plan to use DataHub for that. The data files to be uploaded are currently saved on Compy under /compyfs/wanh895/scidac_share/cflx_2023/

Wan, Hui↗
DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Incorporating the Effect of Large‐Scale Vertical Motion on Convection Through Convective Mass Flux Adjustment in E3SMv2

Abstract Recent observational studies suggest that the large‐scale dynamical forcing (vertical motion) plays important roles in deep convection development. In this study we propose a convective mass flux adjustment (MAdj) approach to represent the dynamical effects of large‐scale vertical motion on convection in the Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2). With MAdj, convection is enhanced (suppressed) when there is large‐scale ascending (descending) motion at the planetary boundary layer top. The coupling of convection with large‐scale circulation significantly improves the simulation of climate variability in E3SMv2 across multiple scales from the diurnal cycle, convectively coupled equatorial waves, to the Madden‐Julian Oscillation (MJO). The standard E3SMv2 tends to simulate overly weak diurnal amplitude of precipitation and overly weak variance of convectively coupled equatorial Kelvin and westward inertio‐gravity (WIG) waves. It also fails to simulate the essential characteristics of the MJO: continuous eastward propagation. With MAdj, the amplitude of diurnal cycle of precipitation is systematically increased and its probability density distribution is much closer to observations. The MAdj can also simulate more realistic eastward propagation of the MJO and much stronger convectively coupled Kelvin and WIG waves. Moreover, the MAdj approach slightly improves the climatology simulations in precipitation, cloud, radiation, circulation, temperature, and moisture fields, with overall root‐mean‐square error (RMSE) of major climatological fields reduced by about 2%. The MAdj approach suppresses excessive grid‐scale precipitation, reducing precipitation wet biases over South China Sea, Philippine Sea, Himalayas, and South Pacific Convergence Zone in western Pacific in summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

In Situ Visualization of the Pinning Effect of Planar Defects on Li Ion Insertion

Longevity of Li ion batteries strongly depends on the interaction of transporting Li ions in electrode crystals with defects. However, detailed interactions between the Li ion flux and structural defects in the host crystal remain obscure due to the transient nature of such interactions. Here, by in situ transmission electron microscopy and density function theory calculations, we reveal how the diffusion pathways and transport kinetics of a Li ion can be affected by planar defects in a tungsten trioxide lattice. We uncover that changes in charge distribution and lattice spacing along the planar defects disrupt the continuity of ion conduction channels and dramatically increase the energy barrier of Li diffusion, thus, arresting Li ions at the defect sites and twisting the lithiation front. In conclusion, the atomic-scale understanding holds critical implications for rational interface design in solid-state batteries and solid oxide fuel cells.

25 ENERGY STORAGE↗
DOE OSTI2023

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e., clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗
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 OSTI2022

Effective radiative forcing of anthropogenic aerosols in E3SMv1: historical changes, causality, decomposition, and parameterization sensitivities

The effective radiative forcing of anthropogenic aerosols (ERF aer ) is an important measure of the anthropogenic aerosol effects simulated by a global climate model. Here we analyze ERF aer simulated by the E3SMv1 atmosphere model using both century-long free-running atmosphere-land simulations and short nudged simulations. We relate the simulated ERF aer to characteristics of the aerosol composition and optical properties, and evaluate the relationships between key aerosol and cloud properties. In terms of historical changes from the year 1870 to 2014, our results show that the global mean anthropogenic aerosol burden and optical depth increase during the simulation period as expected, but the regional averages show large differences in the temporal evolution. The largest regional differences are found in the emission-induced evolution of the burden and optical depth of the sulfate aerosol: a strong decreasing trend is seen in the Northern Hemisphere high-latitude region after around 1970, while a continued increase is simulated in the tropics. Consequently, although the global mean anthropogenic aerosol burden and optical depth increase from 1870 to 2014, the ERF aer magnitude does not increase after around year 1970. The relationships between key aerosol and cloud properties (relative changes between preindustrial and present-day conditions) also show evident changes after 1970, diverging from the linear relationships exhibited for the period from 1870 to 2014. The ERF aer in E3SMv1 is relatively large compared to the recently published multi-model estimates; the primary reason is the large indirect aerosol effect (i.e., through aerosol-cloud interactions). Compared to other models, E3SMv1 features a stronger sensitivity of the cloud droplet effective radius to changes in the cloud droplet number concentration. Large sensitivity is also seen in the liquid cloud optical depth, which is determined by changes in both the effective radius and liquid water path. Aerosol-induced changes in liquid and ice cloud properties in E3SMv1 are found to have a strong correlation, as the evolution of anthropogenic sulfate aerosols affects both the liquid cloud formation and the homogeneous ice nucleation in cirrus clouds. The ERF aer estimates in E3SMv1 for the shortwave and longwave components are sensitive to the parameterization changes in both liquid and ice cloud processes. When the parameterization of ice cloud processes is modified, the top-of-atmosphere forcing changes in the shortwave and longwave components largely offset each other, so the net effect is negligible. This suggests that, to reduce the magnitude of the net ERF aer , it would be more effective to reduce the anthropogenic aerosol effect through liquid or mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗
DOE OSTI2022

HAERO: High-performance aerosols

SAND2021-15059 O Haero is a software library for representing aerosol size distributions, chemical species, and parameterized processes to use with a coupled atmospheric dynamics model or as a standalone tool. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Bosler, Peter↗
DOE OSTI2021

New efforts in improving numerical aspects of atmospheric physics parameterizations

This is an invited article for GEWEX Quarterly, the newsletter of the Global Energy and Water Exchanges (GEWEX) program which is a core project of the World Climate Research Programme (https://www.gewex.org/resources/gewex-news/). The article introduces recent research results and ongoing activities being carried out in a SciDAC project funded by DOE that focuses on identifying and addressing numerical issues in physics parameterizations in global atmospheric general circulation models.

global atmospheric model↗
DOE OSTI2021

Earth System Model Development - Energy Exascale Earth System Model (2020 Principal Investigator Meeting Report)

As part of the U.S. Department of Energy (DOE) Office of Biological and Environmental Research (BER), the Earth System Model Development (ESMD) program area supports significant research in Earth system model development, including major efforts in the Energy Exascale Earth System Model (E3SM) project. Convened virtually on October 26–29, the 2020 ESMD-E3SM Principal Investigator (PI) Meeting provided a forum to survey ESMD sponsored research, recent progress made by the ESMD science team and the broader Earth system modeling community within and outside of DOE, foster interactions across ESMD projects, and discuss research gaps and opportunities for future advancements.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

The challenges of theory-software translation

Background: Software is now ubiquitous within research. In addition to the general challenges common to all software development projects, research software must also represent, manipulate, and provide data for complex theoretical constructs. Ensuring this process of theory-software translation is robust is essential to maintaining the integrity of the science resulting from it, and yet there has been little formal recognition or exploration of the challenges associated with it. Methods: We thematically analyse the outputs of the discussion sessions at the Theory-Software Translation Workshop 2019, where academic researchers and research software engineers from a variety of domains, and with particular expertise in high performance computing, explored the process of translating between scientific theory and software. Results: We identify a wide range of challenges to implementing scientific theory in research software and using the resulting data and models for the advancement of knowledge. We categorise these within the emergent themes of design, infrastructure, and culture, and map them to associated research questions. Conclusions: Systematically investigating how software is constructed and its outputs used within science has the potential to improve the robustness of research software and accelerate progress in its development. We propose that this issue be examined within a new research area of theory-software translation, which would aim to significantly advance both knowledge and scientific practice.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Improving Time Step Convergence in an Atmosphere Model With Simplified Physics: The Impacts of Closure Assumption and Process Coupling

Convergence testing is a common practice in the development of dynamical cores of atmospheric models but is not as often exercised for the parameterization of sub-grid physics. An earlier study revealed that the stratiform cloud parameterizations in several predecessors of the Energy Exascale Earth System Model (E3SM) showed strong time-step sensitivity and slower-than-expected convergence when the model's time step was systematically refined. In this work, a simplified atmosphere model is configured that consists of the spectral-element dynamical core of the E3SM atmosphere model coupled with a large-scale condensation parameterization based on commonly used assumptions. This simplified model also resembles E3SM and its predecessors in the numerical implementation of process coupling and shows poor time-step convergence in short ensemble tests. We present a formal error analysis to reveal the expected time-step convergence rate and the conditions for obtaining such convergence. Numerical experiments are conducted to investigate the root causes of convergence problems. We show that revisions in the process coupling and closure assumption help to improve convergence in short simulations using the simplified model; the same revisions applied to a full atmosphere model lead to significant changes in the simulated long-term climate. This work demonstrates that causes of convergence issues in atmospheric simulations can be understood by combining analyses from physical and mathematical perspectives. Addressing convergence issues can help to obtain a discrete model that is more consistent with the intended representation of the physical phenomena.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Improving Solution Accuracy and Convergence for Stochastic Physics Parameterizations with Colored Noise

Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the Itô correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the Itô correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the Itô correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized Itô correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized Itô correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation.

54 ENVIRONMENTAL SCIENCES↗