Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary layer structure”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulating Real Atmospheric Boundary Layers at Gray-Zone Resolutions: How Do Currently Available Turbulence Parameterizations Perform?

Recent computational and modeling advances have led a diverse modeling community to experiment with atmospheric boundary layer (ABL) simulations at subkilometer horizontal scales. Accurately parameterizing turbulence at these scales is a complex problem. The modeling solutions proposed to date are still in the development phase and remain largely unvalidated. This work assesses the performance of methods currently available in the Weather Research and Forecasting (WRF) model to represent ABL turbulence at a gray-zone grid spacing of 333 m. We consider three one-dimensional boundary layer parameterizations (MYNN, YSU and Shin-Hong) and coarse large-eddy simulations (LES). The reference dataset consists of five real-case simulations performed with WRF-LES nested down to 25 m. Results reveal that users should refrain from coarse LES and favor the scale-aware, Shin-Hong parameterization over traditional one-dimensional schemes. Overall, the spread in model performance is large for the cellular convection regime corresponding to the majority of our cases, with coarse LES overestimating turbulent energy across scales and YSU underestimating it and failing to reproduce its horizontal structure. Despite yielding the best results, the Shin-Hong scheme overestimates the effect of grid dependence on turbulent transport, highlighting the outstanding need for improved solutions to seamlessly parameterize turbulence across scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of Surface and Top Wind Shear on the Spatial Organization of Marine Stratocumulus-Topped Boundary Layers

Typically, the convective nature of Stratocumulus topped boundary layers (STBL) involves the motion of updrafts and downdrafts, driven by surface fluxes and radiative cooling, respectively. The balance between shear and buoyant forcings at the surface can determine the organization of updrafts between cellular and roll structures. We investigate the effect of varying shear at the surface and top of the STBL using Large Eddy Simulations, taking DYCOMS II RF01 as a base case. We focus on spatial identification of the following features: coherent updrafts and downdrafts, and observe how they are affected by varying shear. Stronger surface shear organizes the updrafts in rolls, causes less well-mixed thermodynamic profiles, and decreases cloud fraction and liquid water path (LWP). Stronger top shear also decreases cloud fraction and LWP more than surface shear, by thinning the cloud from the top. Features with stronger top than surface shear are associated with a net downward momentum transport and show early signs of decoupling. Classifying updrafts and downdrafts based on their vertical span and horizontal size confirms the dominance of tall objects spanning the whole STBL. Tall objects occupy 30% of the volume in the STBL, while short ones occupy less than 1%. For updraft and downdraft fluxes, these tall objects explain 65% of the vertical velocity variance and 83% of the buoyancy flux, on average. Stronger top shear also weakens the contribution of downdrafts to the turbulent fluxes and tilts the otherwise vertical development of updrafts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Scientific challenges to characterizing the wind resource in the marine atmospheric boundary layer

