Search NASA⌕ Search

Engineering topics

Zhang, Tao

Publications and source records attributed to Zhang, Tao.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Direct numerical simulations of activation and deactivation in turbulent atmospheric clouds

Significant knowledge gaps remain in our understanding of turbulence–cloud–aerosol interactions in the Earth's atmosphere, and direct numerical simulation (DNS) has increasingly become an indispensable tool to fill such critical knowledge gaps. Here, this study is an extension of our previous DNS model [Gao et al., J. Geophys. Res.: Atmos., 123(4), 2194–2214 (2018)], with a focus on the activation of aerosol particles into cloud droplets and deactivation of cloud droplets into aerosol particles in a microscale cloud environment. The effects of turbulence intensity, particle curvature, and solute, as well as the initial distributions of the aerosol particles (monodisperse or polydisperse) are investigated. The governing equations for the flow of air, temperature, and water vapor mixing ratio are solved numerically in the Eulerian fashion, assuming homogeneous and isotropic turbulence. The dynamics of the aerosol and cloud particles are calculated with the Lagrangian particle tracking method. The results show that the deviations of the thermodynamic variables from their respective means are significantly reduced, the activation process appears to be delayed, and the deactivation process occurs more rapidly, as the turbulence intensity is increased. The inclusion of particle curvature and solute effects, as well as polydispersity, tends to retard the activation of aerosols into cloud droplets. It is also observed that fluctuations in supersaturation broaden the spread of particle radii, and the broadening is followed by a narrowing as turbulent homogenization reduces thermodynamic fluctuations over time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Convergence of Cloud Droplet Spectral Relative Dispersion During Entrainment‐Mixing Based on Particle‐Resolved Direct Numerical Simulations

Entrainment-mixing processes critically impact cloud microphysical properties, but their effects on the relative dispersion (d) of cloud droplet size distributions (CDSDs) remain elusive. A direct numerical simulation model is initialized with different CDSDs to fill the gap. These results show that d decreases for broad CDSDs and increases for narrow ones, ultimately converging to approximately 0.5 regardless of initial CDSDs during the evaporation-dominated mixing stage. The supersaturation fluctuation and the shape of CDSDs jointly influence the convergence behavior of d. Further sensitivity tests show that the initial microphysical/dynamical/thermodynamical conditions exert negligible effects on the final converged value of d but affect the convergence rate (k). The k generally increases with increasing droplet number concentration and dissipation rate, and increases with decreasing liquid water content, relative humidity of entrained air, and mixing fraction of cloudy air. A conceptual model with two timescales is proposed; k and the timescales are negatively correlated, meaning that slow mixing and/or evaporation process results in slow convergence of d. In conclusion, this finding provides an important reference for improving understanding and parameterization of d during the entrainment-mixing processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Causal Directions Matter: How Environmental Factors Drive Convective Cloud Detrainment Heights

This study investigates how environmental factors influence the level of maximum detrainment (LMD) in deep convective clouds. Through a novel application of the Linear Non‐Gaussian Acyclic Model (LiNGAM), we discover causal structures between environmental variables and LMD, observed at six tropical sites operated by the Atmospheric Radiation Measurement (ARM) user facility. LiNGAM effectively identifies causal directions among variables of interest, revealing robust relationships such as those among the lifting condensation level (LCL), level of free convection (LFC), and convective inhibition (CIN), aligning with prior knowledge. Relative humidity is shown to directly influence LMD; however, this relationship exhibits strong nonlinearity and becomes difficult to detect when the contrast between oceanic and continental environments is excluded from the analysis. This study highlights the importance of establishing causal relationships before performing statistical inference.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Search for Solar Boosted Dark Matter Particles at the PandaX-4T Experiment

We present a novel constraint on light dark matter utilizing 1.54 metric ton/year of data acquired from the PandaX-4T dual-phase xenon time projection chamber. This constraint is derived through detecting electronic recoil signals resulting from the interaction with solar-enhanced dark matter flux. Low-mass dark matter particles, lighter than a few MeV / c 2 , can scatter with the thermal electrons in the Sun. Consequently, with higher kinetic energy, the boosted dark matter component becomes detectable via contact scattering with xenon electrons, resulting in a few keV energy deposition that exceeds the threshold of PandaX-4T. We calculate the expected recoil energy in PandaX-4T considering the Sun’s acceleration with heavy mediators and the detection capabilities of the xenon detector. The first experimental search results using the xenon detector yield the most stringent upper limits cross section of 3.51 × 10 − 39 cm 2 at 0.08 MeV / c 2 for a solar boosted dark matter mass ranging from 0.02 to 10 MeV / c 2 , achieving a 23-fold improvement compared with earlier experimental studies. Published by the American Physical Society 2025

