Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesoscale”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical case study of the aerosol–cloud interactions in warm boundary layer clouds over the eastern North Atlantic with an interactive chemistry module

The presence of warm boundary layer stratiform clouds over the eastern North Atlantic (ENA) region is commonly influenced by the Azores High, especially during the summer season. To investigate comprehensive aerosol–cloud interactions, this study employs the Weather Research and Forecasting model coupled with a chemistry component (WRF-Chem), incorporating aerosol chemical components that are relevant to the formation of cloud condensation nuclei (CCN) and accounting for aerosol spatiotemporal variation. This study focuses on aerosol indirect effects, particularly the long-range transport of aerosols, in the ENA region under three different weather regimes: a ridge with a surface high-pressure system, a post-trough with a surface high-pressure system, and a weak trough. The WRF-Chem simulations conducted at a near-large-eddy scale offer valuable insights into the model's performance, especially in terms of its ability to use high spatial resolution to capture mesoscale cloud features across various weather regimes. Our result shows that introducing 5 times more aerosols to either non-precipitating or precipitating clouds significantly increases ambient CCN numbers, resulting in, to varying degrees, higher liquid water path (LWP) values. The substantial aerosol–cloud interaction especially occurs in the precipitating clouds and demonstrates the susceptibility of the LWP to changes in CCN under different regimes. Conversely, thin, non-rain clouds at the edges of a cloud system are prone to evaporation, exhibiting an aerosol drying effect. The aerosols released during this process transition back to the accumulation mode, facilitating future activation. This dynamic behavior is not adequately represented in prescribed-aerosol simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

The interactive role of subsynoptic scale jet sreak and planetary boundary layer adjustments in organizing an apparently isolated convective complex

A mesoscale atmospheric simulation system is described that is being developed in order to improve the simulation of subsynoptic and mesoscale adjustments associated with cyclogenesis, severe storm development, and significant atmospheric transport processes. Present emphasis in model development is in the parameterization of physical processes, time-dependent boundary conditions, sophisticated initialization and analysis procedures, nested grid solutions, and applications software development. Basic characteristics of the system as of March 1982 are listed. In a case study, the Grand Island tornado outbreak of 3 June 1980 is considered in substantial detail. Results of simulations with a mesoscale atmospheric simulation system indicate that over the high plains subtle interactions between existing jet streaks and deep well mixed boundary layers can lead to well organized patterns of mesoscale divergence and pressure falls. The amplitude and positioning of these mesoscale features is a function of the subtle nonlinear interaction between the pre-existing jet-streak and deep well mixed boundary layers. Model results for the case study indicate that the model has the potential for forecasting the precursor mesoscale convective environment.

Kaplan, M. L.↗
NASA NTRS1993

Session on validation of coupled models

The session on validation of coupled models is reviewed. The current use of a mesoscale model with a grid size of 20-km during STORM-FEST in 1992 has proven to be extremely valuable. The availability of forecast products at a much higher temporal and spatial resolution was very helpful for mesoscale forecasting, mission planning, and the guidance of research aircraft. Recent numerical simulation of ocean cyclones and mesoscsle convective systems using nonhydrostatic cloud/mesoscale models with a grid size as small as 2-km have demonstrated the potential of these models for predicting mesoscale convective systems, squall lines, hurricane rainbands, mesoscale gravity waves, and mesoscale frontal structures embedded within an extratropical cyclone. Although mesoscale/cloud scale models have demonstrated strong potential for use in operational forecasting, very limited quantitative evaluation (and verification) of these models were performed. As a result, the accuracy, the systematic biases, and the useful forecasts limits were not properly defined for these models. Also, no serious attempts were made to use these models for operational prediction of mesoscale convective systems.

Kuo, Bill↗
NASA NTRS1998

A Parameterization for the Triggering of Landscape Generated Moist Convection