Abstract. With the increasing level of offshore wind energy investment, it is correspondingly important to be able to accurately characterize the wind resource in terms of energy potential as well as operating conditions affecting wind plant performance, maintenance, and lifespan. Accurate resource assessment at a particular site supports investment decisions. Following construction, accurate wind forecasts are needed to support efficient power markets and integration of wind power with the electrical grid. To optimize the design of wind turbines, it is necessary to accurately describe the environmental characteristics, such as precipitation and waves, that erode turbine surfaces and generate structural loads as a complicated response to the combined impact of shear, atmospheric turbulence, and wave stresses. Despite recent considerable progress both in improvements to numerical weather prediction models and in coupling these models to turbulent flows within wind plants, major challenges remain, especially in the offshore environment. Accurately simulating the interactions among winds, waves, wakes, and their structural interactions with offshore wind turbines requires accounting for spatial (and associated temporal) scales from O(1 m) to O(100 km). Computing capabilities for the foreseeable future will not be able to resolve all of these scales simultaneously, necessitating continuing improvement in subgrid-scale parameterizations within highly nonlinear models. In addition, observations to constrain and validate these models, especially in the rotor-swept area of turbines over the ocean, remains largely absent. Thus, gaining sufficient understanding of the physics of atmospheric flow within and around wind plants remains one of the grand challenges of wind energy, particularly in the offshore environment. This paper provides a review of prominent scientific challenges to characterizing the offshore wind resource using as examples phenomena that occur in the rapidly developing wind energy areas off the United States. Such phenomena include horizontal temperature gradients that lead to strong vertical stratification; consequent features such as low-level jets and internal boundary layers; highly nonstationary conditions, which occur with both extratropical storms (e.g., nor'easters) and tropical storms; air–sea interaction, including deformation of conventional wind profiles by the wave boundary layer; and precipitation with its contributions to leading-edge erosion of wind turbine blades. The paper also describes the current state of modeling and observations in the marine atmospheric boundary layer and provides specific recommendations for filling key current knowledge gaps.

17 WIND ENERGY↗
DOE OSTI2024

Deep-learning based artificial intelligence tool for melt pools and defect segmentation

Accelerating fabrication of additively manufactured components with precise microstructures is important for quality and qualification of built parts, as well as for a fundamental understanding of process improvement. Accomplishing this requires fast and robust characterization of melt pool geometries and structural defects in images. This paper proposes a pragmatic approach based on implementation of deep learning models and self-consistent workflow that enable systematic segmentation of defects and melt pools in optical images. Deep learning is based on an image-to-image translation–conditional generative adversarial neural network architecture. An artificial intelligence (AI) tool based on this deep learning model enables fast and incrementally more accurate predictions of the prevalent geometric features, including melt pool boundaries and printing-induced structural defects. We present statistical analysis of geometric features that is enabled by the AI tool, showing strong spatial correlation of defects and the melt pool boundaries. The correlations of widths and heights of melt pools with dataset processing parameters show the highest sensitivity to thermal influences resulting from laser passes in adjacent and subsequent layer passes. The presented models and tools are demonstrated on the aluminum alloy and datasets produced with different sets of processing parameters. However, they have universal quality and could easily be adapted to different material compositions. The method can be easily generalized to microstructural characterizations other than optical microscopy.

additive manufacturing↗
DOE OSTI2021

Composited structure of non-precipitating shallow cumulus clouds

The normalized distributions of thermodynamic and dynamical variables both within and outside shallow clouds are investigated through a composite algorithm using large-eddy simulations of oceanic and continental cases. The normalized magnitude is maximum near the cloud centre and decreases outwards. While relative humidity (RH) and cloud liquid water ( q l ) decrease smoothly to match the environment, the vertical velocity, virtual potential temperature (θ v ), and potential temperature (θ) perturbations have more complicated behaviour towards the cloud boundary. Below the inversion layer, θ' v becomes negative before the vertical velocity has turned from an updraft to a subsiding shell outside the cloud, indicating the presence of a transition zone where the updraft is negatively buoyant. Due to the downdraft outside the cloud and enhanced horizontal turbulent mixing across the edge, the normalized turbulent kinetic energy (TKE) and horizontal turbulent kinetic energy (HTKE) decrease more slowly from the cloud centre outwards than the thermodynamic variables. The distributions all present asymmetric structures in response to the vertical wind shear, with more negatively buoyant air, stronger downdrafts, and larger TKE on the downshear side. We discuss several implications of the distributions for theoretical models and parameterizations. Positive buoyancy near the cloud base is mostly due to the virtual effect of water vapour, emphasizing the role of moisture in triggering. The mean vertical velocity is found to be approximately half the maximum vertical velocity within each cloud, providing a constraint to achieve possible power-law distributions for some models. Finally, the normalized distributions for different variables are used to estimate the vertical heat and moisture fluxes within clouds. The results suggest that distributions near the cloud edge and variability of maximum perturbations need careful treatment. The fluxes are underestimated in the inversion layer because cloud-top downdrafts cannot be captured well.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Insights into the relation between vertical cloud structure and dynamics of three polar lows: Observations from COMBLE

The vertical structure of three polar lows is described using profiling radar, lidar, and passive remote sensors deployed at a coastal site in northern Norway. The data were collected as part of the Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE). These data are examined in the context of satellite imagery, operational weather radar reflectivity imagery, and output from the AROME-Arctic model. A comparison to satellite images shows that AROME-Arctic captured the polar lows well. All three polar lows form in marine cold-air outbreaks and are surrounded by open-cellular convection typical of such outbreaks. Two of the lows form in a forward shear environment, and one under reverse shear. Initially, the three polar lows are comma-shaped; two of them transition to be more spirali-form at their mature stage and have mostly cloud-free warm cores. The warm cores are the result of a warm-seclusion process. All three lows have stratiform precipitation bands, marked by little cloud liquid water, and rather high surface precipitation rate. The vertical drafts and turbulence in these stratiform clouds are generally weak. In conclusion, all three lows also have convective clouds, which have stronger vertical drafts, stronger turbulence, and pockets of high liquid water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of divertor electrical drifts on particle distribution and detachment near the divertor target plate in DIII-D