Shen, Guofang↗
DOE OSTI2025

A Fortran–Python interface for integrating machine learning parameterization into earth system models

Abstract. Parameterizations in earth system models (ESMs) are subject to biases and uncertainties arising from subjective empirical assumptions and incomplete understanding of the underlying physical processes. Recently, the growing representational capability of machine learning (ML) in solving complex problems has spawned immense interests in climate science applications. Specifically, ML-based parameterizations have been developed to represent convection, radiation, and microphysics processes in ESMs by learning from observations or high-resolution simulations, which have the potential to improve the accuracies and alleviate the uncertainties. Previous works have developed some surrogate models for these processes using ML. These surrogate models need to be coupled with the dynamical core of ESMs to investigate the effectiveness and their performance in a coupled system. In this study, we present a novel Fortran–Python interface designed to seamlessly integrate ML parameterizations into ESMs. This interface showcases high versatility by supporting popular ML frameworks like PyTorch, TensorFlow, and scikit-learn. We demonstrate the interface's modularity and reusability through two cases: an ML trigger function for convection parameterization and an ML wildfire model. We conduct a comprehensive evaluation of memory usage and computational overhead resulting from the integration of Python codes into the Fortran ESMs. By leveraging this flexible interface, ML parameterizations can be effectively developed, tested, and integrated into ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

STFM: Accurate Spatio-Temporal Fusion Model for Weather Forecasting

Meteorological prediction is crucial for various sectors, including agriculture, navigation, daily life, disaster prevention, and scientific research. However, traditional numerical weather prediction (NWP) models are constrained by their high computational resource requirements, while the accuracy of deep learning models remains suboptimal. In response to these challenges, we propose a novel deep learning-based model, the Spatiotemporal Fusion Model (STFM), designed to enhance the accuracy of meteorological predictions. Our model leverages Fifth-Generation ECMWF Reanalysis (ERA5) data and introduces two key components: a spatiotemporal encoder module and a spatiotemporal fusion module. The spatiotemporal encoder integrates the strengths of convolutional neural networks (CNNs) and recurrent neural networks (RNNs), effectively capturing both spatial and temporal dependencies. Meanwhile, the spatiotemporal fusion module employs a dual attention mechanism, decomposing spatial attention into global static attention and channel dynamic attention. This approach ensures comprehensive extraction of spatial features from meteorological data. The combination of these modules significantly improves prediction performance. Experimental results demonstrate that STFM excels in extracting spatiotemporal features from reanalysis data, yielding predictions that closely align with observed values. In comparative studies, STFM outperformed other models, achieving a 7% improvement in ground and high-altitude temperature predictions, a 5% enhancement in the prediction of the u/v components of 10 m wind speed, and an increase in the accuracy of potential height and relative humidity predictions by 3% and 1%, respectively. This enhanced performance highlights STFM’s potential to advance the accuracy and reliability of meteorological forecasting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Wetting mechanism and alteration of nano-sized shale pores: Insights from contrast variation small angle neutron scattering

Wettability of tight shale is crucial for fluid flow and mass transport process in energy geosciences. However, understanding the interfacial chemistry and wetting mechanisms at sub-nano-pore scales remains a formidable challenge. Here, in this study, the Contrast Variation technique of Small Angle Neutron Scattering (CV-SANS) is employed to investigate shale’s interfacial chemistry using reagents that possess a range of different polarities, including water, n-decane, toluene, and dimethyl methanamide. Through five different experimental strategies, we have demonstrated a successful modification of shale wettability, ranging from enhancement, weakening, to reversal. Delving into the mechanisms, we illustrated the crucial role of pre-existing liquid films in these changes, where the uniquely co-existing polar and non-polar functional groups in dimethyl methanamide acted as a conduit for interfacial chemistry adjustments. Furthermore, a solvent immersion led to matrix dilation as well as liberation of residual oil-occupied pores, resulting in altered pore size distributions, with hydrogen bonding playing a significant role in the polar groups. Interestingly, despite shale exhibiting a stronger affinity for oil over water, hydrophilic solvents induced more substantial dilation than lipophilic ones. Collectively, this work elucidates the dynamic change of interfacial chemistry via the configuration of polarity using chemical reagents, and the CV-SANS technique underscores its invaluable utilities in decoding the interfacial wettability traits in nanopore space of shale.

58 GEOSCIENCES↗
DOE OSTI2024

Improving North American Wildfire Prediction by Integrating a Machine-Learning Fire Model in a Land Surface Model

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM) version 2.1. A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Towards accelerating particle-resolved direct numerical simulation with neural operators