A set of relatively high resolution three-dimensional (3D) simulations were produced to investigate the triggering of moist convection by landscape generated mesoscale circulations. The local accumulated rainfall varied monotonically (linearly) with the size of individual landscape patches, demonstrating the need to develop a trigger function that is sensitive to the size of individual patches. A new triggering function that includes the effect of landscapes generated mesoscale circulations over patches of different sizes consists of a parcel's perturbation in vertical velocity (nu(sub 0)), temperature (theta(sub 0)), and moisture (q(sub 0)). Each variable in the triggering function was also sensitive to soil moisture gradients, atmospheric initial conditions, and moist processes. The parcel's vertical velocity, temperature, and moisture perturbation were partitioned into mesoscale and turbulent components. Budget equations were derived for theta(sub 0) and q(sub 0). Of the many terms in this set of budget equations, the turbulent, vertical flux of the mesoscale temperature and moisture contributed most to the triggering of moist convection through the impact of these fluxes on the parcel's temperature and moisture profile. These fluxes needed to be parameterized to obtain theta(sub 0) and q(sub 0). The mesoscale vertical velocity also affected the profile of nu(sub 0). We used similarity theory to parameterize these fluxes as well as the parcel's mesoscale vertical velocity.

Lynn, Barry H.↗
NASA NTRS2007

Mars-GRAM Applications for Mars Science Laboratory Mission Site Selection Processes

An overview is presented of the Mars-Global Reference Atmospheric Model (Mars-GRAM 2005) and its new features. One important new feature is the "auxiliary profile" option, whereby a simple input file is used to replace mean atmospheric values from Mars-GRAM's conventional (General Circulation Model) climatology. An auxiliary profile can be generated from any source of data or alternate model output. Results are presented using auxiliary profiles produced from mesoscale model output (Southwest Research Institute's Mars Regional Atmospheric Modeling System (MRAMS) model and Oregon State University's Mars mesoscale model (MMM5) model) for three candidate Mars Science Laboratory (MSL) landing sites (Terby Crater, Melas Chasma, and Gale Crater). A global Thermal Emission Spectrometer (TES) database has also been generated for purposes of making 'Mars-GRAM auxiliary profiles. This data base contains averages and standard deviations of temperature, density, and thermal wind components, averaged over 5-by-5 degree latitude bins and 15 degree L(sub S) bins, for each of three Mars years of TES nadir data. Comparisons show reasonably good consistency between Mars-GRAM with low dust optical depth and both TES observed and mesoscale model simulated density at the three study sites. Mean winds differ by a more significant degree. Comparisons of mesoscale and TES standard deviations' with conventional Mars-GRAM values, show that Mars-GRAM density perturbations are somewhat conservative (larger than observed variability), while mesoscale-modeled wind variations are larger than Mars-GRAM model estimates. Input parameters rpscale (for density perturbations) and rwscale (for wind perturbations) can be used to "recalibrate" Mars-GRAM perturbation magnitudes to better replicate observed or mesoscale model variability.

Justh, Hilary↗
DOE OSTI2022

Overestimated Eddy Kinetic Energy in the Eddy‐Rich Regions Simulated by Eddy‐Resolving Global Ocean–Sea Ice Models

Abstract The performance of eddy‐resolving global ocean–sea ice models in simulating mesoscale eddies is evaluated using six eddy‐resolving experiments forced by different atmospheric reanalysis products. Interestingly, eddy‐resolving ocean general circulation models (OGCMs) tend to simulate more (less) energetic eddy‐rich (eddy‐poor) regions with a smaller (larger) spatial extent than satellite observation, which finally shows that larger (smaller) mesoscale energy intensity (EI) is simulated in the eddy‐rich (eddy‐poor) regions. Quantitatively, there is an approximately 27%–60% overestimation of EI in the eddy‐rich regions, which are mainly located in the Kuroshio–Oyashio Extension, the Gulf Stream, and the Antarctic Circumpolar Currents regions, although the global mean EI is underestimated by 25%–45%. Apparently, the eddy kinetic energy in the eddy‐poor region is underestimated. Further analyses based on coherent mesoscale eddy properties show that the overestimation in the eddy‐rich regions is mainly attributed to mesoscale eddies’ intensity and is more prominent when mesoscale eddies are in their growth stage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Diagnosis of convective organization and cold pools using ARM datasets and evaluation of a unified convection parameterization (UNICON)

