A Simple, Efficient Solution of Flux–Profile Relationships in the Atmospheric Surface Layer

2006 ◽  
Vol 45 (2) ◽  
pp. 341-347 ◽  
Author(s):  
Jonathan E. Pleim

Abstract This note describes a simple scheme for analytical estimation of the surface-layer similarity functions from state variables. What distinguishes this note from the many previous papers on this topic is that this method is specifically targeted for numerical models in which simplicity and economic execution are critical. In addition, it has been in use in a mesoscale meteorological model for several years. For stable conditions, a very simple scheme is presented that compares well to the iterative solution. The stable scheme includes a very stable regime in which the slope of the stability functions is reduced to permit significant fluxes to occur, which is particularly important for numerical models in which decoupling from the surface can be an important problem. For unstable conditions, simple schemes generalized for varying ratios of aerodynamic roughness to thermal roughness (z0/z0h) are less satisfactory. Therefore, a simple scheme has been empirically derived for a fixed z0/z0h ratio, which represents quasi-laminar sublayer resistance.

2021 ◽  
Author(s):  
Victoria Anne Sinclair ◽  
Jenna Ritvanen ◽  
Gabin Urbancic ◽  
Irina Statnaia ◽  
Yurii Batrak ◽  
...  

Abstract. The planetary boundary layer (BL) height and stratification are key parameters in determining the exchange of heat, momentum, moisture and trace gases between the surface and the free troposphere. Numerous different methods have been used to quantify the BL height and these methods have been applied to a wide variety of observational data sets obtained from different instruments and to numerical model output. We investigate the BL height at the Hyytiälä SMEAR II station in southern Finland diagnosed from radiosonde observations, a microwave radiometer (MWR) and ERA5 reanalysis. Four different algorithms are used to diagnose the BL height from the radiosondes. The diagnosed BL height is sensitive to the method used and the level of agreement, and the sign of systematic bias, between the 4 different methods depends on the surface-layer stability. For example, for very unstable situations, the median BL height diagnosed from the radiosondes varies from 600 m to 1500 m depending on which method is applied. Good agreement between the BL height in ERA5 and diagnosed from the radiosondes using Richardson number-based methods is found for almost all stability classes, suggesting that ERA5 has adequate vertical resolution near the surface to resolve the BL structure. However, ERA5 overestimates the BL height in very stable conditions highlighting the on-going challenge for numerical models to correctly resolve the stable BL. Furthermore, ERA5 BL height differs most from the radiosondes at 18 UTC suggesting ERA5 does not resolve the evening transition correctly. This study has also shown that BL height estimates from the MWR are reliable in unstable situations but often are inaccurate under stable conditions when, in comparison to ERA5 BL heights, they are much deeper. The errors in the MWR BL height estimates originate from the limitations of the manufacturers algorithm for stable conditions and also the mis-identification of the type of BL. A climatology of the annual and diurnal cycle of BL height and observed surface layer stability was created. The shallowest (353 m) monthly median BL height occurs in February and the deepest (576 m) in June. In winter there is no diurnal cycle in BL height, unstable BLs are rare yet so are very stable BLs. The shallowest BLs occur at night in spring and summer and very stable conditions are most common at night in the warm season. Finally, using ERA5 gridded data we determined that the BL height observed at Hyytiälä is representative of most land areas in southern and central Finland. However, the spatial variability of the BL height is largest during daytime in summer reducing the area over which BL height observations from Hyytiälä would be representative of.


2016 ◽  
Vol 38 ◽  
pp. 394
Author(s):  
Janete Gonçalves Lira ◽  
Eduardo Stüker ◽  
Jean Jonathan Schuster ◽  
Cristiano Henrique Schuster ◽  
Daniel Michelon dos Santos ◽  
...  

The proper description of the atmospheric flow in the stable atmospheric boundary layer is one of the most complex tasks to be performed by numerical models of weather and climate prediction. Most of the parameterizations used in these models are based on the statistical theory of turbulence in their conception. However, this theory is valid only for a flow in which turbulence is homogeneous, isotropic and stationary, a conditions that are not commonly found overnight. Basically, the nighttime flow can be classified in two major regimes: very stable, where turbulence is almost entirely suppressed; and weakly stable regime, with intense turbulent mixing. The transition between these regimes is known as atmospheric coupling, and it can happens a lot of times along the same night. In this work, we implemented a single column model with turbulence closure 1.5, in three different configurations and three different turbulence formulations. In general, the model results show that, all the configurations are able to reproduce the average characteristics of the flow in the weakly stable conditions.