In this paper, we present our ongoing work aimed at accelerating a particle-resolved direct numerical simulation model designed to study aerosol–cloud–turbulence interactions. The dynamical model consists of two main components—a set of fluid dynamics equations for air velocity, temperature, and humidity, coupled with a set of equations for particle (i.e., cloud droplet) tracing. Rather than attempting to replace the original numerical solution method in its entirety with a machine learning (ML) method, we consider developing a hybrid approach. We exploit the potential of neural operator learning to yield fast and accurate surrogate models and, in this study, develop such surrogates for the velocity and vorticity fields. We discuss results from numerical experiments designed to assess the performance of ML architectures under consideration as well as their suitability for capturing the behavior of relevant dynamical systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Emulator of PR‐DNS: Accelerating Dynamical Fields With Neural Operators in Particle‐Resolved Direct Numerical Simulation

Abstract Particle‐resolved direct numerical simulations (PR‐DNS) play an increasing role in investigating aerosol‐cloud‐turbulence interactions at the most fundamental level of processes. However, the high computational cost associated with high resolution simulations poses considerable challenges for large domain or long duration simulation using PR‐DNS. To address these issues, here we present an emulator of the complex physics‐based PR‐DNS developed by use of the data‐driven Fourier Neural Operator (FNO) method. The effectiveness of the method is showcased by presenting turbulence and temperature fields in a two‐dimensional space. The results demonstrate high accuracy at various resolutions and the emulator is two orders of magnitude cheaper in terms of computational demand compared to the physics‐based PR‐DNS model. Furthermore, the FNO emulator exhibits strong generalization capabilities for different initial conditions and ultra‐high‐resolution without the need to retrain models. These findings highlight the potential of the FNO method as a promising tool to simulate complex fluid dynamics problems with high accuracy, computational efficiency, and generalization capabilities, enhancing our understanding of the aerosol‐cloud‐precipitation system.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Review of Carbon Support Coordination Environments for Single Metal Atom Electrocatalysts (SACS)

This topical review focuses on the distinct role of carbon support coordination environment of single-atom catalysts (SACs) for electrocatalysis. The article begins with an overview of atomic coordination configurations in SACs, including a discussion of the advanced characterization techniques and simulation used for understanding the active sites. A summary of key electrocatalysis applications is then provided. These processes are oxygen reduction reaction (ORR), oxygen evolution reaction (OER), hydrogen evolution reaction (HER), nitrogen reduction reaction (NRR), and carbon dioxide reduction reaction (CO 2 RR). The review then shifts to modulation of the metal atom-carbon coordination environments, focusing on nitrogen and other non-metal coordination through modulation at the first coordination shell and modulation in the second and higher coordination shells. Representative case studies are provided, starting with the classic four-nitrogen-coordinated single metal atom (M$-$N 4 ) based SACs. Bimetallic coordination models including homo-paired and hetero-paired active sites are also discussed, being categorized as emerging approaches. The theme of the discussions is the correlation between synthesis methods for selective doping, the carbon structure–electron configuration changes associated with the doping, the analytical techniques used to ascertain these changes, and the resultant electrocatalysis performance. In conclusion, critical unanswered questions as well as promising underexplored research directions are identified.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data driven pathway analysis and forecast of global warming and sea level rise

Climate change is a critical issue of our time, and its causes, pathways, and forecasts remain a topic of broader discussion. In this paper, we present a novel data driven pathway analysis framework to identify the key processes behind mean global temperature and sea level rise, and to forecast the magnitude of their increase from the present to 2100. Based on historical data and dynamic statistical modeling alone, we have established the causal pathways that connect increasing greenhouse gas emissions to increasing global mean temperature and sea level, with its intermediate links encompassing humidity, sea ice coverage, and glacier mass, but not for sunspot numbers. Our results indicate that if no action is taken to curb anthropogenic greenhouse gas emissions, the global average temperature would rise to an estimated 3.28 °C (2.46–4.10 °C) above its pre-industrial level while the global sea level would be an estimated 573 mm (474–671 mm) above its 2021 mean by 2100. However, if countries adhere to the greenhouse gas emission regulations outlined in the 2021 United Nations Conference on Climate Change (COP26), the rise in global temperature would lessen to an average increase of 1.88 °C (1.43–2.33 °C) above its pre-industrial level, albeit still higher than the targeted 1.5 °C, while the sea level increase would reduce to 449 mm (389–509 mm) above its 2021 mean by 2100.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Physics-Informed Deep Learning for Reconstruction of Spatial Missing Climate Information in the Antarctic