Tropical thunderstorms often cluster together. Studies have suggested that the degree to which the tropical thunderstorms are clustered impacts Earth's energy balance and water cycle, as well as extreme precipitation events. However, the processes controlling the spatial distribution of the thunderstorms are poorly understood and are not properly represented in most computer models for weather and climate prediction. Under the goals of better understanding how convection organizes at the mesoscale and advancing the representation of mesoscale convective organization in global climate models, we i) objectively quantified the degrees of convective organization and diagnosed cold pool processes using ARM field campaign observations, ii) examined the organization processes in storm-resolving model simulations, and iii) evaluated the impacts of mesoscale convective organization in global model simulations. The project yielded a firm reference against which the global model representation of mesoscale convective organization and cold pools can be evaluated against and shed new light into the role of parameterized convective organization in the global model simulation of the basic state and variability. Our results revealed two distinct phases of convective clustering during the two-day rain episodes (Cheng et al. 2018) and a new mechanism through which vertical wind shear in the low-troposphere can aid convective organization over tropical oceans (Cheng et al. 2020). It was demonstrated that the interactive representation of cold pools and mesoscale convective organization is key for global models to successfully simulate both the mean state and intraseasonal variability in the tropics (Ahn et al. 2019; 2020).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

The measurement of the synoptic scale wind over the ocean

Mesoscale and microscale features of the turbulent winds over the ocean are related to the synoptic scale winds in terms of published spectral forms for the microscale, a mesoscale valley and published values of U*, VAR u', VAR v' and z/L, as defined in the text and as obtained for moderate to gale force winds. The frequencies involved correspond to periods longer than 1 hour and extend to the microscale, which starts at a period near 2 minutes, or so, and continues to the Kolmogorov inertial range. Nondimensional spectra that span both the mesoscale and the microscale are derived as a function of u, f(= n z/u) and z/L, where z is 10 meters, L is the Monin Obukov stability length and u is evaluated at 10 meters. For the same u, different values of z/L produce a range of values of u which in turn result in variations of the eddy structure of the mesoscale and microscale spectra. Both conventional anemometer averages and remotely sensed winds contain a random component of the mesoscale wind in their values. These components are differnces and not errors when winds are compared, and quantitative values for these differences are given. Ways to improve the measurement of the synoptic scale wind by transient ships, data buoys and scatterometers on future spacecraft are described. These ways are loner averaging times for ships and data buoys, depending on the synoptic conditions, and pooling spacecraft to form super observations. Design considerations for future remote sensing systems are given.

Pierson, W. J.↗
NASA NTRS1984

Data analysis and model simulation of the initiation of convection on April 24, 1982

A combined modelling and data analysis is used to examine the initiation and development of a hail-producing mesoscale convective system which developed over the Texas panhandle near Amarillo, Texas, on April 24, 1982, and propagated east into Oklahoma. This event occurred during the NASA AVE-VAS IV experiment, as a dense upper air mesoscale network was collecting data at three-hour intervals. The results show that both synoptic and mesoscale processes separately and in concert help to force the convection. The destabilization of the atmosphere due to cold advection is important in removing any cap to deep convective development. The upper level moistening due to flow from the northwest also helps to increase the convective instability of the atmosphere. The mesoscale convergence due to cloud shading appears dominant over that of topography, but the latter cannot be ignored, especially in its sustained effect. The developing mesoscale circulation also appears to tap a surface supply of moisture feeding the convection.

Mcnider, R. T.↗
NASA NTRS1985