2016 ◽  
Vol 17 (11) ◽  
pp. 2835-2851 ◽  
Author(s):  
Yusri Yusup ◽  
Heping Liu

Abstract Widely used numerical models to estimate turbulent exchange of latent heat flux (LE) and sensible heat flux H across the water–atmosphere interface are based on the bulk transfer relations linked indirectly to atmospheric stability, even though the accurate prediction of the influence of stability on fluxes is uncertain. Here eddy covariance data collected over the water surface of Ross Barnett Reservoir, Mississippi, was analyzed to study how atmospheric stability and other variables (wind speed, vapor pressure gradient, and temperature gradient) in the atmospheric surface layer (ASL) modulated LE and H variations in different stability ranges. LE and H showed right-skewed, bell-shaped distributions as the ASL stability shifted from very unstable to near neutral and then stable conditions. The results demonstrate that the maximum (minimum) LE and H did not necessarily occur under the most unstable (stable) conditions, but rather in the intermediate stability ranges. No individual variables were able to explain the dependence of LE and H variations on stability. The coupling effects of stability, wind speed, and vapor pressure gradient (temperature gradient) on LE (H) primarily caused the observed variations in LE and H in different stability ranges. These results have important implications for improving parameterization schemes to estimate fluxes over water surfaces in numerical models.


2021 ◽  
Vol 11 (9) ◽  
pp. 4136
Author(s):  
Rosario Pecora

Oleo-pneumatic landing gear is a complex mechanical system conceived to efficiently absorb and dissipate an aircraft’s kinetic energy at touchdown, thus reducing the impact load and acceleration transmitted to the airframe. Due to its significant influence on ground loads, this system is generally designed in parallel with the main structural components of the aircraft, such as the fuselage and wings. Robust numerical models for simulating landing gear impact dynamics are essential from the preliminary design stage in order to properly assess aircraft configuration and structural arrangements. Finite element (FE) analysis is a viable solution for supporting the design. However, regarding the oleo-pneumatic struts, FE-based simulation may become unpractical, since detailed models are required to obtain reliable results. Moreover, FE models could not be very versatile for accommodating the many design updates that usually occur at the beginning of the landing gear project or during the layout optimization process. In this work, a numerical method for simulating oleo-pneumatic landing gear drop dynamics is presented. To effectively support both the preliminary and advanced design of landing gear units, the proposed simulation approach rationally balances the level of sophistication of the adopted model with the need for accurate results. Although based on a formulation assuming only four state variables for the description of landing gear dynamics, the approach successfully accounts for all the relevant forces that arise during the drop and their influence on landing gear motion. A set of intercommunicating routines was implemented in MATLAB® environment to integrate the dynamic impact equations, starting from user-defined initial conditions and general parameters related to the geometric and structural configuration of the landing gear. The tool was then used to simulate a drop test of a reference landing gear, and the obtained results were successfully validated against available experimental data.


2012 ◽  
Vol 145 (3) ◽  
pp. 539-550 ◽  
Author(s):  
Hendrik Wouters ◽  
Koen De Ridder ◽  
Nicole P. M. van Lipzig

2020 ◽  
Author(s):  
Bernd Schalge ◽  
Gabriele Baroni ◽  
Barbara Haese ◽  
Daniel Erdal ◽  
Gernot Geppert ◽  
...  