Understanding the influence of the Antarctic on the global climate is crucial for the prediction of global warming. However, due to very few observation sites, it is difficult to reconstruct the rational spatial pattern by filling in the missing values from the limited site observations. To tackle this challenge, regional spatial gap-filling methods, such as Kriging and inverse distance weighted (IDW), are regularly used in geoscience. Nevertheless, the reconstructing credibility of these methods is undesirable when the spatial structure has massive missing pieces. Inspired by image inpainting, we propose a novel deep learning method that demonstrates a good effect by embedding the physics-aware initialization of deep learning methods for rapid learning and capturing the spatial dependence for the high-fidelity imputation of missing areas. We create the benchmark dataset that artificially masks the Antarctic region with ratios of 30%, 50% and 70%. The reconstructing monthly mean surface temperature using the deep learning image inpainting method RFR (Recurrent Feature Reasoning) exhibits an average of 63% and 71% improvement of accuracy over Kriging and IDW under different missing rates. With regard to wind speed, there are still 36% and 50% improvements. In particular, the achieved improvement is even better for the larger missing ratio, such as under the 70% missing rate, where the accuracy of RFR is 68% and 74% higher than Kriging and IDW for temperature and also 38% and 46% higher for wind speed. In addition, the PI-RFR (Physics-Informed Recurrent Feature Reasoning) method we proposed is initialized using the spatial pattern data simulated by the numerical climate model instead of the unified average. Compared with RFR, PI-RFR has an average accuracy improvement of 10% for temperature and 9% for wind speed. When applied to reconstruct the spatial pattern based on the Antarctic site observations, where the missing rate is over 90%, the proposed method exhibits more spatial characteristics than Kriging and IDW.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Tuning the Intermolecular Electrostatic Interaction toward High-Efficiency and Low-Cost Organic Solar Cells

Organic solar cells (OSCs) have achieved much progress with rapidly increasing power conversion efficiencies (PCEs). It should be noted that the top-performance OSCs are generally consisted of active materials with complex chemical structures, resulting in high costs. Here, combining the material design and morphology control, high-efficiency OSCs are fabricated by a low-cost donor: acceptor blend. A completely non-fused electron acceptor named Tz is designed and synthesized via introducing thiazole units on both sides of a bithiophene core, which shows an outstanding PCE of 13.3% with a typical polythiophene donor. More importantly, optimization guidelines are presented to get excellent morphology for low-cost donor:acceptor systems. Three polythiophenes are selected, poly(3-hexylthiophene) and its two derivatives with electron-withdrawing substitutions (PDCBT and PDCBT-2F), as donors to fabricate the cell devices. The computational and experimental data reveal that decreasing the electrostatic interaction between polythiophene and Tz is the key to getting a suppressed miscibility and thus a high phase purity. Finally, this study provides insight into the molecular design and donor:acceptor matching requirements for high-efficiency and low-cost OSCs.

14 SOLAR ENERGY↗
DOE OSTI2023

Trace level of atomic copper in N-doped graphene quantum dots switching the selectivity from C 1 to C 2 products in CO electroreduction

To unravel the relationship between trace Cu on the metal-free catalysts toward CO/CO 2 reduction reaction (CO/CO 2 RR), we investigated the effect of trace Cu loading in N-doped graphene quantum dots (NGQDs) on CO/CO 2 RR. A general trend is that increasing the Cu loading in NGQDs switches the selectivity from C 1 (CH 4 ) to C 2 products in CORR. When 2.5 μg/cm 2 Cu with the atomic size is loaded on NGQDs, the selectivity shifts from 62% Faradaic efficiency (FE) of CH 4 to 52% FE of C 2 products in CORR. Further increasing the atomic Cu loading to 3.8 μg/cm 2 promotes the FE of C 2 products to 78%. CO 2 RR requires one order of magnitude higher Cu loading than CORR to switch the selectivity from C 1 to C 2 products due to the low partial pressure of CO. Finally, this study clarifies the distinct impact of trace (ppm level) Cu on the activity/selectivity between CORR and CO 2 RR.

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

Realization of thousand-second improved confinement plasma with Super I-mode in Tokamak EAST

Mastering nuclear fusion, which is an abundant, safe, and environmentally competitive energy, is a great challenge for humanity. Tokamak represents one of the most promising paths toward controlled fusion. Obtaining a high-performance, steady-state, and long-pulse plasma regime remains a critical issue. Recently, a big breakthrough in steady-state operation was made on the Experimental Advanced Superconducting Tokamak (EAST). A steady-state plasma with a world-record pulse length of 1056 s was obtained, where the density and the divertor peak heat flux were well controlled, with no core impurity accumulation, and a new high-confinement and self-organizing regime (Super I-mode = I-mode + e-ITB) was discovered and demonstrated. These achievements contribute to the integration of fusion plasma technology and physics, which is essential to operate next-step devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