Cooperative VAS program with NASA Marshall Space Flight Center

The importance of mesoscale moisture information in the forecasting of weather events is discussed. A test is being prepared with the March 6, 1982 Vertical Atmospheric Sounder/Atmospheric Variability Experiment (VAS/AVE) case study to be run on the sub-synoptic scale model (SSM). Model intracomparison of three carefully designed simulations should isolate the role of mesoscale information in the initial conditions for both mositure and vertical motion. Three simulations will be made for the period 1200Z 6 March to 0000Z 7 March starting with the regular synoptic time data of 1200Z. The distinction between the three cases arises from data manipulation at 1800Z midway through the forecast period: (1) Control Case: no alterations at 1800Z. (2) Moisture Information Suppression: at 1800Z the model simulation field for water vapor, which by this time contains considerable mesoscale structure, is replaced by a smoothed version from which the mesoscale structure has been eliminated. No other variables are structure altered. (3) Vertical Motion Information Suppression: at 1800Z the model is subjected to a dead start procedure which primarily imposes a nondivergent constraint of the horizontal motion field field suppressing vertical motion. The main effect is to eliminate the mesoscale component. Comparison of results for the 1800Z to 0000Z time period will give the information desired.

Houghton, D. D.↗
NASA NTRS1990

Diabatic forcing and initialization with assimilation of cloud and rain water in a forecast model: Methodology

The focus of this part of the investigation is to find one or more general modeling techniques that will help reduce the time taken by numerical forecast models to initiate or spin-up precipitation processes and enhance storm intensity. If the conventional data base could explain the atmospheric mesoscale flow in detail, then much of our problem would be eliminated. But the data base is primarily synoptic scale, requiring that a solution must be sought either in nonconventional data, in methods to initialize mesoscale circulations, or in ways of retaining between forecasts the model generated mesoscale dynamics and precipitation fields. All three methods are investigated. The initialization and assimilation of explicit cloud and rainwater quantities computed from conservation equations in a mesoscale regional model are examined. The physical processes include condensation, evaporation, autoconversion, accretion, and the removal of rainwater by fallout. The question of how to initialize the explicit liquid water calculations in numerical models and how to retain information about precipitation processes during the 4-D assimilation cycle are important issues that are addressed. The explicit cloud calculations were purposely kept simple so that different initialization techniques can be easily and economically tested. Precipitation spin-up processes associated with three different types of weather phenomena are examined. Our findings show that diabatic initialization, or diabatic initialization in combination with a new diabatic forcing procedure, work effectively to enhance the spin-up of precipitation in a mesoscale numerical weather prediction forecast. Also, the retention of cloud and rain water during the analysis phase of the 4-D data assimilation procedure is shown to be valuable. Without detailed observations, the vertical placement of the diabatic heating remains a critical problem.

Raymond, William H.↗
NASA NTRS1999

Large Eddy Simulation of Cirrus Clouds

The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent heating could be comparable in magnitude in the cloudy layer. This implies that latent heating cannot be neglected in the construction of a cirrus cloud model. The probability density function (PDF) of w was analyzed to assist in the parameterization of cloud-scale velocities in large-scale models. For the more radiatively-driven, thin cirrus case, the PDFs are approximately Gaussian. However, in the interior of the deep, convectively unstable case, the PDFs of w are multi-modal and very broad, indicating that parameterizing cloud-scale motions for such clouds can be very challenging. The results of this research are described in detail in a paper submitted to the Journal of Atmospheric Science (Wu and Cotton, 1999), which is enclosed as Appendix 2. Using soundings extracted from a mesoscale simulation of the November 26, 1991 cirrus event, the radiative effects on vapor deposition/sublimation of ice crystals was studied using a two-dimensional cloud-resolving model (CRM) version of RAMS, coupled to an explicit bin-resolving microphysics. The CRM simulations of the November 26, 1991 cirrus event demonstrate that the radiative impact on the diffusional growth (or sublimation) of ice crystals is significant. In this case, the ice particles experienced radiative warming. Model results show that radiative feedbacks in the diffusional growth of ice particles can be very complex. Radiative warming of an ice particle will restrict the particle's diffusional growth. In the case of radiative warming, ice particles larger than a certain size will experience so much radiative warming that surface ice saturation vapor pressures become large enough to cause sublimation of the larger crystals, while smaller crystals are growing by vapor deposition. However, ice mass production can be enhanced in the case of radiative cooling of an ice particle. For the November 26, 1991 cirrus event, radiative feedback results in significant reduction in the total ice mass, especially in the production of large ice crystals, and consequently, both radiative and dynamic properties of the cirrus cloud are significantly affected. A complete description of this research has been submitted as a paper to the Journal of Atmospheric Science (Wu et al., 1999), and included as Appendix 3.