Strong impacts of drifts on the divertor plasma in–out asymmetry and detachment are demonstrated in DIII-D with an open divertor configuration. For forward toroidal field, B T , i.e., with the ion B × ∇B drift toward the divertor, the particle flux to the inner divertor, as represented by the Langmuir probe measured ion saturation current (J sat ), exhibits a double peak structure, with electron temperature, lower at the inner target. Reversing the B T direction reverses both the radial and poloidal E × B flows, leading to a broad particle flux profile in the outboard scrape-off layer (SOL) with a similar double-peak structure to that observed at the inner target with forward B T . The correlation of a double peak structure with divertor temperature profiles confirms physical coupling between the drift flow and sheath boundary condition and their strong impact on divertor profiles. In addition, under reversed B T conditions, increasing the density flattens the target temperature profile. However, J sat remains high away from the strike point, rendering it difficult to achieve an “effective” detached plasma, i.e., with effective reduction in both peak heat flux and peak temperature (in the far SOL). In contrast, divertor detachment with a cold and flat temperature profile can be achieved at both target plates with the forward B T .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Towards exascale for wind energy simulations

We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.

17 WIND ENERGY↗
DOE OSTI2022

FY2022 Q4: Demonstrate multi-turbine simulation with hybrid-structured/unstructured-moving-grid software stack running primarily on GPUs and propose improvements for successful KPP-2 [Poster]

Milestone accomplishments staged the ExaWind team for successful completion of KPP-2 challenge problem in FY23, which requires the simulation on Frontier of at least four MW-scale turbines in an atmospheric boundary layer with at least 20B gridpoints. The ExaWind project and software stack is many faceted, with team members working on multiple areas, including linear-system solvers (Trilinos, hypre, AMReX), overset meshes, turbulence modeling, and in situ visualization, all with an aim for high fidelity predictions and performance portability. This milestone marks significant improvements on many fronts and provides the team with a pathway to exascale wind farm simulations in FY23.

17 WIND ENERGY↗
DOE OSTI2022

Microstructural examination of interactions between chromia-based refractory and nuclear glass in a melter

This study shows the effects of nuclear waste glass production on Monofrax K-3 refractory corrosion. A continuously fed research-scale melter containing an Fe- and Ni-rich simulated nuclear waste feed with borosilicate glass-forming chemicals was cyclically melted at 1150°C and idled at 1050°C for a total of 11 weeks. Chemical maps using scanning electron microscopy show the interactions between the quenched melt and the refractory. Nanoscale X-ray-computed tomography was used for a three-dimensional visualization of certain parts of the interface. Unreacted K-3 consists of primarily corundum (Al,Cr) 2 O 3 and spinel (Fe 2+ ,Mg)(Al,Cr) 2 O 4 interlocking crystalline phases. Some of the Cr from the refractory interacts with the Ni and Fe from the melt to form a reaction layer comprising (Ni,Fe 2+ )(Cr,Fe 3+ ) 2 O 4 spinel crystals. Simultaneously, melt components (Na,Si) infiltrate into the refractory. This interaction proceeds at the expense of the integrity of the refractory structure. Intact refractory grains (e.g., (Al,Cr) 2 O 3 ) as well as the reaction layer itself can lose mechanical integrity and spall off into the melt, especially near the top of the melter. As the reaction layer can be a protective boundary for the refractory against further melt infiltration, a reduction in the reaction layer thickness allows an increase in refractory corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Wind and Temperature Consensus at Horn Point, HU-Beltsville, Piney Run (Maryland) in support of CoURAGE