Abstract. Coupled numerical models, which simulate water and energy fluxes in the subsurface-land surface-atmosphere system in a physically consistent way are a prerequisite for the analysis and a better understanding of heat and matter exchange fluxes at compartmental boundaries and interdependencies of states across these boundaries. Complete state evolutions generated by such models may be regarded as a proxy of the real world, provided they are run at sufficiently high resolution and incorporate the most important processes. Such a virtual reality can be used to test hypotheses on the functioning of the coupled terrestrial system. Coupled simulation systems, however, face severe problems caused by the vastly different scales of the processes acting in and between the compartments of the terrestrial system, which also hinders comprehensive tests of their realism. We used the Terrestrial Systems Modeling Platform TerrSysMP, which couples the meteorological model COSMO, the land-surface model CLM, and the subsurface model ParFlow, to generate a virtual catchment for a regional terrestrial system mimicking the Neckar catchment in southwest Germany. Simulations for this catchment are made for the period 2007–2015, and at a spatial resolution of 400 m for the land surface and subsurface and 1.1 km for the atmosphere. Among a discussion of modelling challenges, the model performance is evaluated based on real observations covering several variables of the water cycle. We find that the simulated (virtual) catchment behaves in many aspects quite close to observations of the real Neckar catchment, e.g. concerning atmospheric boundary-layer height, precipitation, and runoff. But also discrepancies become apparent, both in the ability of the model to correctly simulate some processes which still need improvement such as overland flow, and in the realism of some observation operators like the satellite based soil moisture sensors. The whole raw dataset is available for interested users. The dataset described here is available via the CERA database (Schalge et al., 2020): https://doi.org/10.26050/WDCC/Neckar_VCS_v1.


2007 ◽  
Vol 135 (10) ◽  
pp. 3474-3483 ◽  
Author(s):  
Kyung-Ja Ha ◽  
Yu-Kyung Hyun ◽  
Hyun-Mi Oh ◽  
Kyung-Eak Kim ◽  
Larry Mahrt

Abstract The Monin–Obukhov similarity theory and a generalized formulation of the mixing length for the stable boundary layer are evaluated using the Cooperative Atmosphere–Surface Exchange Study-1999 (CASES-99) data. The large-scale wind forcing is classified into weak, intermediate, and strong winds. Although the stability parameter, z/L, is inversely dependent on the mean wind speed, the speed of the large-scale flow includes independent influences on the flux–gradient relationship. The dimensionless mean wind shear is found to obey existing stability functions when z/L is less than unity, particularly for the strong and intermediate wind classes. For weak mean winds and/or strong stability (z/L > 1), this similarity theory breaks down. Deviations from similarity theory are examined in terms of intermittency. A case study of a weak-wind night indicates important modulation of the turbulence flux by mesoscale motions of unknown origin.


Complexity ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-15
Author(s):  
Hao Wang ◽  
Shousheng Xie ◽  
Weixuan Wang ◽  
Lei Wang ◽  
Jingbo Peng

The problem of unmeasured parameters estimation for distributed control systems is studied in this paper. The Takagi–Sugeno fuzzy model which can appropriate any nonlinear systems is employed, and based on the model, an observer-based fuzzy H∞ filter which has robustness against time-delay, external noise, and system uncertainties is designed. The sufficient condition for the existence of the desired filter is derived in terms of linear matrix inequalities (LMIs) solutions. Moreover, the underdetermined estimation problem in which the number of sensors available is typically less than the number of state variables to be estimated is specifically addressed. A systematic method is proposed to produce a model tuning parameter vector of appropriate dimension for the estimation of the filter, and the optimal transformation matrix is selected via iterative solution to minimize the estimated error. Finally, a simulation example for turbofan aeroengine is given to illustrate the effectiveness of the proposed method, and the estimated error is less than 2.5%.


2002 ◽  
Vol 2 (3/4) ◽  
pp. 147-155 ◽  
Author(s):  
Ch. Jaedicke ◽  
A. D. Sandvik

Abstract. Blowing snow and snow drifts are common features in the Arctic. Due to sparse vegetation, low temperatures and high wind speeds, the snow is constantly moving. This causes severe problems for transportation and infrastructure in the affected areas. To minimise the effect of drifting snow already in the designing phase of new structures, adequate models have to be developed and tested. In this study, snow distribution in Arctic topography is surveyed in two study areas during the spring of 1999 and 2000. Snow depth is measured by ground penetrating radar and manual methods. The study areas encompass four by four kilometres and are partly glaciated. The results of the surveys show a clear pattern of erosion, accumulation areas and the evolution of the snow cover over time. This high resolution data set is valuable for the validation of numerical models. A simple numerical snow drift model was used to simulate the measured snow distribution in one of the areas for the winter of 1998/1999. The model is a two-level drift model coupled to the wind field, generated by a mesoscale meteorological model. The simulations are based on five wind fields from the dominating wind directions. The model produces a satisfying snow distribution but fails to reproduce the details of the observed snow cover. The results clearly demonstrate the importance of quality field data to detect and analyse errors in numerical simulations.


Sign in / Sign up

Export Citation Format

Share Document