Wu, Ting↗
NASA NTRS1999

A TRMM Rainfall Estimation Method Applicable to Land Areas

Methods developed to estimate rain rate on a footprint scale over land with the satellite-borne multispectral dual-polarization Special Sensor Microwave Imager (SSM/1) radiometer have met with limited success. Variability of surface emissivity on land and beam filling are commonly cited as the weaknesses of these methods. On the contrary, we contend a more significant reason for this lack of success is that the information content of spectral and polarization measurements of the SSM/I is limited. because of significant redundancy. As a result, the complex nature and vertical distribution C, of frozen and melting ice particles of different densities, sizes, and shapes cannot resolved satisfactorily. Extinction in the microwave region due to these complex particles can mask the extinction due to rain drops. Because of these reasons, theoretical models that attempt to retrieve rain rate do not succeed on a footprint scale. To illustrate the weakness of these models, as an example we can consider the brightness temperature measurement made by the radiometer in the 85 GHz channel (T85). Models indicate that T85 should be inversely related to the rain rate, because of scattering. However, rain rate derived from 15-minute rain gauges on land indicate that this is not true in a majority of footprints. This is also supported by the ship-borne radar observations of rain in the Tropical Oceans and Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA-COARE) region over the ocean. Based on these observations. we infer that theoretical models that attempt to retrieve rain rate do not succeed on a footprint scale. We do not follow the above path of rain retrieval on a footprint scale. Instead, we depend on the limited ability of the microwave radiometer to detect the presence of rain. This capability is useful to determine the rain area in a mesoscale region. We find in a given rain event that this rain area is closely related to the mesoscale-average rain rate. Based on this observation, in this study we have developed a method to estimate the mesoscale-average rain rate over land utilizing microwave radiometer data. Because of the high degree of geographic and seasonal variability in the nature and intensity of rain, this method requires some tuning with 15-minute rain gauge data on land. After tuning the method, it can be applied to an independent set of rain events that are close in time and space. We find that the mesoscale rain rates retrieved over the period of a month on land with this method shows a correlation of about 0.85 with respect to the surface rain-gauge observations. This mesoscale-average rain rate estimation method can be useful to extend the spatial and temporal coverage of the rainfall data provided by the Precipitation Radar on board the Tropical Rainfall Measuring Mission (TRMM) satellite.

Prabhakara, C.↗
DOE OSTI2022

Evaluation of the Fitch Wind-Farm Wake Parameterization with Large-Eddy Simulations of Wakes Using the Weather Research and Forecasting Model