The Maryland Department of the Environment (MDE) operates a ground-based atmospheric profiling network consisting of collocated radar wind profilers (RWP) and radio acoustic sounding systems (RASS) as part of its Ambient Air Monitoring Program. This network provides continuous observations of wind and temperature structure in the lower troposphere to support air quality forecasting, regulatory analysis, and atmospheric research. The network currently includes three fixed sites across Maryland: Horn Point (HP, lower eastern shore) [38.587525°,-76.141006°], Howard University-Beltsville (HUB, central Maryland) [39.055277°, -76.878632°], and Piney Run (PR, western Maryland) [39.705950°, -79.012000°] The network is designed to capture regional variability in atmospheric transport and boundary-layer processes. These systems measure vertical profiles of horizontal wind speed and direction using Doppler radar techniques, with observations typically spanning from ~100 m above ground level up to approximately 2.5–4 km. Measurements are derived from the Doppler shift of backscattered electromagnetic signals, enabling retrieval of wind vectors at multiple altitudes with high temporal resolution (e.g., 30-minute averages reported every 6 minutes). Each radar wind profiler is paired with a Radio Acoustic Sounding System (RASS) to provide profiles of virtual temperature in the lower atmosphere (~100–200 m AGL) by measuring the propagation speed of acoustic waves. Together, the RWP/RASS system yields a coupled data set of thermodynamic and kinematic atmospheric structure, including additional parameters such as vertical velocity, radial velocity, signal-to-noise ratio, and spectral width for advanced analysis. There are two types of files for each station: wind data (files with a "w" prefix) and virtual temperature RASS data (files with a "t" prefix). The wind data files are in the format wYYDDD.cns, where YY is the 2-digit year and DDD is the day of the year. The RASS virtual temperature data files are in the format tYYDDD.cns. Each record has the following header structure: Line 1 : Station Name RASS files Line 2 : RASS rev DeTect_2.0, WINDS files Line 2 : WINDS rev ATI 5.1 Line 3 : N latitude, W longitude, and site elevation (m) Line 4 : Date and begin time of consensus: yy mm dd hh mn ss plus # minutes to add to get UTC Line 5 : Consensus averaging time (minutes); number of beams; number of range gates Line 6 : Number of records required to make consensus (num) total number of records (tot) and the consensus window size (m/s) in the format: num:tot (window) RASS files Line 7 : no. of coded cells, no. of spec, pulse width (ns), and inter-pulse period (µs), WINDS files Line 7 : No. of coded cells, no. of spectra, pulse width (ns), and inter-pulse period (µs), each with a pair of values: first value is for oblique beams, second for vertical RASS files Line 8 : Full scale Doppler value (m/s) Delay to first gate (ns) Number of gates Spacing of gates (ns), WINDS files Line 8 : Full scale Doppler velocity (m/s), oblique and vertical Vertical correction applied to oblique beams? (0 = no, 1 = yes) Delay to first gate (ns), oblique and vertical Number of gates, oblique and vertical Spacing of gates (ns), oblique and vertical Line 9 : Azimuth and elevation (9s indicate vertical beam not used) RASS files Line 10, values : HT = Height above ground (km), T = Uncorrected virtual temperature consensus (deg C), Tc = Corrected virtual temperature consensus (deg C), W = Vertical wind consensus (9s indicate vertical beam not used, w-component, positive upward, m/s), CNT = Number of records that made consensus (for the 3 values in same order), SNR = Average signal to noise ratio (dB) of records in consensus (same order) WINDS files Line 10, values : HT = Height above ground (km), SPD = Wind speed (m/s), DIR = Wind direction (deg E of N from N), RAD = Radial velocities for each beam (m/s) in order given in azimuth and elevation line (positive toward radar; 9s indicate vertical beam not used, CNT = Number of records that made consensus, SNR = Average signal to noise ratio (dB) of records in consensus

{"wind speed and direction",temperature}↗
DOE OSTI2022

Electric-field control of the nucleation and motion of isolated three-fold polar vertices

Recently various topological polar structures have been discovered in oxide thin films. Despite the increasing evidence of their switchability under electrical and/or mechanical fields, the dynamic property of isolated ones, which is usually required for applications such as data storage, is still absent. Here, we show the controlled nucleation and motion of isolated three-fold vertices under an applied electric field. At the PbTiO 3 /SrRuO 3 interface, a two-unit-cell thick SrTiO 3 layer provides electrical boundary conditions for the formation of three-fold vertices. Utilizing the SrTiO 3 layer and in situ electrical testing system, we find that isolated three-fold vertices can move in a controllable and reversible manner with a velocity up to ~629 nm s -1 . Microstructural evolution of the nucleation and propagation of isolated three-fold vertices is further revealed by phase-field simulations. This work demonstrates the ability to electrically manipulate isolated three-fold vertices, shedding light on the dynamic property of isolated topological polar structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Additively manufactured β-Ti5553 with laser powder bed fusion: Microstructures and mechanical properties of bulk and lattice parts

Ti5553 (Ti-5Al-5Mo-5V-3Cr wt%) is a titanium alloy widely used for its high strength-to-weight ratio and good formability at elevated temperatures. Unlike Ti-6Al-4V, Ti5553 does not undergo martensitic transformation, preventing cracking of brittle martensite upon rapid cooling. This makes it a strong candidate for additive manufacturing (AM), particularly laser powder bed fusion (L-PBF). L-PBF offers the unique opportunity to make fine lattice structures to reduce component weight. Despite the growing field of AM, there have been limited studies on L-PBF Ti5553 lattices and how their properties differ from the bulk. Here the present work addresses this knowledge gap by investigating microstructures and properties of L-PBF bulk and lattice parts and the effect of post L-PBF heat treatments. Electron microscopy and mechanical testing show that the high dislocation density formed during L-PBF increases bulk part’s yield strength by approximately 100 MPa compared to the conventional alloy. Digital image correlation during compression testing of octet truss lattices reveals a layer-by-layer failure mode. Compared to the bulk, the lattice contains copious ω nanoprecipitation, weaker <001> texture, smaller average grain sizes, and larger content of high-angle grain boundaries. These features elicit differences in Taylor factor distributions for the lattice depending on load direction, underlining challenges in predicting lattice mechanical response based on bulk properties. By examining the processing-structure-property relationships in the bulk and lattice, the present results delineate their microstructural and mechanical differences and establish a benchmark for the future design applications of L-PBF Ti5553.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Performance analysis of functionally graded multifunctional piezoelectric energy harvesting microgyroscopes

Functionally graded materials (FGMs) are composite materials with varying material properties in one or more directions. These materials possess distinct characteristics compared to their constituent components. The ability to control material distribution and compositions in FGMs offers improved constraints over the natural frequencies of the system, making them highly suitable for energy harvesting applications. In this study, we focus on a piezoelectric FGM composed of Platinum (Pt) and Lead Titanate Zirconate (PZT). By utilizing FGMs for energy harvesting, we simplify the system from a multilayer structure to a single-layer system, thereby increasing power density by reducing volume. In this work, a power law distribution is used to model the material variation throughout the thickness of the single-layered beam. The governing equations of motion and boundary conditions for FG energy harvesting microgyroscope are derived using Euler-Bernoulli beam theory, the constitutive piezoelectric principle, and the extended Hamilton's principle. To simulate the nonlinear motion of the FG energy harvester, we employ the differential quadrature method (DQM). Subsequently, an investigation is carried out to examine the effects of various factors including material distribution, platinum percentage, base rotation, DC voltage, and electrical load resistance on the energy harvesting microgyroscope with functionally graded materials. Here, the findings suggest that functionally graded energy harvesting gyroscopes can be adjusted to obtain effective energy harvesting and sensing capabilities by carefully selecting the material distribution, input voltage, and electrical load resistance.

42 ENGINEERING↗
DOE OSTI2024

Atomic-Scale Scanning of Domain Network in the Ferroelectric HfO 2 Thin Film

Ferroelectric HfO 2 -based thin films have attracted much interest in the utilization of ferroelectricity at the nanoscale for next-generation electronic devices. However, the structural origin and stabilization mechanism of the ferroelectric phase are not understood because the film is typically nanocrystalline with active yet stochastic ferroelectric domains. Here, in this study, electron microscopy is used to map the in-plane domain network structures of epitaxially grown ferroelectric Y:HfO 2 films in atomic resolution. The ferroelectricity is confirmed in free-standing Y:HfO 2 films, allowing for investigating the structural origin for their ferroelectricity by 4D-STEM, high-resolution STEM, and iDPC-STEM. At the grain boundaries of <111>-oriented Pca2 1 orthorhombic grains, a high-symmetry mixed-(R3m, Pnm2 1 ) phase is induced, exhibiting enhanced polarization due to in-plane compressive strain. Nanoscale Pca2 1 orthorhombic grains and their grain boundaries with mixed-(R3m, Pnm2 1 ) phases of higher symmetry cooperatively determine the ferroelectricity of the Y:HfO 2 film. It is also found that such ferroelectric domain networks emerge when the film thickness is beyond a finite value. Furthermore, in-plane mapping of oxygen positions overlaid on ferroelectric domains discloses that polarization is suppressed at vertical domain walls, while it is active when domains are aligned horizontally with subangstrom domain walls. In addition, randomly distributed 180° charged domain walls are confined by spacer layers.

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

Metabolic activity and community structure of prokaryotes associated with particles in the twilight zone of the South China Sea