Abstract Wind-farm parameterizations in weather models can be used to predict both the power output and farm effects on the flow; however, their correctness has not been thoroughly assessed. We evaluate the wind-farm parameterization of the Weather Research and Forecasting Model with large-eddy simulations (LES) of the wake performed with the same model. We study the impact on the velocity and turbulence kinetic energy (TKE) of inflow velocity, roughness, resolution, number of turbines (one or two), and inversion height and strength. We compare the mesoscale with the LES by spatially averaging the LES within areas correspondent to the mesoscale horizontal spacing: one covering the turbine area and two downwind. We find an excellent agreement of the velocity within the turbine area between the two types of simulations. However, within the same area, we find the largest TKE discrepancies because in mesoscale simulations, the turbine-added TKE has to be highest at the turbine position to be advected downwind. Within the downwind areas, differences between velocities increase as the wake recovers faster in the LES, whereas for the TKE both types of simulations show similar levels. From the various configurations, the impact of inversion height and strength is small for these heights and inversion levels. The highest impact for the one-turbine simulations appears under the low-speed case due to the higher thrust, whereas the impact of resolution is low for the large-eddy simulations but high for the mesoscale simulations. Our findings demonstrate that higher-fidelity simulations are needed to validate wind-farm parameterizations.

17 WIND ENERGY↗
DOE OSTI2022

Computing diffraction patterns of microstructures from phase-field simulations

The diffraction pattern of a material contains information not only on the crystal structures of its constituting phases, but also on its mesoscale spatial distributions of phases, grains, and ferroelastic, ferroelectric, and ferromagnetic domains. While diffraction patterns from experiments such as X-ray diffraction are presented in the reciprocal or Fourier space, mesoscale microstructure models such as the phase-field method naturally produce real-space images of spatial distribution of chemical composition, structural, and ferroic domains. Although one could rather readily compute the Fourier amplitudes of chemical and structural domain distributions generated by mesoscale simulations, they only contain information about the length scale and alignment of the real-space chemical and structure domains. Therefore, a direct comparison between diffraction experiments and mesoscale microstructure simulations is not possible. Here, we develop a theoretical approach to directly compute the crystal diffraction patterns of microstructures predicted by phase-field simulations. In particular, we consider five representative examples of microstructure patterns involving purely compositional domains, a single pair of tetragonal twin structures, multiple twin variants in a hexagonal system, ferroelectric polar vortices, and polycrystalline grains. The results are compared with previous experimental observations as well as X-ray diffraction experiments performed in the present study. Furthermore, the theoretical framework allows one to directly connect material microstructures and diffraction patterns predicted from phase-field simulations and the corresponding diffraction patterns from experiments, and thus providing guidance to experimental diffraction characterization and interpretation of microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗
DOE OSTI2022

Long‐Term Observations of Turbulence Vertical Velocity Spectra in a Convective Mixed Layer: Dependence on Land‐Surface Forcing in the U.S. Southern Great Plains

Doppler lidar vertical velocity retrievals were analyzed for the scale and structure of mixed-layer turbulence over a 7-year period, on fair-weather warm season days (442 cases) in the U.S. Southern Great Plains. Based on the hypothesized influence of land-surface forcing and updraft size on convective boundary layer clouds, a spectral analysis was performed to quantify effects of surface forcing (surface buoyancy flux, friction velocity, and evaporative fraction) on turbulence scale. Significant (order-of-magnitude) variations in spectral density were found in the energy-production subrange and mesoscale regimes. Integral scale decreased with increasing buoyancy flux (and Monin-Obukhov stability parameter), while spectral density in the energy-production subrange increased, implying a transition to buoyancy-driven cellular structures with narrower updrafts. The influence of stability parameter was limited to the neutral to convective transition, and could not explain the wide variation in spectral density in the mesoscale regime. However, high friction velocity (u* > 0.5 m s –1 ) was associated with larger integral scale (updraft width), and a greater portion of spectral density in the mesoscale regime. Lidar profiles and radar reflectivity for individual cases revealed evidence of roll and wave-like structures contributing to mesoscale variability on high friction velocity days, suggesting shear instabilities on days with wetter land surface conditions and smaller buoyancy flux. The role of friction velocity emphasizes the multivariate nature of surface influences on turbulence, and implies that variance-based turbulence closures may not be adequate for capturing effects of surface forcing on updraft size.

54 ENVIRONMENTAL SCIENCES↗