The twilight zone is an important depth of the ocean where particulate organic matter (POM) remineralization takes place, and prokaryotes contribute to more than 70% of the estimated remineralization. However, little is known about the microbial community and metabolic activity associated with different particles in the twilight zone. The composition and distribution of particle-attached prokaryotes in the twilight zone of the South China Sea (SCS) were investigated using high-throughput sequencing and quantitative PCR, together with the Biolog Ecoplate™ microplates culture to analyze the microbial metabolic activity. We found that α- and γ-Proteobacteria dominating at the lower and upper boundary of the twilight zone, respectively; Methanosarcinales and Halobacteriales of the Euyarchaeota occupied in the larger particles at the upper boundary. Similar microbial community existed between euphotic layer and the upper boundary. Higher amount of shared Operational Taxonomic Units (OTUs) in the larger particles along the water depths, might be due to the fast sinking and major contribution of carbon flux of the larger particles from the euphotic layer. In addition to polymers as the major carbon source, carbohydrates and amino acids were preferentially used by microbial community at the upper and lower boundary, respectively. This could potentially be attributed to the metabolic capabilities of attached microbial groups in different particles, and reflected the initial preference of the carbon source by the natural microbes in the twilight zone as well. The microbial structure and carbon metabolic profiles could be complemented with metatranscriptomic analysis in future studies to augment the understanding of the complex carbon cycling pathways in the twilight zone.

Liu, Hao↗
DOE OSTI2020

Observations of Greenhouse Gas Changes across Summer Frontal Boundaries in the Eastern United States

Our knowledge about synoptic-scale variations of atmospheric CO2 produced by interactions between mid-latitude cyclones and regional-scale surface fluxes remains limited by scarcity of observations. We synthesize observations of greenhouse gases (GHG) with respect to frontal passages to understand how GHG distributions change vertically and horizontally during a synoptic event. We use the airborne in-situ measurements of GHGs collected during the first Atmospheric Carbon and Transport–America (ACT) summer field campaign across three regions of the eastern US. We find that frontal structures are consistently and clearly associated with spatially coherent, significant changes in GHG composition. Warm sector CO2 mole fractions (henceforth, [CO2]) are significantly higher than those in the cold sector. The CO2 frontal differences (5-30 ppm) are largest in the atmospheric boundary layer (ABL) with smaller differences (3-5 ppm) in the free troposphere (FT). Bands of elevated CO2 (5-20 ppm enhancement compared to warm-sector mean [CO2]) repeatedly occurred in the ABL adjacent to the frontal boundary. Analyses of 84 vertical profiles reveal higher [CO2] in the FT than in the ABL in the cold-sector, the opposite pattern in the warm-sector. Average ABL-to-FT [CO2] differences are 12 ppm and -6 ppm in the warm and cold sectors, respectively. The same analysis for CH4 shows higher [CH4] in the ABL than in the FT in both warm and cold sectors with at least 30-ppb higher ABL-to-FT [CH4] differences in the warm-sector than in the cold-sector. The observational analyses presented here define new metrics for model evaluation of GHG transport.

Pal, Sandip↗
DOE OSTI2022

Framework for idealized climate simulations with spatiotemporal stochastic clouds and planetary-scale circulations

In climate predictions, clouds are the leading source of uncertainty. This is partly because, to simulate the fluid dynamics of climate over the entire globe, a large grid spacing must be used, so clouds are a subgrid-scale parametrization rather than a resolved feature. Here, a framework is investigated with finer grid spacing of O(1) or O(10)km so that some clouds are not subgrid-scale; instead, clouds evolve on the numerical grid. This cloud evolution is achieved using stochastic modeling. Hence, the framework is idealized in the sense that the full fluid dynamics of cloud circulations is still not resolved, and simplified vertical structures are used. Nevertheless, the fluid dynamics model includes evolving clouds that interactively adjust in size, shape, lifetime, and regional coverage. In addition, different cloud types are included with different roles in the climate system, including deep convective clouds and also boundary-layer clouds such as shallow cumulus and stratocumulus clouds. Other basic aspects of the idealized climate system are planetary-scale circulations (e.g., Walker circulation) and radiation. With these ingredients (evolving clouds, planetary-scale circulations, and radiation), the framework has the potential for idealized investigations of climate change with interactive cloud–radiative feedback of individual clouds. Here, the formulation of the model equations is presented, and numerical simulations are shown to illustrate the model dynamics and climate change.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